<div dir="ltr"><div><div><div><div>Dear gradsusr mailing list:<br></div>Sorry to ask multiple questions in the same day, but I think that this might be relatively easy fix.  I am running a 70 year time series of near-surface temp. poleward of 60N.  I have been struggling with one final part of generating the time series.  The final year exhibits a great deal of noise, reflecting seasonality, likely due to an issue with averaging over a data set that is ending.  I have attached both code code and figure.<br>
<br></div><div>For file 1:<br>define meanmonth=(aave(temp.1,lon=-180,lon=180,lat=60,lat=90))<br>*set time jan2020 dec2089<br>define mean=ave(meanmonth,t+0,t+11)<br>d tloop(mean)<br></div><div>From the corresponding .ctl file:<br>
SET /geo1/RCP4.5/RCP4.5/GISS/Temps/tas_Amon_GISS-E2-R_rcp45_r1i1p1_202001-208912.nc<br>DTYPE netcdf<br>TITLE Yearly Weighted Area Avg. Temp.(K)_N.Hem. lat. &gt; 60<br>UNDEF  1.e+20<br>XDEF lon 144 linear 1.25 2.5<br>YDEF lat 90 linear -89 2<br>
TDEF time 840 linear 01jan2020 1mo<br>VARS 1 <br>tas=&gt;temp 0 t,y,x temperature<br>ENDVARS<br><br><br></div><div><br></div>My question is this: where do I redefine or reset t or time in each of the time series on the plot?  Also, I am opening using data from ten different files.  Each data set spans my time interval, but some begin in 2006 rather than 2020.<br>
<br></div>Thanks again,<br></div>Corey<br><div><div><div><div><div><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Jul 11, 2013 at 3:47 PM,  <span dir="ltr">&lt;<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Send gradsusr mailing list submissions to<br>
        <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
or, via email, send a message with subject or body &#39;help&#39; to<br>
        <a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:gradsusr-owner@gradsusr.org">gradsusr-owner@gradsusr.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than &quot;Re: Contents of gradsusr digest...&quot;<br>
<br>
<br>
Today&#39;s Topics:<br>
<br>
   1. Re: help - NetCDF file (Celso Bandeira de Melo Ribeiro)<br>
   2. Re: Least-Squares Regression and Climatology (Time Series)<br>
      (Hosmay Lopez)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Thu, 11 Jul 2013 16:26:00 -0300<br>
From: Celso Bandeira de Melo Ribeiro &lt;<a href="mailto:celso.bandeira@ufjf.edu.br">celso.bandeira@ufjf.edu.br</a>&gt;<br>
Subject: Re: [gradsusr] help - NetCDF file<br>
To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
Message-ID: &lt;1122f292023cec649d99f1fda15546f8@localhost&gt;<br>
Content-Type: text/plain; charset=UTF-8<br>
<br>
Hi James,<br>
Thank you for attention and help again.<br>
I included the line &quot;OPTIONS template&quot;, in amazonia.ctl, as you sugest,<br>
show below:<br>
<br>
DSET ^amazonia.clm2.h0.2002-01-%<a href="http://m2-00000.nc" target="_blank">m2-00000.nc</a><br>
OPTIONS template<br>
OPTIONS 365_day_calendar<br>
tdef time 4 linear 01jan2002 1dy<br>
<br>
The file open perfectily now, see the message below:<br>
<br>
ga-&gt; xdfopen amazonia.ctl<br>
Scanning Descriptor File:  amazonia.ctl<br>
SDF file amazonia.clm2.h0.2002-01-%<a href="http://m2-00000.nc" target="_blank">m2-00000.nc</a> is open as file 1<br>
LON set to 279.062 311.875<br>
LAT set to -20.2999 6.45372<br>
LEV set to 0.00710064 0.00710064<br>
Time values set: 2002:1:1:0 2002:1:1:0<br>
E set to 1 1<br>
ga-&gt;<br>
<br>
But I still could not see the time values ??of four, the following message<br>
appears below:<br>
<br>
ga-&gt; set t 1 4<br>
Time values set: 2002:1:1:0 2002:1:4:0<br>
ga-&gt; d wind<br>
Contouring: 0.5 to 5.5 interval 0.5<br>
NetCDF Error (gancrow, nc_get_vara_double): NetCDF: Index exceeds<br>
dimension bound<br>
Data Request Error:  Error for variable &#39;wind&#39;<br>
  Error ocurred at column 1<br>
DISPLAY error:  Invalid expression<br>
  Expression = wind<br>
