defining dimensions

From: Micah Sklut <micahs2005_at_nyahnyahspammersnyahnyah>
Date: Fri, 16 Jun 2006 16:36:01 -0400

Hello,

I am making sst contour maps, and my contours values seem to be slightly off
as far
as their lat/lon locations. As you can see in my script below, I take only a
portion of the whole dataset. I recieve memory errors, when I try to read
in the whole large 1/12 degree global dataset, and I'm only using the
Atlantic values anyway. When I don't declare lat and lon variables, the
output is exactly the same - so I'm thinking this might be where the problem
is. The basic script is attached below:
Thanks for any input you can give.

begin

f = addfile("sst.grb","r")
sst = f->TMP_173_SFC(300:1500,2700:3700)
;lat = f->lat_173(300:1500)
;lon = f->lon_173(2700:3700)

sst = (sst-273) * 1.8 + 32

res_at_mpMinLatF = minlat
res_at_mpMinLonF = minlon
res_at_mpMaxLatF = maxlat
res_at_mpMaxLonF = maxlon

wks = gsn_open_wks("ps",name_out)
gsn_define_colormap(wks,"sstgraph ") ; choose colormap

res = True ; plot mods desired
res_at_mpDataBaseVersion = "Ncarg4_1"
res_at_mpOutlineBoundarySets = "GeophysicalandUSStates"
res_at_mpFillDrawOrder = "postDraw"
res_at_gsnAddCyclic = False

res_at_mpMinLatF = minlat
res_at_mpMinLonF = minlon
res_at_mpMaxLatF = maxlat
res_at_mpMaxLonF = maxlon

res_at_cnLevelSelectionMode = "ManualLevels"
res_at_cnLevelSpacingF = 1.0
res_at_cnMinLevelValF = 33.0
res_at_cnMaxLevelValF = 92.0

res_at_gsnSpreadColors = True
res_at_gsnSpreadColorStart = 2
res_at_gsnSpreadColorEnd = 62

contour = gsn_csm_contour_map_ce(wks,sst,res)

delete(res)

end

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Jun 16 2006 - 14:36:01 MDT

This archive was generated by hypermail 2.2.0 : Mon Jun 19 2006 - 10:30:07 MDT