#!/bin/sh

m=ERA5

for iyear in {2019..2019} ; do

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

paste ancors.$m/ans.$n.dat DJETS1 > .tmp
echo $iyear $(gawk '$3==1{x=x+$2;n=n+1}END{print x/n}' .tmp)

done