ga-&gt;<br>
<br>
Do you know how solve this problem?<br>
<br>
Tank you very much!<br>
Celso<br>
<br>
<br>
On Thu, 11 Jul 2013 07:33:38 -1000, &quot;James T. Potemra&quot; &lt;<a href="mailto:jimp@hawaii.edu">jimp@hawaii.edu</a>&gt;<br>
wrote:<br>
&gt; Hi Celso:<br>
&gt;<br>
&gt; You also have to add the option &quot;template&quot;:<br>
&gt; OPTIONS template<br>
&gt;<br>
&gt; Jim<br>
&gt; On 7/11/13 6:57 AM, Celso Bandeira de Melo Ribeiro wrote:<br>
&gt;&gt; Hi James,<br>
&gt;&gt; Thank you for attention and help.<br>
&gt;&gt; I modified the file data to diaries as show below:<br>
&gt;&gt;<br>
&gt;&gt; DSET ^amazonia.clm2.h0.2002-01-%<a href="http://m2-00000.nc" target="_blank">m2-00000.nc</a><br>
&gt;&gt; OPTIONS 365_day_calendar<br>
&gt;&gt; tdef time 4 linear 01jan2000 1dy<br>
&gt;&gt;<br>
&gt;&gt; But it did not work, appeared the following error message:<br>
&gt;&gt;<br>
&gt;&gt; ga-&gt; xdfopen amazonia.ctl<br>
&gt;&gt; Scanning Descriptor File:  amazonia.ctl<br>
&gt;&gt; Error: nc_open failed to open file<br>
amazonia.clm2.h0.2002-01-%<a href="http://m2-00000.nc" target="_blank">m2-00000.nc</a><br>
&gt;&gt; No such file or directory<br>
&gt;&gt; gadsdf: Couldn&#39;t ingest SDF metadata.<br>
&gt;&gt; SDF Descriptor file amazonia.ctl was not successfully opened &amp; parsed.<br>
&gt;&gt; ga-&gt;<br>
&gt;&gt;<br>
&gt;&gt; All the best!<br>
&gt;&gt;<br>
&gt;&gt; Celso<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; On Wed, 10 Jul 2013 17:06:21 -1000, &quot;James T. Potemra&quot;<br>
&lt;<a href="mailto:jimp@hawaii.edu">jimp@hawaii.edu</a>&gt;<br>
&gt;&gt; wrote:<br>
&gt;&gt;&gt; Celso:<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; Your descriptor file lists four times (months), but it is specifically<br>
&gt;&gt;&gt; reading only one file (<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a>).  You<br>
&gt;&gt;&gt; likely want to read all four, so you&#39;ll have to add an option<br>
&gt;&gt;&gt; &quot;template&quot;.  Also, it&#39;s not clear from the file names if they<br>
reference<br>
&gt;&gt;&gt; month-day or day-month, e.g., is <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a><br>
&gt;&gt;&gt; for January 4, 2002 or April 1, 2002?  In your descriptor file you<br>
&gt;&gt;&gt; specify the time increment as monthly (1mo).  If the files have single<br>
&gt;&gt;&gt; time values (one monthly value), and the format of the file name is<br>
&gt;&gt;&gt; day-month, your DSET line should be something like:<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; DSET ^amazonia.clm2.h0.2002-01-%<a href="http://m2-00000.nc" target="_blank">m2-00000.nc</a><br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; The %m2 indicating a substitution template of 2-digit month.<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; Jim<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; On 7/10/13 3:58 PM, Celso Bandeira de Melo Ribeiro wrote:<br>
&gt;&gt;&gt;&gt; Hi Matt,<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; Thank you very much for your attention and help.<br>
&gt;&gt;&gt;&gt; Worked perfectly I created the following file:<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; DSET ^<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a><br>
&gt;&gt;&gt;&gt; OPTIONS 365_day_calendar<br>
&gt;&gt;&gt;&gt; tdef time 4 linear 01jan2000 1mo<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; I am working with 4 files (Daily), whose names are:<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; <a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a><br>
&gt;&gt;&gt;&gt; <a href="http://amazonia.clm2.h0.2002-01-02-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-02-00000.nc</a><br>
&gt;&gt;&gt;&gt; <a href="http://amazonia.clm2.h0.2002-01-03-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-03-00000.nc</a><br>
&gt;&gt;&gt;&gt; <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a><br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; The file opens perfectly in GrADS:<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; ga-&gt; xdfopen amazonia.ctl<br>
&gt;&gt;&gt;&gt; Scanning Descriptor File:  amazonia.ctl<br>
&gt;&gt;&gt;&gt; SDF file <a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a> is open as file 1<br>
&gt;&gt;&gt;&gt; LON set to 279.062 311.875<br>
&gt;&gt;&gt;&gt; LAT set to -20.2999 6.45372<br>
&gt;&gt;&gt;&gt; LEV set to 0.00710064 0.00710064<br>
&gt;&gt;&gt;&gt; Time values set: 2000:1:1:0 2000:1:1:0<br>
&gt;&gt;&gt;&gt; E set to 1 1<br>
&gt;&gt;&gt;&gt; ga-&gt; q file<br>
&gt;&gt;&gt;&gt; File 1 : CLM History file information<br>
&gt;&gt;&gt;&gt;     Descriptor: amazonia.ctl<br>
&gt;&gt;&gt;&gt;     Binary: <a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a><br>
&gt;&gt;&gt;&gt;     Type = Gridded<br>
&gt;&gt;&gt;&gt;     Xsize = 106  Ysize = 115  Zsize = 15  Tsize = 4  Esize = 1<br>
&gt;&gt;&gt;&gt;     Number of Variables = 175<br>
&gt;&gt;&gt;&gt;       ...<br>
&gt;&gt;&gt;&gt;        tsoi  15  t,z,y,x  soil temperature<br>
&gt;&gt;&gt;&gt;        tsoi_10cm  0  t,y,x  soil temperature in top 10cm of soil<br>
&gt;&gt;&gt;&gt;        tv  0  t,y,x  vegetation temperature<br>
&gt;&gt;&gt;&gt;        u10  0  t,y,x  10-m wind<br>
&gt;&gt;&gt;&gt;        urban_ac  0  t,y,x  urban air conditioning flux<br>
&gt;&gt;&gt;&gt;    ...<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; But when I try to open in GrADS putting time as a variable (set 1 t<br>
4)<br>
&gt;&gt;&gt;&gt; appears the following error message:<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; ga-&gt; set t 1 4<br>
&gt;&gt;&gt;&gt; Time values set: 2000:1:1:0 2000:4:1:0<br>
&gt;&gt;&gt;&gt; ga-&gt; d tsoi<br>
&gt;&gt;&gt;&gt; Contouring: 273 to 279 interval 0.5<br>
&gt;&gt;&gt;&gt; gancgrid error: nc_get_vara_double failed; NetCDF: Index exceeds<br>
&gt;&gt;&gt;&gt; dimension<br>
&gt;&gt;&gt;&gt; bound<br>
&gt;&gt;&gt;&gt; Data Request Error:  Error for variable &#39;tsoi&#39;<br>
&gt;&gt;&gt;&gt;     Error ocurred at column 1<br>
&gt;&gt;&gt;&gt; DISPLAY error:  Invalid expression<br>
&gt;&gt;&gt;&gt;     Expression = tsoi<br>
&gt;&gt;&gt;&gt; ga-&gt;<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; You know how to help me change the time?<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; All the best!<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; Celso<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; On Wed, 10 Jul 2013 12:43:38 -0700, Matt Masarik<br>
&gt;&gt;&gt;&gt; &lt;<a href="mailto:mmasarik@atmos.ucla.edu">mmasarik@atmos.ucla.edu</a>&gt;<br>
&gt;&gt;&gt;&gt; wrote:<br>
&gt;&gt;&gt;&gt;&gt; Hi Celso,<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; you can solve this by using a small (3 line) control file<br>
&gt;&gt;&gt;&gt;&gt; that has the 365_day_calendar option specified.  Copy<br>
&gt;&gt;&gt;&gt;&gt; the following text into a file, say, amazonia.ctl<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; DSET ^<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a><br>
&gt;&gt;&gt;&gt;&gt; OPTIONS 365_day_calendar<br>
&gt;&gt;&gt;&gt;&gt; tdef time 3084 linear 16jan1850 1mo<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; I put your file name in the DSET line already.  You will<br>
&gt;&gt;&gt;&gt;&gt; need to change the 3084 to the total number of time<br>
&gt;&gt;&gt;&gt;&gt; records in the file.  You can find this out by issuing<br>
&gt;&gt;&gt;&gt;&gt; the following at a command prompt:<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; ncdump -c <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a><br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; you will also need to change 16jan1850 to the start<br>
&gt;&gt;&gt;&gt;&gt; date of your data set, again use ncdump.<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; Once you make those changes you should be able to open<br>
&gt;&gt;&gt;&gt;&gt; the file by using<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; xdfopen amazonia.ctl<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; Matt<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; On 07/10/2013 11:50 AM, Celso Bandeira de Melo Ribeiro wrote:<br>
&gt;&gt;&gt;&gt;&gt;&gt; Hello,<br>
&gt;&gt;&gt;&gt;&gt;&gt; I&#39;m working with GrADS Version 2.0.a9 in Linux Ubuntu 64-bit and I<br>
am<br>
&gt;&gt;&gt;&gt; not<br>
&gt;&gt;&gt;&gt;&gt;&gt; able to open NetCDF file.<br>
&gt;&gt;&gt;&gt;&gt;&gt; I am using the function &quot;sdfopen&quot; but an error message appears as<br>
&gt;&gt;&gt;&gt; sending<br>
&gt;&gt;&gt;&gt;&gt;&gt; the attachment.<br>
&gt;&gt;&gt;&gt;&gt;&gt; The message, as follows:<br>
&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt; ga-&gt; sdfopen <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a><br>
&gt;&gt;&gt;&gt;&gt;&gt; Scanning self-describing file:<br>
<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a><br>
&gt;&gt;&gt;&gt;&gt;&gt; SDF Error: 365 day calendars are no longer supported by sdfopen.<br>
&gt;&gt;&gt;&gt;&gt;&gt;      To open this file with GrADS, use a descriptor file with<br>
&gt;&gt;&gt;&gt;&gt;&gt;      a complete TDEF entry and OPTIONS 365_day_calendar.<br>
&gt;&gt;&gt;&gt;&gt;&gt;      Documentation is at<br>
&gt;&gt;&gt;&gt; <a href="http://iges.org/grads/gadoc/SDFdescriptorfile.html" target="_blank">http://iges.org/grads/gadoc/SDFdescriptorfile.html</a><br>
&gt;&gt;&gt;&gt;&gt;&gt; ga-&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt; I would like to know how solve this problem.<br>
&gt;&gt;&gt;&gt;&gt;&gt; Thank you!<br>
&gt;&gt;&gt;&gt;&gt;&gt; Celso<br>
&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; _______________________________________________<br>
&gt;&gt;&gt;&gt;&gt; gradsusr mailing list<br>
&gt;&gt;&gt;&gt;&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt;&gt;&gt;&gt;&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt;&gt;&gt; _______________________________________________<br>
&gt;&gt;&gt; gradsusr mailing list<br>
&gt;&gt;&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt;&gt;&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br>
--<br>
____________________________________________<br>
Prof. Celso Bandeira de Melo Ribeiro<br>
Universidade Federal de Juiz de Fora - UFJF<br>
Faculdade de Engenharia, Plataforma 4<br>
Departamento de Engenharia Sanit?ria e Ambiental - ESA<br>
Bairro Martelos, Juiz de Fora (MG) CEP. 36036-900<br>
tel. (32) 2102-3419, Ramal 30<br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Thu, 11 Jul 2013 19:39:13 +0000<br>
From: Hosmay Lopez &lt;<a href="mailto:hlopez@rsmas.miami.edu">hlopez@rsmas.miami.edu</a>&gt;<br>
Subject: Re: [gradsusr] Least-Squares Regression and Climatology (Time<br>
        Series)<br>
To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
Message-ID:<br>
        &lt;<a href="mailto:8A902EB84795C94C8726F11DC8C574A81D26BA7B@EXCHANGE1.amazon.rsmas.miami.edu">8A902EB84795C94C8726F11DC8C574A81D26BA7B@EXCHANGE1.amazon.rsmas.miami.edu</a>&gt;<br>
