[Thread Prev][Thread Next][Index]

[ferret_users] Unable to write to the next vertical level

Hi Ferreters,

   I have been trying to work with 8 ensemble outputs, which I have only at the surface level.
I added a z-axis to the data using a pseudo-dataset which has existing vertical levels.
Now my problem is though the vertical axis is created andI couldi write at k=1, but 
I could not write data from k=2:8. 
   I am attaching my script. Does any of you have some suggestions?

use Ensemble_1.cdf  ! Ensemble run 1
use test.cdf        ! A dataset which has z-axis (znum)

! Now add the z-axis (znum) to the existing grid

def grid/x=I/y=J263_289/z=znum/t=TIME tgrid  ! I, J263_289,TIME are axes names in Ensemble_1.cdf

!Create a new variable having a zaxis

let sst = temp[d=1,g=tgrid]+0*z[gz=znum,z=1]
set var/title="Ensembles"/units=ensemble sst

!save the new variable as k=1 in a file having 8 vertical levels
save/file=Outfile.cdf/KLIMITS=1:8/k=1/clobber sst
cancel data/all

! Similarly with the 2nd ensemble member
use Ensemble_2.cdf
let sst = temp[g=tgrid]+0*z[gz=znum,z=2]
save/append/file=Outfile.cdf/k=2 sst

!Though the Outfile.cdf has 8 vertical levels yet it has data only in k=1,
!but do not write the data of Ensemble_2.cdf in k=2

Thanking you,
Graduate student
Department of Oceanography
Physical Oceanography
Texas A&M University

[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce / NOAA / OAR / PMEL / Ferret

Privacy Policy | Disclaimer | Accessibility Statement