361 |
|
|
362 |
C Laplacian and bi-harmonic terms, Zonal Fluxes -> fZon |
C Laplacian and bi-harmonic terms, Zonal Fluxes -> fZon |
363 |
CALL MOM_U_XVISCFLUX(bi,bj,k,uFld,v4F,fZon, |
CALL MOM_U_XVISCFLUX(bi,bj,k,uFld,v4F,fZon, |
364 |
I viscAh_Z,viscAh_D,viscA4_Z,viscA4_D,myThid) |
I viscAh_D,viscA4_D,myThid) |
365 |
|
|
366 |
C Laplacian and bi-harmonic termis, Merid Fluxes -> fMer |
C Laplacian and bi-harmonic termis, Merid Fluxes -> fMer |
367 |
CALL MOM_U_YVISCFLUX(bi,bj,k,uFld,v4F,hFacZ,fMer, |
CALL MOM_U_YVISCFLUX(bi,bj,k,uFld,v4F,hFacZ,fMer, |
368 |
I viscAh_Z,viscAh_D,viscA4_Z,viscA4_D,myThid) |
I viscAh_Z,viscA4_Z,myThid) |
369 |
|
|
370 |
C Eddy component of vertical flux (interior component only) -> fVrUp & fVrDw |
C Eddy component of vertical flux (interior component only) -> fVrUp & fVrDw |
371 |
IF (.NOT.implicitViscosity) THEN |
IF (.NOT.implicitViscosity) THEN |
403 |
C-- No-slip and drag BCs appear as body forces in cell abutting topography |
C-- No-slip and drag BCs appear as body forces in cell abutting topography |
404 |
IF (no_slip_sides) THEN |
IF (no_slip_sides) THEN |
405 |
C- No-slip BCs impose a drag at walls... |
C- No-slip BCs impose a drag at walls... |
406 |
CALL MOM_U_SIDEDRAG(bi,bj,k,uFld,v4F,hFacZ,vF,myThid) |
CALL MOM_U_SIDEDRAG( |
407 |
|
I bi,bj,k, |
408 |
|
I uFld, v4f, hFacZ, |
409 |
|
I viscAh_Z,viscA4_Z, |
410 |
|
I harmonic,biharmonic,useVariableViscosity, |
411 |
|
O vF, |
412 |
|
I myThid) |
413 |
DO j=jMin,jMax |
DO j=jMin,jMax |
414 |
DO i=iMin,iMax |
DO i=iMin,iMax |
415 |
gUdiss(i,j) = gUdiss(i,j) + vF(i,j) |
gUdiss(i,j) = gUdiss(i,j) + vF(i,j) |
548 |
|
|
549 |
C Laplacian and bi-harmonic terms, Zonal Fluxes -> fZon |
C Laplacian and bi-harmonic terms, Zonal Fluxes -> fZon |
550 |
CALL MOM_V_XVISCFLUX(bi,bj,k,vFld,v4f,hFacZ,fZon, |
CALL MOM_V_XVISCFLUX(bi,bj,k,vFld,v4f,hFacZ,fZon, |
551 |
I viscAh_Z,viscAh_D,viscA4_Z,viscA4_D,myThid) |
I viscAh_Z,viscA4_Z,myThid) |
552 |
|
|
553 |
C Laplacian and bi-harmonic termis, Merid Fluxes -> fMer |
C Laplacian and bi-harmonic termis, Merid Fluxes -> fMer |
554 |
CALL MOM_V_YVISCFLUX(bi,bj,k,vFld,v4f,fMer, |
CALL MOM_V_YVISCFLUX(bi,bj,k,vFld,v4f,fMer, |
555 |
I viscAh_Z,viscAh_D,viscA4_Z,viscA4_D,myThid) |
I viscAh_D,viscA4_D,myThid) |
556 |
|
|
557 |
C Eddy component of vertical flux (interior component only) -> fVrUp & fVrDw |
C Eddy component of vertical flux (interior component only) -> fVrUp & fVrDw |
558 |
IF (.NOT.implicitViscosity) THEN |
IF (.NOT.implicitViscosity) THEN |
590 |
C-- No-slip and drag BCs appear as body forces in cell abutting topography |
C-- No-slip and drag BCs appear as body forces in cell abutting topography |
591 |
IF (no_slip_sides) THEN |
IF (no_slip_sides) THEN |
592 |
C- No-slip BCs impose a drag at walls... |
C- No-slip BCs impose a drag at walls... |
593 |
CALL MOM_V_SIDEDRAG(bi,bj,k,vFld,v4F,hFacZ,vF,myThid) |
CALL MOM_V_SIDEDRAG( |
594 |
|
I bi,bj,k, |
595 |
|
I vFld, v4f, hFacZ, |
596 |
|
I viscAh_Z,viscA4_Z, |
597 |
|
I harmonic,biharmonic,useVariableViscosity, |
598 |
|
O vF, |
599 |
|
I myThid) |
600 |
DO j=jMin,jMax |
DO j=jMin,jMax |
601 |
DO i=iMin,iMax |
DO i=iMin,iMax |
602 |
gvDiss(i,j) = gvDiss(i,j) + vF(i,j) |
gvDiss(i,j) = gvDiss(i,j) + vF(i,j) |