Re: [ferret_users] Define data mask valid in all time space

["Ansley C. Manke" <ansley.b.manke@xxxxxxxx>]

Hi,

You don't need to use a loop. You can define a mask in XY that shows where all the time data id present. One method would be to define a test based on the number of valid points in the T direction:

let nt = `var,return=lsize`
let var_mask = if var[L=@ngd] EQ nt then 1

On 6/18/2018 11:27 AM, Ge Peng - NOAA Affiliate wrote:

Dear Ferreters,

 

I have arrays of x, y, and t of a regular grid (idim,jdim,ldim).

 

For each array, I need to define a mask for only grid cells that contain valid data in all t, i.e., common valid data grid cells though out t dimension.

 

I am thinking about a repeat loop to go through t dimension. Just wonder if anyone has a ferret script to do that - perhaps in a more efficient or elegant way, or have a suggestion.

 

Thanks a lot,

 

--- Peng

--

Ge Peng, PhD
Research Scholar
Cooperative Institute for Climate and Satellites - NC (CICS-NC)/NCSU at

NOAA’s National Centers for Environmental Information (NCEI)

Center for Weather and Climate (CWC)

151 Patton Ave, Asheville, NC 28801
+1 828 257 3009; ge.peng@xxxxxxxx

ORCID: http://orcid.org/0000-0002-1986-9115

Following CICS-NC on Facebook