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


I'm working on some southern hemisphere polar plots for both a global
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

Is there a way to do this? Thank you in advance for your help.


Anji Seth, Geography, University of Connecticut

ncl-talk mailing list
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