I am trying to convert my NCL routine to plot 3-D variables from a
netcdf file as a zonal average. I have tried to use the routine
"dim_avg" to do this.
When I try to run the routine I get the following
error:
Variable: u
Type: float
Total Size: 1437696 bytes
359424 values
Number of Dimensions: 4
Dimensions and sizes: [time | 1] x [lev_p | 26] x [lat | 96] x [lon
| 144]
Coordinates:
time: [ 95.. 95]
lev_p: [91500..0.03]
lat: [ -90..89.99999999999999]
lon: [ 0..357.5]
Number Of Attributes: 4
units : m/s
long_name : Zonal wind
cell_method : time: mean
_FillValue : -9999
(0) gsn_csm_pres_hgt: Fatal: The first dimension of the input
data must
(0) have a coordinate variable called 'lev.'
(0) Cannot create plot.
fatal:Illegal right-hand side type for assignment
fatal:Execute: Error occurred at or near line 82 in file zonalav.ncl
This is the routine I am trying to modify to produce
the zonal average:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
; input file
diri = "./"
fili = "cosfoura.cam2.h0.0000-12-05-00000.nc"
fi = addfile (diri+fili, "r")
; output file
diro = "./"
filo = "helen.nc"
system ("/bin/rm -f "+diro+filo) ; remove any pre-exist file
fo = addfile (diro+filo, "c")
; add any file attributes
fo_at_title = fi_at_title + ": Selected Variables at pressure levels"
fo_at_history = systemfunc ("date")
fo_at_source = fi_at_source
fo_at_case = fi_at_case
fo_at_Conventions = fi_at_Conventions
Var = (/ "T" , "Q", "U", "V"/) ; select variable to be
interpolated.
nVar = dimsizes (Var)
; desired output levels
;lev_p = (/ 91500., 90000., 87500., 80000., 67500., 50000., 30000.,
20000., 10000., 5000., 3000., 1500., 750., 390., 200., 100., 50.,25.,
20., 15.,10., 8., 5.,3.,1. /)
;lev_p = (/ 91500., 90000., 87500., 80000., 67500., 50000., 30000.,
20000., 10000., 5000., 3000., 1500., 750., 390., 200., 100., 50.,25.,
14., 7., 3.5, 1.9, 0.95, 0.4, 0.1, 0.01 /)
lev_p = (/ 91500., 90000., 87500., 80000., 67500., 50000., 30000.,
20000., 10000., 5000., 3000., 1500., 750., 390., 200., 100., 50.,25.,
14., 7., 3.5, 1.9, 0.95, 0.4, 0.1, 0.03 /)
;lev_p = (/ 91500., 90000., 87500., 80000., 67500., 50000., 30000.,
20000., 10000., 5000., 3000., 1500., 750., 390., 200., 100., 50.,25.,
14., 7., 3.5, 2.0 /)
lev_p!0 = "lev_p" ; variable and dimension name
the same
lev_p&lev_p = lev_p ; create coordinate variable
lev_p_at_long_name = "pressure" ; attach some attributes
lev_p_at_units = "hPa"
lev_p_at_positive = "down"
hyam = fi->hyam ; read hybrid info
hybm = fi->hybm
PS = fi->PS
P0mb = 0.01*fi->P0
do n=0,nVar-1 ; loop over the variables
X = fi->$Var(n)$
Xp = vinth2p (X, hyam, hybm, lev_p ,PS, 1, P0mb, 2, True)
copy_VarAtts(X, Xp)
fo->$Var(n)$ = Xp ; write to netCDF file
print (Var(n)+": interpolated and written to netCDF")
end do
end
;***********************
; zonal.ncl
;***********************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;***********************
begin
f = addfile ("helen.nc","r")
print(f)
u = f->U
printVarSummary( u )
dimu = dimsizes( u )
ntim = dimu(0)
klvl = dimu(1)
nlon = dimu(2)
mlon = dimu(3)
uavg = dim_avg(u)
;***********************
; Create Plot
;***********************
wks = gsn_open_wks ("pdf", "cosfoura001205lon40" ) ;
open workstation
gsn_define_colormap(wks,"rainbow") ; choose colormap
res = True
res_at_cnFillOn = True
res_at_lbLabelAutoStride = True
res_at_gsnMaximize = True ; if [ps, eps, pdf] make large
res_at_gsnSpreadColors = True ; span color map
do nt=0,ntim-1
res_at_gsnCenterString = "cosfoura001205lon40"
; plot = gsn_csm_pres_hgt(wks, u(nt,:,:,40), res ) ; (lev,lat)
plot = gsn_csm_pres_hgt(wks, uavg(nt,:,:), res ) ; (lev,lat)
res_at_trYReverse = True
; plot = gsn_csm_contour (wks, u(nt,:,:,0), res ) ; (lev,lat)
; plot = gsn_csm_contour (wks, u(nt,:,0,:), res ) ; (lev,lon)
kl = 5
res_at_mpFillOn = False
res_at_mpGridAndLimbOn = True
res_at_mpGridLineDashPattern = 2
res_at_mpOutlineBoundarySets = "NoBoundaries"
res_at_mpCenterLonF = 180.
res_at_gsnCenterString = res_at_gsnCenterString+" p="+u&lev_p(kl)
; plot = gsn_csm_contour_map_ce (wks, u(nt,kl,:,:), res )
end do
end
Any idea what the problem is?
Thanks,
Helen.
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Aug 28 2008 - 17:30:24 MDT
This archive was generated by hypermail 2.2.0 : Fri Aug 29 2008 - 07:39:21 MDT