#!/bin/sh

n=$(gawk '$1=='2019'{print $2}' MAXTIMES)

mlist="ERA5 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"

for m in $mlist ; do
  gawk '$2>0.5{print 1950+int((FNR-1)/92),1}$2<=0.5{print 1950+int((FNR-1)/92),0}' ancors.$m/ane.$n.dat > EVENTS/events.$m
  rm EVENTS/pers.$m
  for y in {1950..2022} ; do
    gawk '$1=='$y'&&$2==1{p=p+1}$1=='$y'&&$2==0&&p>0{print '$y',p;p=0}$1=='$y'&&p>0&&$1==2022{p=p+1;q=q+p;n=n+1;print '$y',p}' EVENTS/events.$m >> EVENTS/pers.$m
#   gawk '$1=='$y'&&$2==1{p=p+1}$1=='$y'&&$2==0&&p>0{q=q+p;n=n+1;p=0}END{if(p>0){p=p+1;q=q+p;n=n+1};print '$y',q/(n+1e-8)}' EVENTS/events.$m >> EVENTS/pers.$m
  done
  gawk '{print $2}' covariates/SGSAT1950-2022.$m.dat > .tmp
  echo $m $(paste .tmp EVENTS/pers.$m | gawk '{print $1,$3}' | gawk -f caltrendp.awk)
done