<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p style="margin-top:0;margin-bottom:0">Hey Jennifer,</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Thanks a ton for you time. I'm almost there. I have the models plotted, now I just need to get the observations live plotted.</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">You sent a piece of code to write the ascii data as a binary time series. Does this does need to replace anything in the plume script? I am unsure what to do with it.<br>
</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0"></p>
<div>'set gxout fwrite'<br class="">
t=0<br class="">
while (1)<br class="">
res=read(ascfile)<br class="">
line1=sublin(res,1)<br class="">
line2=sublin(res,2)<br class="">
rc1 = subwrd(line1,1)<br class="">
if (rc1); break; endif<br class="">
ck=subwrd(line2,1)<br class="">
if (ck!='YR')<br class="">
t=t+1<br class="">
year=ck<br class="">
m=subwrd(line2,2)<br class="">
val=subwrd(line2,10)<br class="">
'd 'val<br class="">
if (t=1) <br class="">
iyear=year<br class="">
imon=subwrd(mons,m)<br class="">
endif<br class="">
endif<br class="">
endwhile<br class="">
'disable fwrite'</div>
<br>
<p></p>
<br>
<br>
<div style="color: rgb(0, 0, 0);">
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> gradsusr <gradsusr-bounces@gradsusr.org> on behalf of Jennifer M Adams <jadams21@gmu.edu><br>
<b>Sent:</b> Wednesday, June 27, 2018 12:17 PM<br>
<b>To:</b> GrADS Users Forum<br>
<b>Subject:</b> Re: [gradsusr] ENSO Forecast Plumes Script</font>
<div> </div>
</div>
<meta content="text/html; charset=utf-8">
<div class="" style="word-wrap:break-word; line-break:after-white-space">Use this :
<div class=""><br class="">
</div>
<div class=""><font class="" size="2" face="Courier">* define a grid of Nino3.4 index values in the T,E dimensions<br class="">
'set dfile 'fn<br class="">
'set t 1 'tmax<br class="">
'set x 1'<br class="">
'set y 1'<br class="">
'set e 1 'numens.m<br class="">
'grid'mused'=eloop(tloop(aave(tmpsfc,lon=-170,lon=-120,lat=-5,lat=5)))’</font><br class="">
<br class="">
instead of this:</div>
<div class=""><br class="">
</div>
<div class="">
<div class=""><font class="" size="2" face="Courier">* get the range of times</font></div>
<div class=""><font class="" size="2" face="Courier"> 'set dfile 'fn</font></div>
<div class=""><font class="" size="2" face="Courier"> 'set t 1 'tmax</font></div>
<div class=""><font class="" size="2" face="Courier"> 'q time'</font></div>
<div class=""><font class="" size="2" face="Courier"> time1=subwrd(result,3)</font></div>
<div class=""><font class="" size="2" face="Courier"> time2=subwrd(result,5)</font></div>
<div class=""><font class="" size="2" face="Courier"><br class="">
</font></div>
<div class=""><font class="" size="2" face="Courier">* build a server-side analysis URL</font></div>
<div class=""><font class="" size="2" face="Courier"> dset='nmme/'year%mm'/tmpsfc_'model.m</font></div>
<div class=""><font class="" size="2" face="Courier"> expr='eloop(tloop(aave(tmpsfc,lon=-170,lon=-120,lat=-5,lat=5)))' ;* NINO3.4 definition</font></div>
<div class=""><font class="" size="2" face="Courier"> dims='0:0,-90:-90,0:0,'time1':'time2',01:'numens.m</font></div>
<div class=""><font class="" size="2" face="Courier"> URL=server'_expr_{'dset'}{'expr'}{'dims'}'</font></div>
<div class=""><font class="" size="2" face="Courier"> 'close 'fn</font></div>
<div class=""><font class="" size="2" face="Courier">* This may take some time if the calculation result is not cached</font></div>
<div class=""><font class="" size="2" face="Courier"> 'sdfopen 'URL</font></div>
<div class=""><font class="" size="2" face="Courier"><br class="">
</font></div>
<div class=""><font class="" size="2" face="Courier">* define a grid of Nino3.4 index values in the T,E dimensions</font></div>
<div class=""><font class="" size="2" face="Courier"> 'set dfile 'fn</font></div>
<div class=""><font class="" size="2" face="Courier"> 'set t 1 'tmax</font></div>
<div class=""><font class="" size="2" face="Courier"> 'set x 1'</font></div>
<div class=""><font class="" size="2" face="Courier"> 'set y 1'</font></div>
<div class=""><font class="" size="2" face="Courier"> 'set e 1 'numens.m</font></div>
<div class=""><font class="" size="2" face="Courier"> 'grid'mused'=result'</font></div>
<div class=""><br class="">
</div>
—Jennifer<br class="">
<br class="">
<blockquote type="cite" class="">On Jun 27, 2018, at 10:28 AM, James Spinks <<a href="mailto:jspin881@hotmail.com" class="OWAAutoLink" id="LPlnk920057" previewremoved="true">jspin881@hotmail.com</a>> wrote:<br class="">
<br class="">
Hey Jennifer,<br class="">
<br class="">
I have successfully created a shell script to download the data with no problem and I have set up my directory to mimic some of the defined paths in the plume script. However, I am having an issue on getting past the "build a server-side analysis URL" section.<br class="">
<br class="">
This is what my section looks like:<br class="">
<br class="">
* build a server-side analysis URL<br class="">
dset='NMME/data/'year%mm'/'model.m'.tmpsfc.'year''mm'.anom.ctl'<br class="">
expr='eloop(tloop(aave(tmpsfc,lon=-170,lon=-120,lat=-5,lat=5)))' ;* NINO3.4 definition<br class="">
dims='0:0,-90:-90,0:0,'time1':'time2',01:'numens.m<br class="">
URL=server'_expr_{'dset'}{'expr'}{'dims'}'<br class="">
'close 'fn<br class="">
* This may take some time if the calculation result is not cached<br class="">
'open 'URL<br class="">
<br class="">
I modified the dset variable, but I "can't open the description" on the open URL part. I am assuming since this was server-side analysis, this particular section needs to be modified since I now download the data? What would be the best way to go about modifying
this section? Does the expr and dims variable remain the same? What would I do about the URL variable?<br class="">
<br class="">
Thanks,<br class="">
James<br class="">
<br class="">
<br class="">
From: gradsusr <<a href="mailto:gradsusr-bounces@gradsusr.org" class="OWAAutoLink" id="LPlnk412451" previewremoved="true">gradsusr-bounces@gradsusr.org</a>> on behalf of Jennifer M Adams <jadams21@gmu.edu><br class="">
Sent: Tuesday, June 26, 2018 6:25 PM<br class="">
To: GrADS Users Forum<br class="">
Subject: Re: [gradsusr] ENSO Forecast Plumes Script<br class="">
<br class="">
I use the ‘dtype netcdf’ descriptor to make sure the dimensions are mapping properly, set the ensemble names with leading zeros, and alias the variable name so the plume script will work with new and old data. <br class="">
<br class="">
dset ^NASA_GEOS5v2.tmpsfc.201806.anom.nc<br class="">
dtype netcdf<br class="">
options yrev<br class="">
undef 9.999E+20<br class="">
title NASA_GEOS5v2 anomalies for tmpsfc initialized Jun 2018<br class="">
ydef 181 linear -90 1<br class="">
xdef 360 linear 0 1<br class="">
zdef 1 linear 1 1<br class="">
tdef 9 linear 01Jun2018 1mo<br class="">
edef 10 names 01 02 03 04 05 06 07 08 09 10 <br class="">
vars 1<br class="">
fcst=>tmpsfc 0 e,t,y,x surface temperature [K]<br class="">
endvars<br class="">
<br class="">
—Jennifer<br class="">
<br class="">
<br class="">
<blockquote type="cite" class="">On Jun 26, 2018, at 3:59 PM, James Spinks <jspin881@hotmail.com> wrote:<br class="">
<br class="">
Hey Jennifer,<br class="">
<br class="">
Thanks. You're right. After downloading the data in netcdf format, the files opened with no problem. What would I use the descriptor file if the nc files open without it? Just curious.<br class="">
<br class="">
I'll give your suggestions a try. Thanks for the help.<br class="">
From: gradsusr <gradsusr-bounces@gradsusr.org> on behalf of Jennifer M Adams <jadams21@gmu.edu><br class="">
Sent: Tuesday, June 26, 2018 3:00 PM<br class="">
To: GrADS Users Forum<br class="">
Subject: Re: [gradsusr] ENSO Forecast Plumes Script<br class="">
<br class="">
Hi, James — <br class="">
The data I use to draw the plumes is freely available from NOAA:<br class="">
<br class="">
ftp://ftp.cpc.ncep.noaa.gov/wd52dg/data/indices/sstoi.indices (for the obs)<br class="">
ftp://ftp.cpc.ncep.noaa.gov/NMME/realtime_anom (for the model forecasts)<br class="">
<br class="">
You will need to download the data and create your own descriptor files, but that is a very straightforward task — no special tricks. Recently the format of the model data changed -- now it is NetCDF instead of GRIB — which makes it even easier. Below is the
bit of script code for writing the ascii data as a binary time series: <br class="">
<br class="">
'set gxout fwrite'<br class="">
t=0<br class="">
while (1)<br class="">
res=read(ascfile)<br class="">
line1=sublin(res,1)<br class="">
line2=sublin(res,2)<br class="">
rc1 = subwrd(line1,1)<br class="">
if (rc1); break; endif<br class="">
ck=subwrd(line2,1)<br class="">
if (ck!='YR')<br class="">
t=t+1<br class="">
year=ck<br class="">
m=subwrd(line2,2)<br class="">
val=subwrd(line2,10)<br class="">
'd 'val<br class="">
if (t=1) <br class="">
iyear=year<br class="">
imon=subwrd(mons,m)<br class="">
endif<br class="">
endif<br class="">
endwhile<br class="">
'disable fwrite'<br class="">
<br class="">
<br class="">
—Jennifer<br class="">
<br class="">
<blockquote type="cite" class="">On Jun 26, 2018, at 10:40 AM, James Spinks <jspin881@hotmail.com> wrote:<br class="">
<br class="">
Hi all,<br class="">
<br class="">
I was wondering if there was an updated link for the ENSO Forecast Plumes script since the URL http://monsoondata.org:9090/dods/ has been decommissioned.<br class="">
<br class="">
The script is hosted here http://cola.gmu.edu/enso_forecast_plumes.html and looks to still be running just fine, but the script does not contain an updated link. If there is a replacement link, can some please provide the new link. <br class="">
<br class="">
Thank you,<br class="">
James<br class="">
_______________________________________________<br class="">
gradsusr mailing list<br class="">
gradsusr@gradsusr.org<br class="">
http://secure-web.cisco.com/10cV7uvP1gMb6aGhsXoFmP9FcAbE16taY_R_CEn2DeLZcOVmzB2z8hEltkCJAzIqJc222mfbwJD1BXPCTrrkdCW-92A_3qXj4-zk3pcAGXRXDsfDVGyoNujgAffAPmVw613uXORG6J4-KHEuak38h5fRR3-mpQvMbpqwaSo_paZU-ZzrCz_sA6HXlZ89camd3E1ugdUdARvcxwvWRbX8rH9ZZSlsIGM92_lQYuTp69dFfq4B4T7GhQ2WibEwTgESteR32Zlwkwd9pQsrMU03xMITiJuRmVXIhsvqg3q1eaSud5BN6-nYbJf_z6FattWnWt2f6TGx0vl_nGZBb226xlQb0_sb5y_UnBh_mJFCN-KWM1VE66DaDGW1VIHUe48-EOTQxnnp7IvtlxkSenu6GSQOcOu9k87zqLTajNFICvw6HGFTbrE9OGDNxb-fyIKk2/http%3A%2F%2Fgradsusr.org%2Fmailman%2Flistinfo%2Fgradsusr<br class="">
</blockquote>
<br class="">
--<br class="">
Jennifer Miletta Adams<br class="">
Center for Ocean-Land-Atmosphere Studies (COLA)<br class="">
George Mason University<br class="">
<br class="">
<br class="">
<br class="">
_______________________________________________<br class="">
gradsusr mailing list<br class="">
gradsusr@gradsusr.org<br class="">
http://secure-web.cisco.com/1j6CyJMJVEv8urj2N4H7sGUrV2dRpugRfTlhcHDzRNMticAKWRHmaBmIKjMfaYXw_RSaJEyzOaYOdgU5QvfUuUZptVu_wCu14zwse0TMWu0S4yWZ0gzglRu0F_tkiti4HwylPGmUISs-lIYnhLCsshVTep0mJw9VE0TERJR9CFMAJ-1G2VsIDikZMcWPNXIrAbmBsKydCKQQpK9x7g4YRsnvjWp7yk1lmqXglJ-2b6LBMHdxUkQ9JCgoavyfykgMfTIblwe7xVmaK3oFAaZwnVisWGEI_OlJtAz-c6vjI9MX_RgzCBRgE1IlPeWTKuzbxoS1UjfDMu0_YDiD6ZaR_pdbTEJ_C3rD4UsZUvNQhwcy5nw52Kvpnrm7UJ3ZZ4i0ONpBtZ5FDQ4GiEU2JJBqxFXIzHQJhFUmu4zUaw2sGoDeawoG63un6BdkxPNnhtvI7/http%3A%2F%2Fgradsusr.org%2Fmailman%2Flistinfo%2Fgradsusr<br class="">
</blockquote>
<br class="">
--<br class="">
Jennifer Miletta Adams<br class="">
Center for Ocean-Land-Atmosphere Studies (COLA)<br class="">
George Mason University<br class="">
<br class="">
<br class="">
<br class="">
_______________________________________________<br class="">
gradsusr mailing list<br class="">
gradsusr@gradsusr.org<br class="">
http://secure-web.cisco.com/1RRyoo5uzU0u6Cqodtm93v7Z_GBwgPJOQu83b-zumZdQC5VOws-H5L740i07_Hy8WYFWEzu1F6oH6Kyxrh_mA-mmnb3EE-lWBcCy16poXhELNxFlIX4RFAhs3IWps7ka3807b_bSTuRxwwB8HEQVRR1DEG6PBYHJ7hR5FdGBNqorJmkj9NUjRzKVbrwjqeKZGiwVOewxzy8gXRuGylmSdm6YxY3B_ZaxDGGIf6HnyaekLzxKaV3NuvWDlPq4c5L6TYFmk2_996AzRbm-YfeU5ZP28aZ6o8jkPJnTvYr9KuRrtCzZknrkCca9AJzLZjf0vmI-ig0R39PfSvufdrjiYD_PFHu4P7cRtNHrWAtOyNepLTD-OvMdCAIPyAKoKRZwtTburkfmfAJncuuC0enNmOzzQgKilI52q94xFndpP_M8chyeVXO4Cxkqup6f7Zc93/http%3A%2F%2Fgradsusr.org%2Fmailman%2Flistinfo%2Fgradsusr<br class="">
</blockquote>
<br class="">
<div class="">--<br class="">
Jennifer Miletta Adams<br class="">
Center for Ocean-Land-Atmosphere Studies (COLA)<br class="">
George Mason University<br class="">
<br class="">
<br class="">
</div>
<br class="">
</div>
</div>
</div>
</div>
</body>
</html>