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_at_gsnMaximize = True
res_at_gsnSpreadColors = True
res_at_cnFillOn = True
res_at_cnLinesOn = False
res_at_cnLineLabelsOn = False
res_at_gsnAddCyclic = False
; use WRF_contributed procedure to set map resources
WRF_map_c (fin,res,0)
; set True for native projection (faster)
res_at_tfDoNDCOverlay = True
; turn on lat/lon labeling
res_at_pmTickMarkDisplayMode = "Always"
; make the plot
;res_at_mpProjection = "mercator"
;res_at_mpLimitMode = "Corners"
;res_at_mpCenterLonF = 180.0
;res_at_mpLeftCornerLatF = -19.5815
;res_at_mpLeftCornerLonF = 40.52136
;res_at_mpRightCornerLatF = 59.65268
;res_at_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
This archive was generated by hypermail 2.2.0 : Wed Mar 29 2006 - 11:53:03 MST