question about wrf_user_intrp3d, number of vertical layers

From: xiaoming Hu <yuanfangcan_at_nyahnyahspammersnyahnyah>
Date: Thu Aug 19 2010 - 11:28:35 MDT

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