Regridding 3 dimensional data

[Davide Sacchetti]

at first sight there is a missing "disable fwrite" in the script: you
should put this statement before the script end.
Bye bye
Davide


On Mon, 2007-06-18 at 13:24 -0600, John Guhin wrote:
> Hi,
>      I am trying to regrid a 4 dimensional data  (time,  height, lat,
> lon). I understand that the regrid function works for 2D data. This is
> the script I wrote, but only the first time step is written. Please
> help....
>
> 'reinit'
> 'c'
> *regridding 3D data to 0.5 degree resolution
> 'open u_mon_96to00.ctl'
> 'set gxout fwrite'
> 'set fwrite u_0.5x0.5.dat'
> tt=1
> zz=1
> while(tt<61)
> 'set t 'tt
> while(zz<25)
> 'set z 'zz
> 'set x 1 1110' ; 'set lat -30 45'
> 'd regrid2(u,0.5,0.5,ba_p1_ma000, 0.0,-30.0)'
> zz=zz+1
> endwhile
> tt=tt+1
> endwhile
> 'close 1'
>
> How can I write the data for all time steps?
>
> thanks.
>
> --
> John Guhin
--
Sacchetti Davide
ARPAL UO3 Centro Meteo Idrologico Regione Liguria - Dir. Scientifica
V.le Brigare Partigiane 2 16121 Genova (I)
tel: +39 010 6437535
mail: davide.sacchetti at arpal.org   web: www.meteoliguria.it