map projection

From: Michael Notaro <mnotaro_at_nyahnyahspammersnyahnyah>
Date: Tue, 12 Feb 2008 10:35:24 -0600

I am trying to plot data with the following projection information.
However, the data doesn't seem to correctly plot in the correct
position on the map. A clip of the related code is included below.
Any idea what might be wrong? Since I have 4 points specified
in the projection data, is mpLimitMode="Corners" correct to use,
since it only allows 2 points? I am also not entirely sure about
mpLambertMeridianF.

Thanks, Mike

Lower left lat : -2.417 deg
Lower left lon : -143.913 deg
Upper left lat : 54.626 deg
Upper left lon : 157.412 deg
Lower right lat : -2.382 deg
Lower right lon : -62.143 deg
Upper right lat : 54.686 deg
Upper right lon : -3.476 deg
projection = ALBERS Conical Equal-area projection
                            uses the clark ellipsoid
first standard parallel = 20
second standard parallel = 60
lat0 = 45.0
lon0 = -103.0
cell height = 8 km
cell width = 8 km
xsize = 1280
ysize = 1024

--------------------------

res_at_gsnAddCyclic = False
res_at_mpLimitMode = "Corners"
res_at_mpLeftCornerLatF = lat(0,0)
res_at_mpLeftCornerLonF = lon(0,0)
res_at_mpRightCornerLatF = lat(1023,1279)
res_at_mpRightCornerLonF = lon(1023,1279)
res_at_mpProjection = "LambertConformal"
res_at_mpLambertParallel1F = 20.
res_at_mpLambertParallel2F = 60.
res_at_mpLambertMeridianF = -106.
res_at_tfDoNDCOverlay = True
res_at_trYReverse=True
plot = gsn_csm_contour_map(wks,var,res)

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Feb 12 2008 - 09:35:24 MST

This archive was generated by hypermail 2.2.0 : Tue Feb 12 2008 - 14:45:27 MST