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-talkReceived 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