2 |
C $Name$ |
C $Name$ |
3 |
|
|
4 |
#include "MOM_COMMON_OPTIONS.h" |
#include "MOM_COMMON_OPTIONS.h" |
|
#define SIDEDRAG_OLD_VERSION |
|
5 |
|
|
6 |
CBOP |
CBOP |
7 |
C !ROUTINE: MOM_V_SIDEDRAG |
C !ROUTINE: MOM_V_SIDEDRAG |
53 |
C i,j :: loop indices |
C i,j :: loop indices |
54 |
C hFacZClosedE :: fractional open water to east |
C hFacZClosedE :: fractional open water to east |
55 |
C hFacZClosedW :: fractional open water to west |
C hFacZClosedW :: fractional open water to west |
56 |
INTEGER I,J |
INTEGER i,j |
57 |
_RL hFacZClosedE,hFacZClosedW |
_RL hFacZClosedE,hFacZClosedW |
58 |
CEOP |
CEOP |
|
|
|
|
#ifdef SIDEDRAG_OLD_VERSION |
|
59 |
_RL Ahtmp,A4tmp |
_RL Ahtmp,A4tmp |
60 |
|
|
61 |
C - Laplacian and bi-harmonic terms |
|
62 |
|
IF ( sideDragFactor.LE.0. ) THEN |
63 |
|
C-- Laplacian and bi-harmonic terms: variable-Viscosity coeff. from |
64 |
|
C MOM_CALC_VISC are not used here (corresponds to the "old" version). |
65 |
|
|
66 |
DO j=2-Oly,sNy+Oly-1 |
DO j=2-Oly,sNy+Oly-1 |
67 |
DO i=2-Olx,sNx+Olx-1 |
DO i=2-Olx,sNx+Olx-1 |
68 |
#ifdef NONLIN_FRSURF |
#ifdef NONLIN_FRSURF |
92 |
& +hFacZClosedE*_dyU(i+1,j,bi,bj) |
& +hFacZClosedE*_dyU(i+1,j,bi,bj) |
93 |
& *_recip_dxV(i+1,j,bi,bj) ) |
& *_recip_dxV(i+1,j,bi,bj) ) |
94 |
& *drF(k)*2.*( |
& *drF(k)*2.*( |
95 |
& Ahtmp*vFld(i,j)*cosFacV(J,bi,bj) |
& Ahtmp*vFld(i,j)*cosFacV(j,bi,bj) |
96 |
#ifdef COSINEMETH_III |
#ifdef COSINEMETH_III |
97 |
& -viscA4*del2v(i,j)*sqcosFacV(J,bi,bj) |
& -viscA4*del2v(i,j)*sqcosFacV(j,bi,bj) |
98 |
#else |
#else |
99 |
& -A4tmp*del2v(i,j)*cosFacV(J,bi,bj) |
& -A4tmp*del2v(i,j)*cosFacV(j,bi,bj) |
100 |
#endif |
#endif |
101 |
& ) |
& ) |
102 |
ENDDO |
ENDDO |
103 |
ENDDO |
ENDDO |
104 |
|
|
105 |
#else /* SIDEDRAG_OLD_VERSION */ |
ELSE |
106 |
|
C-- Laplacian and bi-harmonic terms: using variable-Viscosity coeff. |
107 |
|
C from MOM_CALC_VISC, consistent with dissipation in the interior |
108 |
|
|
|
C - Laplacian and bi-harmonic terms |
|
109 |
DO j=2-Oly,sNy+Oly-1 |
DO j=2-Oly,sNy+Oly-1 |
110 |
DO i=2-Olx,sNx+Olx-1 |
DO i=2-Olx,sNx+Olx-1 |
111 |
C- this will not give any side-drag along thin wall. |
C- this will not give any side-drag along thin wall. |
123 |
& -_recip_hFacS(i,j,k,bi,bj) |
& -_recip_hFacS(i,j,k,bi,bj) |
124 |
& *recip_drF(k)*recip_rAs(i,j,bi,bj) |
& *recip_drF(k)*recip_rAs(i,j,bi,bj) |
125 |
& *( hFacZClosedW*_dyU( i ,j,bi,bj)*_recip_dxV( i ,j,bi,bj) |
& *( hFacZClosedW*_dyU( i ,j,bi,bj)*_recip_dxV( i ,j,bi,bj) |
126 |
& *(viscAh_Z(i ,j)*vFld(i,j)*cosFacV(J,bi,bj) |
& *( viscAh_Z(i ,j)*vFld(i,j)*cosFacV(j,bi,bj) |
127 |
#ifdef COSINEMETH_III |
#ifdef COSINEMETH_III |
128 |
& -viscA4_Z(i ,j)*del2v(i,j)*sqcosFacV(J,bi,bj)) |
& -viscA4_Z(i ,j)*del2v(i,j)*sqcosFacV(j,bi,bj) ) |
129 |
#else |
#else |
130 |
& -viscA4_Z(i ,j)*del2v(i,j)*cosFacV(J,bi,bj)) |
& -viscA4_Z(i ,j)*del2v(i,j)*cosFacV(j,bi,bj) ) |
131 |
#endif |
#endif |
132 |
& +hFacZClosedE*_dyU(i+1,j,bi,bj)*_recip_dxV(i+1,j,bi,bj) |
& +hFacZClosedE*_dyU(i+1,j,bi,bj)*_recip_dxV(i+1,j,bi,bj) |
133 |
& *(viscAh_Z(i+1,j)*vFld(i,j)*cosFacV(J,bi,bj) |
& *( viscAh_Z(i+1,j)*vFld(i,j)*cosFacV(j,bi,bj) |
134 |
#ifdef COSINEMETH_III |
#ifdef COSINEMETH_III |
135 |
& -viscA4_Z(i+1,j)*del2v(i,j)*sqcosFacV(J,bi,bj)) |
& -viscA4_Z(i+1,j)*del2v(i,j)*sqcosFacV(j,bi,bj) ) |
136 |
#else |
#else |
137 |
& -viscA4_Z(i+1,j)*del2v(i,j)*cosFacV(J,bi,bj)) |
& -viscA4_Z(i+1,j)*del2v(i,j)*cosFacV(j,bi,bj) ) |
138 |
#endif |
#endif |
139 |
& ) |
& )*drF(k)*sideDragFactor |
|
& *drF(k)*2. |
|
140 |
ENDDO |
ENDDO |
141 |
ENDDO |
ENDDO |
142 |
|
|
143 |
#endif /* (not) SIDEDRAG_OLD_VERSION */ |
C-- end old-version / new-version IF blocks |
144 |
|
ENDIF |
145 |
|
|
146 |
#ifdef ALLOW_DIAGNOSTICS |
#ifdef ALLOW_DIAGNOSTICS |
147 |
IF (useDiagnostics) THEN |
IF (useDiagnostics) THEN |