<br>
Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
<br>
Dear Corey,<br>
<br>
There is a problem with how you define the seasonal cycle on your code, specifically line 3. You have:<br>
define nino34clim=ave(nino34,time=jan2020,time=dec2089,12)<br>
<br>
where the averaging should be:<br>
define nino34clim=ave(nino34,t+0,time=dec2089,12)<br>
<br>
Please note the 2nd argument on the average function. It should average through all months, skipping by 12, but should include all times from jan2020 to dec2089.<br>
Hope this helps!!!<br>
<br>
Hosmay Lopez<br>
Rosenstiel School of Marine and Atmospheric Science<br>
University of Miami<br>
4600 Rickenbacker Causeway Miami, FL 33149<br>
Office: MSC337<br>
Phone: (305) 298-8006<br>
________________________________<br>
From: <a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a> [<a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>] on behalf of Corey Gabriel [<a href="mailto:cjgabriel7@gmail.com">cjgabriel7@gmail.com</a>]<br>

Sent: Thursday, July 11, 2013 2:51 PM<br>
To: gradsusr<br>
Subject: [gradsusr] Least-Squares Regression and Climatology (Time Series)<br>
<br>
***Identical to previous message.  Subject line included this time***<br>
<br>
Dear gradsusr mailing list:<br>
<br>
I am attempting to perform a least-squares regression of a 3-month running<br>
mean time series of sea surface temperatures over the Nino 3.4 region over<br>
the 70 year period 2020 to 2089.  My goal is to diagnose significant events<br>
over that 70 year period by performing the least-squares regression.  I<br>
need to remove the seasonal cycle, and then compare the time series to the<br>
regression line.  A couple problems:<br>
1.) GrADS will not allow me to use time as the ind. var. or dep. var. in<br>
the tregr() function or when defining the coeff.(b).<br>
2.) I attempted to invoke climatology, based on the documentation provided<br>
here:<a href="http://grads.wikidot.com/algorithms" target="_blank">http://grads.wikidot.com/algorithms</a>.  However, I cannot seem to remove<br>
the seasonality.<br>
<br>
  Please find the relevant portion of my code copied below.  An alternate<br>
version of my code is also copied below, using a different variable.  I<br>
have also attached the anomaly figure that code2(see below) generates.<br>
<br>
Any assistance would be tremendously appreciated!<br>
Corey<br>
<br>
*Code1: *&#39;et time jan2020 dec2020<br>
define nino34=(aave((sst-273.15),lon=<br>
120,lon=170,lat=-5,lat=5))<br>
define nino34clim=ave(nino34,time=jan2020,time=dec2089,12)<br>
modify nino34clim seasonal<br>
define b=tregr(t,nino34,t=1,t=840)<br>
define linearnino34clim= ( b*t + 420) + nino34clim<br>
d linearnino34clim<br>
c<br>
set lat 1<br>
set lon 1<br>
set grads off<br>
set cstyle 3<br>
set ccolor 3<br>
set cmark 0<br>
d nino34-linearnino34clim&#39;<br>
<br>
*Code2(figure attached)<br>
*&#39;set time jan2020 dec2089<br>
set t 1 840<br>
define nino34=(aave((sst-273.15),lon=120,lon=170,lat=-5,lat=5))<br>
define nino34clim=ave(nino34,time=jan2020,time=dec2089,12)<br>
modify nino34clim seasonal<br>
define b=tregr(nino34,nino34clim,t=1,t=840)<br>
define linearnino34clim=(b*nino34+nino34clim)<br>
d linearnino34clim<br>
c<br>
set grads off<br>
set lat 1<br>
set lon 1<br>
set cstyle 3<br>
set ccolor 3<br>
set cmark 0<br>
d nino34-linearnino34clim&#39;<br>
*<br>
<br>
<br>
On Thu, Jul 11, 2013 at 2:22 PM, &lt;<a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a>&gt;&gt; wrote:<br>

Send gradsusr mailing list submissions to<br>
        <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
or, via email, send a message with subject or body &#39;help&#39; to<br>
        <a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a>&gt;<br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:gradsusr-owner@gradsusr.org">gradsusr-owner@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr-owner@gradsusr.org">gradsusr-owner@gradsusr.org</a>&gt;<br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than &quot;Re: Contents of gradsusr digest...&quot;<br>
<br>
<br>
Today&#39;s Topics:<br>
<br>
   1. Re: gradsusr Digest, Vol 41, Issue 19 (Corey Gabriel)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Thu, 11 Jul 2013 14:22:07 -0400<br>
From: Corey Gabriel &lt;<a href="mailto:cjgabriel7@gmail.com">cjgabriel7@gmail.com</a>&lt;mailto:<a href="mailto:cjgabriel7@gmail.com">cjgabriel7@gmail.com</a>&gt;&gt;<br>
Subject: Re: [gradsusr] gradsusr Digest, Vol 41, Issue 19<br>
To: <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
Message-ID:<br>
        &lt;CAEhnGEru0gWKX=<a href="mailto:tSMgLyWCm3iFgMthwhxkEdnnY6yTTHANZA1A@mail.gmail.com">tSMgLyWCm3iFgMthwhxkEdnnY6yTTHANZA1A@mail.gmail.com</a>&lt;mailto:<a href="mailto:tSMgLyWCm3iFgMthwhxkEdnnY6yTTHANZA1A@mail.gmail.com">tSMgLyWCm3iFgMthwhxkEdnnY6yTTHANZA1A@mail.gmail.com</a>&gt;&gt;<br>

Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
<br>
Dear gradsusr mailing list:<br>
<br>
I am attempting to perform a least-squares regression of a 3-month running<br>
mean time series of sea surface temperatures over the Nino 3.4 region over<br>
the 70 year period 2020 to 2089.  My goal is to diagnose significant events<br>
over that 70 year period by performing the least-squares regression.  I<br>
need to remove the seasonal cycle, and then compare the time series to the<br>
regression line.  A couple problems:<br>
1.) GrADS will not allow me to use time as the ind. var. or dep. var. in<br>
the tregr() function or when defining the coeff.(b).<br>
2.) I attempted to invoke climatology, based on the documentation provided<br>
here:<a href="http://grads.wikidot.com/algorithms" target="_blank">http://grads.wikidot.com/algorithms</a>.  However, I cannot seem to remove<br>
the seasonality.<br>
<br>
  Please find the relevant portion of my code copied below.  An alternate<br>
version of my code is also copied below, using a different variable.  I<br>
have also attached the anomaly figure that code2(see below) generates.<br>
<br>
Any assistance would be tremendously appreciated!<br>
Corey<br>
<br>
*Code1: *&#39;et time jan2020 dec2020<br>
define nino34=(aave((sst-273.15),lon=120,lon=170,lat=-5,lat=5))<br>
define nino34clim=ave(nino34,time=jan2020,time=dec2089,12)<br>
modify nino34clim seasonal<br>
define b=tregr(t,nino34,t=1,t=840)<br>
define linearnino34clim= ( b*t + 420) + nino34clim<br>
d linearnino34clim<br>
c<br>
set lat 1<br>
set lon 1<br>
set grads off<br>
set cstyle 3<br>
set ccolor 3<br>
set cmark 0<br>
d nino34-linearnino34clim&#39;<br>
<br>
*Code2(figure attached)<br>
*&#39;set time jan2020 dec2089<br>
set t 1 840<br>
define nino34=(aave((sst-273.15),lon=120,lon=170,lat=-5,lat=5))<br>
define nino34clim=ave(nino34,time=jan2020,time=dec2089,12)<br>
modify nino34clim seasonal<br>
define b=tregr(nino34,nino34clim,t=1,t=840)<br>
define linearnino34clim=(b*nino34+nino34clim)<br>
d linearnino34clim<br>
c<br>
set grads off<br>
set lat 1<br>
set lon 1<br>
set cstyle 3<br>
set ccolor 3<br>
set cmark 0<br>
d nino34-linearnino34clim&#39;<br>
*<br>
*<br>
<br>
<br>
<br>
<br>
On Thu, Jul 11, 2013 at 1:40 PM, &lt;<a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a>&gt;&gt; wrote:<br>

