Hello
I would like to extract a N-S cross section.
The number of layers of WRF output is 49. Why the extracted N-S cross section has 100 vertical layers? What is the vertical coordinate of the extracted cross section?
I pasted my code below:
pvo = wrf_user_getvar(f,"pvo",0)
printVarSummary(pvo)
z = wrf_user_getvar(f, "z",0) ; grid point height
printVarSummary(z)
opt = True
opt@returnInt = False ; Return real values (set to True for integer values - integer is default)
; loc_barrow = wrf_user_ll_to_ij (f, -156.788, 71.290, opt)
loc_barrow = wrf_user_ll_to_ij (f, -156.788, 60, opt)
loc_polar = wrf_user_ll_to_ij (f, -156.788, 90, opt)
print(loc_barrow)
print(loc_polar)
dimspvo = dimsizes(pvo)
plane = (/ loc_barrow(0),loc_barrow(1), loc_polar(0),loc_polar(1) /) ; pivot point is center of domain
pvo_plane1 = wrf_user_intrp3d(pvo,z,"v",plane,0.0,True)
printVarSummary(pvo_plane1)
Thanks
Xiaoming
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Aug 19 11:28:44 2010
This archive was generated by hypermail 2.1.8 : Mon Sep 20 2010 - 15:45:07 MDT