[gradsusr] Using maskout function to apply TWO masks

[Roberto Mera]

Grads crew:

I would like to get the area average for the temperature for the California Central Valley. I would also like to mask out the ocean so I don't get those temperatures when I take the average. Right now I have a script that masks out the ocean and I have also adapted it for other purposes.

My question is this: How do I apply two masks?

I want to mask out the ocean and also, say >700 m altitude.

Right now my script loops through ensemble members and calculates the time average within each ensemble member and it gives me a number for the area average:

*mskgrd=maskout(ht/ht,ht-0.1) ****this masks out the ocean
count = 1
while (count < 18)
'set e 'count
'areal=aave(maskout(ave(field88,t=1,t=6),mskgrd(t=1)),x=147,x=250,y=94.1815,y=182.057)'
'd areal'
areal=subwrd(result,4)
say areal
if (rc != 0) ; break ; endif
count = count +1
endwhile

Robert



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140603/e28a2163/attachment.html