pkg/generic_advdiff/gad_c2_adv_x.F

31566b6684 Alis*0001 #include "GAD_OPTIONS.h"
                0002
aba3ddc7b7 Alis*0003 CBOP
                0004 C !ROUTINE: GAD_C2_ADV_X
                0005
                0006 C !INTERFACE: ==========================================================
31566b6684 Alis*0007       SUBROUTINE GAD_C2_ADV_X(
                0008      I           bi,bj,k,
                0009      I           uTrans,
                0010      I           tracer,
                0011      O           uT,
                0012      I           myThid )
                0013
aba3ddc7b7 Alis*0014 C !DESCRIPTION:
527a84022c Alis*0015 C Calculates the area integrated zonal flux due to advection of a tracer using
aba3ddc7b7 Alis*0016 C centered second-order interpolation:
                0017 C \begin{equation*}
527a84022c Alis*0018 C F^x_{adv} = U \overline{\theta}^i
aba3ddc7b7 Alis*0019 C \end{equation*}
                0020
                0021 C !USES: ===============================================================
                0022       IMPLICIT NONE
31566b6684 Alis*0023 #include "SIZE.h"
                0024 #include "GRID.h"
                0025
aba3ddc7b7 Alis*0026 C !INPUT PARAMETERS: ===================================================
                0027 C  bi,bj                :: tile indices
                0028 C  k                    :: vertical level
                0029 C  uTrans               :: zonal volume transport
                0030 C  tracer               :: tracer field
                0031 C  myThid               :: thread number
31566b6684 Alis*0032       INTEGER bi,bj,k
                0033       _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
                0034       _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
                0035       INTEGER myThid
                0036
aba3ddc7b7 Alis*0037 C !OUTPUT PARAMETERS: ==================================================
                0038 C  uT                   :: zonal advective flux
                0039       _RL uT    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
                0040
                0041 C !LOCAL VARIABLES: ====================================================
                0042 C  i,j                  :: loop indices
31566b6684 Alis*0043       INTEGER i,j
aba3ddc7b7 Alis*0044 CEOP
31566b6684 Alis*0045
                0046       DO j=1-Oly,sNy+Oly
                0047        uT(1-Olx,j)=0.
360ad14abb Mart*0048       ENDDO
                0049       DO j=1-Oly,sNy+Oly
31566b6684 Alis*0050        DO i=1-Olx+1,sNx+Olx
                0051         uT(i,j) =
                0052      &   uTrans(i,j)*(Tracer(i,j)+Tracer(i-1,j))*0.5 _d 0
                0053        ENDDO
                0054       ENDDO
                0055
                0056       RETURN
                0057       END