Fwd: Contour plot on Lambert Conformal Projection (Extremely Slow!)

Hi,

This is Ling. I am trying to make a contour plot on Lambert Conformal
Projection. The data is a two dimensional array with 872 x 764 grids.
However, the problem is it is extremely slow to generate a contour plot,
sometimes more than a hour. I am wondering is there something wrong with my
code or is there any way speed up the plotting? Following is my code.

Thanks

Ling

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin

start_time = get_cpu_time()
filename = "isop_200505.nc"
latfile = "lat.csv"
lonfile = "lon.csv"

f = addfile(filename,"r")
lat_1d = asciiread(latfile,-1,"float")
lon_1d = asciiread(lonfile,-1,"float")
lat2d = onedtond(lat_1d,(/872,764/))
lon2d = onedtond(lon_1d,(/872,764/))
dimlc = dimsizes(lat2d)
nlat = dimlc(0)
nlon = dimlc(1)
isop_dailyavg = f->ISOP_DAILYAVG
isop_dailyavg = where(ismissing(isop_dailyavg),0,isop_dailyavg)
isop_dailyavg@lat2d = lat2d
isop_dailyavg@lon2d = lon2d

;******************create plots*******************
fon = "May_lambert"
wks = gsn_open_wks("pdf",fon)

res = True
res@gsnDraw = False
res@gsnFrame = False
res@cnFillOn = True
res@cnLinesOn = False
res@cnFillPalette = "dayisop"
res@mpGridAndLimbOn = True
res@pmTickMarkDisplayMode = "Always"
res@tmXTOn = False
res@gsnAddCyclic = False

res@mpLimitMode = "Corners"
res@mpLeftCornerLatF = lat2d(0,0)
res@mpLeftCornerLonF = lon2d(0,0)
res@mpRightCornerLatF = lat2d(nlat-1,nlon-1)
res@mpRightCornerLonF = lon2d(nlat-1,nlon-1)
res@tfDoNDCOverlay = True
res@mpProjection = "LambertConformal"
res@mpLambertParallel1F = 33
res@mpLambertParallel2F = 45
res@mpLambertMeridianF = -97

res@mpFillOn = False

res@mpOutlineBoundarySets = "AllBoundaries"
res@mpDataBaseVersion = "Ncarg4_1"
res@mpDataSetName = "Earth..4"
res@mpUSStateLineThicknessF = 1.5
res@mpUSStateLineColor = "Black"
res@mpNationalLineThicknessF = 1.5
res@mpNationalLineColor = "Black"
res@mpCountyLineThicknessF = 0.5
res@mpCountyLineColor = "grey70"

res@mpGridAndLimbOn = False
res@mpPerimOn = True

res@cnLevelSelectionMode = "ExplicitLevels"
res@cnLevels =
(/0.004,0.008,0.012,0.016,0.020,0.024,0.028,0.032/)

res@gsnCenterString = "isoprene daily average
emssion(tons/km2/day)"
plot = gsn_csm_contour_map(wks,isop_dailyavg,res)

draw(plot)
frame(wks)

end_time = get_cpu_time()
print(end_time - start_time)

end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jul 2 10:33:35 2013