#!/bin/bash
#SBATCH --account=nn5002k
#SBATCH --job-name=MyJob
#SBATCH --time=1-0:0:0

Y1=2016
Y2=2100

JY=${Y1}

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

    cdo div NAA10KM_5d_${JY}0101_${JY}1231_morten_grid_U3D.nc e3u.nc dum_u_0.nc
    cdo div dum_u_0.nc e2u.nc dum_u_1.nc 
    cdo shiftx,1,cyclic dum_u_1.nc s_dum_u_1.nc 
    cdo add dum_u_1.nc s_dum_u_1.nc dum_u.nc
    cdo divc,2. dum_u.nc u_apoint.nc
    mv u_apoint.nc TPOINT_NAA10KM_5d_${JY}0101_${JY}1231_morten_grid_U3D.nc
    rm -rf *dum*.nc

    cdo div NAA10KM_5d_${JY}0101_${JY}1231_morten_grid_V3D.nc e3v.nc dum_v_0.nc
    cdo div dum_v_0.nc e1v.nc dum_v_1.nc
    cdo shifty,1,cyclic dum_v_1.nc s_dum_v_1.nc
    cdo add dum_v_1.nc s_dum_v_1.nc dum_v.nc
    cdo divc,2. dum_v.nc v_apoint.nc
    mv v_apoint.nc TPOINT_NAA10KM_5d_${JY}0101_${JY}1231_morten_grid_V3D.nc
    rm -rf *dum*.nc

    cdo mul TPOINT_NAA10KM_5d_${JY}0101_${JY}1231_morten_grid_U3D.nc gcost.nc dum_u_0.nc
    cdo mul TPOINT_NAA10KM_5d_${JY}0101_${JY}1231_morten_grid_V3D.nc gsint.nc dum_u_1.nc
    cdo sub dum_u_0.nc dum_u_1.nc u_eastward.nc
    mv u_eastward.nc TPOINT_NAA10KM_5d_${JY}0101_${JY}1231_morten_grid_U_EAST.nc
    rm -rf *dum*.nc

    cdo mul TPOINT_NAA10KM_5d_${JY}0101_${JY}1231_morten_grid_U3D.nc gsint.nc dum_v_0.nc
    cdo mul TPOINT_NAA10KM_5d_${JY}0101_${JY}1231_morten_grid_V3D.nc gcost.nc dum_v_1.nc
    cdo add dum_v_0.nc dum_v_1.nc v_northward.nc
    mv v_northward.nc TPOINT_NAA10KM_5d_${JY}0101_${JY}1231_morten_grid_V_NORTH.nc
    rm -rf *dum*.nc
    
    rm -rf TPOINT*grid_U3D.nc TPOINT*grid_V3D.nc

    JY=`expr ${JY} + 1`

done

exit $?
#