map shifting

From: ahmed lasheen <ahmed4kernel_at_nyahnyahspammersnyahnyah>
Date: Wed Jan 05 2011 - 10:44:54 MST

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
Received on Wed Jan 5 10:45:41 2011

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