<br>
&gt; Send gradsusr mailing list submissions to<br>
&gt;         <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
&gt;<br>
&gt; To subscribe or unsubscribe via the World Wide Web, visit<br>
&gt;         <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt; or, via email, send a message with subject or body &#39;help&#39; to<br>
&gt;         <a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a>&gt;<br>
&gt;<br>
&gt; You can reach the person managing the list at<br>
&gt;         <a href="mailto:gradsusr-owner@gradsusr.org">gradsusr-owner@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr-owner@gradsusr.org">gradsusr-owner@gradsusr.org</a>&gt;<br>
&gt;<br>
&gt; When replying, please edit your Subject line so it is more specific<br>
&gt; than &quot;Re: Contents of gradsusr digest...&quot;<br>
&gt;<br>
&gt;<br>
&gt; Today&#39;s Topics:<br>
&gt;<br>
&gt;    1. Re: help - NetCDF file (Celso Bandeira de Melo Ribeiro)<br>
&gt;    2. Re: Shift in the lattitude while plotting in grads<br>
&gt;       mercator/latlon projction (Mohsen Soltani)<br>
&gt;    3. Re: Shift in the lattitude while plotting in grads<br>
&gt;       mercator/latlon projction (itesh dash)<br>
&gt;    4. Re: help - NetCDF file (James T. Potemra)<br>
&gt;    5. Re: Shift in the lattitude while plotting in grads<br>
&gt;       mercator/latlon projction (James T. Potemra)<br>
&gt;<br>
&gt;<br>
&gt; ----------------------------------------------------------------------<br>
&gt;<br>
&gt; Message: 1<br>
&gt; Date: Thu, 11 Jul 2013 13:57:21 -0300<br>
&gt; From: Celso Bandeira de Melo Ribeiro &lt;<a href="mailto:celso.bandeira@ufjf.edu.br">celso.bandeira@ufjf.edu.br</a>&lt;mailto:<a href="mailto:celso.bandeira@ufjf.edu.br">celso.bandeira@ufjf.edu.br</a>&gt;&gt;<br>
&gt; Subject: Re: [gradsusr] help - NetCDF file<br>
&gt; To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;&gt;<br>
&gt; Message-ID: &lt;34fb6a6e970b1fff4dc36f87278e432e@localhost&gt;<br>
&gt; Content-Type: text/plain; charset=UTF-8<br>
&gt;<br>
&gt; Hi James,<br>
&gt; Thank you for attention and help.<br>
&gt; I modified the file data to diaries as show below:<br>
&gt;<br>
&gt; DSET ^amazonia.clm2.h0.2002-01-%<a href="http://m2-00000.nc" target="_blank">m2-00000.nc</a>&lt;<a href="http://m2-00000.nc" target="_blank">http://m2-00000.nc</a>&gt;<br>
&gt; OPTIONS 365_day_calendar<br>
&gt; tdef time 4 linear 01jan2000 1dy<br>
&gt;<br>
&gt; But it did not work, appeared the following error message:<br>
&gt;<br>
&gt; ga-&gt; xdfopen amazonia.ctl<br>
&gt; Scanning Descriptor File:  amazonia.ctl<br>
&gt; Error: nc_open failed to open file amazonia.clm2.h0.2002-01-%<a href="http://m2-00000.nc" target="_blank">m2-00000.nc</a>&lt;<a href="http://m2-00000.nc" target="_blank">http://m2-00000.nc</a>&gt;<br>
&gt; No such file or directory<br>
&gt; gadsdf: Couldn&#39;t ingest SDF metadata.<br>
&gt; SDF Descriptor file amazonia.ctl was not successfully opened &amp; parsed.<br>
&gt; ga-&gt;<br>
&gt;<br>
&gt; All the best!<br>
&gt;<br>
&gt; Celso<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; On Wed, 10 Jul 2013 17:06:21 -1000, &quot;James T. Potemra&quot; &lt;<a href="mailto:jimp@hawaii.edu">jimp@hawaii.edu</a>&lt;mailto:<a href="mailto:jimp@hawaii.edu">jimp@hawaii.edu</a>&gt;&gt;<br>
&gt; wrote:<br>
&gt; &gt; Celso:<br>
&gt; &gt;<br>
&gt; &gt; Your descriptor file lists four times (months), but it is specifically<br>
&gt; &gt; reading only one file (<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-01-00000.nc</a>&gt;).  You<br>

&gt; &gt; likely want to read all four, so you&#39;ll have to add an option<br>
&gt; &gt; &quot;template&quot;.  Also, it&#39;s not clear from the file names if they reference<br>
&gt; &gt; month-day or day-month, e.g., is <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-04-00000.nc</a>&gt;<br>

&gt; &gt; for January 4, 2002 or April 1, 2002?  In your descriptor file you<br>
&gt; &gt; specify the time increment as monthly (1mo).  If the files have single<br>
&gt; &gt; time values (one monthly value), and the format of the file name is<br>
&gt; &gt; day-month, your DSET line should be something like:<br>
&gt; &gt;<br>
&gt; &gt; DSET ^amazonia.clm2.h0.2002-01-%<a href="http://m2-00000.nc" target="_blank">m2-00000.nc</a>&lt;<a href="http://m2-00000.nc" target="_blank">http://m2-00000.nc</a>&gt;<br>
&gt; &gt;<br>
&gt; &gt; The %m2 indicating a substitution template of 2-digit month.<br>
&gt; &gt;<br>
&gt; &gt; Jim<br>
&gt; &gt;<br>
&gt; &gt; On 7/10/13 3:58 PM, Celso Bandeira de Melo Ribeiro wrote:<br>
&gt; &gt;&gt; Hi Matt,<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; Thank you very much for your attention and help.<br>
&gt; &gt;&gt; Worked perfectly I created the following file:<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; DSET ^<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-01-00000.nc</a>&gt;<br>

&gt; &gt;&gt; OPTIONS 365_day_calendar<br>
&gt; &gt;&gt; tdef time 4 linear 01jan2000 1mo<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; I am working with 4 files (Daily), whose names are:<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; <a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-01-00000.nc</a>&gt;<br>

&gt; &gt;&gt; <a href="http://amazonia.clm2.h0.2002-01-02-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-02-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-02-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-02-00000.nc</a>&gt;<br>

&gt; &gt;&gt; <a href="http://amazonia.clm2.h0.2002-01-03-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-03-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-03-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-03-00000.nc</a>&gt;<br>

&gt; &gt;&gt; <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-04-00000.nc</a>&gt;<br>

&gt; &gt;&gt;<br>
&gt; &gt;&gt; The file opens perfectly in GrADS:<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; ga-&gt; xdfopen amazonia.ctl<br>
&gt; &gt;&gt; Scanning Descriptor File:  amazonia.ctl<br>
&gt; &gt;&gt; SDF file <a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-01-00000.nc</a>&gt; is open as file 1<br>

&gt; &gt;&gt; LON set to 279.062 311.875<br>
&gt; &gt;&gt; LAT set to -20.2999 6.45372<br>
&gt; &gt;&gt; LEV set to 0.00710064 0.00710064<br>
&gt; &gt;&gt; Time values set: 2000:1:1:0 2000:1:1:0<br>
&gt; &gt;&gt; E set to 1 1<br>
&gt; &gt;&gt; ga-&gt; q file<br>
&gt; &gt;&gt; File 1 : CLM History file information<br>
&gt; &gt;&gt;    Descriptor: amazonia.ctl<br>
&gt; &gt;&gt;    Binary: <a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-01-00000.nc</a>&gt;<br>

&gt; &gt;&gt;    Type = Gridded<br>
&gt; &gt;&gt;    Xsize = 106  Ysize = 115  Zsize = 15  Tsize = 4  Esize = 1<br>
&gt; &gt;&gt;    Number of Variables = 175<br>
&gt; &gt;&gt;      ...<br>
&gt; &gt;&gt;       tsoi  15  t,z,y,x  soil temperature<br>
&gt; &gt;&gt;       tsoi_10cm  0  t,y,x  soil temperature in top 10cm of soil<br>
&gt; &gt;&gt;       tv  0  t,y,x  vegetation temperature<br>
&gt; &gt;&gt;       u10  0  t,y,x  10-m wind<br>
&gt; &gt;&gt;       urban_ac  0  t,y,x  urban air conditioning flux<br>
&gt; &gt;&gt;   ...<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; But when I try to open in GrADS putting time as a variable (set 1 t 4)<br>
&gt; &gt;&gt; appears the following error message:<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; ga-&gt; set t 1 4<br>
&gt; &gt;&gt; Time values set: 2000:1:1:0 2000:4:1:0<br>
&gt; &gt;&gt; ga-&gt; d tsoi<br>
&gt; &gt;&gt; Contouring: 273 to 279 interval 0.5<br>
&gt; &gt;&gt; gancgrid error: nc_get_vara_double failed; NetCDF: Index exceeds<br>
&gt; &gt;&gt; dimension<br>
&gt; &gt;&gt; bound<br>
&gt; &gt;&gt; Data Request Error:  Error for variable &#39;tsoi&#39;<br>
&gt; &gt;&gt;    Error ocurred at column 1<br>
&gt; &gt;&gt; DISPLAY error:  Invalid expression<br>
&gt; &gt;&gt;    Expression = tsoi<br>
&gt; &gt;&gt; ga-&gt;<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; You know how to help me change the time?<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; All the best!<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; Celso<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; On Wed, 10 Jul 2013 12:43:38 -0700, Matt Masarik<br>
&gt; &gt;&gt; &lt;<a href="mailto:mmasarik@atmos.ucla.edu">mmasarik@atmos.ucla.edu</a>&lt;mailto:<a href="mailto:mmasarik@atmos.ucla.edu">mmasarik@atmos.ucla.edu</a>&gt;&gt;<br>
&gt; &gt;&gt; wrote:<br>
&gt; &gt;&gt;&gt; Hi Celso,<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; you can solve this by using a small (3 line) control file<br>
&gt; &gt;&gt;&gt; that has the 365_day_calendar option specified.  Copy<br>
&gt; &gt;&gt;&gt; the following text into a file, say, amazonia.ctl<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; DSET ^<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-04-00000.nc</a>&gt;<br>

