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