fwrite problem of 4D data

[John Guhin]

Thanks Heiner. Now it works.

On Nov 27, 2007 12:06 PM, Heiner Körnich <heiner at misu.su.se> wrote:

> Hi John,
>
> my first reply did not work. But I think you forgot the z1=1 in front of
> the z-loop. After the first loop z was 19, so all other loops were
> skipped.
>
> regards,
> Heiner
> On Tue, 2007-11-27 at 11:54 -0500, John Guhin wrote:
> > Thanks for the quick reply. Your method worked, except I am not sure
> > how the control file will look!
> > So I split the output in two files, one for U and the other for V.
> > The following is just minor modification to your script, may be
> > helpful to others who are following. But, I am still not sure why my
> > original script did not work.
> > ---------------------------------------------------
> > 'reinit'
> > 'c'
> > 'sdfopen uwnd.2002.nc'
> > 'sdfopen vwnd.2002.nc'
> >
> > 'set gxout fwrite'
> > 'set fwrite u_anom_ncep.grd'
> > t1=100
> > while(t1<161)
> > 'set t 't1
> > z1=1
> > while(z1<18)
> > 'set z 'z1
> > 'set x 1 144' ; 'set y 1 73'
> > 'define ubar=ave(uwnd,time=10apr02,time=10jun02)'
> > 'd uwnd-ubar'
> >
> > z1=z1+1
> > endwhile
> > t1=t1+1
> > endwhile
> > 'disable fwrite'
> >
> > 'set gxout fwrite'
> > 'set fwrite v_anom_ncep.grd'
> > t1=100
> > while(t1<161)
> > 'set t 't1
> > z1=1
> > while(z1<18)
> > 'set z 'z1
> > 'set x 1 144' ; 'set y 1 73'
> > 'define vbar=ave(vwnd.2,time=10apr02,time=10jun02)'
> > 'd vwnd.2-vbar'
> > z1=z1+1
> > endwhile
> > t1=t1+1
> > endwhile
> > 'disable fwrite'
> >
> ---------------------------------------------------------------------------
> >
> > On Nov 27, 2007 11:02 AM, Mary Jo Nath < Mary-Jo.Nath at noaa.gov> wrote:
> >         Try the following. (This assumes that the U and V files have
> >         the same time coordinates.) You might also want to post the
> >         associated CTL file for "u_v_anom_ncep.grd"
> >          if you continue to have difficulties. - MJ
> >
> >         'reinit'
> >         'c'
> >
> >         'sdfopen uwnd.2002.nc'
> >         'sdfopen vwnd.2002.nc'
> >
> >         'set gxout fwrite'
> >         'set fwrite u_v_anom_ncep.grd'
> >
> >         t1=100
> >
> >
> >         while(t1<161)
> >
> >         'set t 't1
> >
> >         * write U values at time "t1"  for all Z levels
> >         z1=1
> >         while(z1<18)
> >         'set z 'z1
> >         'set x 1 144' ; 'set y 1 73'
> >         'define ubar=ave(uwnd,time=10apr02,time=10jun02)'
> >
> >         'd uwnd-ubar'
> >         z1=z1+1
> >         endwhile
> >
> >         * Then write V values at time "t1" for all Z levels
> >         z1=1
> >         while(z1<18)
> >         'set z 'z1
> >         'set x 1 144' ; 'set y 1 73'
> >
> >         'define vbar=ave(vwnd.2,time=10apr02,time=10jun02)'
> >
> >         'd vwnd.2-vbar'
> >         z1=z1+1
> >         endwhile
> >
> >         t1=t1+1
> >         endwhile
> >         'disable fwrite'
> >
> >
> >
> >
> >
> >         John Guhin wrote:
> >         > Hi,
> >         >       I am trying to write a 4D  (time, height, lat,
> >         >  lon). This is the script I wrote, but only the first time
> >         step is written.
> >         > Any help?
> >         >
> >         > 'reinit'
> >         > 'c'
> >         >
> >         > 'sdfopen uwnd.2002.nc'
> >         > 'sdfopen vwnd.2002.nc'
> >         >
> >         > 'set gxout fwrite'
> >         > 'set fwrite u_v_anom_ncep.grd'
> >         >
> >         > t1=100
> >         > z1=1
> >         > while(t1<161)
> >         > 'set t 't1
> >         > while(z1<18)
> >         > 'set z 'z1
> >         >
> >         > 'set x 1 144' ; 'set y 1 73'
> >         >
> >         > 'define ubar=ave(uwnd,time=10apr02,time=10jun02)'
> >         > 'define vbar=ave(vwnd.2,time=10apr02,time=10jun02)'
> >         >
> >         > 'd uwnd-ubar'
> >         > 'd vwnd.2-vbar'
> >         >
> >         > z1=z1+1
> >         > endwhile
> >         > t1=t1+1
> >         > endwhile
> >         > 'disable fwrite'
> >         >
> >         >
> >         >  How can I write the data for all time steps?
> >         >
> >         >  thanks.
> >         >
> >         >
> >         >  John Guhin
> >         > -
> >         >
> >         >
> >
> >
> >
> >
> > --
> > John Guhin
> --
> Heiner Körnich
> Dept. of Meteorology            Tel:  +46 8 164333
> Stockholms University, SE-106 91 Stockholm, Sweden
> Email: heiner at misu.su.se   www.misu.su.se/~heiner/<http://www.misu.su.se/%7Eheiner/>
>



-- 
John Guhin
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20071127/ffd38b6b/attachment.html