[gradsusr] FW: ensemble windchill

[Jeffrey Duda]

Eric,
What exactly is not going right in your images that you need help with?

Jeff Duda

On Wed, Dec 15, 2010 at 5:08 AM, eric delaere <karimhamid2 at hotmail.com>wrote:

>  hello
> I still don't found a solution to fix this problem.
> I'm able to make ensembles with some simple parameters (see
> http://users.fulladsl.be/spb4195/gfs.htm) but with parameters, derived
> from formula (like the wind chill), it doesn't work...
> Anybody?
>
> ------------------------------
> From: karimhamid2 at hotmail.com
> To: gradsusr at gradsusr.org
> Date: Fri, 10 Dec 2010 14:05:07 +0100
> Subject: [gradsusr] ensemble windchill
>
> Hi
>
> I have problems with printing the ensemble meteogram for the wind chill
> (with  whiskers & bars). No problem to make a ensemble for eg. T850 or so,
> but to calculate the wind chill, I need to use some formula & also some
> maskout & const stuff...
> Therefore , the min, max,..etc aren't correct in the tloop formula...
>
> Obviously I do something wrong in my script philosophy but what??
> Thx
>
> (I use grads 2.0)
> ***********
> function main( args )
> 'reinit'
> 'sdfopen
> http://nomads.ncep.noaa.gov:9090/dods/gens_bc/gens'datum''mon''day'/gep_all_'run'z<http://nomads.ncep.noaa.gov:9090/dods/gens_bc/gens%27datum%27%27mon%27%27day%27/gep_all_%27run%27z>
> '
>
> * Calculate the ensemble mean
> * ---------------------------
>   'set lon 4'
>   'set lat 51'
>   'set t 1 37'
>
> 'define wind = (mag(ugrd10m,vgrd10m))'
> 'define windH = const(const(maskout(wind,wind-1),1),0,-u)'
> 'define windL = const(const(maskout(wind,0.999-wind),1),0,-u)'
> 'define W =  wind*2.2369363'
> 'define T =  (9/5)*(tmp2m-273.15)+32'
> 'define chillf=
> 35.74+(0.6215*T)-(35.75*pow(W,0.16))+(0.4275*T*pow(W,0.16))'
> 'define chillC =  (chillf-32)*(5/9)'
> 'define chill2 = (windL*(tmp2m-273.15))+(windH*chillC)'
>
> 'define ensmean=ave(chill2,e=1,e=21)'
>
>  'set cmark 0'
>  'set vrange -20 10'
>  'set ylint 5'
>  'set cstyle 2'
>  'set ccolor 4'
>  'set cthick 3'
>  'display ensmean'
>
> * Calculate the variance
> * ----------------------
>   diffsq = 'pow(chill2-ensmean,2)'
>
> variance = 'ave('diffsq',e=1,e=21)'
> 'define stddev=sqrt('variance')'
>
> * Calculate the min/max
> * ---------------------
>  'define ensmin=tloop(min(chill2,e=1,e=21))'
>
>  'set cmark 0'
>  'set vrange -20 10'
>  'set ylint 5'
>  'set cstyle 1'
>  'set ccolor 4'
>  'set cthick 3'
>  'd ensmin'
>
>   'define ensmax=tloop(max(chill2,e=1,e=21))'
>
>  'set cmark 0'
>  'set vrange -20 10'
>  'set ylint 5'
>  'set cstyle 1'
>  'set ccolor 2'
>  'set cthick 3'
>  'd ensmax'
>
> * Plot the results
> * ----------------
>   'set t 1 37'
>   'set lon 4'
>   'set lat 51'
>   'set vrange -20 10'
>   'set ylint 5'
>
> * Draw error bars for min/max
> * ---------------------------
>   'set gxout errbar'
>   'set bargap 70'
>   'set rgb 16 183 183 183'
>   'set ccolor 16'
>   'set cthick 3'
>   'd ensmin;ensmax'
>
> * Draw bars for +/- standard deviation
> * ------------------------------------
>   plus  = '(ensmean+stddev)'
>   minus = '(ensmean-stddev)'
>
>   'set gxout bar'
>   'set bargap 50'
>   'set baropts filled'
>   'set rgb 16 183 183 183'
>   'set ccolor 16'
>   'd 'minus';'plus
>
> * Draw line for Ensemble mean
> * ---------------------------
>   'set gxout line'
>   'set cmark 0'
>   'set cthick 6'
>   'set digsiz 0.05'
>   'set ccolor 2'
>   'd ensmean'
>
> 'printim c:\model\upload\test\ensemble.png'
>
> *clear
> 'quit'
>
> _______________________________________________ gradsusr mailing list
> gradsusr at gradsusr.org http://gradsusr.org/mailman/listinfo/gradsusr
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>


-- 
Jeff Duda
Iowa State University
Meteorology Graduate Student
3134 Agronomy Hall
www.meteor.iastate.edu/~jdduda
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101215/9c44719d/attachment-0003.html