I can draw a netCDF file in GrADS

Hi,

Thank you for previous help for solving NCL problem.
But I have a problem again.
I made a netCDF file interpolated from fixed grid to gaussian grid.

*************************************************
; read and calcualte of multiple files
;*************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
;*************************************************
begin
;*************************************************
; open file and read in data
;*************************************************
  filename = "stream.daily.grid.1978.nc"
  fu = addfile (filename,"r")
  sf = fu->sf(:,9,:,:)
  time = fu->time
  lat = fu->lat
  lon = fu->lon
  dimu = getfilevardimsizes(fu,"sf")
  nt = dimu(0)
  nlev = dimu(1)
  nlat = dimu(2)
  nlon = dimu(3)
;************************************************
;interpolates a scalar quantity from a fixed grid to a Gaussian grid
;************************************************
  jlat=32
  ilon = 64
  sfT21 = f2gsh_Wrap (sf, (/jlat,ilon/), 21)
  lonT21 = lonGlobeF (ilon, "lon", "longitude", "degrees_east")
  latT21 = latGau (jlat, "lat", "latitude" , "degrees_north")
;************************************************
; outfile of interpolated stream function
;************************************************
  fout = addfile("stream.daily.trun.1978.nc","c")
  setfileoption(fout,"DefineMode",True)
  fileAtt = True
  fileAtt_at_creation_date = systemfunc("date")
  fileAtt_at_title = "efficient write netcdf data "
  fileattdef(fout,fileAtt)
  dimNames = (/"time","lat","lon"/)
  dimSizes = (/nt,jlat,ilon/)
  dimUnlim = (/False,False, False/)
  filedimdef(fout,dimNames,dimSizes,dimUnlim)
  filevardef (fout,"time",typeof(time),getvardims(time))
  filevardef (fout,"lat" ,typeof(lat) ,getvardims(latT21))
  filevardef (fout,"lon" ,typeof(lon) ,getvardims(lonT21))
  filevardef (fout,"sfT21" ,typeof(sfT21) ,getvardims(sfT21))
  filevarattdef(fout,"time",time)
  filevarattdef(fout,"lat",latT21)
  filevarattdef(fout,"lon",lonT21)
  filevarattdef(fout,"sfT21",sfT21)
  fout->time = (/time/)
  fout->latT21 = (/latT21/)
  fout->lonT21 = (/lonT21/)
  fout->sfT21 = (/sfT21/)
end

but I have a error message during drawing in GrADS.
=====================================
sdfopen stream.daily.trun.1978.nc
Scanning self-describing file: stream.daily.trun.1978.nc
SDF file stream.daily.trun.1978.nc is open as file 1
LON set to 0 360
LAT set to 9.96921e+36 9.96921e+36
LEV set to 0 0
Time values set: 1978:1:1:0 1978:1:1:0
ga-> d sft21
Map Projection Error: Internal Logic check
Cannot contour grid - all undefined values
==============================
I don't know why?
Thank you for help in advance!

Hee-Jeong
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jul 26 2006 - 18:26:14 MDT