&gt; &gt;&gt;&gt; OPTIONS 365_day_calendar<br>
&gt; &gt;&gt;&gt; tdef time 3084 linear 16jan1850 1mo<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; I put your file name in the DSET line already.  You will<br>
&gt; &gt;&gt;&gt; need to change the 3084 to the total number of time<br>
&gt; &gt;&gt;&gt; records in the file.  You can find this out by issuing<br>
&gt; &gt;&gt;&gt; the following at a command prompt:<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; ncdump -c <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-04-00000.nc</a>&gt;<br>

&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; you will also need to change 16jan1850 to the start<br>
&gt; &gt;&gt;&gt; date of your data set, again use ncdump.<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; Once you make those changes you should be able to open<br>
&gt; &gt;&gt;&gt; the file by using<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; xdfopen amazonia.ctl<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; Matt<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; On 07/10/2013 11:50 AM, Celso Bandeira de Melo Ribeiro wrote:<br>
&gt; &gt;&gt;&gt;&gt; Hello,<br>
&gt; &gt;&gt;&gt;&gt; I&#39;m working with GrADS Version 2.0.a9 in Linux Ubuntu 64-bit and I am<br>
&gt; &gt;&gt; not<br>
&gt; &gt;&gt;&gt;&gt; able to open NetCDF file.<br>
&gt; &gt;&gt;&gt;&gt; I am using the function &quot;sdfopen&quot; but an error message appears as<br>
&gt; &gt;&gt; sending<br>
&gt; &gt;&gt;&gt;&gt; the attachment.<br>
&gt; &gt;&gt;&gt;&gt; The message, as follows:<br>
&gt; &gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt; ga-&gt; sdfopen <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-04-00000.nc</a>&gt;<br>

&gt; &gt;&gt;&gt;&gt; Scanning self-describing file:  <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-04-00000.nc</a>&gt;<br>

&gt; &gt;&gt;&gt;&gt; SDF Error: 365 day calendars are no longer supported by sdfopen.<br>
&gt; &gt;&gt;&gt;&gt;     To open this file with GrADS, use a descriptor file with<br>
&gt; &gt;&gt;&gt;&gt;     a complete TDEF entry and OPTIONS 365_day_calendar.<br>
&gt; &gt;&gt;&gt;&gt;     Documentation is at<br>
&gt; &gt;&gt; <a href="http://iges.org/grads/gadoc/SDFdescriptorfile.html" target="_blank">http://iges.org/grads/gadoc/SDFdescriptorfile.html</a><br>
&gt; &gt;&gt;&gt;&gt; ga-&gt;<br>
&gt; &gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt; I would like to know how solve this problem.<br>
&gt; &gt;&gt;&gt;&gt; Thank you!<br>
&gt; &gt;&gt;&gt;&gt; Celso<br>
&gt; &gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; _______________________________________________<br>
&gt; &gt;&gt;&gt; gradsusr mailing list<br>
&gt; &gt;&gt;&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
&gt; &gt;&gt;&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt; &gt;<br>
&gt; &gt; _______________________________________________<br>
&gt; &gt; gradsusr mailing list<br>
&gt; &gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
&gt; &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt;<br>
&gt; --<br>
&gt; ____________________________________________<br>
&gt; Prof. Celso Bandeira de Melo Ribeiro<br>
&gt; Universidade Federal de Juiz de Fora - UFJF<br>
&gt; Faculdade de Engenharia, Plataforma 4<br>
&gt; Departamento de Engenharia Sanit?ria e Ambiental - ESA<br>
&gt; Bairro Martelos, Juiz de Fora (MG) CEP. 36036-900<br>
&gt; tel. (32) 2102-3419, Ramal 30<br>
&gt;<br>
&gt;<br>
&gt; ------------------------------<br>
&gt;<br>
&gt; Message: 2<br>
&gt; Date: Thu, 11 Jul 2013 09:58:42 -0700<br>
&gt; From: Mohsen Soltani &lt;<a href="mailto:soltani.clima@gmail.com">soltani.clima@gmail.com</a>&lt;mailto:<a href="mailto:soltani.clima@gmail.com">soltani.clima@gmail.com</a>&gt;&gt;<br>
&gt; Subject: Re: [gradsusr] Shift in the lattitude while plotting in grads<br>
&gt;         mercator/latlon projction<br>
&gt; To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;&gt;<br>
&gt; Message-ID:<br>
&gt;         &lt;<br>
&gt; <a href="mailto:CAPoCDiM50XFQjKOHUN86kCNNdB3kp69YFaJYb5emWgMyYLpCOA@mail.gmail.com">CAPoCDiM50XFQjKOHUN86kCNNdB3kp69YFaJYb5emWgMyYLpCOA@mail.gmail.com</a>&lt;mailto:<a href="mailto:CAPoCDiM50XFQjKOHUN86kCNNdB3kp69YFaJYb5emWgMyYLpCOA@mail.gmail.com">CAPoCDiM50XFQjKOHUN86kCNNdB3kp69YFaJYb5emWgMyYLpCOA@mail.gmail.com</a>&gt;&gt;<br>

&gt; Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
&gt;<br>
&gt; itesh,<br>
&gt;<br>
&gt; Make sure Lat/Lon in your .ctl file are the same as your script!<br>
&gt;<br>
&gt; Mohsen<br>
&gt;<br>
&gt; --<br>
&gt; April showers bring May flowers!<br>
&gt; --<br>
&gt; *Mohsen Soltani, M.Sc.*<br>
&gt; Climatology Grad Student,<br>
&gt; Department of Climatology, Faculty of Geography,<br>
&gt; University of Tehran, Tehran, Iran<br>
&gt; voice: <a href="tel:%28%2B98%29%209119772934" value="+989119772934">(+98) 9119772934</a>&lt;tel:%28%2B98%29%209119772934&gt;<br>
&gt; e-mail: <a href="mailto:soltani.clima@gmail.com">soltani.clima@gmail.com</a>&lt;mailto:<a href="mailto:soltani.clima@gmail.com">soltani.clima@gmail.com</a>&gt;<br>
&gt;<br>
&gt;<br>
&gt; On Thu, Jul 11, 2013 at 9:05 AM, itesh dash &lt;<a href="mailto:itesh@rimes.int">itesh@rimes.int</a>&lt;mailto:<a href="mailto:itesh@rimes.int">itesh@rimes.int</a>&gt;&gt; wrote:<br>
&gt;<br>
&gt; &gt; Dear Sir,<br>
&gt; &gt;<br>
&gt; &gt; I am working with the WRF and tried to plot the output in grads. I found<br>
&gt; &gt; there is nearly 1 degree shift in the projection and the plot.<br>
&gt; &gt;<br>
&gt; &gt; Please guide how to correct this.<br>
&gt; &gt;<br>
&gt; &gt; I am using mercator projection for the WRF model. Please find attached<br>
&gt; &gt; image.<br>
&gt; &gt;<br>
&gt; &gt; Seeking your advice on the same.<br>
&gt; &gt;<br>
&gt; &gt; regards,<br>
&gt; &gt; itesh<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; _______________________________________________<br>
&gt; &gt; gradsusr mailing list<br>
&gt; &gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
&gt; &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; -------------- next part --------------<br>
&gt; An HTML attachment was scrubbed...<br>
&gt; URL:<br>
&gt; <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20130711/7cf013d5/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20130711/7cf013d5/attachment-0001.html</a><br>
&gt;<br>
&gt; ------------------------------<br>
&gt;<br>
&gt; Message: 3<br>
&gt; Date: Fri, 12 Jul 2013 00:31:29 +0700<br>
&gt; From: itesh dash &lt;<a href="mailto:itesh@rimes.int">itesh@rimes.int</a>&lt;mailto:<a href="mailto:itesh@rimes.int">itesh@rimes.int</a>&gt;&gt;<br>
&gt; Subject: Re: [gradsusr] Shift in the lattitude while plotting in grads<br>
&gt;         mercator/latlon projction<br>
&gt; To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;&gt;<br>
&gt; Message-ID:<br>
&gt;         &lt;CAFTUJN465mGcZ8k=<br>
&gt; <a href="mailto:JDB1t6UvJRG9jO%2B1L%2BfAD20Gd4G8X_Qayg@mail.gmail.com">JDB1t6UvJRG9jO+1L+fAD20Gd4G8X_Qayg@mail.gmail.com</a>&lt;mailto:<a href="mailto:JDB1t6UvJRG9jO%252B1L%252BfAD20Gd4G8X_Qayg@mail.gmail.com">JDB1t6UvJRG9jO%2B1L%2BfAD20Gd4G8X_Qayg@mail.gmail.com</a>&gt;&gt;<br>

