<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>
Your code has a minor mistake, This is how I would remove the seasonal cycle:<br>
<br>
<span dir="ltr"><b>set t 1 840</b><b><br>
define nino34=(aave((sst-273.15),lon=120,lon=170,lat=-5,lat=5))</b><b><br>
set t 1 12</b><b><br>
define nino34clim=ave(nino34,t+0,t=840,12)</b><b><br>
modify nino34clim seasonal</b><b><br>
set t 1 840</b><b><br>
d nino34 - nino34clim<br>
<br>
</b>Note, there is a minor difference in how you and I defined the second argument from the average function. You need to loop through all times, that is the reason why you need the second argument as t+0 instead of t=1.<br>
<br>
Hope this help!!!<br>
</span>
<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="divRpF433415"><font color="#000000" size="2" face="Tahoma"><b>From:</b> gradsusr-bounces@gradsusr.org [gradsusr-bounces@gradsusr.org] on behalf of Corey Gabriel [cjgabriel7@gmail.com]<br>
<b>Sent:</b> Friday, July 12, 2013 10:49 AM<br>
<b>To:</b> gradsusr@gradsusr.org<br>
<b>Subject:</b> [gradsusr] Monthly Clim. and least Squares_2020-2089 time series (montly data)<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">
<div>
<div>
<div>
<div>Hosmay and James:<br>
</div>
Thank you so much for taking the time out of your day to answer my question.<br>
</div>
Dear All: As per Jim and Hosmay's instructions, I invoked the following code to produce climatology in the attached figure. Their instructions correspond with the documentation provided here:
<a href="http://grads.wikidot.com/algorithms" target="_blank">http://grads.wikidot.com/algorithms</a>. There are still two problems:<br>
</div>
1.) I thought that the code copied below would remove the seasonal cycle from my anomaly map. It still looks like the anomaly is centered around an annual average, rather than making a comparison between Jan clim. and Jan actual sst, February, March and so
on. here is the code:<br>
</div>
<b>**Bolded portion produces the attached figure.<br>
</b>
<div><b>set t 1 840<br>
define nino34=(aave((sst-273.15),lon=120,lon=170,lat=-5,lat=5))<br>
set t 1 12<br>
define nino34clim=ave(nino34,t=1,t=840,12)<br>
modify nino34clim seasonal<br>
set t 1 840<br>
d nino34 - nino34clim</b><br>
</div>
<div>***This is the attempt to perform the least-squares regression.<br>
</div>
<div>c<br>
set cstyle 3<br>
set ccolor 1<br>
set cmark 0<br>
define cff=tregr(nino34,nino34clim,t=1,t=840)<br>
define linearnino34clim=cff*nino34+nino34clim<br>
d nino34-linearnino34clim<br>
</div>
<div>2.) After successfully creating the anomaly from climatology, what I'd like to do is draw a regression line throughout the time series, and then compare the actual sst to the best fit for each month in the time series. The departure from the trend will
signal ENSO events and allow me to look deeper into the system.<br>
<br>
</div>
<div>Again, any help would be greatly appreciated!<br>
</div>
<div>Corey<br>
</div>
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">On Thu, Jul 11, 2013 at 3:47 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: 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 <<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: <1122f292023cec649d99f1fda15546f8@localhost><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 "OPTIONS template", 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-> 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-><br>
<br>
But I still could not see the time values ??of four, the following message<br>
appears below:<br>
<br>
ga-> set t 1 4<br>
Time values set: 2002:1:1:0 2002:1:4:0<br>
ga-> 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 'wind'<br>
Error ocurred at column 1<br>
DISPLAY error: Invalid expression<br>
Expression = wind<br>
ga-><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, "James T. Potemra" <<a href="mailto:jimp@hawaii.edu" target="_blank">jimp@hawaii.edu</a>><br>
wrote:<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<br>
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"<br>
<<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<br>
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<br>
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<br>
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:<br>
<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>
> 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 19:39:13 +0000<br>
From: Hosmay Lopez <<a href="mailto:hlopez@rsmas.miami.edu" target="_blank">hlopez@rsmas.miami.edu</a>><br>
Subject: Re: [gradsusr] Least-Squares Regression and Climatology (Time<br>
Series)<br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>
Message-ID:<br>
<<a href="mailto:8A902EB84795C94C8726F11DC8C574A81D26BA7B@EXCHANGE1.amazon.rsmas.miami.edu" target="_blank">8A902EB84795C94C8726F11DC8C574A81D26BA7B@EXCHANGE1.amazon.rsmas.miami.edu</a>><br>
<br>
Content-Type: text/plain; charset="iso-8859-1"<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" target="_blank">gradsusr-bounces@gradsusr.org</a> [<a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a>] on behalf of Corey Gabriel [<a href="mailto:cjgabriel7@gmail.com" target="_blank">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: *'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'<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>
On Thu, Jul 11, 2013 at 2:22 PM, <<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a><mailto:<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a>>> wrote:<br>
Send gradsusr mailing list submissions to<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>
<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><mailto:<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><mailto:<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><mailto:<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><mailto:<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><mailto:<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><mailto:<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><mailto:<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><mailto:<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><mailto:<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><mailto:<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><mailto:<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><<a href="http://m2-00000.nc" target="_blank">http://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><<a href="http://m2-00000.nc" target="_blank">http://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><mailto:<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><<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://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><<a href="http://m2-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-02-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-03-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://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><mailto:<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><<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://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><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>
> > 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>
> 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><mailto:<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><mailto:<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><mailto:<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><tel:%28%2B98%29%209119772934><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><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 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><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>
> -------------- 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><mailto:<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><mailto:<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><mailto:<a href="mailto:JDB1t6UvJRG9jO%252B1L%252BfAD20Gd4G8X_Qayg@mail.gmail.com" target="_blank">JDB1t6UvJRG9jO%2B1L%2BfAD20Gd4G8X_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><mailto:<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><tel:%28%2B98%29%209119772934><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><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 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><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>
> -------------- 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><mailto:<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><mailto:<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><mailto:<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><<a href="http://m2-00000.nc" target="_blank">http://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><<a href="http://m2-00000.nc" target="_blank">http://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><mailto:<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><<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://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><<a href="http://m2-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-02-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-03-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-01-00000.nc" target="_blank">http://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><mailto:<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><<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://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><<a href="http://amazonia.clm2.h0.2002-01-04-00000.nc" target="_blank">http://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><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>
> >> 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>
> 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><mailto:<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><mailto:<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><mailto:<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>> <mailto:<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><tel:%28%2B98%29%209119772934><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>> <mailto:<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><mailto:<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><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>> <mailto:<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>> <mailto:<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><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>
> -------------- 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><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>
> 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><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>
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" 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 23<br>
****************************************<br>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div>
</body>
</html>