linint1_Wrap

From: Tummuri, Spandana <spandana.tummuri_at_nyahnyahspammersnyahnyah>
Date: Tue, 13 Nov 2007 23:51:44 -0600

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