&gt; Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
&gt;<br>
&gt; Dear Mohsen,<br>
&gt;<br>
&gt; thanks for the response. But I am not even trying to use any script. Just<br>
&gt; by opening and plotting the ctl file.<br>
&gt;<br>
&gt; Here is my ctl file<br>
&gt;<br>
&gt; dset ^Jul1013_12.dat<br>
&gt; undef 1.e35<br>
&gt; xdef  639 linear   61.1815   0.0935<br>
&gt; ydef  459 linear   -3.1006   0.0873<br>
&gt; zdef    3 levels<br>
&gt;   850.00000<br>
&gt;   500.00000<br>
&gt;   200.00000<br>
&gt; tdef         8 linear 12z10jul2013  6hr<br>
&gt; vars 1<br>
&gt; HGT 0 0 Terrain Height<br>
&gt; endvars<br>
&gt;<br>
&gt; The wrf model ran with mercator projection 640 x 460 grid points 9 km<br>
&gt; resolution.<br>
&gt;<br>
&gt; I could not figure out why the shift is there. Is there a way i can balance<br>
&gt; the shift.<br>
&gt;<br>
&gt; regards,<br>
&gt; itesh<br>
&gt;<br>
&gt;<br>
&gt; On Thu, Jul 11, 2013 at 11:58 PM, Mohsen Soltani &lt;<a href="mailto:soltani.clima@gmail.com">soltani.clima@gmail.com</a>&lt;mailto:<a href="mailto:soltani.clima@gmail.com">soltani.clima@gmail.com</a>&gt;<br>
&gt; &gt;wrote:<br>
&gt;<br>
&gt; &gt; itesh,<br>
&gt; &gt;<br>
&gt; &gt; Make sure Lat/Lon in your .ctl file are the same as your script!<br>
&gt; &gt;<br>
&gt; &gt; Mohsen<br>
&gt; &gt;<br>
&gt; &gt; --<br>
&gt; &gt; April showers bring May flowers!<br>
&gt; &gt; --<br>
&gt; &gt; *Mohsen Soltani, M.Sc.*<br>
&gt; &gt; Climatology Grad Student,<br>
&gt; &gt; Department of Climatology, Faculty of Geography,<br>
&gt; &gt; University of Tehran, Tehran, Iran<br>
&gt; &gt; voice: <a href="tel:%28%2B98%29%209119772934" value="+989119772934">(+98) 9119772934</a>&lt;tel:%28%2B98%29%209119772934&gt;<br>
&gt; &gt; e-mail: <a href="mailto:soltani.clima@gmail.com">soltani.clima@gmail.com</a>&lt;mailto:<a href="mailto:soltani.clima@gmail.com">soltani.clima@gmail.com</a>&gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; On Thu, Jul 11, 2013 at 9:05 AM, itesh dash &lt;<a href="mailto:itesh@rimes.int">itesh@rimes.int</a>&lt;mailto:<a href="mailto:itesh@rimes.int">itesh@rimes.int</a>&gt;&gt; wrote:<br>
&gt; &gt;<br>
&gt; &gt;&gt; Dear Sir,<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; I am working with the WRF and tried to plot the output in grads. I found<br>
&gt; &gt;&gt; there is nearly 1 degree shift in the projection and the plot.<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; Please guide how to correct this.<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; I am using mercator projection for the WRF model. Please find attached<br>
&gt; &gt;&gt; image.<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; Seeking your advice on the same.<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; regards,<br>
&gt; &gt;&gt; itesh<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; _______________________________________________<br>
&gt; &gt;&gt; gradsusr mailing list<br>
&gt; &gt;&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
&gt; &gt;&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt;<br>
&gt; &gt;<br>
&gt; &gt; _______________________________________________<br>
&gt; &gt; gradsusr mailing list<br>
&gt; &gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
&gt; &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; -------------- next part --------------<br>
&gt; An HTML attachment was scrubbed...<br>
&gt; URL:<br>
&gt; <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20130712/9a256213/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20130712/9a256213/attachment-0001.html</a><br>
&gt;<br>
&gt; ------------------------------<br>
&gt;<br>
&gt; Message: 4<br>
&gt; Date: Thu, 11 Jul 2013 07:33:38 -1000<br>
&gt; From: &quot;James T. Potemra&quot; &lt;<a href="mailto:jimp@hawaii.edu">jimp@hawaii.edu</a>&lt;mailto:<a href="mailto:jimp@hawaii.edu">jimp@hawaii.edu</a>&gt;&gt;<br>
&gt; Subject: Re: [gradsusr] help - NetCDF file<br>
&gt; To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;&gt;<br>
&gt; Message-ID: &lt;<a href="mailto:51DEEC72.8050809@hawaii.edu">51DEEC72.8050809@hawaii.edu</a>&lt;mailto:<a href="mailto:51DEEC72.8050809@hawaii.edu">51DEEC72.8050809@hawaii.edu</a>&gt;&gt;<br>
&gt; Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>
&gt;<br>
&gt; Hi Celso:<br>
&gt;<br>
&gt; You also have to add the option &quot;template&quot;:<br>
&gt; OPTIONS template<br>
&gt;<br>
&gt; Jim<br>
&gt; On 7/11/13 6:57 AM, Celso Bandeira de Melo Ribeiro wrote:<br>
&gt; &gt; Hi James,<br>
&gt; &gt; Thank you for attention and help.<br>
&gt; &gt; I modified the file data to diaries as show below:<br>
&gt; &gt;<br>
&gt; &gt; DSET ^amazonia.clm2.h0.2002-01-%<a href="http://m2-00000.nc" target="_blank">m2-00000.nc</a>&lt;<a href="http://m2-00000.nc" target="_blank">http://m2-00000.nc</a>&gt;<br>
&gt; &gt; OPTIONS 365_day_calendar<br>
&gt; &gt; tdef time 4 linear 01jan2000 1dy<br>
&gt; &gt;<br>
&gt; &gt; But it did not work, appeared the following error message:<br>
&gt; &gt;<br>
&gt; &gt; ga-&gt; xdfopen amazonia.ctl<br>
&gt; &gt; Scanning Descriptor File:  amazonia.ctl<br>
&gt; &gt; Error: nc_open failed to open file amazonia.clm2.h0.2002-01-%<a href="http://m2-00000.nc" target="_blank">m2-00000.nc</a>&lt;<a href="http://m2-00000.nc" target="_blank">http://m2-00000.nc</a>&gt;<br>
&gt; &gt; No such file or directory<br>
&gt; &gt; gadsdf: Couldn&#39;t ingest SDF metadata.<br>
&gt; &gt; SDF Descriptor file amazonia.ctl was not successfully opened &amp; parsed.<br>
&gt; &gt; ga-&gt;<br>
&gt; &gt;<br>
&gt; &gt; All the best!<br>
&gt; &gt;<br>
&gt; &gt; Celso<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; On Wed, 10 Jul 2013 17:06:21 -1000, &quot;James T. Potemra&quot; &lt;<a href="mailto:jimp@hawaii.edu">jimp@hawaii.edu</a>&lt;mailto:<a href="mailto:jimp@hawaii.edu">jimp@hawaii.edu</a>&gt;&gt;<br>
&gt; &gt; wrote:<br>
&gt; &gt;&gt; Celso:<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; Your descriptor file lists four times (months), but it is specifically<br>
&gt; &gt;&gt; reading only one file (<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-01-00000.nc</a>&gt;).  You<br>

&gt; &gt;&gt; likely want to read all four, so you&#39;ll have to add an option<br>
&gt; &gt;&gt; &quot;template&quot;.  Also, it&#39;s not clear from the file names if they reference<br>
&gt; &gt;&gt; month-day or day-month, e.g., is <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-04-00000.nc</a>&gt;<br>

