Re: about the output data

From: Saji Njarackalazhikam Hameed <saji_at_nyahnyahspammersnyahnyah>
Date: Wed, 15 Mar 2006 21:39:20 -1000

Alice,

The Netcdf file may not have written properly. The variable "airsea"
seems to be of type short, but is written out as float. I guess, you can
read the data once again in NCL, "unpack" it and then write it out
properly. Check out the following code for an example.

; ----------------------------------------------
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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"

procedure add_attributes(var,att_array)
begin
 ndims=dimsizes(att_array)
  do i = 0,ndims-1
   var!i=att_array(i)
  end do
end

begin
fin=addfile("air.season.mean5.nc","r")
asea=fin->airsea
ndims=dimsizes(asea)
asea2=new( ndims, "float")
asea2=asea*asea_at_scale_factor+asea_at_add_offset
add_attributes(asea2,(/"time","level","lat","lon"/))
asea2&time=asea&time
asea2&level=asea&level
asea2&lat=asea&lat
asea2&lon=asea&lon

fout=addfile("test.nc","c")
fout->airsea=asea2
end

; -------------------------------------------

* ?$BId?(B???? <fjl_at_pku.edu.cn> [2006-03-16 14:53:55 +0800]:

> Hi$B!'(B
> Attached is the output data by ncl, when I run Gradsnc to plot a latitude-height section figure by using the data, it can't complished. And the error information given by Gradsnc is that:
> SDF data variable has unsupported packed data type.
> Data Request Error: Error for variable 'airsea'
> Error ocurred at column 1
> DISPLAY error: Invalid expression
> Expression = airsea

> Can somebody tell me the reason. Thanks a lot!

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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"

procedure add_attributes(var,att_array)
begin
 ndims=dimsizes(att_array)
  do i = 0,ndims-1
   var!i=att_array(i)
  end do
end

begin
fin=addfile("air.season.mean5.nc","r")
asea=fin->airsea
ndims=dimsizes(asea)
print(ndims)
printVarSummary(asea)
asea2=new( ndims, "float")
asea2=asea*asea_at_scale_factor+asea_at_add_offset
add_attributes(asea2,(/"time","level","lat","lon"/))
asea2&time=asea&time
asea2&level=asea&level
asea2&lat=asea&lat
asea2&lon=asea&lon

fout=addfile("test.nc","c")
fout->airsea=asea2
end
> Best Wishes!

> Alice
> $B!!!!!!!!!!!!!!!!(Bfjl_at_pku.edu.cn
> $B!!!!!!!!!!!!!!!!!!!!(B2006-03-16

> mutt_octet_view: file type is unsupported

> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
Saji N. Hameed,
Assistant Researcher,			Tel:808 9569534
International Pacific Research Center,	email:saji_at_hawaii.edu
2525 Correa Rd, Honolulu, HI 96822	http://iprc.soest.hawaii.edu/~saji
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Mar 16 2006 - 00:39:20 MST

This archive was generated by hypermail 2.2.0 : Wed Mar 29 2006 - 11:53:03 MST