global coordinate for wrf regional modeling ?

From: Yan Bao <ybao2009_at_nyahnyahspammersnyahnyah>
Date: Wed, 22 Apr 2009 19:01:54 -0700

Hi, Dear all, My script are like the following to check if I can get the
correct met data to drive wrf model. As a result, I got something with
grobal coordinate(-90,90; -180,180) for a regional simulation( Like
California, Lambert projection).

My ncl script is as following:I don't think the pink lines work for the
scripts. The wrf_map_c line always some error like : WRF_resLamCon_c:
unexpected lat/lon rank: rank=1;
Also, whether I chose the True or False ( In this case, I have removed the
wrf_map_c line), the coordinate is global one.
Any suggestions?

      Thanks very much for your help.

      Yan

;*************************************************
; WRF: plot data with "missing_values"
;************************************************
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/wrf/WRF_contributed.ncl"

begin
;************************************************
; open file
; (1) Read SST and SMOIS at specified time/level
; (2) Read character variable Times; Convert to string for plots
;************************************************
  f = addfile ("met_em.d01.1996-01-01_00:00:00.nc", "r")
  fmap = addfile ("geo_em.d01.nc", "r")
  sm = f->SM000010(0,:,:) ; (Time, soil_layers_stag, south_north,
west_east )
  sm@_FillValue = 1.0 ; manually set _FillValue
  sst = f->SKINTEMP(0,:,:) ; (Time, south_north, west_east )
  sst@_FillValue = 0.0 ; manually set _FillValue

  times = chartostring(f->Times) ; convert to type string
  ntim = dimsizes(times) ; # time steps

;************************************************
; create plots
;************************************************
  wks = gsn_open_wks("ps" ,"WRF_test") ; ps,pdf,x11,ncgm,eps
  gsn_define_colormap(wks,"BlAqGrYeOrReVi200") ; select color map

  res = True ; plot mods desired
;;res_at_gsnMaximize = True ; uncomment to maximize size
  res_at_gsnSpreadColors = True ; use full range of colormap
  res_at_cnFillOn = True ; color plot desired
  res_at_cnLinesOn = False ; turn off contour lines
  res_at_cnLineLabelsOn = False ; turn off contour labels
  res_at_lbLabelAutoStride = True ; let NCL figure lb stride
  res_at_lbOrientation = "Vertical" ; default is horizontal
;************************************************
; Use WRF_contributed procedure to set map resources
;************************************************
; WRF_map_c(fmap, res, 0) ; reads info from file

;************************************************
; if appropriate, set True for native (direct) mapping (faster)
;************************************************
  res_at_tfDoNDCOverlay = False
; res_at_tfDoNDCOverlay = True

;************************************************
; associate the 2-dimensional coordinates to the variable for plotting
;************************************************
  if (.not.res_at_tfDoNDCOverlay) then
; lat2d = fmap->XLAT_U(0,50:100,50:100)
 ; lon2d = fmap->XLONG_U(0,50:100,50:100)
      lat2d = fmap->XLAT_U(0,:,:)
     lon2d = fmap->XLONG_U(0,:,:)
      sm_at_lat2d = lat2d
      sm_at_lon2d = lon2d
      sst_at_lat2d = lat2d
      sst_at_lon2d = lon2d
  end if

;************************************************
; For individual plots
;************************************************
  res_at_gsnDraw = False ; do not draw
  res_at_gsnFrame = False ; do not advance 'frame'
;************************************************
; Allocate array to store plots: specify time step
;************************************************
  plts = new (2 , "graphic") ; 1d array to hold plots

;************************************************
; NCL contouring does not like _FillValue=0.0
;************************************************
  sst@_FillValue = -999. ; change to -999

;************************************************
; create panel: panel plots have their own set of resources
; for demo: use Raster fo soil moisture
;************************************************
  resP = True ; modify the panel plot
  resP_at_gsnMaximize = True ; maximize panel area

  res_at_cnFillMode = "RasterFill" ; turn raster on
  plts(0) = gsn_csm_contour_map(wks,sm,res)
  delete(res_at_cnFillMode) ; delete raster mode
  plts(1) = gsn_csm_contour_map(wks,sst,res)

  resP_at_txString = times(0)
  gsn_panel(wks,plts,(/2,1/),resP) ; now draw as one plot
end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Apr 22 2009 - 20:01:54 MDT

This archive was generated by hypermail 2.2.0 : Fri Apr 24 2009 - 14:12:40 MDT