&gt; &gt;&gt; for January 4, 2002 or April 1, 2002?  In your descriptor file you<br>
&gt; &gt;&gt; specify the time increment as monthly (1mo).  If the files have single<br>
&gt; &gt;&gt; time values (one monthly value), and the format of the file name is<br>
&gt; &gt;&gt; day-month, your DSET line should be something like:<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; DSET ^amazonia.clm2.h0.2002-01-%<a href="http://m2-00000.nc" target="_blank">m2-00000.nc</a>&lt;<a href="http://m2-00000.nc" target="_blank">http://m2-00000.nc</a>&gt;<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; The %m2 indicating a substitution template of 2-digit month.<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; Jim<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; On 7/10/13 3:58 PM, Celso Bandeira de Melo Ribeiro wrote:<br>
&gt; &gt;&gt;&gt; Hi Matt,<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; Thank you very much for your attention and help.<br>
&gt; &gt;&gt;&gt; Worked perfectly I created the following file:<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; DSET ^<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-01-00000.nc</a>&gt;<br>

&gt; &gt;&gt;&gt; OPTIONS 365_day_calendar<br>
&gt; &gt;&gt;&gt; tdef time 4 linear 01jan2000 1mo<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; I am working with 4 files (Daily), whose names are:<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; <a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-01-00000.nc</a>&gt;<br>

&gt; &gt;&gt;&gt; <a href="http://amazonia.clm2.h0.2002-01-02-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-02-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-02-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-02-00000.nc</a>&gt;<br>

&gt; &gt;&gt;&gt; <a href="http://amazonia.clm2.h0.2002-01-03-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-03-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-03-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-03-00000.nc</a>&gt;<br>

&gt; &gt;&gt;&gt; <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-04-00000.nc</a>&gt;<br>

&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; The file opens perfectly in GrADS:<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; ga-&gt; xdfopen amazonia.ctl<br>
&gt; &gt;&gt;&gt; Scanning Descriptor File:  amazonia.ctl<br>
&gt; &gt;&gt;&gt; SDF file <a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-01-00000.nc</a>&gt; is open as file 1<br>

&gt; &gt;&gt;&gt; LON set to 279.062 311.875<br>
&gt; &gt;&gt;&gt; LAT set to -20.2999 6.45372<br>
&gt; &gt;&gt;&gt; LEV set to 0.00710064 0.00710064<br>
&gt; &gt;&gt;&gt; Time values set: 2000:1:1:0 2000:1:1:0<br>
&gt; &gt;&gt;&gt; E set to 1 1<br>
&gt; &gt;&gt;&gt; ga-&gt; q file<br>
&gt; &gt;&gt;&gt; File 1 : CLM History file information<br>
&gt; &gt;&gt;&gt;     Descriptor: amazonia.ctl<br>
&gt; &gt;&gt;&gt;     Binary: <a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-01-00000.nc</a>&gt;<br>

&gt; &gt;&gt;&gt;     Type = Gridded<br>
&gt; &gt;&gt;&gt;     Xsize = 106  Ysize = 115  Zsize = 15  Tsize = 4  Esize = 1<br>
&gt; &gt;&gt;&gt;     Number of Variables = 175<br>
&gt; &gt;&gt;&gt;       ...<br>
&gt; &gt;&gt;&gt;        tsoi  15  t,z,y,x  soil temperature<br>
&gt; &gt;&gt;&gt;        tsoi_10cm  0  t,y,x  soil temperature in top 10cm of soil<br>
&gt; &gt;&gt;&gt;        tv  0  t,y,x  vegetation temperature<br>
&gt; &gt;&gt;&gt;        u10  0  t,y,x  10-m wind<br>
&gt; &gt;&gt;&gt;        urban_ac  0  t,y,x  urban air conditioning flux<br>
&gt; &gt;&gt;&gt;    ...<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; But when I try to open in GrADS putting time as a variable (set 1 t 4)<br>
&gt; &gt;&gt;&gt; appears the following error message:<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; ga-&gt; set t 1 4<br>
&gt; &gt;&gt;&gt; Time values set: 2000:1:1:0 2000:4:1:0<br>
&gt; &gt;&gt;&gt; ga-&gt; d tsoi<br>
&gt; &gt;&gt;&gt; Contouring: 273 to 279 interval 0.5<br>
&gt; &gt;&gt;&gt; gancgrid error: nc_get_vara_double failed; NetCDF: Index exceeds<br>
&gt; &gt;&gt;&gt; dimension<br>
&gt; &gt;&gt;&gt; bound<br>
&gt; &gt;&gt;&gt; Data Request Error:  Error for variable &#39;tsoi&#39;<br>
&gt; &gt;&gt;&gt;     Error ocurred at column 1<br>
&gt; &gt;&gt;&gt; DISPLAY error:  Invalid expression<br>
&gt; &gt;&gt;&gt;     Expression = tsoi<br>
&gt; &gt;&gt;&gt; ga-&gt;<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; You know how to help me change the time?<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; All the best!<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; Celso<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt; On Wed, 10 Jul 2013 12:43:38 -0700, Matt Masarik<br>
&gt; &gt;&gt;&gt; &lt;<a href="mailto:mmasarik@atmos.ucla.edu">mmasarik@atmos.ucla.edu</a>&lt;mailto:<a href="mailto:mmasarik@atmos.ucla.edu">mmasarik@atmos.ucla.edu</a>&gt;&gt;<br>
&gt; &gt;&gt;&gt; wrote:<br>
&gt; &gt;&gt;&gt;&gt; Hi Celso,<br>
&gt; &gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt; you can solve this by using a small (3 line) control file<br>
&gt; &gt;&gt;&gt;&gt; that has the 365_day_calendar option specified.  Copy<br>
&gt; &gt;&gt;&gt;&gt; the following text into a file, say, amazonia.ctl<br>
&gt; &gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt; DSET ^<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-04-00000.nc</a>&gt;<br>

&gt; &gt;&gt;&gt;&gt; OPTIONS 365_day_calendar<br>
&gt; &gt;&gt;&gt;&gt; tdef time 3084 linear 16jan1850 1mo<br>
&gt; &gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt; I put your file name in the DSET line already.  You will<br>
&gt; &gt;&gt;&gt;&gt; need to change the 3084 to the total number of time<br>
&gt; &gt;&gt;&gt;&gt; records in the file.  You can find this out by issuing<br>
&gt; &gt;&gt;&gt;&gt; the following at a command prompt:<br>
&gt; &gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt; ncdump -c <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-04-00000.nc</a>&gt;<br>

&gt; &gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt; you will also need to change 16jan1850 to the start<br>
&gt; &gt;&gt;&gt;&gt; date of your data set, again use ncdump.<br>
&gt; &gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt; Once you make those changes you should be able to open<br>
&gt; &gt;&gt;&gt;&gt; the file by using<br>
&gt; &gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt; xdfopen amazonia.ctl<br>
&gt; &gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt; Matt<br>
&gt; &gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt; On 07/10/2013 11:50 AM, Celso Bandeira de Melo Ribeiro wrote:<br>
&gt; &gt;&gt;&gt;&gt;&gt; Hello,<br>
&gt; &gt;&gt;&gt;&gt;&gt; I&#39;m working with GrADS Version 2.0.a9 in Linux Ubuntu 64-bit and I am<br>
&gt; &gt;&gt;&gt; not<br>
&gt; &gt;&gt;&gt;&gt;&gt; able to open NetCDF file.<br>
&gt; &gt;&gt;&gt;&gt;&gt; I am using the function &quot;sdfopen&quot; but an error message appears as<br>
&gt; &gt;&gt;&gt; sending<br>
&gt; &gt;&gt;&gt;&gt;&gt; the attachment.<br>
&gt; &gt;&gt;&gt;&gt;&gt; The message, as follows:<br>
&gt; &gt;&gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt;&gt; ga-&gt; sdfopen <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-04-00000.nc</a>&gt;<br>

&gt; &gt;&gt;&gt;&gt;&gt; Scanning self-describing file:  <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a>&lt;<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://amazonia.clm2.h0.2002-01-04-00000.nc</a>&gt;<br>

