<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">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>
<div><br>
<div style="font-family:Tahoma; font-size:13px">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>
</div>
</div>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div style="direction: ltr;" id="divRpF85"><font color="#000000" face="Tahoma" size="2"><b>From:</b> gradsusr-bounces@gradsusr.org [gradsusr-bounces@gradsusr.org] on behalf of Corey Gabriel [cjgabriel7@gmail.com]<br>
<b>Sent:</b> Thursday, July 11, 2013 2:51 PM<br>
<b>To:</b> gradsusr<br>
<b>Subject:</b> [gradsusr] Least-Squares Regression and Climatology (Time Series)<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">***Identical to previous message. Subject line included this time***<br>
<div><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: *'et time jan2020 dec2020<br>
define nino34=(aave((sst-273.15),lon=
<div id=":3y">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'<br>
<br>
*Code2(figure attached)<br>
*'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'<br>
*</div>
</div>
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">On Thu, Jul 11, 2013 at 2:22 PM, <span dir="ltr"><<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex; border-left:1px #ccc solid; padding-left:1ex">
Send gradsusr mailing list submissions to<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">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 'help' to<br>
<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:gradsusr-owner@gradsusr.org" target="_blank">gradsusr-owner@gradsusr.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of gradsusr digest..."<br>
<br>
<br>
Today'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 <<a href="mailto:cjgabriel7@gmail.com" target="_blank">cjgabriel7@gmail.com</a>><br>
Subject: Re: [gradsusr] gradsusr Digest, Vol 41, Issue 19<br>
To: <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
Message-ID:<br>
<CAEhnGEru0gWKX=<a href="mailto:tSMgLyWCm3iFgMthwhxkEdnnY6yTTHANZA1A@mail.gmail.com" target="_blank">tSMgLyWCm3iFgMthwhxkEdnnY6yTTHANZA1A@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<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: *'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'<br>
<br>
*Code2(figure attached)<br>
*'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'<br>
*<br>
*<br>
<br>
<br>
<br>
<br>
On Thu, Jul 11, 2013 at 1:40 PM, <<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a>> wrote:<br>
<br>
> Send gradsusr mailing list submissions to<br>
> <a href="mailto:gradsusr@gradsusr.org" target="_blank">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 'help' to<br>
> <a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a><br>
><br>
> You can reach the person managing the list at<br>
> <a href="mailto:gradsusr-owner@gradsusr.org" target="_blank">gradsusr-owner@gradsusr.org</a><br>
><br>
> When replying, please edit your Subject line so it is more specific<br>
> than "Re: Contents of gradsusr digest..."<br>
><br>
><br>
> Today's Topics:<br>
><br>
> 1. Re: help - NetCDF file (Celso Bandeira de Melo Ribeiro)<br>
> 2. Re: Shift in the lattitude while plotting in grads<br>
> mercator/latlon projction (Mohsen Soltani)<br>
> 3. Re: Shift in the lattitude while plotting in grads<br>
> mercator/latlon projction (itesh dash)<br>
> 4. Re: help - NetCDF file (James T. Potemra)<br>
> 5. Re: Shift in the lattitude while plotting in grads<br>
> mercator/latlon projction (James T. Potemra)<br>
><br>
><br>
> ----------------------------------------------------------------------<br>
><br>
> Message: 1<br>
> Date: Thu, 11 Jul 2013 13:57:21 -0300<br>
> From: Celso Bandeira de Melo Ribeiro <<a href="mailto:celso.bandeira@ufjf.edu.br" target="_blank">celso.bandeira@ufjf.edu.br</a>><br>
> Subject: Re: [gradsusr] help - NetCDF file<br>
> To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>
> Message-ID: <34fb6a6e970b1fff4dc36f87278e432e@localhost><br>
> Content-Type: text/plain; charset=UTF-8<br>
><br>
> Hi James,<br>
> Thank you for attention and help.<br>
> I modified the file data to diaries as show below:<br>
><br>
> DSET ^amazonia.clm2.h0.2002-01-%<a href="http://m2-00000.nc" target="_blank">m2-00000.nc</a><br>
> OPTIONS 365_day_calendar<br>
> tdef time 4 linear 01jan2000 1dy<br>
><br>
> But it did not work, appeared the following error message:<br>
><br>
> ga-> xdfopen amazonia.ctl<br>
> Scanning Descriptor File: amazonia.ctl<br>
> Error: nc_open failed to open file amazonia.clm2.h0.2002-01-%<a href="http://m2-00000.nc" target="_blank">m2-00000.nc</a><br>
> No such file or directory<br>
> gadsdf: Couldn't ingest SDF metadata.<br>
> SDF Descriptor file amazonia.ctl was not successfully opened & parsed.<br>
> ga-><br>
><br>
> All the best!<br>
><br>
> Celso<br>
><br>
><br>
><br>
> On Wed, 10 Jul 2013 17:06:21 -1000, "James T. Potemra" <<a href="mailto:jimp@hawaii.edu" target="_blank">jimp@hawaii.edu</a>><br>
> wrote:<br>
> > Celso:<br>
> ><br>
> > Your descriptor file lists four times (months), but it is specifically<br>
> > 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>
> > likely want to read all four, so you'll have to add an option<br>
> > "template". Also, it's not clear from the file names if they reference<br>
> > 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>
> > for January 4, 2002 or April 1, 2002? In your descriptor file you<br>
> > specify the time increment as monthly (1mo). If the files have single<br>
> > time values (one monthly value), and the format of the file name is<br>
> > day-month, your DSET line should be something like:<br>
> ><br>
> > DSET ^amazonia.clm2.h0.2002-01-%<a href="http://m2-00000.nc" target="_blank">m2-00000.nc</a><br>
> ><br>
> > The %m2 indicating a substitution template of 2-digit month.<br>
> ><br>
> > Jim<br>
> ><br>
> > On 7/10/13 3:58 PM, Celso Bandeira de Melo Ribeiro wrote:<br>
> >> Hi Matt,<br>
> >><br>
> >> Thank you very much for your attention and help.<br>
> >> Worked perfectly I created the following file:<br>
> >><br>
> >> DSET ^<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a><br>
> >> OPTIONS 365_day_calendar<br>
> >> tdef time 4 linear 01jan2000 1mo<br>
> >><br>
> >> I am working with 4 files (Daily), whose names are:<br>
> >><br>
> >> <a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a><br>
> >> <a href="http://amazonia.clm2.h0.2002-01-02-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-02-00000.nc</a><br>
> >> <a href="http://amazonia.clm2.h0.2002-01-03-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-03-00000.nc</a><br>
> >> <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a><br>
> >><br>
> >> The file opens perfectly in GrADS:<br>
> >><br>
> >> ga-> xdfopen amazonia.ctl<br>
> >> Scanning Descriptor File: amazonia.ctl<br>
> >> 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>
> >> 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: 2000:1:1:0 2000:1:1:0<br>
> >> E set to 1 1<br>
> >> ga-> q file<br>
> >> File 1 : CLM History file information<br>
> >> Descriptor: amazonia.ctl<br>
> >> Binary: <a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">
amazonia.clm2.h0.2002-01-01-00000.nc</a><br>
> >> Type = Gridded<br>
> >> Xsize = 106 Ysize = 115 Zsize = 15 Tsize = 4 Esize = 1<br>
> >> Number of Variables = 175<br>
> >> ...<br>
> >> tsoi 15 t,z,y,x soil temperature<br>
> >> tsoi_10cm 0 t,y,x soil temperature in top 10cm of soil<br>
> >> tv 0 t,y,x vegetation temperature<br>
> >> u10 0 t,y,x 10-m wind<br>
> >> urban_ac 0 t,y,x urban air conditioning flux<br>
> >> ...<br>
> >><br>
> >> But when I try to open in GrADS putting time as a variable (set 1 t 4)<br>
> >> appears the following error message:<br>
> >><br>
> >> ga-> set t 1 4<br>
> >> Time values set: 2000:1:1:0 2000:4:1:0<br>
> >> ga-> d tsoi<br>
> >> Contouring: 273 to 279 interval 0.5<br>
> >> gancgrid error: nc_get_vara_double failed; NetCDF: Index exceeds<br>
> >> dimension<br>
> >> bound<br>
> >> Data Request Error: Error for variable 'tsoi'<br>
> >> Error ocurred at column 1<br>
> >> DISPLAY error: Invalid expression<br>
> >> Expression = tsoi<br>
> >> ga-><br>
> >><br>
> >> You know how to help me change the time?<br>
> >><br>
> >> All the best!<br>
> >><br>
> >> Celso<br>
> >><br>
> >><br>
> >><br>
> >> On Wed, 10 Jul 2013 12:43:38 -0700, Matt Masarik<br>
> >> <<a href="mailto:mmasarik@atmos.ucla.edu" target="_blank">mmasarik@atmos.ucla.edu</a>><br>
> >> wrote:<br>
> >>> Hi Celso,<br>
> >>><br>
> >>> you can solve this by using a small (3 line) control file<br>
> >>> that has the 365_day_calendar option specified. Copy<br>
> >>> the following text into a file, say, amazonia.ctl<br>
> >>><br>
> >>> DSET ^<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a><br>
> >>> OPTIONS 365_day_calendar<br>
> >>> tdef time 3084 linear 16jan1850 1mo<br>
> >>><br>
> >>> I put your file name in the DSET line already. You will<br>
> >>> need to change the 3084 to the total number of time<br>
> >>> records in the file. You can find this out by issuing<br>
> >>> the following at a command prompt:<br>
> >>><br>
> >>> 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>
> >>><br>
> >>> you will also need to change 16jan1850 to the start<br>
> >>> date of your data set, again use ncdump.<br>
> >>><br>
> >>> Once you make those changes you should be able to open<br>
> >>> the file by using<br>
> >>><br>
> >>> xdfopen amazonia.ctl<br>
> >>><br>
> >>><br>
> >>> Matt<br>
> >>><br>
> >>><br>
> >>> On 07/10/2013 11:50 AM, Celso Bandeira de Melo Ribeiro wrote:<br>
> >>>> Hello,<br>
> >>>> I'm working with GrADS Version 2.0.a9 in Linux Ubuntu 64-bit and I am<br>
> >> not<br>
> >>>> able to open NetCDF file.<br>
> >>>> I am using the function "sdfopen" but an error message appears as<br>
> >> sending<br>
> >>>> the attachment.<br>
> >>>> The message, as follows:<br>
> >>>><br>
> >>>> ga-> sdfopen <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">
amazonia.clm2.h0.2002-01-04-00000.nc</a><br>
> >>>> 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><br>
> >>>> SDF Error: 365 day calendars are no longer supported by sdfopen.<br>
> >>>> To open this file with GrADS, use a descriptor file with<br>
> >>>> a complete TDEF entry and OPTIONS 365_day_calendar.<br>
> >>>> Documentation is at<br>
> >> <a href="http://iges.org/grads/gadoc/SDFdescriptorfile.html" target="_blank">
http://iges.org/grads/gadoc/SDFdescriptorfile.html</a><br>
> >>>> ga-><br>
> >>>><br>
> >>>> I would like to know how solve this problem.<br>
> >>>> Thank you!<br>
> >>>> Celso<br>
> >>>><br>
> >>> _______________________________________________<br>
> >>> gradsusr mailing list<br>
> >>> <a href="mailto:gradsusr@gradsusr.org" target="_blank">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>
> > gradsusr mailing list<br>
> > <a href="mailto:gradsusr@gradsusr.org" target="_blank">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>
> ____________________________________________<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 09:58:42 -0700<br>
> From: Mohsen Soltani <<a href="mailto:soltani.clima@gmail.com" target="_blank">soltani.clima@gmail.com</a>><br>
> Subject: Re: [gradsusr] Shift in the lattitude while plotting in grads<br>
> mercator/latlon projction<br>
> To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>
> Message-ID:<br>
> <<br>
> <a href="mailto:CAPoCDiM50XFQjKOHUN86kCNNdB3kp69YFaJYb5emWgMyYLpCOA@mail.gmail.com" target="_blank">
CAPoCDiM50XFQjKOHUN86kCNNdB3kp69YFaJYb5emWgMyYLpCOA@mail.gmail.com</a>><br>
> Content-Type: text/plain; charset="iso-8859-1"<br>
><br>
> itesh,<br>
><br>
> Make sure Lat/Lon in your .ctl file are the same as your script!<br>
><br>
> Mohsen<br>
><br>
> --<br>
> April showers bring May flowers!<br>
> --<br>
> *Mohsen Soltani, M.Sc.*<br>
> Climatology Grad Student,<br>
> Department of Climatology, Faculty of Geography,<br>
> University of Tehran, Tehran, Iran<br>
> voice: <a href="tel:%28%2B98%29%209119772934" value="+989119772934" target="_blank">
(+98) 9119772934</a><br>
> e-mail: <a href="mailto:soltani.clima@gmail.com" target="_blank">soltani.clima@gmail.com</a><br>
><br>
><br>
> On Thu, Jul 11, 2013 at 9:05 AM, itesh dash <<a href="mailto:itesh@rimes.int" target="_blank">itesh@rimes.int</a>> wrote:<br>
><br>
> > Dear Sir,<br>
> ><br>
> > I am working with the WRF and tried to plot the output in grads. I found<br>
> > there is nearly 1 degree shift in the projection and the plot.<br>
> ><br>
> > Please guide how to correct this.<br>
> ><br>
> > I am using mercator projection for the WRF model. Please find attached<br>
> > image.<br>
> ><br>
> > Seeking your advice on the same.<br>
> ><br>
> > regards,<br>
> > itesh<br>
> ><br>
> ><br>
> > _______________________________________________<br>
> > gradsusr mailing list<br>
> > <a href="mailto:gradsusr@gradsusr.org" target="_blank">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>
> -------------- next part --------------<br>
> An HTML attachment was scrubbed...<br>
> URL:<br>
> <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>
><br>
> ------------------------------<br>
><br>
> Message: 3<br>
> Date: Fri, 12 Jul 2013 00:31:29 +0700<br>
> From: itesh dash <<a href="mailto:itesh@rimes.int" target="_blank">itesh@rimes.int</a>><br>
> Subject: Re: [gradsusr] Shift in the lattitude while plotting in grads<br>
> mercator/latlon projction<br>
> To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>
> Message-ID:<br>
> <CAFTUJN465mGcZ8k=<br>
> <a href="mailto:JDB1t6UvJRG9jO%2B1L%2BfAD20Gd4G8X_Qayg@mail.gmail.com" target="_blank">
JDB1t6UvJRG9jO+1L+fAD20Gd4G8X_Qayg@mail.gmail.com</a>><br>
> Content-Type: text/plain; charset="iso-8859-1"<br>
><br>
> Dear Mohsen,<br>
><br>
> thanks for the response. But I am not even trying to use any script. Just<br>
> by opening and plotting the ctl file.<br>
><br>
> Here is my ctl file<br>
><br>
> dset ^Jul1013_12.dat<br>
> undef 1.e35<br>
> xdef 639 linear 61.1815 0.0935<br>
> ydef 459 linear -3.1006 0.0873<br>
> zdef 3 levels<br>
> 850.00000<br>
> 500.00000<br>
> 200.00000<br>
> tdef 8 linear 12z10jul2013 6hr<br>
> vars 1<br>
> HGT 0 0 Terrain Height<br>
> endvars<br>
><br>
> The wrf model ran with mercator projection 640 x 460 grid points 9 km<br>
> resolution.<br>
><br>
> I could not figure out why the shift is there. Is there a way i can balance<br>
> the shift.<br>
><br>
> regards,<br>
> itesh<br>
><br>
><br>
> On Thu, Jul 11, 2013 at 11:58 PM, Mohsen Soltani <<a href="mailto:soltani.clima@gmail.com" target="_blank">soltani.clima@gmail.com</a><br>
> >wrote:<br>
><br>
> > itesh,<br>
> ><br>
> > Make sure Lat/Lon in your .ctl file are the same as your script!<br>
> ><br>
> > Mohsen<br>
> ><br>
> > --<br>
> > April showers bring May flowers!<br>
> > --<br>
> > *Mohsen Soltani, M.Sc.*<br>
> > Climatology Grad Student,<br>
> > Department of Climatology, Faculty of Geography,<br>
> > University of Tehran, Tehran, Iran<br>
> > voice: <a href="tel:%28%2B98%29%209119772934" value="+989119772934" target="_blank">
(+98) 9119772934</a><br>
> > e-mail: <a href="mailto:soltani.clima@gmail.com" target="_blank">soltani.clima@gmail.com</a><br>
> ><br>
> ><br>
> > On Thu, Jul 11, 2013 at 9:05 AM, itesh dash <<a href="mailto:itesh@rimes.int" target="_blank">itesh@rimes.int</a>> wrote:<br>
> ><br>
> >> Dear Sir,<br>
> >><br>
> >> I am working with the WRF and tried to plot the output in grads. I found<br>
> >> there is nearly 1 degree shift in the projection and the plot.<br>
> >><br>
> >> Please guide how to correct this.<br>
> >><br>
> >> I am using mercator projection for the WRF model. Please find attached<br>
> >> image.<br>
> >><br>
> >> Seeking your advice on the same.<br>
> >><br>
> >> regards,<br>
> >> itesh<br>
> >><br>
> >><br>
> >> _______________________________________________<br>
> >> gradsusr mailing list<br>
> >> <a href="mailto:gradsusr@gradsusr.org" target="_blank">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>
> ><br>
> > _______________________________________________<br>
> > gradsusr mailing list<br>
> > <a href="mailto:gradsusr@gradsusr.org" target="_blank">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>
> -------------- next part --------------<br>
> An HTML attachment was scrubbed...<br>
> URL:<br>
> <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>
><br>
> ------------------------------<br>
><br>
> Message: 4<br>
> Date: Thu, 11 Jul 2013 07:33:38 -1000<br>
> From: "James T. Potemra" <<a href="mailto:jimp@hawaii.edu" target="_blank">jimp@hawaii.edu</a>><br>
> Subject: Re: [gradsusr] help - NetCDF file<br>
> To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>
> Message-ID: <<a href="mailto:51DEEC72.8050809@hawaii.edu" target="_blank">51DEEC72.8050809@hawaii.edu</a>><br>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>
><br>
> Hi Celso:<br>
><br>
> You also have to add the option "template":<br>
> OPTIONS template<br>
><br>
> Jim<br>
> On 7/11/13 6:57 AM, Celso Bandeira de Melo Ribeiro wrote:<br>
> > Hi James,<br>
> > Thank you for attention and help.<br>
> > I modified the file data to diaries as show below:<br>
> ><br>
> > DSET ^amazonia.clm2.h0.2002-01-%<a href="http://m2-00000.nc" target="_blank">m2-00000.nc</a><br>
> > OPTIONS 365_day_calendar<br>
> > tdef time 4 linear 01jan2000 1dy<br>
> ><br>
> > But it did not work, appeared the following error message:<br>
> ><br>
> > ga-> xdfopen amazonia.ctl<br>
> > Scanning Descriptor File: amazonia.ctl<br>
> > Error: nc_open failed to open file amazonia.clm2.h0.2002-01-%<a href="http://m2-00000.nc" target="_blank">m2-00000.nc</a><br>
> > No such file or directory<br>
> > gadsdf: Couldn't ingest SDF metadata.<br>
> > SDF Descriptor file amazonia.ctl was not successfully opened & parsed.<br>
> > ga-><br>
> ><br>
> > All the best!<br>
> ><br>
> > Celso<br>
> ><br>
> ><br>
> ><br>
> > On Wed, 10 Jul 2013 17:06:21 -1000, "James T. Potemra" <<a href="mailto:jimp@hawaii.edu" target="_blank">jimp@hawaii.edu</a>><br>
> > wrote:<br>
> >> Celso:<br>
> >><br>
> >> Your descriptor file lists four times (months), but it is specifically<br>
> >> 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>
> >> likely want to read all four, so you'll have to add an option<br>
> >> "template". Also, it's not clear from the file names if they reference<br>
> >> 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>
> >> for January 4, 2002 or April 1, 2002? In your descriptor file you<br>
> >> specify the time increment as monthly (1mo). If the files have single<br>
> >> time values (one monthly value), and the format of the file name is<br>
> >> day-month, your DSET line should be something like:<br>
> >><br>
> >> DSET ^amazonia.clm2.h0.2002-01-%<a href="http://m2-00000.nc" target="_blank">m2-00000.nc</a><br>
> >><br>
> >> The %m2 indicating a substitution template of 2-digit month.<br>
> >><br>
> >> Jim<br>
> >><br>
> >> On 7/10/13 3:58 PM, Celso Bandeira de Melo Ribeiro wrote:<br>
> >>> Hi Matt,<br>
> >>><br>
> >>> Thank you very much for your attention and help.<br>
> >>> Worked perfectly I created the following file:<br>
> >>><br>
> >>> DSET ^<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a><br>
> >>> OPTIONS 365_day_calendar<br>
> >>> tdef time 4 linear 01jan2000 1mo<br>
> >>><br>
> >>> I am working with 4 files (Daily), whose names are:<br>
> >>><br>
> >>> <a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-01-00000.nc</a><br>
> >>> <a href="http://amazonia.clm2.h0.2002-01-02-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-02-00000.nc</a><br>
> >>> <a href="http://amazonia.clm2.h0.2002-01-03-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-03-00000.nc</a><br>
> >>> <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a><br>
> >>><br>
> >>> The file opens perfectly in GrADS:<br>
> >>><br>
> >>> ga-> xdfopen amazonia.ctl<br>
> >>> Scanning Descriptor File: amazonia.ctl<br>
> >>> 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>
> >>> 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: 2000:1:1:0 2000:1:1:0<br>
> >>> E set to 1 1<br>
> >>> ga-> q file<br>
> >>> File 1 : CLM History file information<br>
> >>> Descriptor: amazonia.ctl<br>
> >>> Binary: <a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">
amazonia.clm2.h0.2002-01-01-00000.nc</a><br>
> >>> Type = Gridded<br>
> >>> Xsize = 106 Ysize = 115 Zsize = 15 Tsize = 4 Esize = 1<br>
> >>> Number of Variables = 175<br>
> >>> ...<br>
> >>> tsoi 15 t,z,y,x soil temperature<br>
> >>> tsoi_10cm 0 t,y,x soil temperature in top 10cm of soil<br>
> >>> tv 0 t,y,x vegetation temperature<br>
> >>> u10 0 t,y,x 10-m wind<br>
> >>> urban_ac 0 t,y,x urban air conditioning flux<br>
> >>> ...<br>
> >>><br>
> >>> But when I try to open in GrADS putting time as a variable (set 1 t 4)<br>
> >>> appears the following error message:<br>
> >>><br>
> >>> ga-> set t 1 4<br>
> >>> Time values set: 2000:1:1:0 2000:4:1:0<br>
> >>> ga-> d tsoi<br>
> >>> Contouring: 273 to 279 interval 0.5<br>
> >>> gancgrid error: nc_get_vara_double failed; NetCDF: Index exceeds<br>
> >>> dimension<br>
> >>> bound<br>
> >>> Data Request Error: Error for variable 'tsoi'<br>
> >>> Error ocurred at column 1<br>
> >>> DISPLAY error: Invalid expression<br>
> >>> Expression = tsoi<br>
> >>> ga-><br>
> >>><br>
> >>> You know how to help me change the time?<br>
> >>><br>
> >>> All the best!<br>
> >>><br>
> >>> Celso<br>
> >>><br>
> >>><br>
> >>><br>
> >>> On Wed, 10 Jul 2013 12:43:38 -0700, Matt Masarik<br>
> >>> <<a href="mailto:mmasarik@atmos.ucla.edu" target="_blank">mmasarik@atmos.ucla.edu</a>><br>
> >>> wrote:<br>
> >>>> Hi Celso,<br>
> >>>><br>
> >>>> you can solve this by using a small (3 line) control file<br>
> >>>> that has the 365_day_calendar option specified. Copy<br>
> >>>> the following text into a file, say, amazonia.ctl<br>
> >>>><br>
> >>>> DSET ^<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">amazonia.clm2.h0.2002-01-04-00000.nc</a><br>
> >>>> OPTIONS 365_day_calendar<br>
> >>>> tdef time 3084 linear 16jan1850 1mo<br>
> >>>><br>
> >>>> I put your file name in the DSET line already. You will<br>
> >>>> need to change the 3084 to the total number of time<br>
> >>>> records in the file. You can find this out by issuing<br>
> >>>> the following at a command prompt:<br>
> >>>><br>
> >>>> 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>
> >>>><br>
> >>>> you will also need to change 16jan1850 to the start<br>
> >>>> date of your data set, again use ncdump.<br>
> >>>><br>
> >>>> Once you make those changes you should be able to open<br>
> >>>> the file by using<br>
> >>>><br>
> >>>> xdfopen amazonia.ctl<br>
> >>>><br>
> >>>><br>
> >>>> Matt<br>
> >>>><br>
> >>>><br>
> >>>> On 07/10/2013 11:50 AM, Celso Bandeira de Melo Ribeiro wrote:<br>
> >>>>> Hello,<br>
> >>>>> I'm working with GrADS Version 2.0.a9 in Linux Ubuntu 64-bit and I am<br>
> >>> not<br>
> >>>>> able to open NetCDF file.<br>
> >>>>> I am using the function "sdfopen" but an error message appears as<br>
> >>> sending<br>
> >>>>> the attachment.<br>
> >>>>> The message, as follows:<br>
> >>>>><br>
> >>>>> ga-> sdfopen <a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">
amazonia.clm2.h0.2002-01-04-00000.nc</a><br>
> >>>>> 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><br>
> >>>>> SDF Error: 365 day calendars are no longer supported by sdfopen.<br>
> >>>>> To open this file with GrADS, use a descriptor file with<br>
> >>>>> a complete TDEF entry and OPTIONS 365_day_calendar.<br>
> >>>>> Documentation is at<br>
> >>> <a href="http://iges.org/grads/gadoc/SDFdescriptorfile.html" target="_blank">
http://iges.org/grads/gadoc/SDFdescriptorfile.html</a><br>
> >>>>> ga-><br>
> >>>>><br>
> >>>>> I would like to know how solve this problem.<br>
> >>>>> Thank you!<br>
> >>>>> Celso<br>
> >>>>><br>
> >>>> _______________________________________________<br>
> >>>> gradsusr mailing list<br>
> >>>> <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
> >>>> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
> >> _______________________________________________<br>
> >> gradsusr mailing list<br>
> >> <a href="mailto:gradsusr@gradsusr.org" target="_blank">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>
><br>
> ------------------------------<br>
><br>
> Message: 5<br>
> Date: Thu, 11 Jul 2013 07:40:30 -1000<br>
> From: "James T. Potemra" <<a href="mailto:jimp@hawaii.edu" target="_blank">jimp@hawaii.edu</a>><br>
> Subject: Re: [gradsusr] Shift in the lattitude while plotting in grads<br>
> mercator/latlon projction<br>
> To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>
> Message-ID: <<a href="mailto:51DEEE0E.8030905@hawaii.edu" target="_blank">51DEEE0E.8030905@hawaii.edu</a>><br>
> Content-Type: text/plain; charset="iso-8859-1"<br>
><br>
> Itesh:<br>
><br>
> Your model output is 640 x 460, but your descriptor file has<br>
> xdef = 639 and ydef = 459. These should be 640 and 460,<br>
> respectively. Also, I think what Mohsen meant was to make<br>
> sure your starting lat/lon is the same in the descriptor file<br>
> and model grid, i.e., you have 61.1815 as the starting lon<br>
> and -3.1006 as the starting lat. Are these the same as the<br>
> model grid?<br>
><br>
> Jim<br>
><br>
> On 7/11/13 7:31 AM, itesh dash wrote:<br>
> > Dear Mohsen,<br>
> ><br>
> > thanks for the response. But I am not even trying to use any script.<br>
> > Just by opening and plotting the ctl file.<br>
> ><br>
> > Here is my ctl file<br>
> ><br>
> > dset ^Jul1013_12.dat<br>
> > undef 1.e35<br>
> > xdef 639 linear 61.1815 0.0935<br>
> > ydef 459 linear -3.1006 0.0873<br>
> > zdef 3 levels<br>
> > 850.00000<br>
> > 500.00000<br>
> > 200.00000<br>
> > tdef 8 linear 12z10jul2013 6hr<br>
> > vars 1<br>
> > HGT 0 0 Terrain Height<br>
> > endvars<br>
> ><br>
> > The wrf model ran with mercator projection 640 x 460 grid points 9 km<br>
> > resolution.<br>
> ><br>
> > I could not figure out why the shift is there. Is there a way i can<br>
> > balance the shift.<br>
> ><br>
> > regards,<br>
> > itesh<br>
> ><br>
> ><br>
> > On Thu, Jul 11, 2013 at 11:58 PM, Mohsen Soltani<br>
> > <<a href="mailto:soltani.clima@gmail.com" target="_blank">soltani.clima@gmail.com</a> <mailto:<a href="mailto:soltani.clima@gmail.com" target="_blank">soltani.clima@gmail.com</a>>> wrote:<br>
> ><br>
> > itesh,<br>
> ><br>
> > Make sure Lat/Lon in your .ctl file are the same as your script!<br>
> ><br>
> > Mohsen<br>
> ><br>
> > --<br>
> > April showers bring May flowers!<br>
> > --<br>
> > *Mohsen Soltani, M.Sc.*<br>
> > Climatology Grad Student,<br>
> > Department of Climatology, Faculty of Geography,<br>
> > University of Tehran, Tehran, Iran<br>
> > voice: <a href="tel:%28%2B98%29%209119772934" value="+989119772934" target="_blank">
(+98) 9119772934</a><br>
> > e-mail: <a href="mailto:soltani.clima@gmail.com" target="_blank">soltani.clima@gmail.com</a> <mailto:<a href="mailto:soltani.clima@gmail.com" target="_blank">soltani.clima@gmail.com</a>><br>
> ><br>
> ><br>
> > On Thu, Jul 11, 2013 at 9:05 AM, itesh dash <<a href="mailto:itesh@rimes.int" target="_blank">itesh@rimes.int</a><br>
> > <mailto:<a href="mailto:itesh@rimes.int" target="_blank">itesh@rimes.int</a>>> wrote:<br>
> ><br>
> > Dear Sir,<br>
> ><br>
> > I am working with the WRF and tried to plot the output in<br>
> > grads. I found there is nearly 1 degree shift in the<br>
> > projection and the plot.<br>
> ><br>
> > Please guide how to correct this.<br>
> ><br>
> > I am using mercator projection for the WRF model. Please find<br>
> > attached image.<br>
> ><br>
> > Seeking your advice on the same.<br>
> ><br>
> > regards,<br>
> > itesh<br>
> ><br>
> ><br>
> > _______________________________________________<br>
> > gradsusr mailing list<br>
> > <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a> <mailto:<a href="mailto:gradsusr@gradsusr.org" target="_blank">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>
> ><br>
> > _______________________________________________<br>
> > gradsusr mailing list<br>
> > <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a> <mailto:<a href="mailto:gradsusr@gradsusr.org" target="_blank">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>
> ><br>
> ><br>
> > _______________________________________________<br>
> > gradsusr mailing list<br>
> > <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
> > <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
><br>
> -------------- next part --------------<br>
> An HTML attachment was scrubbed...<br>
> URL:<br>
> <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>
><br>
> ------------------------------<br>
><br>
> _______________________________________________<br>
> gradsusr mailing list<br>
> <a href="mailto:gradsusr@gradsusr.org" target="_blank">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 19<br>
> ****************************************<br>
><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" target="_blank">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 20<br>
****************************************<br>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div>
</body>
</html>