#!/bin/bash
#SBATCH --account=nn5002k --job-name=diag
#SBATCH --partition=bigmem
#SBATCH --time=0-2:0:0
#SBATCH --nodes=1 --ntasks-per-node=1 --cpus-per-task=1
#SBATCH --mem-per-cpu=36G

Y1=2015
Y2=2050

JY=${Y1}


while [ ${JY} -le ${Y2} ]; do

    ncks -d x,143,773 -d y,13,513  -v  vocetr_eff NAA10KM_1m_${JY}0101_${JY}1231_grid_V.nc dum.nc
    ncra dum.nc AMOC_NAA10KM_1y_${JY}0101_${JY}1231_grid_V.nc
    ncap2 -s 'mflux=vocetr_eff.total($x)' AMOC_NAA10KM_1y_${JY}0101_${JY}1231_grid_V.nc  out.nc
    ncks -v mflux out.nc MFLUX_AMOC_NAA10KM_1y_${JY}0101_${JY}1231_grid_V.nc 
    rm -rf dum.nc out.nc

    JY=`expr ${JY} + 1`

done