&gt; &gt;&gt;&gt;&gt;&gt; SDF Error: 365 day calendars are no longer supported by sdfopen.<br>
&gt; &gt;&gt;&gt;&gt;&gt;      To open this file with GrADS, use a descriptor file with<br>
&gt; &gt;&gt;&gt;&gt;&gt;      a complete TDEF entry and OPTIONS 365_day_calendar.<br>
&gt; &gt;&gt;&gt;&gt;&gt;      Documentation is at<br>
&gt; &gt;&gt;&gt; <a href="http://iges.org/grads/gadoc/SDFdescriptorfile.html" target="_blank">http://iges.org/grads/gadoc/SDFdescriptorfile.html</a><br>
&gt; &gt;&gt;&gt;&gt;&gt; ga-&gt;<br>
&gt; &gt;&gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt;&gt; I would like to know how solve this problem.<br>
&gt; &gt;&gt;&gt;&gt;&gt; Thank you!<br>
&gt; &gt;&gt;&gt;&gt;&gt; Celso<br>
&gt; &gt;&gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt; _______________________________________________<br>
&gt; &gt;&gt;&gt;&gt; gradsusr mailing list<br>
&gt; &gt;&gt;&gt;&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
&gt; &gt;&gt;&gt;&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt; &gt;&gt; _______________________________________________<br>
&gt; &gt;&gt; gradsusr mailing list<br>
&gt; &gt;&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
&gt; &gt;&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; ------------------------------<br>
&gt;<br>
&gt; Message: 5<br>
&gt; Date: Thu, 11 Jul 2013 07:40:30 -1000<br>
&gt; From: &quot;James T. Potemra&quot; &lt;<a href="mailto:jimp@hawaii.edu">jimp@hawaii.edu</a>&lt;mailto:<a href="mailto:jimp@hawaii.edu">jimp@hawaii.edu</a>&gt;&gt;<br>
&gt; Subject: Re: [gradsusr] Shift in the lattitude while plotting in grads<br>
&gt;         mercator/latlon projction<br>
&gt; To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;&gt;<br>
&gt; Message-ID: &lt;<a href="mailto:51DEEE0E.8030905@hawaii.edu">51DEEE0E.8030905@hawaii.edu</a>&lt;mailto:<a href="mailto:51DEEE0E.8030905@hawaii.edu">51DEEE0E.8030905@hawaii.edu</a>&gt;&gt;<br>
&gt; Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
&gt;<br>
&gt; Itesh:<br>
&gt;<br>
&gt; Your model output is 640 x 460, but your descriptor file has<br>
&gt; xdef = 639 and ydef = 459.  These should be 640 and 460,<br>
&gt; respectively.  Also, I think what Mohsen meant was to make<br>
&gt; sure your starting lat/lon is the same in the descriptor file<br>
&gt; and model grid, i.e., you have 61.1815 as the starting lon<br>
&gt; and -3.1006 as the starting lat.  Are these the same as the<br>
&gt; model grid?<br>
&gt;<br>
&gt; Jim<br>
&gt;<br>
&gt; On 7/11/13 7:31 AM, itesh dash wrote:<br>
&gt; &gt; Dear Mohsen,<br>
&gt; &gt;<br>
&gt; &gt; thanks for the response. But I am not even trying to use any script.<br>
&gt; &gt; Just by opening and plotting the ctl file.<br>
&gt; &gt;<br>
&gt; &gt; Here is my ctl file<br>
&gt; &gt;<br>
&gt; &gt; dset ^Jul1013_12.dat<br>
&gt; &gt; undef 1.e35<br>
&gt; &gt; xdef  639 linear   61.1815   0.0935<br>
&gt; &gt; ydef  459 linear   -3.1006   0.0873<br>
&gt; &gt; zdef    3 levels<br>
&gt; &gt;   850.00000<br>
&gt; &gt;   500.00000<br>
&gt; &gt;   200.00000<br>
&gt; &gt; tdef         8 linear 12z10jul2013  6hr<br>
&gt; &gt; vars 1<br>
&gt; &gt; HGT 0 0 Terrain Height<br>
&gt; &gt; endvars<br>
&gt; &gt;<br>
&gt; &gt; The wrf model ran with mercator projection 640 x 460 grid points 9 km<br>
&gt; &gt; resolution.<br>
&gt; &gt;<br>
&gt; &gt; I could not figure out why the shift is there. Is there a way i can<br>
&gt; &gt; balance the shift.<br>
&gt; &gt;<br>
&gt; &gt; regards,<br>
&gt; &gt; itesh<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; On Thu, Jul 11, 2013 at 11:58 PM, Mohsen Soltani<br>
&gt; &gt; &lt;<a href="mailto:soltani.clima@gmail.com">soltani.clima@gmail.com</a>&lt;mailto:<a href="mailto:soltani.clima@gmail.com">soltani.clima@gmail.com</a>&gt; &lt;mailto:<a href="mailto:soltani.clima@gmail.com">soltani.clima@gmail.com</a>&lt;mailto:<a href="mailto:soltani.clima@gmail.com">soltani.clima@gmail.com</a>&gt;&gt;&gt; wrote:<br>

&gt; &gt;<br>
&gt; &gt;     itesh,<br>
&gt; &gt;<br>
&gt; &gt;     Make sure Lat/Lon in your .ctl file are the same as your script!<br>
&gt; &gt;<br>
&gt; &gt;     Mohsen<br>
&gt; &gt;<br>
&gt; &gt;     --<br>
&gt; &gt;     April showers bring May flowers!<br>
&gt; &gt;     --<br>
&gt; &gt;     *Mohsen Soltani, M.Sc.*<br>
&gt; &gt;     Climatology Grad Student,<br>
&gt; &gt;     Department of Climatology, Faculty of Geography,<br>
&gt; &gt;     University of Tehran, Tehran, Iran<br>
&gt; &gt;     voice: <a href="tel:%28%2B98%29%209119772934" value="+989119772934">(+98) 9119772934</a>&lt;tel:%28%2B98%29%209119772934&gt;<br>
&gt; &gt;     e-mail: <a href="mailto:soltani.clima@gmail.com">soltani.clima@gmail.com</a>&lt;mailto:<a href="mailto:soltani.clima@gmail.com">soltani.clima@gmail.com</a>&gt; &lt;mailto:<a href="mailto:soltani.clima@gmail.com">soltani.clima@gmail.com</a>&lt;mailto:<a href="mailto:soltani.clima@gmail.com">soltani.clima@gmail.com</a>&gt;&gt;<br>

&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;     On Thu, Jul 11, 2013 at 9:05 AM, itesh dash &lt;<a href="mailto:itesh@rimes.int">itesh@rimes.int</a>&lt;mailto:<a href="mailto:itesh@rimes.int">itesh@rimes.int</a>&gt;<br>
&gt; &gt;     &lt;mailto:<a href="mailto:itesh@rimes.int">itesh@rimes.int</a>&lt;mailto:<a href="mailto:itesh@rimes.int">itesh@rimes.int</a>&gt;&gt;&gt; wrote:<br>
&gt; &gt;<br>
&gt; &gt;         Dear Sir,<br>
&gt; &gt;<br>
&gt; &gt;         I am working with the WRF and tried to plot the output in<br>
&gt; &gt;         grads. I found there is nearly 1 degree shift in the<br>
&gt; &gt;         projection and the plot.<br>
&gt; &gt;<br>
&gt; &gt;         Please guide how to correct this.<br>
&gt; &gt;<br>
&gt; &gt;         I am using mercator projection for the WRF model. Please find<br>
&gt; &gt;         attached image.<br>
&gt; &gt;<br>
&gt; &gt;         Seeking your advice on the same.<br>
&gt; &gt;<br>
&gt; &gt;         regards,<br>
&gt; &gt;         itesh<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;         _______________________________________________<br>
&gt; &gt;         gradsusr mailing list<br>
&gt; &gt;         <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt; &lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;&gt;<br>

&gt; &gt;         <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;     _______________________________________________<br>
&gt; &gt;     gradsusr mailing list<br>
&gt; &gt;     <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt; &lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;&gt;<br>

&gt; &gt;     <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; _______________________________________________<br>
&gt; &gt; gradsusr mailing list<br>
&gt; &gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
&gt; &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt;<br>
&gt; -------------- next part --------------<br>
&gt; An HTML attachment was scrubbed...<br>
&gt; URL:<br>
&gt; <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20130711/44403521/attachment.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20130711/44403521/attachment.html</a><br>
&gt;<br>
&gt; ------------------------------<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt;<br>
&gt;<br>
&gt; End of gradsusr Digest, Vol 41, Issue 19<br>
&gt; ****************************************<br>
&gt;<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20130711/721238f0/attachment.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20130711/721238f0/attachment.html</a><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: NINO34GISSR1G4.jpeg<br>
Type: image/jpeg<br>
Size: 15210 bytes<br>
Desc: not available<br>
Url : <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20130711/721238f0/attachment.jpeg" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20130711/721238f0/attachment.jpeg</a><br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br>
<br>
End of gradsusr Digest, Vol 41, Issue 20<br>
****************************************<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20130711/0c843818/attachment.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20130711/0c843818/attachment.html</a><br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br>
<br>
End of gradsusr Digest, Vol 41, Issue 23<br>
****************************************<br>
</blockquote></div><br></div></div></div></div></div></div></div>