display two times

[Stegert]

Dear Rafael,

what is the time-resolution?
you say: 'each files cotain a one hour of my run...'
The ctl names look as if you have 1 timestep in each
ctl (i.e. 5 timesteps alltogether), is it right?

as a possibility I'd suggest, you define
a dummy ctl/grads file with x/y/z resolution  = 1
and the sum of timesteps (here 5) for one variable
(say 'var') with all values set to '1'.
Then, recursively collect the values from your ctl's and
set them on the new variable:

open var.ctl (will be dfile 1)
open wrfnmm_d01_00.ctl (dfile 2)
open wrfnmm_d01_06.ctl (dfile 3)
open wrfnmm_d01_12.ctl (dfile 4)
open wrfnmm_d01_18.ctl (dfile 5)
open wrfnmm_d01_24.ctl (dfile 6)

itime0 = 1
itime1 = 5 (your end time)

while (itime0 <= itime1)
*first take the value from your ctl
   'set dfile 'itime0+1
   set lat 9.5
   set lon -66.5
   set lev 850
   set t 1
   'd 'apcpsfc
   vartmp = subwrd(result,4)
*then put it on 'var'
   'set dfile 1'
   'set x 1'
   'set y 1'
   'set z 1'
   'set t 'itime0
   'set defval var 1 1 ' vartmp
   itime0 = itime0+1
endwhile

*afterwards you should be able to display the new variable
'set dfile 1'
'set x 1'
'set y 1'
'set z 1'
'set t 1 5'
'd var'

Implement this into your meteogram skript (I didn't have
a closer look on that...)
You can use 'var' for each other varialbe of your ctl, of course.

Hope this works for you.
Christoph

Rafael Mundaray schrieb:
> Hi my dears GrADS USERS, i have files of WRF NMM in format .ctl, then i
> can read these with GrADS, but each files cotain a one hour of my run...
> For example:
> wrfnmm_d01_00.ctl
> wrfnmm_d01_06.ctl
> wrfnmm_d01_12.ctl
> wrfnmm_d01_18.ctl
> wrfnmm_d01_24.ctl
> Then if i want to make a meteorogram for example of the variable
> PRECIPITATION (in my .ctl it is named apcpsfc), then i make this
> open wrfnmm_d01_00.ctl
> open wrfnmm_d01_06.ctl
> open wrfnmm_d01_12.ctl
> open wrfnmm_d01_18.ctl
> open wrfnmm_d01_24.ctl
> set lat 9.5
> set lon -66.5
> set lev 850
> set t (this i don't know as put it)
> d apcpsfc (this i don't know as put it, too)
>
> I want a meteorogram with this five time of my five files .ctl... How do
> i make  this???
>
> Please... Thanks...
>
> I add my script for make meteorogram...
> You can see it...
>
> --
> Ing. Rafael E. Mundaray M.

--
Christoph Stegert
Institute for Oceanography - Ecological Modelling group
University of Hamburg, GER - (ZMK/ZMAW alliance member)
Bundesstr.53 20146 Hamburg - +49-40/42838-7486 Room 348