Re: To overlay data on cylindrical equidistant projection grid map

The GRIB file you are reading looks something like:

    float A_PCP_GDS0_SFC_acc3h ( g0_lat_0, g0_lon_1 )
    float g0_lat_0 ( g0_lat_0 )
    float g0_lon_1 ( g0_lon_1 )

You can plot this directly ....

   apcp = f->A_PCP_GDS0_SFC_acc3h
           :
           :
  plot=gsn_csm_contour_map_ce(wks,apcp,rescn)

Jin-Young Kim wrote:
> Greetings.
>
> I have problem to overlay data on cylindrical equidistant projection
> grid.
>
> When i used following option,
> Map looks strange.. Would you look at this one to make correct.?
>
>
> lat = afile->g0_lat_0
> lon = afile->g0_lon_1
>
> ;=============================
> ; create lat and long coordinate variables
> ;============================
> lon!0 = "lon"
> lon_at_long_name = "lon"
> lon_at_units = "degrees-east"
> lon&lon = lon
>
> lat!0 = "lat"
> lat_at_long_name = "lat"
> lat_at_units = "degrees_north"
> lat&lat = lat
>
> printVarSummary(lat)
>
> print(max(lat))
> print(min(lat))
> print(max(lon))
> print(min(lon))
>
> ;************************************************
> wks=gsn_open_wks("ps",ofim(i))
>
> rescn = True
> rescnFillOn = True ; turn on color
> rescn_at_cnLinesOn = True ; no contour lines
> rescn_at_gsnSpreadColors = True ; use full
> color map
> rescn_at_lbLabelAutoStride = True ; every
> other label
>
> rescn_at_tiMainString = "Precipitation (mm) " + titl(i) + "h
> Fcst.";title
> rescn_at_pmTickMarkDisplayMode = "Always" ; turn on
> tick marks
>
> rescn_at_gsnAddCyclic = False ; regional data, don't
> add pt
> ; rescn_at_mpLimitMode = NhlMAXIMALAREA
> rescn_at_mpMinLatF = min(lat)
> rescn_at_mpMinLonF = min(lon)
> rescn_at_mpMaxLatF = max(lat)
> rescn_at_mpMaxLonF = max(lon)
>
>
> plot=gsn_csm_contour_map_ce(wks,apcp,rescn)
>
> After run the program..
>
> Variable: lat
> Type: float
> Total Size: 532 bytes
> 133 values
> Number of Dimensions: 1
> Dimensions and sizes: [lat | 133]
> Coordinates:
> lat: [12.436..4.516]
> Number Of Attributes: 9
> La1 : 12.436
> Lo1 : -76.428
> La2 : 20.356
> Lo2 : -67.068
> Di : 0.06
> Dj : 0.06
> units : degrees_north
> GridType : Cylindrical Equidistant Projection Grid
> long_name : lat
> (0) 12.436
> (0) 4.516
> (0) -67.068
> (0) -76.428
>
> --
> Ph. D. Student
> Jin-Young Kim (jin_at_climate.pknu.ac.kr <mailto:jin_at_climate.pknu.ac.kr>)
> Integrated Climate System Modeling Lab.
> Department of Environmental Atmospheric Sciences
> Pukyong National University, Busan 608-737, South Korea
> Web : http:climate.pknu.ac.kr/~jin <http://climate.pknu.ac.kr/%7Ejin>
>
> Visiting Scholar
> Jin-Young Kim (Jin-Young.Kim_at_noaa.gov <mailto:Jin-Young.Kim_at_noaa.gov>)
> Global Systems Division
> Earth System Research Laboratory
> R/GSD7
> 325 Broadway
> Boulder, CO80305-3328
>
> +1-303- 497-4422 (Office)
> +1-303- 497-7262 (FAX)
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
======================================================
Dennis J. Shea                  tel: 303-497-1361    |
P.O. Box 3000                   fax: 303-497-1333    |
Climate Analysis Section                             |
Climate & Global Dynamics Div.                       |
National Center for Atmospheric Research             |
Boulder, CO  80307                                   |
USA                        email: shea 'at' ucar.edu |
======================================================
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk