<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&#39;t represent the locations correctly. For example on y-axis &quot;40N&quot; 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">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=&#39;LAMCON&#39;) ! 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 &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl&quot;<br><br><br> begin<br><br> fin  = addfile(&quot;DOMAIN90km.nc&quot;,&quot;r&quot;)            ; Open a netcdf file<br> ht   = fin-&gt;ht                                 ; Terrain Elevation<br> lu   = fin-&gt;lu                                 ; Land Use Type<br> xlon = fin-&gt;xlon                              ; Cross Grid Longitude<br> xlat = fin-&gt;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=&quot;eps&quot;                                                  ; for ps/pdf plots<br> wkstyp@wkOrientation = &quot;Portrait&quot;                             ; for ps/pdf plots<br> wks = gsn_open_wks(wkstyp,&quot;Domains_90km&quot;)               ; Open an X11 workstation.<br><br><br> res1                     = True                                                  ; Indicate you want to set some resources.<br> res1@gsnMaximize         = True<br> res1@gsnPaperOrientation = &quot;Portrait&quot;              ;for ps/pdf plots<br> res1@cnFillMode          = &quot;CellFill&quot;<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 = &quot;ExplicitLevels&quot;<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,&quot;OceanLakeLandSnow&quot;)<br> res1@gsnSpreadColors    = True<br> res1@gsnSpreadColorEnd  = -2<br> res1@gsnCenterString    = &quot; &quot;<br><br> res1@tiMainString    = &quot;&quot;<br> res1@gsnLeftString   = &quot;&quot;<br> res1@gsnRightString  = &quot;&quot;<br><br> res1@pmLabelBarDisplayMode = &quot;Always&quot;<br> res1@lbPerimOn             = False<br> res1@lbLabelAlignment      = &quot;ExternalEdges&quot;<br> res1@lbLabelStride         = 10<br> res1@lbOrientation         = &quot;Vertical&quot;<br> res1@tmBorderThicknessF    = 1.<br><br> res1@pmTickMarkDisplayMode  = &quot;Always&quot;<br><br> res1@mpDataSetName     = &quot;Earth..4&quot;<br> res1@mpDataBaseVersion = &quot;MediumRes&quot;<br> res1@mpOutlineBoundarySets = &quot;AllBoundaries&quot;<br><br> res1@mpPerimOn             = True<br> res1@mpFillDrawOrder = &quot;Postdraw&quot;<br> res1@mpLandFillColor = &quot;Transparent&quot;<br> res1@mpOceanFillColor = &quot;light blue&quot;<br> res1@mpInlandWaterFillColor = &quot;light blue&quot;<br><br> res1@mpLimitMode       = &quot;Corners&quot;<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        = &quot;LambertConformal&quot;<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>