NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
Tom, The problem exists because the AVN thinned grids provide 6 hour accumulations: 135 010920/1200F006 0 NONE P06M 209 010920/1200F012 0 NONE P06M 283 010920/1200F018 0 NONE P06M 357 010920/1200F024 0 NONE P06M 431 010920/1200F030 0 NONE P06M 505 010920/1200F036 0 NONE P06M 579 010920/1200F042 0 NONE P06M 653 010920/1200F048 0 NONE P06M 727 010920/1200F060 0 NONE P06M 801 010920/1200F072 0 NONE P06M When you ask to plot P12I, GEMPAK knows how to compute this using the 6 hourly accumulations. However, to plot a P12I for F060, you need the P06I at f054 and f060. And, to plot P12I for F072, you need the P06I at f066 and F072. The problem boils down to not having AVN grids at f054 and f066 when trying to compute 12 hour accumulations. You can plot the 6 hour accumulations because they exist in the grid file. In the case of the ETA 211 grids, the following grids exist: 305 010920/1200F006 0 NONE P06M 614 010920/1200F018 0 NONE P06M 925 010920/1200F030 0 NONE P06M 1234 010920/1200F042 0 NONE P06M 1543 010920/1200F054 0 NONE P06M 459 010920/1200F012 0 NONE P12M 770 010920/1200F024 0 NONE P12M 1079 010920/1200F036 0 NONE P12M 1388 010920/1200F048 0 NONE P12M 1699 010920/1200F060 0 NONE P12M The 12 hour accumulations are at the times evenly divisible by 12. The 6 hourly accumulations are at the intervening times. GEMPAK can compute the 6 hourly accumulations at all times from f006 to f060 given the above fields, eg: gdattim = f006-f060 glevel = 0 gvcord = none gfunc = p06i or the 12 hourly accumulations: gdatttim = f012-f060-12 or f012-f060-6 gfunc = p12i Since EMSL is available at every time step f000-f060, you should be able to produce a plot for that parameter at all forecast times (providing you received the grid). You can verify the times you have EMSL grid for using gdinfo. gdattim = all gfunc = emsl Steve Chiswell Unidata User Support >From: Tom Priddy <address@hidden> >Organization: University of Kentucky >Keywords: 200109181935.f8IJZj114835 > >--------------9A400F51D4051C8FCF1729F2 >Content-Type: text/plain; charset=us-ascii >Content-Transfer-Encoding: 7bit > >From: Tom Priddy > >I recently started plotting surface forecast maps for the additional >times (60 and 72 hrs) using the avn_thin.gem files. But I can't seem to >get the 12 qpf >plotted on these maps. The 12 hr qpf for 60 and 72 hours should be >P12I... >but all I can plot is the 6 hr forecast accumulation (P06I). I have a >similar >problem is the 60 hour eta surface map for surface pressure. Can you >take >a look at my script for problems. >---- >#!/bin/csh -f ># ># This script will produce the forecasted surface conditions when it is ># called upon by the script metforecast.csh. That script is called upon ># by several other scripts after the eta, ngm and nhmrf models are run. ># The eta and ngm are run twice a day and the nhmrf once a day. ># ># Approximate run time: ># >cd /export/home/gempak/scripts/LDM >source ~gempak/.cshrc >source ~gempak/GEMsetup > ># ># ># Check arguments ># >if ( $#argv == 4 ) then > set model = $argv[1] > set run = $argv[2] > set valid = $argv[3] > set YYMMDD = `date -u "+%y%m%d"` > set ddate2 = `date "+%D"` > set level = $argv[4] >else if ( $#argv == 5 ) then > set model = $argv[1] > set run = $argv[2] > set valid = $argv[3] > set YYMMDD = $argv[4] > set ddate2 = $argv[5] > set level = all >else > echo "" > echo USAGE: metforecast.csh model run time _date_ YYMMDD GMT format > echo "" > exit >endif ># ># Check model file ># >set model_file = `find /ldmdata/gempak/hds/${YYMMDD}${run}_${model}.gem >-print` > >if ($#model_file == 0) then > echo ${YYMMDD}${run}_${model} not in - maps not made | mail >address@hidden > exit >endif > >set model_size = `ls -l /ldmdata/gempak/hds/${YYMMDD}${run}_${model}.gem >| awk '{print $5}'` > >if ($model_size <= 400000) then > exit >endif > >#if ($model_size <= 6000000 && $model == "ngm" || $model_size <= 6000000 >&& $model == "eta") then ># exit >#endif ># ># Check model choice and set pressure variable. Establish filter for ># wind barbs for gdwind. ># >set skip = "/-4" >if (${model} == ngm||${model} == thin) then > set presfield = "PMSL" > else if (${model} == nhmrf||${model} == ruc) then > set presfield = "PMSL" > set skip = "0" > else if (${model} == eta) then > set presfield = "EMSL" > else > echo "" > echo INVALID MODEL: ngm eta nhmrf > echo "" > exit >endif ># ># Set title information and select correct precipitation variable based ># on the valid forecast time. ># >set titleinfo = "${valid} Hr Forecast" > >if (${valid} == 6 || ${valid} == 18 || ${valid} == 30 || ${valid} == 42) >then > set precip = "P06I " >else > set precip = "P12I" >endif > ># ># Begin section for surface data. ># >set fname = "${YYMMDD}${run}_${model}.${valid}_sfc" >set output = "${fname}.gif" >set device = "gfw|${output}|792;612" > >rm -f ${output} > >gdcntr << EOF > \$MAPFIL = hipowo.cia > \$RESPOND = no > GDATTIM = f${valid} > GLEVEL = 0 > GVCORD = none > GFUNC = ${precip} > GDFILE = /ldmdata/gempak/hds/${YYMMDD}${run}_${model}.gem > CINT = .05;.1;.25;.5;.75;1;1.5;2;3;4;5 > LINE = 23/2 > MAP = 10 > TITLE = 0 > DEVICE = ${device} > SATFIL = 0 > RADFIL = 0 > PROJ = lcc > GAREA = uslcc > CLEAR = yes > PANEL = 0 > TEXT = 1/21//hw > SCALE = 0 > LATLON = 0 > HILO = 0 > HLSYM = 0 > CLRBAR = 24/H/LL//.75;.02/-1 > CONTUR = 2 > SKIP = 0 > FINT = .05;.1;.25;.5;.75;1;1.5;2;3;4;5 > FLINE = 101;20-25;16-21 > CTYPE = f > LUTFIL = 0 > STNPLT = 0 > run > > exit >EOF > >gdcntr << EOF > CLEAR = no > GLEVEL = 1000 > GVCORD = pres > GFUNC = tmpf > CINT = 10 > LINE = 2/3 > TITLE = 24/-3/${ddate2} ${run}Z ${model} - ${titleinfo} - Surface >Conditions > CTYPE = c > LUTFIL = 0 > STNPLT = 0 > run > > exit >EOF > >gdcntr << EOF > GLEVEL = 0 > GVCORD = none > GFUNC = ${presfield} > CINT = 4 > LINE = 4/1/2 > HILO = 4;2/H#;L#///1;1/ > HLSYM = 2;1/2/2//HW > CTYPE = c > LUTFIL = 0 > STNPLT = 0 > run > > exit >EOF > >gpend > >#/export/home/gempak/scripts/add_logo >/export/home/gempak/scripts/LDM/${fname} ${model}_${valid}_sfc logo.gif > >/usr/local/bin/combine -compose over -colors 256 logo.gif ${fname}.gif >${model}_${valid}_sfc.gif >/usr/local/bin/giftool -B -i -C ${model}_${valid}_sfc.gif >mv -f ${model}_${valid}_sfc.gif /ldmdata/images >rm -f *.nts ${output} > > >-- >============================================================ >Tom Priddy address@hidden >Extension Ag. Meteorologist ph: 859-257-3000 ext 245 >Ag. Weather Center fax: 859-257-5671 >243 Charles E. Barnhart Bldg. >Dept. of Biosystems and Ag. Engr. >University of Kentucky >Lexington, KY 40546-0276 >HOMEPAGE http://wwwagwx.ca.uky.edu >============================================================ > > >--------------9A400F51D4051C8FCF1729F2 >Content-Type: text/html; charset=us-ascii >Content-Transfer-Encoding: 7bit > ><!doctype html public "-//w3c//dtd html 4.0 transitional//en"> ><html> >From: Tom Priddy ><p>I recently started plotting surface forecast maps for the additional >times (60 and 72 hrs) using the avn_thin.gem files. But I can't seem to >get the 12 qpf ><br>plotted on these maps. The 12 hr qpf for 60 and 72 hours should be >P12I... ><br>but all I can plot is the 6 hr forecast accumulation (P06I). I have >a similar ><br>problem is the 60 hour eta surface map for surface pressure. Can you >take ><br>a look at my script for problems. ><br>---- ><br><tt><font size=-1>#!/bin/csh -f</font></tt> ><br><tt><font size=-1>#</font></tt> ><br><tt><font size=-1># This script will produce the forecasted surface >conditions when it is</font></tt> ><br><tt><font size=-1># called upon by the script metforecast.csh. That >script is called upon</font></tt> ><br><tt><font size=-1># by several other scripts after the eta, ngm and >nhmrf models are run.</font></tt> ><br><tt><font size=-1># The eta and ngm are run twice a day and the nhmrf >once a day.</font></tt> ><br><tt><font size=-1>#</font></tt> ><br><tt><font size=-1># Approximate run time:</font></tt> ><br><tt><font size=-1>#</font></tt> ><br><tt><font size=-1>cd /export/home/gempak/scripts/LDM</font></tt> ><br><tt><font size=-1>source ~gempak/.cshrc</font></tt> ><br><tt><font size=-1>source ~gempak/GEMsetup</font></tt><tt><font size=-1></f > ont></tt> ><p><tt><font size=-1>#</font></tt> ><br><tt><font size=-1>#</font></tt> ><br><tt><font size=-1># Check arguments</font></tt> ><br><tt><font size=-1>#</font></tt> ><br><tt><font size=-1>if ( $#argv == 4 ) then</font></tt> ><br><tt><font size=-1> set model = $argv[1]</font></tt> ><br><tt><font size=-1> set run = $argv[2]</font></tt> ><br><tt><font size=-1> set valid = $argv[3]</font></tt> ><br><tt><font size=-1> set YYMMDD = `date -u "+%y%m%d"`</font></tt> ><br><tt><font size=-1> set ddate2 = `date "+%D"`</font></tt> ><br><tt><font size=-1> set level = $argv[4]</font></tt> ><br><tt><font size=-1>else if ( $#argv == 5 ) then</font></tt> ><br><tt><font size=-1> set model = $argv[1]</font></tt> ><br><tt><font size=-1> set run = $argv[2]</font></tt> ><br><tt><font size=-1> set valid = $argv[3]</font></tt> ><br><tt><font size=-1> set YYMMDD = $argv[4]</font></tt> ><br><tt><font size=-1> set ddate2 = $argv[5]</font></tt> ><br><tt><font size=-1> set level = all</font></tt> ><br><tt><font size=-1>else</font></tt> ><br><tt><font size=-1> echo ""</font></tt> ><br><tt><font size=-1> echo USAGE: metforecast.csh model run time >_date_ YYMMDD GMT format</font></tt> ><br><tt><font size=-1> echo ""</font></tt> ><br><tt><font size=-1> exit</font></tt> ><br><tt><font size=-1>endif</font></tt> ><br><tt><font size=-1>#</font></tt> ><br><tt><font size=-1># Check model file</font></tt> ><br><tt><font size=-1>#</font></tt> ><br><tt><font size=-1>set model_file = `find /ldmdata/gempak/hds/${YYMMDD}${ru > n}_${model}.gem >-print`</font></tt><tt><font size=-1></font></tt> ><p><tt><font size=-1>if ($#model_file == 0) then</font></tt> ><br><tt><font size=-1> echo ${YYMMDD}${run}_${model} not in >- maps not made | mail address@hidden</font></tt> ><br><tt><font size=-1> exit</font></tt> ><br><tt><font size=-1>endif</font></tt> ><br><tt><font size=-1> </font></tt> ><br><tt><font size=-1>set model_size = `ls -l /ldmdata/gempak/hds/${YYMMDD}${r > un}_${model}.gem >| awk '{print $5}'`</font></tt><tt><font size=-1></font></tt> ><p><tt><font size=-1>if ($model_size <= 400000) then</font></tt> ><br><tt><font size=-1> exit</font></tt> ><br><tt><font size=-1>endif</font></tt> ><br><tt><font size=-1> </font></tt> ><br><tt><font size=-1>#if ($model_size <= 6000000 && $model >== "ngm" || $model_size <= 6000000 && $model == "eta") then</font>< > /tt> ><br><tt><font size=-1># exit</font></tt> ><br><tt><font size=-1>#endif</font></tt> ><br><tt><font size=-1>#</font></tt> ><br><tt><font size=-1># Check model choice and set pressure variable. Establis > h >filter for</font></tt> ><br><tt><font size=-1># wind barbs for gdwind.</font></tt> ><br><tt><font size=-1>#</font></tt> ><br><tt><font size=-1>set skip = "/-4"</font></tt> ><br><tt><font size=-1>if (${model} == ngm||${model} == thin) then</font></tt> ><br><tt><font size=-1> set presfield = "PMSL"</font></tt> ><br><tt><font size=-1> else if (${model} == nhmrf||${model} == ruc) >then</font></tt> ><br><tt><font size=-1> set presfield = "PMSL"</font></tt> ><br><tt><font size=-1> set skip = "0"</font></tt> ><br><tt><font size=-1> else if (${model} == eta) then</font></tt> ><br><tt><font size=-1> set presfield = "EMSL"</font></tt> ><br><tt><font size=-1> else</font></tt> ><br><tt><font size=-1> echo ""</font></tt> ><br><tt><font size=-1> echo INVALID MODEL: ngm eta nhmrf</font></t > t> ><br><tt><font size=-1> echo ""</font></tt> ><br><tt><font size=-1> exit</font></tt> ><br><tt><font size=-1>endif</font></tt> ><br><tt><font size=-1>#</font></tt> ><br><tt><font size=-1># Set title information and select correct precipitation >variable based</font></tt> ><br><tt><font size=-1># on the valid forecast time.</font></tt> ><br><tt><font size=-1>#</font></tt> ><br><tt><font size=-1>set titleinfo = "${valid} Hr Forecast"</font></tt><tt><f > ont size=-1></font></tt> ><p><tt><font size=-1>if (${valid} == 6 || ${valid} == 18 || ${valid} == >30 || ${valid} == 42) then</font></tt> ><br><tt><font size=-1> set precip = "P06I "</font></tt> ><br><tt><font size=-1>else</font></tt> ><br><tt><font size=-1> set precip = "P12I"</font></tt> ><br><tt><font size=-1>endif</font></tt><tt><font size=-1></font></tt> ><p><tt><font size=-1>#</font></tt> ><br><tt><font size=-1># Begin section for surface data.</font></tt> ><br><tt><font size=-1>#</font></tt> ><br><tt><font size=-1>set fname = "${YYMMDD}${run}_${model}.${valid}_sfc"</fon > t></tt> ><br><tt><font size=-1>set output = "${fname}.gif"</font></tt> ><br><tt><font size=-1>set device = "gfw|${output}|792;612"</font></tt><tt><fon > t size=-1></font></tt> ><p><tt><font size=-1>rm -f ${output}</font></tt><tt><font size=-1></font></tt> ><p><tt><font size=-1>gdcntr << EOF</font></tt> ><br><tt><font size=-1> \$MAPFIL = hipowo.cia</font></tt> ><br><tt><font size=-1> \$RESPOND = no</font></tt> ><br><tt><font size=-1> GDATTIM = f${valid}</font></tt> ><br><tt><font size=-1> GLEVEL = 0</font></tt> ><br><tt><font size=-1> GVCORD = none</font></tt> ><br><tt><font size=-1> GFUNC = ${precip}</font></ > tt> ><br><tt><font size=-1> GDFILE = /ldmdata/gempak/hds/${Y > YMMDD}${run}_${model}.gem</font></tt> ><br><tt><font size=-1> CINT = .05;.1;.25;.5 > ;.75;1;1.5;2;3;4;5</font></tt> ><br><tt><font size=-1> LINE = 23/2</font></ > tt> ><br><tt><font size=-1> MAP = 10</font > ></tt> ><br><tt><font size=-1> TITLE = 0</font></tt> ><br><tt><font size=-1> DEVICE = ${device}</font></tt> ><br><tt><font size=-1> SATFIL = 0</font></tt> ><br><tt><font size=-1> RADFIL = 0</font></tt> ><br><tt><font size=-1> PROJ = lcc</font></t > t> ><br><tt><font size=-1> GAREA = uslcc</font></tt> ><br><tt><font size=-1> CLEAR = yes</font></tt> ><br><tt><font size=-1> PANEL = 0</font></tt> ><br><tt><font size=-1> TEXT = 1/21//hw</fon > t></tt> ><br><tt><font size=-1> SCALE = 0</font></tt> ><br><tt><font size=-1> LATLON = 0</font></tt> ><br><tt><font size=-1> HILO = 0</font></tt> ><br><tt><font size=-1> HLSYM = 0</font></tt> ><br><tt><font size=-1> CLRBAR = 24/H/LL//.75;.02/-1</fo > nt></tt> ><br><tt><font size=-1> CONTUR = 2</font></tt> ><br><tt><font size=-1> SKIP = 0</font></tt> ><br><tt><font size=-1> FINT = .05;.1;.25;.5 > ;.75;1;1.5;2;3;4;5</font></tt> ><br><tt><font size=-1> FLINE = 101;20-25;16-21</f > ont></tt> ><br><tt><font size=-1> CTYPE = f</font></tt> ><br><tt><font size=-1> LUTFIL = 0</font></tt> ><br><tt><font size=-1> STNPLT = 0</font></tt> ><br><tt><font size=-1> run</font></tt> ><br><tt><font size=-1> </font></tt> ><br><tt><font size=-1> exit</font></tt> ><br><tt><font size=-1>EOF</font></tt><tt><font size=-1></font></tt> ><p><tt><font size=-1>gdcntr << EOF</font></tt> ><br><tt><font size=-1> CLEAR = no</font></tt> ><br><tt><font size=-1> GLEVEL = 1000</font></tt> ><br><tt><font size=-1> GVCORD = pres</font></tt> ><br><tt><font size=-1> GFUNC = tmpf</font></tt> ><br><tt><font size=-1> CINT = 10</font></tt> ><br><tt><font size=-1> LINE = 2/3</font></tt> ><br><tt><font size=-1> TITLE = 24/-3/${ddate2} ${run}Z >${model} - ${titleinfo} - Surface Conditions</font></tt> ><br><tt><font size=-1> CTYPE = c</font></tt> ><br><tt><font size=-1> LUTFIL = 0</font></tt> ><br><tt><font size=-1> STNPLT = 0</font></tt> ><br><tt><font size=-1> run</font></tt><tt><font size=-1></font></tt> ><p><tt><font size=-1> exit</font></tt> ><br><tt><font size=-1>EOF</font></tt><tt><font size=-1></font></tt> ><p><tt><font size=-1>gdcntr << EOF</font></tt> ><br><tt><font size=-1> GLEVEL = 0</font></tt> ><br><tt><font size=-1> GVCORD = none</font></tt> ><br><tt><font size=-1> GFUNC = ${presfield}</font></tt> ><br><tt><font size=-1> CINT = 4</font></tt> ><br><tt><font size=-1> LINE = 4/1/2</font></tt> ><br><tt><font size=-1> HILO = 4;2/H#;L#///1;1/</font></ > tt> ><br><tt><font size=-1> HLSYM = 2;1/2/2//HW</font></tt> ><br><tt><font size=-1> CTYPE = c</font></tt> ><br><tt><font size=-1> LUTFIL = 0</font></tt> ><br><tt><font size=-1> STNPLT = 0</font></tt> ><br><tt><font size=-1> run</font></tt><tt><font size=-1></font></tt> ><p><tt><font size=-1> exit</font></tt> ><br><tt><font size=-1>EOF</font></tt> ><br><tt><font size=-1> </font></tt> ><br><tt><font size=-1>gpend</font></tt><tt><font size=-1></font></tt> ><p><tt><font size=-1>#/export/home/gempak/scripts/add_logo /export/home/gempak > /scripts/LDM/${fname} >${model}_${valid}_sfc logo.gif</font></tt><tt><font size=-1></font></tt> ><p><tt><font size=-1>/usr/local/bin/combine -compose over -colors 256 logo.gif >${fname}.gif ${model}_${valid}_sfc.gif</font></tt> ><br><tt><font size=-1>/usr/local/bin/giftool -B -i -C ${model}_${valid}_sfc.gi > f</font></tt> ><br><tt><font size=-1>mv -f ${model}_${valid}_sfc.gif /ldmdata/images</font></ > tt> ><br><tt><font size=-1>rm -f *.nts ${output}</font></tt> ><br><tt><font size=-1></font></tt> ><p>-- ><br>============================================================ ><br>Tom Priddy &nbs > p; &nb > sp; >address@hidden ><br>Extension Ag. Meteorologist   > ; >ph: 859-257-3000 ext 245 ><br>Ag. Weather Center &n > bsp; >fax: 859-257-5671 ><br>243 Charles E. Barnhart Bldg. ><br>Dept. of Biosystems and Ag. Engr. ><br>University of Kentucky ><br>Lexington, KY 40546-0276 ><br>HOMEPAGE <A HREF="http://wwwagwx.ca.uky.edu">http://wwwagwx.ca > .uky.edu</A> ><br>============================================================ ><br> </html> > >--------------9A400F51D4051C8FCF1729F2-- >