Re: help, the problem in drawing MM5 data

> I have writen my MM5 data to the format that can be read by function
> ,nclbinrecrd(). It can work smoothly.
> However, I want to draw the data with a map (lambert), and the lat and lon
> are known at every grid. Anyone can give me some suggestion how to do this?
----------------------------------------------------------

Do you know the necessary map information?

-----------------------------------------------------------

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin
   nlat = ...
   mlon = ...
   fName = "....."
   x = fbinrecread(fName, (/nlat,mlon/), "float" )
;************************************************
; create plots
;************************************************
   wks = gsn_open_wks("ps" ,"MM5_lc") ; ps,pdf,x11,ncgm,eps
   res = True ; plot mods desired
;;res@gsnMaximize = True ; uncomment to maximize size
   res@gsnAddCyclic = False ; regional data: not cyclic
   res@tfDoNDCOverlay = True ; set True for native projection

   res@mpProjection = "LambertConformal"
   res@mpLambertParallel1F = ............
   res@mpLambertParallel2F = ............
   res@mpLambertMeridianF = ............
   res@mpLimitMode = "Corners"
   res@mpLeftCornerLatF = ..........
   res@mpLeftCornerLonF = ..........
   res@mpRightCornerLatF = ..........
   res@mpRightCornerLonF = .......... )
   res@mpOutlineDrawOrder = "PostDraw" ; draw continental outline last
  ;res@mpOutlineBoundarySets = "GeophysicalAndUSStates" ; state boundaries
  ;res@mpFillOn = False ; turn off map fill
  ;res@mpFillDrawOrder = "PreDraw"

   res@tiMainString = "your title"
   plot = gsn_csm_contour_map(wks, x, res)
end
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Mar 20 2006 - 19:17:27 MST