Re: map shifting

From: David Brown <dbrown_at_nyahnyahspammersnyahnyah>
Date: Wed Jan 05 2011 - 11:41:37 MST

Hi Ahmed,

Where did you get the values for the center lat and lon? Are they part of the metadata for the file? It would help if you can upload your data file to our ftp server. See the directions at
http://www.ncl.ucar.edu/report_bug.shtml. Thanks,
 -dave

On Jan 5, 2011, at 10:44 AM, ahmed lasheen wrote:

> Hello
> I trying to plot the terrain of the NMM-WRF ,but the map is being shifted from the terrain (or reverse) , the projection of data is rotated lat lon .
> i have used the CylindricalEquidistant projection ,as i read from the ncl-talk archive that is like the rotated lat lon.
> here is the code i wrote
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> ;******************************************************
> begin
> ;****************** Open file and get variables *************
> a = addfile("../../data/wrf_10111200//geo_nmm_nest.l01.nc","r")
> ter = a->HGT_M (0,:,:)
> lat2d = a->XLAT_M (0,:,:)
> lon2d = a->XLONG_M(0,:,:)
> ;*******************make the graphics ***********************
> wks = gsn_open_wks("pdf" ,"terrain") ; ps,pdf,x11,ncgm,eps
> gsn_define_colormap(wks,"gui_default") ; select color map
>
> res = True ; plot mods desired
> res@cnFillOn = True ; turn on color
> res@cnLinesOn = False ; turn contour lines
> res@mpMaxLatF = max(lat2d)
> res@mpMinLatF = min(lat2d)
> res@mpMaxLonF = max(lon2d)
> res@mpMinLonF = min(lon2d)
> res@gsnSpreadColors = True
> res@lbLabelStride = 2
> res@mpLimitMode = "latlon" ; use lat/lon coordinates to limit area
> res@mpProjection = "CylindricalEquidistant"
> res@mpGreatCircleLinesOn = True
> res@gsnCenterString = res@mpProjection+" Projection"
> res@tfDoNDCOverlay = True
> res@mpOutlineBoundarySets = "National"
> res@mpCenterLatF = 37
> res@mpCenterLonF = 22
> plot = gsn_csm_contour_map_ce(wks,ter(:,:),res)
> end
>
> thanks in advance
>
> --
> ===============
> Ahmed Lasheen
> Egyptian Meteorological Authority(EMA)
> Cairo,Egypt
> ===============
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jan 5 11:41:44 2011

This archive was generated by hypermail 2.1.8 : Thu Jan 13 2011 - 09:24:21 MST