Re: Native-NCL and MM5

I am trying to plot MM5 output on its native grid using example
 ncl scripts. I have all the relevant informataion (see below)
 required and I do get a plot(attached). However, the coastline map
 does not match the terrain elevation data. Whatelse do I need to
 do to correct this problem.

 I have a 130x130 grid on MM5 Polar stereographic, center lat 68
 N,center lon. -80 (west), true lat1 60 true lat2 0; pole position
 90 latitude.

 Here are my relevant ncl settings

 ; Input data from MM5 model
 ;***************************************
  hgt = asciiread("TERRAIN",(/130,130/),"float")
  hgt_at_units = "m"
  hgt_at_long_name = "MM5 Terrain height "
 ;***************************************
 ; Input lat lon data
 ;***************************************
  lat = asciiread("LATD1",(/130,130/),"float")
 lon = asciiread("LOND1",(/130,130/),"float")

 ; the following resources are REQUIRED to plot this projection correctly
 res_at_mpProjection = "Stereographic" ; projection>
res_at_mpLimitMode = "Corners" ; method to zoom
 res_at_mpLeftCornerLatF = lat(0,0)
 res_at_mpLeftCornerLonF = lon(0,0)
 res_at_mpRightCornerLatF = lat(129,129)
 res_at_mpRightCornerLonF = lon(129,129)
 res_at_mpRelativeCenterLon = True ; set a center lon
  res_at_mpCenterLonF = -80.0 ; Coarse domain center longitude
 res_at_mpRelativeCenterLat = True ; set a center lat
 res_at_mpCenterLatF = 68. ; center lat
res_at_tfDoNDCOverlay = True ; do not transform data

 plot = gsn_csm_contour_map (wks,hgt,res)

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Jul 21 2006 - 12:01:21 MDT