I found the answer ( path) in the landsea_mask documentation.
Thanks anyways,
Ioana
From: "Ioana Colfescu" <colfescu@cola.iges.org>
To: "Dennis Shea" <shea@ucar.edu>
Sent: Friday, February 11, 2011 5:26:40 PM
Subject: about land sea mask
I still have a question - when I mask the continents the contour around them is very rough.I would want the mask approximating better the contours of the continents but is kinda far.
I'm reading this power point you wrote - NCL Data Processing and you mention here about the utility called landsea_mask.
I also see in mask5.ncl ( in the online documentation) there is :
begin
a = addfile("ts.t85.nc","r") ; read in sample data file (resolution=T85)
t85 = a->TS(0,:, ; grab the first timestep of the TS field
b = addfile("landsea.nc","r") ; read in land sea mask basemap file
lsm_t85 = landsea_mask(b->LSMASK,t85&lat,t85&lon) ; read in land sea mask, and pass it and
; the t85 lat/lon arrays into landsea_mask
t85 = mask(t85,lsm_t85.eq.0,False) ; mask out all ocean points from t85
Is this landsea.nc file somewhere I can use and get is ? I don't have actually a 1*1 grid in the files I'm using they have both more than 1 degree.
Thank you,
Ioana
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Feb 11 16:58:58 2011
This archive was generated by hypermail 2.1.8 : Tue Feb 15 2011 - 09:43:19 MST