Eric,<br>What exactly is not going right in your images that you need help with?<br><br>Jeff Duda<br><br><div class="gmail_quote">On Wed, Dec 15, 2010 at 5:08 AM, eric delaere <span dir="ltr"><<a href="mailto:karimhamid2@hotmail.com">karimhamid2@hotmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div>
hello<br>I still don't found a solution to fix this problem.<br>I'm able to make ensembles with some simple parameters (see <a href="http://users.fulladsl.be/spb4195/gfs.htm" target="_blank">http://users.fulladsl.be/spb4195/gfs.htm</a>) but with parameters, derived from formula (like the wind chill), it doesn't work...<br>
Anybody?<br><br><hr>From: <a href="mailto:karimhamid2@hotmail.com" target="_blank">karimhamid2@hotmail.com</a><br>To: <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>Date: Fri, 10 Dec 2010 14:05:07 +0100<br>
Subject: [gradsusr] ensemble windchill<br><br>
Hi<br><br>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...<br>
Therefore , the min, max,..etc aren't correct in the tloop formula...<br><br>Obviously I do something wrong in my script philosophy but what??<br>Thx<br><br>(I use grads 2.0)<br>***********<br>function main( args )<br>
'reinit'<br>'sdfopen <a href="http://nomads.ncep.noaa.gov:9090/dods/gens_bc/gens%27datum%27%27mon%27%27day%27/gep_all_%27run%27z" target="_blank">http://nomads.ncep.noaa.gov:9090/dods/gens_bc/gens'datum''mon''day'/gep_all_'run'z</a>'<br>
<br>* Calculate the ensemble mean <br>* ---------------------------<br> 'set lon 4'<br> 'set lat 51'<br> 'set t 1 37'<br><br>'define wind = (mag(ugrd10m,vgrd10m))'<br>'define windH = const(const(maskout(wind,wind-1),1),0,-u)'<br>
'define windL = const(const(maskout(wind,0.999-wind),1),0,-u)'<br>'define W = wind*2.2369363'<br>'define T = (9/5)*(tmp2m-273.15)+32'<br>'define chillf= 35.74+(0.6215*T)-(35.75*pow(W,0.16))+(0.4275*T*pow(W,0.16))'<br>
'define chillC = (chillf-32)*(5/9)'<br>'define chill2 = (windL*(tmp2m-273.15))+(windH*chillC)'<br><br>'define ensmean=ave(chill2,e=1,e=21)'<br><br> 'set cmark 0'<br> 'set vrange -20 10'<br>
'set ylint 5'<br> 'set cstyle 2'<br> 'set ccolor 4'<br> 'set cthick 3'<br> 'display ensmean'<br> <br>* Calculate the variance<br>* ----------------------<br> diffsq = 'pow(chill2-ensmean,2)'<br>
<br>variance = 'ave('diffsq',e=1,e=21)'<br>'define stddev=sqrt('variance')'<br><br>* Calculate the min/max<br>* ---------------------<br> 'define ensmin=tloop(min(chill2,e=1,e=21))'<br>
<br> 'set cmark 0'<br> 'set vrange -20 10'<br> 'set ylint 5'<br> 'set cstyle 1'<br> 'set ccolor 4'<br> 'set cthick 3'<br> 'd ensmin'<br><br> 'define ensmax=tloop(max(chill2,e=1,e=21))'<br>
<br> 'set cmark 0'<br> 'set vrange -20 10'<br> 'set ylint 5'<br> 'set cstyle 1'<br> 'set ccolor 2'<br> 'set cthick 3'<br> 'd ensmax'<br> <br>* Plot the results<br>* ----------------<br>
'set t 1 37'<br> 'set lon 4'<br> 'set lat 51'<br> 'set vrange -20 10'<br> 'set ylint 5'<br><br>* Draw error bars for min/max<br>* ---------------------------<br> 'set gxout errbar'<br>
'set bargap 70'<br> 'set rgb 16 183 183 183'<br> 'set ccolor 16'<br> 'set cthick 3'<br> 'd ensmin;ensmax'<br><br>* Draw bars for +/- standard deviation<br>* ------------------------------------<br>
plus = '(ensmean+stddev)'<br> minus = '(ensmean-stddev)'<br><br> 'set gxout bar'<br> 'set bargap 50'<br> 'set baropts filled'<br> 'set rgb 16 183 183 183'<br> 'set ccolor 16'<br>
'd 'minus';'plus<br><br>* Draw line for Ensemble mean <br>* ---------------------------<br> 'set gxout line'<br> 'set cmark 0'<br> 'set cthick 6'<br> 'set digsiz 0.05'<br>
'set ccolor 2'<br> 'd ensmean'<br><br>'printim c:\model\upload\test\ensemble.png'<br><br>*clear<br>'quit'<br>
<br>_______________________________________________
gradsusr mailing list
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a> </div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>Jeff Duda<br>Iowa State University<br>Meteorology Graduate Student<br>3134 Agronomy Hall<br><a href="http://www.meteor.iastate.edu/~jdduda">www.meteor.iastate.edu/~jdduda</a><br>