! This PyFerret script creates a pipe mask for the wMask UVic code.

! Written by M. Juerchott (mjuerchott@geomar.de) on the 25.10.2023.



! ***  Only run it in -nodisplay -nojnl mode ***


! Save current state

        define REGION/DEFAULT "save"
        set data/save
        set grid/save

! Clean everything

        cancel data/all
        cancel variable/all
        cancel region/all



! Script function

	use  /Users/mjuerchott/runs/2023/2023_10_23_FeDyn_wLand_noAU_restart_files_1765_2100_COU85/tavg.01775.01.01.nc	! [d=1] 



! pipe area

	let OMZ = O_O2[k=5,l=33]*1000
	let mask0 = if (G_KMT gt 4.5) and (y ge -60) and (y le 60) and OMZ ge 50 then 5
	let mask1 = mask0 + 0*o_po4 	! 2D to 3D conversion
	let mask2 = mask1

! plot routine

	cancel windows/all
	set windows/new
	set text/font='Times New Roman'
	set text/isiz=4


	go landscape2x2.jnl
	define view/xlimits=0.0,0.5/ylimits=0.05,0.55 ll4
	define view/xlimits=0.0,0.5/ylimits=0.5,1 ul4
	define view/xlimits=0.5,1/ylimits=0.05,0.55 lr4
        define view/xlimits=0.1,0.9/ylimits=0.1,0.9 ur4


	set view ur4
	shade/nolabel/pal=ocean_blue/level=(2.5,7.5,1)/k=1/l=1 mask1
	go land

	label/nouser 3.5,5.8,0,0,0.15 "Fe_Dyn, Pipe Mask, 60˚N-S, no OMZ, K >= 5" 	! OMZ < 50 mmol m-3 at k=5 




! rename mask to O_mask.nc and copy it into UVic data folder

frame/file = "S1_FIG.pdf"

save/clobber/file=tmp.nc mask2[l=1]
can va mask1, mask2
can da 1
use tmp.nc
let/title="AU_mask" mask=mask2
save/clobber/file=make_mask_60deg_NS_500m_no_OMZ.nc MASK




! Cancel everything

        cancel data/all
        cancel variable/all
        cancel region/all

! Restore current state

        set data/restore
        set grid/restore