#!/bin/sh

n=6377
m=6585

rm EVENTS/PERS* EVENTS/CASES* EVENTS/DAYS*

gawk '$2>0.5{y=1950+int((FNR-1)/92);print y,1}$2<=0.5{y=1950+int((FNR-1)/92);print y,0}' ancors.ERA5/ans.$n.dat > EVENTS/EVENTS1
for y in {1950..2022} ; do
  gawk '$1=='$y'&&$2==1{p=p+1}$1=='$y'&&$2==0&&p>0{q=q+p;n=n+1;p=0}p>0&&FNR==6716{p=p+1;q=q+p;n=n+1}END{print '$y',n+0}' EVENTS/EVENTS1 >> EVENTS/CASES1
  gawk '$1=='$y'&&$2==1{p=p+1}$1=='$y'&&$2==0&&p>0{q=q+p;n=n+1;p=0}p>0&&FNR==6716{p=p+1;q=q+p;n=n+1}END{print '$y',q/(n+1e-8)}' EVENTS/EVENTS1 >> EVENTS/PERS1
  gawk '$1=='$y'{p=p+$2}END{print '$y',p+0}' EVENTS/EVENTS1 >> EVENTS/DAYS1
done

gawk '$2>0.5{y=1950+int((FNR-1)/92);print y,1}$2<=0.5{y=1950+int((FNR-1)/92);print y,0}' ancors.ERA5/ans.$m.dat > EVENTS/EVENTS2
for y in {1950..2022} ; do
  gawk '$1=='$y'&&$2==1{p=p+1}$1=='$y'&&$2==0&&p>0{q=q+p;n=n+1;p=0}p>0&&FNR==6716{p=p+1;q=q+p;n=n+1}END{print '$y',n+0}' EVENTS/EVENTS2 >> EVENTS/CASES2
  gawk '$1=='$y'&&$2==1{p=p+1}$1=='$y'&&$2==0&&p>0{q=q+p;n=n+1;p=0}p>0&&FNR==6716{p=p+1;q=q+p;n=n+1}END{print '$y',q/(n+1e-8)}' EVENTS/EVENTS2 >> EVENTS/PERS2
  gawk '$1=='$y'{p=p+$2}END{print '$y',p+0}' EVENTS/EVENTS2 >> EVENTS/DAYS2
done

echo CASES EVENTS1 : $(gawk -f calatrend.awk EVENTS/CASES1)
echo PERS. EVENTS1 : $(gawk -f calatrend.awk EVENTS/PERS1)
echo DAYS. EVENTS1 : $(gawk -f calatrend.awk EVENTS/DAYS1)

echo CASES EVENTS2 : $(gawk -f calatrend.awk EVENTS/CASES2)
echo PERS. EVENTS2 : $(gawk -f calatrend.awk EVENTS/PERS2)
echo DAYS. EVENTS2 : $(gawk -f calatrend.awk EVENTS/DAYS2)