[Thread Prev][Thread Next][Index]
padding a data set
Hi,
A simple question but one that's stumping me. I have data on a lat/lev
grid where lat goes from -2.5:90:2.5spacing (38 pts). For an existing
code, I need to put this onto a -90:90:2.5spacing grid (73 points). The
numbers that are stored in the -90:-5 section (the new section) don't
matter, i just need the data to have the proper dimensions. I tried to
define a new global axis and regrid to that.
yes? define axis/y=90s:90n:2.5 newy
yes? let ztend=hgttend[d=1,gy=newy@asn]*0 !pad the new section with zeroes
yes? save/file=test.nc/jlimits=1:73/j=1:35 ztend !write the new section to
a file
!this part works fine
yes? let ztend=hgttend[d=1,gy=newy@asn] !now define the real data
yes? save/file=test.nc/append/i=1/l=1/jlimits=1:73/j=36:73 ztend
This approach doesn't work though. I understand why it doesn't work but I
don't know how to shift the axis so that the proper data is written for
the block -2.5:90. I need to @asn the good data to the last 38 points on
the new axis but how?
Any help is appreciated!
Cheers,
Brent
--
Brent A. McDaniel
Dept of Earth and Atmospheric Sciences
Georgia Institute of Technology
Atlanta, Ga. USA
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement