Map shift

Hello,

     I noticed that Andrea and Leilei have posted similar questions
about coastline and map shift but I did not see the final solutions
to their questions. I am having the similar problem. My WRF model
terrain was shifted to the north and west by roughly 20 degrees in
both directions when I plotted it using NCL scripts listed at the end
of this email. My model domain is rather large, 252 x 87 grid points
with 111 km resolution. Latitudes range from -19.5815 to 59.65268
degrees N. Longitudes range from 40.52136 to -40.52136 degrees E.
Mercator map projection was used. Your suggestions on correctly
plotting the model terrain are greatly appreciated.

Thanks,
Yongxin

---------------------------------------------------------------
;
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

fin = addfile("wrfout_d01_2003-03-01_00:00:00.nc","r")

; read character variable Times and convert to string for plots

times = chartostring(fin->Times)

; read model terrain

var = fin->HGT(0,:,:)

; set some resources

wks = gsn_open_wks("ps","HGT")
gsn_define_colormap(wks,"BlAqGrYeOrReVi200")

res = True
res@gsnMaximize = True
res@gsnSpreadColors = True
res@cnFillOn = True
res@cnLinesOn = False
res@cnLineLabelsOn = False
res@gsnAddCyclic = False

; use WRF_contributed procedure to set map resources

WRF_map_c (fin,res,0)

; set True for native projection (faster)

res@tfDoNDCOverlay = True

; turn on lat/lon labeling

res@pmTickMarkDisplayMode = "Always"

; make the plot

;res@mpProjection = "mercator"
;res@mpLimitMode = "Corners"
;res@mpCenterLonF = 180.0
;res@mpLeftCornerLatF = -19.5815
;res@mpLeftCornerLonF = 40.52136
;res@mpRightCornerLatF = 59.65268
;res@mpRightCornerLonF = -40.52136

plot = gsn_csm_contour_map(wks,var,res)

end

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Mar 16 2006 - 14:26:13 MST