#!/bin/sh
# Computes significance for dynamical trends (assumes zanalogues.f90 has been applied)
# Check consistency between xanalogues and zanalogues
# Basis for Supp Table 1

mlist="ERA5 EOBS ACCESS-CM2_r1i1p1f1.hssp585 ACCESS-ESM1-5_r2i1p1f1.hssp585 CAMS-CSM1-0_r2i1p1f1.hssp585 CanESM5_r1i1p1f1.hssp585 CMCC-ESM2_r1i1p1f1.hssp585 CNRM-CM6-1-HR_r1i1p1f2.hssp585 CNRM-CM6-1_r1i1p1f2.hssp585 CNRM-ESM2-1_r1i1p1f2.hssp585 EC-Earth3_r1i1p1f1.hssp585 EC-Earth3-CC_r1i1p1f1.hssp585 EC-Earth3-Veg_r1i1p1f1.hssp585 EC-Earth3-Veg-LR_r1i1p1f1.hssp585 FGOALS-g3_r1i1p1f1.hssp585 GFDL-CM4_r1i1p1f1.hssp585 GISS-E2-1-G_r1i1p1f2.hssp585 HadGEM3-GC31-LL_r1i1p1f3.hssp585 HadGEM3-GC31-MM_r1i1p1f3.hssp585 INM-CM4-8_r1i1p1f1.hssp585 INM-CM5-0_r1i1p1f1.hssp585 IPSL-CM6A-LR_r1i1p1f1.hssp585 KACE-1-0-G_r1i1p1f1.hssp585 KIOST-ESM_r1i1p1f1.hssp585 MIROC-ES2L_r1i1p1f2.hssp585 MIROC6_r1i1p1f1.hssp585 MPI-ESM1-2-HR_r1i1p1f1.hssp585 MPI-ESM1-2-LR_r1i1p1f1.hssp585 MRI-ESM2-0_r1i1p1f1.hssp585 NESM3_r1i1p1f1.hssp585 NorESM2-LM_r1i1p1f1.hssp585 NorESM2-MM_r1i1p1f1.hssp585 TaiESM1_r1i1p1f1.hssp585 UKESM1-0-LL_r1i1p1f2.hssp585"
mlist="ERA5"

yy=$1

for m in $mlist; do

f=TRENDS/txmx.$m.yearly.dat
g=PATTRENDS/txmx.$m.yearly.$yy.dat

gawk '{print $3}' $f | cdo -s -f nc input,GRID_TAS tmp/t$m.nc ; cdo -s setvrange,-100,100 tmp/t$m.nc TRENDS/txm.$m.nc
gawk '{print $4}' $f | cdo -s -f nc input,GRID_TAS tmp/t$m.nc ; cdo -s setvrange,-100,100 tmp/t$m.nc TRENDS/txx.$m.nc
gawk '{print $5}' $f | cdo -s -f nc input,GRID_TAS tmp/t$m.nc ; cdo -s setvrange,-100,100 tmp/t$m.nc TRENDS/txmd.$m.nc
gawk '{print $6}' $f | cdo -s -f nc input,GRID_TAS tmp/t$m.nc ; cdo -s setvrange,-100,100 tmp/t$m.nc TRENDS/txxd.$m.nc
gawk '{print $3}' $g | cdo -s -f nc input,GRID_TAS tmp/t$m.nc ; cdo -s setvrange,-100,100 tmp/t$m.nc TRENDS/txmdex.$m.nc
gawk '{print $4}' $g | cdo -s -f nc input,GRID_TAS tmp/t$m.nc ; cdo -s setvrange,-100,100 tmp/t$m.nc TRENDS/txxdex.$m.nc
gawk '{print $5}' $g | cdo -s -f nc input,GRID_TAS tmp/t$m.nc ; cdo -s setvrange,-100,100 tmp/t$m.nc TRENDS/txmdin.$m.nc
gawk '{print $6}' $g | cdo -s -f nc input,GRID_TAS tmp/t$m.nc ; cdo -s setvrange,-100,100 tmp/t$m.nc TRENDS/txxdin.$m.nc

cdo -s outputtab,value -fldmean -mul box_mask.nc -remapbil,box_mask.nc TRENDS/txmd.$m.nc | gawk 'FNR>1{print 1948+FNR,$1}' > txm.dyn.box.dat
cdo -s outputtab,value -fldmean -mul box_mask.nc -remapbil,box_mask.nc TRENDS/txxd.$m.nc | gawk 'FNR>1{print 1948+FNR,$1}' > txx.dyn.box.dat
cdo -s outputtab,value -fldmean -mul box_mask.nc -remapbil,box_mask.nc TRENDS/txm.$m.nc | gawk 'FNR>1{print 1948+FNR,$1}' > txm.box.dat
cdo -s outputtab,value -fldmean -mul box_mask.nc -remapbil,box_mask.nc TRENDS/txx.$m.nc | gawk 'FNR>1{print 1948+FNR,$1}' > txx.box.dat
cdo -s outputtab,value -fldmean -mul box_mask.nc -remapbil,box_mask.nc TRENDS/txmdex.$m.nc | gawk 'FNR>1{print 1948+FNR,$1}' > txm.dyn.box.ex.dat
cdo -s outputtab,value -fldmean -mul box_mask.nc -remapbil,box_mask.nc TRENDS/txxdex.$m.nc | gawk 'FNR>1{print 1948+FNR,$1}' > txx.dyn.box.ex.dat
cdo -s outputtab,value -fldmean -mul box_mask.nc -remapbil,box_mask.nc TRENDS/txmdin.$m.nc | gawk 'FNR>1&&$1>-100{print 1948+FNR,$1}' > txm.dyn.box.in.dat
cdo -s outputtab,value -fldmean -mul box_mask.nc -remapbil,box_mask.nc TRENDS/txxdin.$m.nc | gawk 'FNR>1&&$1>-100{print 1948+FNR,$1}' > txx.dyn.box.in.dat

done

xmg -param exinpat.par txx.box.dat txx.dyn.box.dat txx.dyn.box.ex.dat txx.dyn.box.in.dat &
# xmg txm.box.dat txm.dyn.box.dat txm.dyn.box.ex.dat txm.dyn.box.in.dat &