<div dir="ltr"><div><div><div><div>Hi,<br><br></div>I think this is an artifact of the distortion of the particular projection being used. To illustrate what's going on, try adding these resources:<br><br>res1@mpGridAndLimbOn = True<br>res1@mpGridAndLimbDrawOrder = "PostDraw"<br><br></div>I think the resulting plot will explain fairly well what is going on, but if not please write back to the group.<br><br></div>Hope that helps...<br></div>Rick<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Mar 25, 2015 at 11:18 AM, Adam Will <span dir="ltr"><<a href="mailto:ashleywats@gmail.com" target="_blank">ashleywats@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Dear NCL Users,<br></div><br> I need help correctly plotting model ( ICTP-RegCM) topography. Below is the script and input netcdf file (195 KB). The problem is that the geographical coordinates on the axes doesn't represent the locations correctly. For example on y-axis "40N" appears at location where it should be actually 47N. I used the map projection information from the domain setting file (also provided below).<br></div><div><br>Thanks in advance<br></div><div><br>Regards<br></div><div>Ashley<br><br><br><a href="https://www.dropbox.com/s/gf30luj63mqckus/DOMAIN90km.nc?dl=0" target="_blank">https://www.dropbox.com/s/gf30luj63mqckus/DOMAIN90km.nc?dl=0</a> (Input Netcdf file)<br><br><a href="https://www.dropbox.com/s/j6p4pktbdyb5438/Domains_90km.eps?dl=0" target="_blank">https://www.dropbox.com/s/j6p4pktbdyb5438/Domains_90km.eps?dl=0</a> (Output eps file)<br><br>Information from Domain setting file.<br><br>parameter(iproj='LAMCON') ! Map projection (LAMCON=Lambert conformal)<br>parameter(clat= 47.00) ! Central latitude (North positive)<br>parameter(clon= 10.50) ! Central longitude (East positive)<br>parameter(truelatL= 30.) ! LAMCON true latitude (low latitude side)<br>parameter(truelatH= 60.) ! LAMCON true latitude (high latitude side)<br><br><br><br><br></div><div><div><div>-------- NCL Script -------<br><br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"<br><br><br> begin<br><br> fin = addfile("DOMAIN90km.nc","r") ; Open a netcdf file<br> ht = fin->ht ; Terrain Elevation<br> lu = fin->lu ; Land Use Type<br> xlon = fin->xlon ; Cross Grid Longitude<br> xlat = fin->xlat ; Cross Grid Latitude<br><br> dims = dimsizes(xlon)<br> nlat=dims(0)<br> nlon =dims(1)<br><br> ht= where(ht.lt.0 , 0 , ht)<br><br> wkstyp="eps" ; for ps/pdf plots<br> wkstyp@wkOrientation = "Portrait" ; for ps/pdf plots<br> wks = gsn_open_wks(wkstyp,"Domains_90km") ; Open an X11 workstation.<br><br><br> res1 = True ; Indicate you want to set some resources.<br> res1@gsnMaximize = True<br> res1@gsnPaperOrientation = "Portrait" ;for ps/pdf plots<br> res1@cnFillMode = "CellFill"<br> res1@cnFillOn = True<br> res1@cnMonoFillPattern = True<br> res1@cnMonoFillColor = False<br> res1@cnLineLabelsOn = False<br> res1@cnInfoLabelOn = False<br> res1@cnLinesOn = False<br> res1@cnLevelSelectionMode = "ExplicitLevels"<br> res1@cnLevels = (/0.,5.,10.,15.,20.,25.,30.,35.,40.,45.,50.,\<br> 95.,140.,185.,230.,275.,320.,365.,410.,455.,500.,\<br> 570.,640.,710.,780.,850.,920.,990.,1060.,1130.,1200.,\<br> 1250.,1300.,1350.,1400.,1450.,1500.,1550.,1600.,1650.,1700.,\<br> 1750.,1800.,1850.,1900.,1950.,2000.,2050.,2100.,2150.,2200.,\<br> 2250.,2300.,2350.,2400.,2450.,2500.,2550.,2600.,2650.,2700./)<br><br> gsn_define_colormap(wks,"OceanLakeLandSnow")<br> res1@gsnSpreadColors = True<br> res1@gsnSpreadColorEnd = -2<br> res1@gsnCenterString = " "<br><br> res1@tiMainString = ""<br> res1@gsnLeftString = ""<br> res1@gsnRightString = ""<br><br> res1@pmLabelBarDisplayMode = "Always"<br> res1@lbPerimOn = False<br> res1@lbLabelAlignment = "ExternalEdges"<br> res1@lbLabelStride = 10<br> res1@lbOrientation = "Vertical"<br> res1@tmBorderThicknessF = 1.<br><br> res1@pmTickMarkDisplayMode = "Always"<br><br> res1@mpDataSetName = "Earth..4"<br> res1@mpDataBaseVersion = "MediumRes"<br> res1@mpOutlineBoundarySets = "AllBoundaries"<br><br> res1@mpPerimOn = True<br> res1@mpFillDrawOrder = "Postdraw"<br> res1@mpLandFillColor = "Transparent"<br> res1@mpOceanFillColor = "light blue"<br> res1@mpInlandWaterFillColor = "light blue"<br><br> res1@mpLimitMode = "Corners"<br> res1@mpLeftCornerLatF = xlat(0,0)<br> res1@mpLeftCornerLonF = xlon(0,0)<br> res1@mpRightCornerLatF = xlat(nlat-2,nlon-2)<br> res1@mpRightCornerLonF = xlon(nlat-2,nlon-2)<br><br> res1@mpProjection = "LambertConformal"<br> res1@mpLambertParallel1F = 30.<br> res1@mpLambertParallel2F = 60.<br> res1@mpLambertMeridianF = 10.5<br><br> res1@gsnAddCyclic = False<br> res1@tfDoNDCOverlay = True<br><br> map = gsn_csm_contour_map(wks,ht,res1)<br>end<br><br></div></div></div></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>