gsn_csm_contour_map_polar masked for land_only?

From: Anji Seth <anji.seth_at_nyahnyahspammersnyahnyah>
Date: Tue, 5 Dec 2006 16:10:32 -0500

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