> 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_at_gsnMaximize = True ; uncomment to maximize size
res_at_gsnAddCyclic = False ; regional data: not cyclic
res_at_tfDoNDCOverlay = True ; set True for native projection
res_at_mpProjection = "LambertConformal"
res_at_mpLambertParallel1F = ............
res_at_mpLambertParallel2F = ............
res_at_mpLambertMeridianF = ............
res_at_mpLimitMode = "Corners"
res_at_mpLeftCornerLatF = ..........
res_at_mpLeftCornerLonF = ..........
res_at_mpRightCornerLatF = ..........
res_at_mpRightCornerLonF = .......... )
res_at_mpOutlineDrawOrder = "PostDraw" ; draw continental outline last
;res_at_mpOutlineBoundarySets = "GeophysicalAndUSStates" ; state boundaries
;res_at_mpFillOn = False ; turn off map fill
;res_at_mpFillDrawOrder = "PreDraw"
res_at_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
This archive was generated by hypermail 2.2.0 : Wed Mar 29 2006 - 11:53:04 MST