Re: Map shift

Hi Yongxin,

Can you upload your data file to ftp.cgd.ucar.edu?
We would like to take a look at this problem.

...
ftp ftp.cgd.ucar.edu
Name: anonymous
Password: <your email address>
cd incoming
bin
put <file>

Let us know when the file has been uploaded and its name.
We can get the file only if we know the name.
  -dave

On Mar 16, 2006, at 2:26 PM, Yongxin Zhang wrote:

> 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

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Mar 16 2006 - 15:34:40 MST