Commit 9839dce4 authored by Willi Rath's avatar Willi Rath

Fix land-mask creation and use transparent land overlay

parent ab52522e
...@@ -11,10 +11,10 @@ cancel data climatological_axes ...@@ -11,10 +11,10 @@ cancel data climatological_axes
se region/y=30S:30N/x=60W:20E se region/y=30S:30N/x=60W:20E
let ocean_mask = if LSM[d=3, l=@max] == 0 then 1 let ocean_mask = if LSM[d=3, l=@max] eq 0 then 1
let taux = IEWS[d=1] * LSM[gxy = IEWS[d=1]@asn] let taux = IEWS[d=1] * ocean_mask
let tauy = INSS[d=2] * LSM[gxy = INSS[d=2]@asn] let tauy = INSS[d=2] * ocean_mask
let taux_seasonal_cycle = taux[gt=month_gregorian@mod] let taux_seasonal_cycle = taux[gt=month_gregorian@mod]
let tauy_seasonal_cycle = tauy[gt=month_gregorian@mod] let tauy_seasonal_cycle = tauy[gt=month_gregorian@mod]
...@@ -27,11 +27,11 @@ se window/asp=2/siz=2 ...@@ -27,11 +27,11 @@ se window/asp=2/siz=2
! plot mean wind stress ! plot mean wind stress
se viewp upper se viewp upper
vector/yskip=2/title="mean wind stress" taux[l=@ave], tauy[l=@ave] vector/yskip=2/title="mean wind stress" taux[l=@ave], tauy[l=@ave]
go fland go land
! plot wind stress anomalies (w.r.t. seasonal cycle) for a given day ! plot wind stress anomalies (w.r.t. seasonal cycle) for a given day
set viewp lower set viewp lower
vector/length=3e-2/title="wind stress an. May 1982"/t="15-Mar-1982" taux_deseasoned, tauy_deseasoned vector/length=3e-2/title="wind stress an. May 1982"/t="15-Mar-1982" taux_deseasoned, tauy_deseasoned
go fland go land
frame/file="example_08_ERAIN_SFC12_MM_deseasoned_Tropical_Atlantic_wind_stress.png" frame/file="example_08_ERAIN_SFC12_MM_deseasoned_Tropical_Atlantic_wind_stress.png"
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment