Re: global coordinate for wrf regional modeling ?

The WRF_contributed.ncl library is designed to
handle files created by the WRF model.

You are reading two non-WRF files.

You are reading two

Yan Bao wrote:
> 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 <http://00.nc>",
> "r")
> fmap = addfile ("geo_em.d01.nc <http://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
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Apr 22 2009 - 20:17:01 MDT