Hello again,
I sorted out the problem, in the output I'm using, land=2.
Sorry for the spam. -Anji
Hello,
I'm working on some southern hemisphere polar plots for both a global
model
and a regional model (using a native grid) and I'd like to mask
the model output to show land only results,
I tried the following ( in both global and regional model scripts):
;************************************************
; use mask function to mask out land and then ocean data
; ocean=0,land=1,sea_ice=2
;************************************************
land_only = ts ; ts is a time average
ocean_only = ts
land_only = mask(ts,oro,1) ; oro looks good from a printVarSummary
ocean_only = mask(ts,oro,0)
; land_only is ok in printVarSummary
and used this in ( for the global model):
plot = gsn_csm_contour_map_polar(wks,land_only({-50.:-90.},:),res)
and (for the regional model:
plot = gsn_csm_contour_map(wks,land_only({-50.:-90.},:),res)
but I get a plot that shows no data and states " NO CONTOUR DATA"
from both scripts ( global and regional) with an error:
warning:ContourPlotInitialize: no valid values in scalar field;
ContourPlot not possi
ble:[errno=1101]
Is there a way to do this? Thank you in advance for your help.
-Anji
Anji Seth, Geography, University of Connecticut
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Dec 05 2006 - 14:10:32 MST
This archive was generated by hypermail 2.2.0 : Thu Dec 07 2006 - 16:52:36 MST