Hello,
I wrote the following program to interpolate levels. I used linint1_Wrap but it still does not retain Metadata. Is there any error in my program?
Thanks in advance.
Spandana Tummuri
lat = in1->lat
lon = in1->lon
time = in1->time
levsoi = in1->levsoi
TSOI = in1->TSOI
nlat = dimsizes(lat)
mlon = dimsizes(lon)
ntime = dimsizes(time)
nlevsoi = dimsizes(levsoi)
lNew = (/ 10, 40, 100, 200 /)
lNew_at_units = "centimeters"
lNew!0 = "levsoi"
x = linint1_Wrap (levsoi, TSOI(time|:,lat|:,lon|:,levsoi|:), False, lNew, 0)
x!0 = "time"
x!1 = "lat"
x!2 = "lon"
x!3 = "levsoi"
x_at_units = "K"
filo="/home/stummuri/CCSM-INPUTS/s01-int.nc"
fo=addfile(filo,"c")
fo->TSOI = x
fo->lat = lat
fo->lon = lon
fo->time = time
fo->levsoi = (/ lNew /)
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Nov 13 2007 - 22:51:44 MST
This archive was generated by hypermail 2.2.0 : Thu Nov 15 2007 - 10:01:54 MST