[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
20010918: trouble with precip on avn_thin
- Subject: 20010918: trouble with precip on avn_thin
- Date: Thu, 20 Sep 2001 16:38:28 -0600
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--
>