Re: mm52ncdf plot

From: Benjamin Lamptey (lamptey AT essc.psu.edu)
Date: Wed Nov 03 2004 - 12:06:14 MST

  • Next message: Derrick Snowden: "Minor details on plots"

    Hi Dennis,
    Could you please elaborate on the method you specified below?

    Some info
    ---------
    1) My file contains the subsect I wish to plot
    2)
    (0) ==========> printVarInfo: lat2d

    Variable: x (parameter)
    Type: float
    Total Size: 11776 bytes
                2944 values
                Number of Dimensions: 2
                Dimensions and sizes: [MLON | 64] x [NLAT | 46]
                Coordinates:
                Number Of Attributes: 2
                  units : degrees_north
                    long_name : array of latitudes
                    (0) Minimum: 34.8703 Maximum: 50.3709

    3)
    (0) ==========> printVarInfo: lon2d

    Variable: x (parameter)
    Type: float
    Total Size: 11776 bytes
                2944 values
                Number of Dimensions: 2
                Dimensions and sizes: [MLON | 64] x [NLAT | 46]
                Coordinates:
                Number Of Attributes: 2
                  units : degrees_east
                    long_name : array of longitudes
    4) I shall plot variables like
    (0) ==========> printVarInfo: raintot1A_DJF

    Variable: x (parameter)
    Type: float
    Total Size: 11776 bytes
                2944 values
                Number of Dimensions: 2
                Dimensions and sizes: [MLON | 64] x [NLAT | 46]
                Coordinates:
                Number Of Attributes: 5
                  lon2d : <ARRAY>
                    lat2d : <ARRAY>
                      NMO : 0
                        _FillValue : -999
                          average_op_ncl : dim_avg over dimension: time
                          (0) Minimum: 0 Maximum: 7.29543
    5) My NCL script will read in two files; one with the normal lat/lon
    and the second with lat2d/lon2d and make a panel plot.

    Thanks
    Ben5) My NCL script will read in two files; one with the normal lat/lon
    and the second with lat2d/lon2d and make a panel plot.

    Thanks
    Ben5) My NCL script will read in two files; one with the normal lat/lon
    and the second with lat2d/lon2d and make a panel plot.

    Thanks
    Ben5) My NCL script will read in two files; one with the normal lat/lon
    and the second with lat2d/lon2d and make a panel plot.

    Thanks
    Ben5) My NCL script will read in two files; one with the normal lat/lon
    and the second with lat2d/lon2d and make a panel plot.

    Thanks
    Ben
    On Wed, Nov 03, 2004 at 10:24:17AM -0700, Dennis Shea wrote:
    > The way to access a subset of (say) a Lambert conformal
    > projection is via standard subscripting.
    >
    > There is no magic bullet. NCL's syntax [ie the {..}]
    > for accessing coordinates is not possible for arrays
    > with two-dimensional coordinate arrays. This is not
    > supported by the underlying netCDF software.
    >
    > The user must manually specify the subscripts.
    >
    > iLowerLeft = 30
    > jLowerLeft = 24
    > iUpperLeft = 72
    > jUpperLeft = 26
    > [etc]
    >
    > This could, of course, be automated.
    >
    > D
    > >
    > >On Tue, 2 Nov 2004, Benjamin Lamptey wrote:
    > >
    > >> Hello,
    > >> The examples on the following website which illustrates how to plot
    > >> data with 2D lat/lon does not have an example on plotting a subset
    > >>
    > >> http://www.cgd.ucar.edu/csm/support/CSM_Graphics/popscal.shtml
    > >>
    > >> I tried
    > >> res@mpMinLatF = minlat
    > >> res@mpMaxLatF = maxlat
    > >> res@mpMinLonF = minlon
    > >> res@mpMaxLonF = maxlon
    > >>
    > >> where minlat,maxlat,minlon,maxlon are defined. But did not get a
    > >> satisfactory plot (see plot attached).
    > >>
    > >> I do not want to use something like
    > >> res@mpLimitMode = "Corners" ; choose range of map
    > >> res@mpLeftCornerLatF = lat2d(0,0)
    > >> res@mpLeftCornerLonF = lon2d(0,0)
    > >> res@mpRightCornerLatF = lat2d(nlon-1,nlat-1)
    > >> res@mpRightCornerLonF = lon2d(nlon-1,nlat-1)
    > >>
    > >> used for making plots of Lambert Confromal projection that I used when
    > >> I read my data from a binary file because the purpose of converting
    > >> the binary data files to netcdf will be defeated.
    > >>
    > >> The reason why I am converting the binary files to netcdf is because I
    > >> shall have to use one NCL script to make a panel plot of data from
    > >> different files; one file is netcdf, the 2nd file has Lambert Conformal
    > >> Project (binary file) and the third file has Mercator projection (binary).
    > >>
    > >> I shall appreciate any ideas.
    > >>
    > >> I also want to avoid having to set different plotting resources for the
    > >> different datasets if possible. If not, then that will be my last resort
    > >> in which case I do not need to convert my binary files into netcdf,
    > >> correct?
    > >>
    > >> 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
    > >>
    > >>
    > >>
    > >>
    > >
    > >
    > >
    > >******************************************************
    > >Sylvia Murphy EML: murphys AT ucar.edu
    > >NCAR CGD/CSEG PHN: 303-497-1720
    > >1850 Table Mesa Drive FAX: 303-497-1333
    > >Boulder CO 80305
    > >
    > >
    > >WEB: http://www.cgd.ucar.edu/csm/support/
    > > http://www.cgd.ucar.edu/csm/support/CSM_Graphics/
    > >******************************************************
    > >
    > >
    > >_______________________________________________
    > >ncl-talk mailing list
    > >ncl-talk AT ucar.edu
    > >http://mailman.ucar.edu/mailman/listinfo/ncl-talk
    >

    -- 
    

    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 AT ucar.edu http://mailman.ucar.edu/mailman/listinfo/ncl-talk



    This archive was generated by hypermail 2b29 : Thu Nov 04 2004 - 08:05:50 MST