Hi Steeven,
Another option would be to make a mask from the original data, which is missing in the locations that you want to keep missing
yes? let mask = if temp then 1
Use MISSING to fill in your computed variable, and then apply this mask.
Paul brings up using SHADE. Often the best plot is a SHADE, overlaid with a FILL/NOLABEL.
-Ansley
On 2/3/2015 1:44 AM, Young, Paul wrote:
Hi Steeven,
One way to do this would be to define a land/ocean mask using the etopo files that come with ferret. E.g.
use etopo40let isOcean = if rose[d=1] le 0 then 1
You could then combine this with your own calculation below just to get the ocean region you want. I would use “shade” to fill in the missing data here.
Hope this helps,
Paul
On 3 Feb 2015, at 09:11, Steeven Paul Yerraguntla <steevenpaul@xxxxxxxxx> wrote:
<HCR.gif>Dear ferret users,
I calculated Heat content ratio(HCR) between top 50 meters heat content(HC50) and HC above 20 degrees isotherm (HC_D20). HCR is infinity ( when HC_D20 is zero or not available ie,. 20 deg isotherm reaches surface) which causes big white patches in the fill plot. I want to fill those white patche (data gap) with a value 1 without filling other missing data in coastal regions. MISSING funtion is useful but it is filling other missing data also (for example land region). Please let me know how can i make a logic to achieve it. Ferret script lines which i used for plotting are as follows. Please see attached example figure.
use temp.cdf
let d20 = temp[z=@loc:20]
let mask_d20 = IF z[g=temp] LT d20 THEN 1
let temp_d20 = mask_d20*temp
let HCR= HC50/HC_D20
fill/l=1/levels=(-inf)(0.2,1,0.1)(inf) HCR