landsea_mask problem - answer found

From: Ioana Colfescu <colfescu_at_nyahnyahspammersnyahnyah>
Date: Fri Feb 11 2011 - 16:58:53 MST

I found the answer ( path) in the landsea_mask documentation.

Thanks anyways,


From: "Ioana Colfescu" <>
To: "Dennis Shea" <>
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 :

  a = addfile("","r") ; read in sample data file (resolution=T85)
  t85 = a->TS(0,:, ; grab the first timestep of the TS field
  b = addfile("","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 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,


ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
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