Hi Adam,
My GCM data is in a netcdf file while my regional model data is
binary.
For the regional model (using the Lambert Conformal projection), I
use resources like
mres@mpLimitMode = "Corners" ; choose range of map
mres@mpLeftCornerLatF = lat2d(0,0)
mres@mpLeftCornerLonF = lon2d(0,0)
mres@mpRightCornerLatF = lat2d(nlon-1,nlat-1)
mres@mpRightCornerLonF = lon2d(nlon-1,nlat-1)
mres@mpProjection = "LambertConformal"
mres@mpLambertParallel1F = (max(lat2d)+min(lat2d))*0.5
mres@mpLambertParallel2F = (max(lat2d)+min(lat2d))*0.5
etc and the correct geographical region is plotted.
Thanks
Ben
On Tue, Jan 18, 2005 at 02:31:57PM -0700, Adam Phillips wrote:
> Hi Ben,
>
> It looks like your regional model has 2D latitudes and longitudes, and thus is
> on a curvilinear grid. NCL can handle this, but you have to assign some
> attributes to your data array "x", using the 2 arrays you sent us printVarInfo
> info about:
>
> x@lon2d = lon2d
> x@lat2d = lat2d
>
> When you plot x, NCL will then recognize that you have 2D lats/lons.
>
> >4) Is the correct data for the GCM extracted? Since the lon is
> >from 0 to 360 for the data, is NCL able to identify the lon value
> >corresponding to -90.0246 (minlon) and -57.9754 (maxlon)?
>
> Yes, NCL is able to recognize the negative longitudes when plotting, even though
> your GCM longitudes go from 1->360...
>
> Adam
>
> >Delivered-To: asphilli@ucar.edu
> >Delivered-To: ncl-talk@ucar.edu
> >Date: Tue, 18 Jan 2005 15:54:31 -0500
> >From: Benjamin Lamptey <lamptey@essc.psu.edu>
> >To: ncl-talk@ucar.edu
> >Mime-Version: 1.0
> >Content-Disposition: inline
> >User-Agent: Mutt/1.4.2.1i
> >Subject: Plots
> >X-BeenThere: ncl-talk@ucar.edu
> >X-Mailman-Version: 2.1.1
> >List-Id: NCAR Command Language User Group <ncl-talk.ucar.edu>
> >List-Unsubscribe: <http://mailman.ucar.edu/mailman/listinfo/ncl-talk>,
> <mailto:ncl-talk-request@ucar.edu?subject=unsubscribe>
> >List-Post: <mailto:ncl-talk@ucar.edu>
> >List-Help: <mailto:ncl-talk-request@ucar.edu?subject=help>
> >List-Subscribe: <http://mailman.ucar.edu/mailman/listinfo/ncl-talk>,
> <mailto:ncl-talk-request@ucar.edu?subject=subscribe>
> >
> >Hello,
> >I am comparing a GCM and a regional model. The longitudes in the
> >go from 1 to 360.
> >
> >2) I have the following from the regional model
> >
> >==========> printVarInfo: lat2d
> >
> >Variable: x (parameter)
> >Type: float
> >Total Size: 11776 bytes
> > 2944 values
> > Number of Dimensions: 2
> > Dimensions and sizes: [64] x [46]
> > Coordinates:
> > (0) Minimum: 34.8703 Maximum: 50.3709
> >
> >
> >(0) ==========> printVarInfo: lon2d
> >
> >
> >Variable: x (parameter)
> >Type: float
> >Total Size: 11776 bytes
> > 2944 values
> > Number of Dimensions: 2
> > Dimensions and sizes: [64] x [46]
> > Coordinates:
> > (0) Minimum: -90.0246 Maximum: -57.9754
> >
> >3) So for the GCM resources, I do the following
> >
> > gres@mpMinLatF = minlat
> > gres@mpMaxLatF = maxlat
> > gres@mpMinLonF = minlon
> > gres@mpMaxLonF = maxlon
> >
> >The tickmarks for both the GCM and the regional model indicate that
> >the same geographical region is being plotted.
> >
> >4) Is the correct data for the GCM extracted? Since the lon is
> >from 0 to 360 for the data, is NCL able to identify the lon value
> >corresponding to -90.0246 (minlon) and -57.9754 (maxlon)?
> >
> >The plots look too different. I want to be sure I am not making a
> >mistake in plotting.
> >
> >Thanks
> >Ben
> >--
> >
> >
> >Benjamin L. Lamptey Phone: (814) 865-9911 (office)
> >Pennstate (EMS Environment Institute) (814) 237-8193 (home)
> >2217 Earth-Engineering Sciences Bldg Fax : (814) 865-3191
> >University Park, PA 16802 WWW:http://www.essc.psu.edu/~lamptey
> >
> >
> >
> >
> >----- End forwarded message -----
> >
> >--
> >
> >
> >Benjamin L. Lamptey Phone: (814) 865-9911 (office)
> >Pennstate (EMS Environment Institute) (814) 237-8193 (home)
> >2217 Earth-Engineering Sciences Bldg Fax : (814) 865-3191
> >University Park, PA 16802 WWW:http://www.essc.psu.edu/~lamptey
> >
> >
> >
> >_______________________________________________
> >ncl-talk mailing list
> >ncl-talk@ucar.edu
> >http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
> -------------------------------------------------------------
> Adam Phillips email: asphilli@ucar.edu
> Climate and Global Dynamics Division tel: (303) 497-1726
> National Center for Atmospheric Research fax: (303) 497-1333
> P.O. Box 3000
> Boulder, CO 80307-3000 http://www.cgd.ucar.edu/~asphilli
>
--Benjamin L. Lamptey Phone: (814) 865-9911 (office) Pennstate (EMS Environment Institute) (814) 237-8193 (home) 2217 Earth-Engineering Sciences Bldg Fax : (814) 865-3191 University Park, PA 16802 WWW:http://www.essc.psu.edu/~lamptey
_______________________________________________ ncl-talk mailing list ncl-talk@ucar.edu http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2b29 : Tue Jan 18 2005 - 15:09:33 MST