Thanks Stegert... I am will see this...<br><br><div><span class="gmail_quote">2007/6/1, Stegert <<a href="mailto:stegert@ifm.uni-hamburg.de">stegert@ifm.uni-hamburg.de</a>>:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Dear Rafael,<br><br>what is the time-resolution?<br>you say: 'each files cotain a one hour of my run...'<br>The ctl names look as if you have 1 timestep in each<br>ctl (i.e. 5 timesteps alltogether), is it right?<br>
<br>as a possibility I'd suggest, you define<br>a dummy ctl/grads file with x/y/z resolution = 1<br>and the sum of timesteps (here 5) for one variable<br>(say 'var') with all values set to '1'.<br>Then, recursively collect the values from your ctl's and
<br>set them on the new variable:<br><br>open var.ctl (will be dfile 1)<br>open wrfnmm_d01_00.ctl (dfile 2)<br>open wrfnmm_d01_06.ctl (dfile 3)<br>open wrfnmm_d01_12.ctl (dfile 4)<br>open wrfnmm_d01_18.ctl (dfile 5)<br>open wrfnmm_d01_24.ctl (dfile 6)
<br><br>itime0 = 1<br>itime1 = 5 (your end time)<br><br>while (itime0 <= itime1)<br>*first take the value from your ctl<br> 'set dfile 'itime0+1<br> set lat 9.5<br> set lon -66.5<br> set lev 850<br> set t 1
<br> 'd 'apcpsfc<br> vartmp = subwrd(result,4)<br>*then put it on 'var'<br> 'set dfile 1'<br> 'set x 1'<br> 'set y 1'<br> 'set z 1'<br> 'set t 'itime0<br>
'set defval var 1 1 ' vartmp<br> itime0 = itime0+1<br>endwhile<br><br>*afterwards you should be able to display the new variable<br>'set dfile 1'<br>'set x 1'<br>'set y 1'<br>'set z 1'
<br>'set t 1 5'<br>'d var'<br><br>Implement this into your meteogram skript (I didn't have<br>a closer look on that...)<br>You can use 'var' for each other varialbe of your ctl, of course.<br><br>
Hope this works for you.<br>Christoph<br><br>Rafael Mundaray schrieb:<br>> Hi my dears GrADS USERS, i have files of WRF NMM in format .ctl, then i<br>> can read these with GrADS, but each files cotain a one hour of my run...
<br>> For example:<br>> wrfnmm_d01_00.ctl<br>> wrfnmm_d01_06.ctl<br>> wrfnmm_d01_12.ctl<br>> wrfnmm_d01_18.ctl<br>> wrfnmm_d01_24.ctl<br>> Then if i want to make a meteorogram for example of the variable
<br>> PRECIPITATION (in my .ctl it is named apcpsfc), then i make this<br>> open wrfnmm_d01_00.ctl<br>> open wrfnmm_d01_06.ctl<br>> open wrfnmm_d01_12.ctl<br>> open wrfnmm_d01_18.ctl<br>> open wrfnmm_d01_24.ctl
<br>> set lat 9.5<br>> set lon -66.5<br>> set lev 850<br>> set t (this i don't know as put it)<br>> d apcpsfc (this i don't know as put it, too)<br>><br>> I want a meteorogram with this five time of my five files .ctl... How do
<br>> i make this???<br>><br>> Please... Thanks...<br>><br>> I add my script for make meteorogram...<br>> You can see it...<br>><br>> --<br>> Ing. Rafael E. Mundaray M.<br><br>--<br>Christoph Stegert
<br>Institute for Oceanography - Ecological Modelling group<br>University of Hamburg, GER - (ZMK/ZMAW alliance member)<br>Bundesstr.53 20146 Hamburg - +49-40/42838-7486 Room 348<br></blockquote></div><br><br clear="all"><br>
-- <br>Ing. Rafael E. Mundaray M.