Hi.
I'm just want to make a contour plot of Land surface temperature from a MODIS-HDF file and I'm getting this error message:
warning:ScalarFieldSetValues: coordinate array sfXArray requires 15475147 elements: defaulting
I don't have any idea how to solve this problem. Could you give me some suggestions? You can see the code below.
Thank you very much.
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
hdf=addfile(dir+"MYD11C2.A2009073.005.2009090205119.hdf","r")
n=hdf->LST_Night_CMG*0.02
n@_FillValue=0
n@_FillValue=-999
n_at_units="(K)"
r=True
r_at_tiMainString = "Temperature (K)"
r_at_tiXAxisString = "Longitude"
r_at_tiYAxisString = "Latitude"
r_at_gsnMaximize = True
r_at_cnFillOn = True
r_at_cnMonoFillPattern = True
r_at_cnMonoFillColor = False
r_at_cnMonoLineColor = False
r_at_cnLineLabelsOn = False
r_at_cnInfoLabelOn = False
r_at_pmLabelBarDisplayMode = "Always"
r_at_lbPerimOn = False
r_at_sfXArray = lon
r_at_sfYArray = lat
r_at_mpMinLatF = 40.0 ; map area
r_at_mpMaxLatF = 41.6 ; latitudes
r_at_mpMinLonF = -75.0 ; and
r_at_mpMaxLonF = -72.0 ; longitudes
nlat = 3600
nlon=7200
lat = latGlobeFo(nlat, "lat", "latitude", "degrees_north")
lat = lat(::-1)
lon = lonGlobeFo(nlon, "lon", "longitude", "degrees_east")
n!0 = "lat"
n!1 = "lon"
n&lat = lat
n&lon = lon
n=lonFlip(n)
printVarSummary(n)
wks = gsn_open_wks("x11","h3")
setvalues NhlGetWorkspaceObjectId()
"wsMaximumSize" : 110000000
end setvalues
gsn_define_colormap(wks,"BlueDarkRed18")
plot = gsn_csm_contour_map(wks,n,r)
end
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Apr 16 2009 - 01:46:10 MDT
This archive was generated by hypermail 2.2.0 : Thu Apr 16 2009 - 10:05:49 MDT