#!/bin/bash


rm -f *.dat
for YYYY in `seq 1940 2023` ; do

 echo ${YYYY}

ferret<<EOF
use MARv3.14-monthly-ERA5-${YYYY}.nc
set mem/size=10000
let SMB3= if(MSK gt 2.5 and SMBcorr ge -100000) then SMBcorr else 0
let SF3 = if(MSK gt 2.5) then SF else 0
let RF3 = if(MSK gt 2.5) then RF else 0
let RU3 = if(MSK gt 2.5 and MEcorr ge -100000) then RUcorr else 0
let ME3 = if(MSK gt 2.5 and MEcorr ge -100000) then MEcorr else 0
list/nohead/format=(f9.1)/file=smb.dat/append SMB3[l=@sum]
list/nohead/format=(f9.1)/file=sf.dat/append sf3[l=@sum]
list/nohead/format=(f9.1)/file=rf.dat/append RF3[l=@sum]
list/nohead/format=(f9.1)/file=ru.dat/append RU3[l=@sum]
list/nohead/format=(f9.1)/file=me.dat/append ME3[l=@sum]
quit
EOF

done

ferret <<EOF
use MARv3.14-monthly-ERA5-1940.nc
set LIST/NCFORMAT=4/DEFLATE=1/SHUFFLE
define axis/t=1940:2023:1 time2
define grid/x=lat/y=lat/t=time2 gg
set mem/size=20000
file/var=SMB/col=1/grid=gg smb.dat
file/var=SF/col=1/grid=gg sf.dat
file/var=RF/col=1/grid=gg rf.dat
file/var=RU/col=1/grid=gg ru.dat
file/var=ME/col=1/grid=gg me.dat
save/file=MARv3.14-ERA5-1940-2023.nc/clobber/d=1 lon,lat,msk,srf
save/file=MARv3.14-ERA5-1940-2023.nc/append SMB[d=2]
save/file=MARv3.14-ERA5-1940-2023.nc/append SF[d=3]
save/file=MARv3.14-ERA5-1940-2023.nc/append RF[d=4]
save/file=MARv3.14-ERA5-1940-2023.nc/append RU[d=5]
save/file=MARv3.14-ERA5-1940-2023.nc/append ME[d=6]
quit
EOF