Coloring land black in oceanic plot

Hi,

I am plotting oceanic data in a lat vs level plot using isolines. I don't
know how to color the land black. Kindly help. Below is my NCL script.

Thanks and regards,
Vinit Asher.

;*******************moc_iso_lines2.ncl***********************
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

;===== user defined params

casename = "eofc2240.01"
datelab_mn = "0950-0999"
filename_mn="/ptmp/vasher/csm/"+casename+"/ocn/tseries/eofc2240.01.pop.h.TVT
S.0950_cat_0999.nc"
;filename_mn =
"/ptmp/nagaraj/data/b30.031/ocn/b30.031.pop.h.TVTS.0830_cat_0879.nc"
plotname_mn =
"/ptmp/vasher/plots/petm/"+casename+"/ocn/timeseries/MOC/"+casename+"_moc_bw
_0950_0999_moc"

;==== read in files and variables
fmn = addfile(filename_mn,"r")

;--test case
time = fmn->time
time = time/365.
;print("Time = "+time)
itime= floattointeger(time)
ntimes=dimsizes(itime)

lev = fmn->lev ;meters
nlev= dimsizes(lev)

lat = fmn->lat
nlat= dimsizes(lat)

gmsf_mn = fmn->gmsf ;meridional stream function (time,lev,lat)

print(casename+ " Data read")

gmsf_avg = dim_avg(gmsf_mn(lev|:,lat|:,time|:))

print("Averaging complete")

wks = gsn_open_wks("ps",plotname_mn)
gsn_define_colormap(wks,"BlAqGrYeOrReVi200")

;define y axis for test case
ilev = floattointeger(lev)
;ylabels = new(nlev,string)
ylabels = new(6,string)
ylabels(0) = "0"
ylabels(1) = "1000"
ylabels(2) = "2000"
ylabels(3) = "3000"
ylabels(4) = "4000"
ylabels(5) = "5000"
plot_mn = new(1,graphic)

; mos
res = True
res_at_gsnFrame = False
res_at_gsnDraw = False
res_at_vpWidthF = 0.8
res_at_vpHeightF = 0.35
res_at_gsnPaperOrientation = "portrait"
res_at_cnFillOn = False
res_at_cnLinesOn = True
res_at_lbLabelStride = 3
res_at_cnInfoLabelOn = True
res_at_cnLineLabelsOn = True
res_at_txFontHeightF = 0.017
res_at_tmXBOn = True
res_at_tmYLMode = "Explicit"
res_at_tmYLPrecision = 2
res_at_tmYLLabelFontHeightF = 0.015
res_at_tiYAxisFontHeightF = 0.02
res_at_tiMainFontHeightF = 0.02
res_at_trYReverse = True
res_at_gsnYAxisIrregular2Linear = True
res_at_lbBoxMinorExtentF = 0.15
res_at_lbLabelFontHeightF = 0.013
res_at_cnLevelSelectionMode = "AutomaticLevels"
;res_at_cnLevelSelectionMode = "ExplicitLevels"
res_at_cnLevelSelectionMode = "ManualLevels"
res_at_cnMaxLevelValF = 40
res_at_cnMinLevelValF = -30
res_at_cnLevelSpacingF = 2
;res_at_lbLabelStride = 2

res_at_cnLevels =
(/-40.,-30.,-25.,-20.,-15.,-10.,-5.,0.,5.,10.,15.,20.,25.,30.,40./)

;test 0 - 5000m
res_at_gsnCenterString = casename + " Global MOS - Annual Mean [yrs " +
datelab_mn + "]"
res_at_sfYArray = lev
res_at_sfXArray = lat
;res_at_tmYLValues = lev
res_at_tmXBLabelsOn = True ; do not draw bottom labels
res_at_tmXBOn = True ; no bottom tickmarks
res_at_gsnYAxisIrregular2Linear = True
res_at_trYMaxF = 5000.
res_at_trYMinF = 0.
res_at_lbLabelBarOn = True
res_at_tiYAxisString = "depth (m)"
res_at_lbLabelBarOn = True
plot_mn0 = gsn_csm_contour(wks,gmsf_avg,res)
plot_mn(0) = plot_mn0

resP = True
resP_at_gsnMaximize = True
gsn_panel(wks,plot_mn,(/1,1/),resP)

system("convert -crop 0x0 -rotate 270 " + plotname_mn + ".ps" + " " +
plotname_mn + ".png")

print("MOS plotted")

end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Feb 24 2009 - 10:35:18 MST