Commit f8db550d authored by Willi Rath's avatar Willi Rath

Add ERAIN_SFC12_MM Ferret example by @joke-luebbecke

parent ab9df771
cancel data/all
cancel var/all
use "/data/c2/TMdata/git_geomar_de_data/ERAIN_SFC12_MM/v2.x.x/data/nc/ERAIN_SFC12_MM_IEWS_229_1-456.nc"
use "/data/c2/TMdata/git_geomar_de_data/ERAIN_SFC12_MM/v2.x.x/data/nc/ERAIN_SFC12_MM_INSS_230_1-456.nc"
! load climatological axes
use climatological_axes
cancel data climatological_axes
se region/y=30S:30N/x=60W:20E
let taux = IEWS[d=1]
let tauy = INSS[d=2]
let taux_seasonal_cycle = taux[gt=month_gregorian@mod]
let tauy_seasonal_cycle = tauy[gt=month_gregorian@mod]
let taux_deseasoned = taux - taux_seasonal_cycle
let tauy_deseasoned = tauy - tauy_seasonal_cycle
se window/asp=2/siz=2
! plot mean wind stress
se viewp upper
vector/yskip=2/title="mean wind stress" taux[l=@ave], tauy[l=@ave]
go fland
! plot wind stress anomalies (w.r.t. seasonal cycle) for a given day
set viewp lower
vector/length=3e-2/title="wind stress an. May 1982"/t="15-Mar-1982" taux_deseasoned, tauy_deseasoned
go fland
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