Re: map shifting

From: ahmed lasheen <ahmed4kernel_at_nyahnyahspammersnyahnyah>
Date: Thu Jan 06 2011 - 01:46:27 MST

Hello
I have uploaded the data file to the ftp server. about the center values
i got them from the wrf namelist ref_lat and ref_lon which indicate the
center in which the domain is rotated.
thanks in advance
ahmed

On Wed, Jan 5, 2011 at 1:41 PM, David Brown <dbrown@ucar.edu> wrote:

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

-- 
===============
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 Thu Jan 6 01:47:19 2011

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