/[MITgcm]/MITgcm/pkg/mom_fluxform/mom_fluxform.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_fluxform/mom_fluxform.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.26 by jmc, Fri Sep 23 15:19:38 2005 UTC revision 1.27 by baylor, Mon Sep 26 15:27:11 2005 UTC
# Line 361  C     Bi-harmonic term del^2 U -> v4F Line 361  C     Bi-harmonic term del^2 U -> v4F
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
# Line 403  C--   Tendency is minus divergence of th Line 403  C--   Tendency is minus divergence of th
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)
# Line 542  C     Bi-harmonic term del^2 V -> v4F Line 548  C     Bi-harmonic term del^2 V -> v4F
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
# Line 584  C--   Tendency is minus divergence of th Line 590  C--   Tendency is minus divergence of th
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)
Legend:
Removed from v.1.26  
changed lines
  Added in v.1.27
  ViewVC Help
Powered by ViewVC 1.1.22