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) 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 Visiting Scholar Jin-Young Kim (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
This archive was generated by hypermail 2.2.0 : Wed Oct 29 2008 - 09:22:41 MDT