coordinate problem

From: JinQ Liu <liujq_0912_at_nyahnyahspammersnyahnyah>
Date: Wed Jan 11 2012 - 20:18:57 MST

Hi,
I have add the "degrees_north" attributes on the variable, but still have this error :
       fatal:Dimension (west_east) of (u_plane1) does not have an associated coordinate variable

Any suggestions will be appreciate.Thanks.
JinQ

--
JinQ,                          M.S. 
Ocean University of China
--
 
ftp> put wrfoutd01_origin
scripts:
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/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl" ; WRF_Times2Udunits_c
;**************************************************
begin
   xp = new((/3,6/),float)            ; allocate memory
   yp = new((/3,6/),float)
   yp(0,:) = (/ 23.2, 24.3, 25.6, 26.8, 27.8, 29.0 /)
   xp(0,:) = (/ 119.0,118.7,118.7,118.7,119.2,119.9  /)
   yp(1,:) = (/ 23.382,  24.1339,   25.2679, 26.1043, 26.8851, 27.8487/)
   xp(1,:) = (/ 119.,  118.945,  118.838,   118.786, 118.884, 119.472/)
   yp(2,:) = (/ 23.382,  24.1336, 25.402,  26.4258, 27.421, 28.3855/)
   xp(2,:) = (/ 119. ,  118.975,  118.929,  119.03,   119.192, 119.571/)
res=True
f = addfile("wrfoutd01_origin.nc","r")
loc=wrf_user_ll_to_ij(f,xp(1,:),yp(1,:),res)
print(loc) 
lon1=loc(0,:)+11.11           
lon2=loc(0,:)-11.11
lat1=loc(1,:)+11.11
lat2=loc(1,:)-11.11
afterloc1=wrf_user_ij_to_ll(f,lon1,lat1,res)
afterloc2=wrf_user_ij_to_ll(f,lon1,lat2,res)
afterloc3=wrf_user_ij_to_ll(f,lon2,lat1,res)
afterloc4=wrf_user_ij_to_ll(f,lon2,lat2,res)
 wks = gsn_open_wks("x11" ,"radius")
write_matrix(afterloc1,"6f10.2",False)
;;;;;;;;;;;;;;;;;;;;;;;;;;;
a=afterloc2(1,0)
b=afterloc1(1,0)
c=afterloc2(0,0)
d=afterloc1(0,0)
print(a)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
  times  = wrf_user_list_times(f)  ; get times in the file
  pressure_levels = (/ 850., 200./)   ; pressure levels to plot
  nlevels         = dimsizes(pressure_levels)           ; number of pressure levels
do it=3,4
    print("Working on time: " + times(it) )
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;   
    u  = wrf_user_getvar(f,"ua",it)        ; u averaged to mass points
    v  = wrf_user_getvar(f,"va",it)        ; v averaged to mass points
    p  = wrf_user_getvar(f, "pressure",it) ; pressure is our vertical coordinate
    z  = wrf_user_getvar(f, "z",it)        ; grid point height
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
      pressure0 = pressure_levels(0)
      pressure1 = pressure_levels(1)
      u_plane0  = wrf_user_intrp3d( u,p,"h",pressure0,0.,False)
      v_plane0  = wrf_user_intrp3d( v,p,"h",pressure0,0.,False)
      u_plane1  = wrf_user_intrp3d( u,p,"h",pressure1,0.,False)
      v_plane1  = wrf_user_intrp3d( v,p,"h",pressure1,0.,False)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
lat2d=f->XLAT(0,:,:)
lon2d=f->XLONG(0,:,:)
lat2d@units="degrees_north"
lon2d@units="degrees_east"
u_plane1@lat2d=lat2d
u_plane1@lon2d=lon2d
v_plane1@lat2d=lat2d
v_plane1@lon2d=lon2d
u_plane0@lat2d=lat2d
u_plane0@lon2d=lon2d
v_plane0@lat2d=lat2d
v_plane0@lon2d=lon2d
printVarSummary(u_plane1)
u_planeb=wgt_areaave(u_plane1({a:b},{c:d}),1.,1.,0)-wgt_areaave(u_plane0({a:b},{c:d}),1.,1.,0)
v_planeb=wgt_areaave(v_plane1({a:b},{c:d}),1.,1.,0)-wgt_areaave(v_plane0({a:b},{c:d}),1.,1.,0)
      u_plane = u_plane*1.94386     ; kts
      v_plane = v_plane*1.94386     ; kts
      u_plane@units = "kts"
      v_plane@units = "kts"
print(u_plane)
wind=((u_plane)^2+(v_plane)^2)^0.5
print(wind)
end do
end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jan 11 20:19:12 2012

This archive was generated by hypermail 2.1.8 : Wed Jan 18 2012 - 09:21:55 MST