/[MITgcm]/MITgcm_contrib/natl_12/code/mom_u_sidedrag.F
ViewVC logotype

Diff of /MITgcm_contrib/natl_12/code/mom_u_sidedrag.F

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

revision 1.1 by cnh, Tue Aug 5 21:22:44 2003 UTC revision 1.2 by cnh, Thu Aug 7 13:03:31 2003 UTC
# Line 49  C  hFacZClosedN         :: fractional op Line 49  C  hFacZClosedN         :: fractional op
49  C  hFacZClosedS         :: fractional open water to south  C  hFacZClosedS         :: fractional open water to south
50        INTEGER I,J        INTEGER I,J
51        _RS hFacZClosedS,hFacZClosedN        _RS hFacZClosedS,hFacZClosedN
52          _RS gridScalingAh
53        _RS gridScalingA4        _RS gridScalingA4
54  CEOP  CEOP
55    
# Line 58  C     - Laplacian  and bi-harmonic terms Line 59  C     - Laplacian  and bi-harmonic terms
59          hFacZClosedS = _hFacW(i,j,k,bi,bj) - hFacZ(i,j)          hFacZClosedS = _hFacW(i,j,k,bi,bj) - hFacZ(i,j)
60          hFacZClosedN = _hFacW(i,j,k,bi,bj) - hFacZ(i,j+1)          hFacZClosedN = _hFacW(i,j,k,bi,bj) - hFacZ(i,j+1)
61          gridScalingA4 = (_dyU(i, j ,bi,bj)**4)/(5000.**4)   ! Realy should split this loop for a more precise formulation          gridScalingA4 = (_dyU(i, j ,bi,bj)**4)/(5000.**4)   ! Realy should split this loop for a more precise formulation
62            gridScalingAh = (_dyU(i, j ,bi,bj)**2)/(5000.**2)   ! Realy should split this loop for a more precise formulation
63          uDragTerms(i,j) =          uDragTerms(i,j) =
64       &   -_recip_hFacW(i,j,k,bi,bj)       &   -_recip_hFacW(i,j,k,bi,bj)
65       &   *recip_drF(k)*recip_rAw(i,j,bi,bj)       &   *recip_drF(k)*recip_rAw(i,j,bi,bj)
# Line 74  C     - Laplacian  and bi-harmonic terms Line 76  C     - Laplacian  and bi-harmonic terms
76       &               -viscA4*del2u(i,j)*cosFacU(J,bi,bj)       &               -viscA4*del2u(i,j)*cosFacU(J,bi,bj)
77  #endif  #endif
78  #else  #else
79       &                viscAh*uFld(i,j)       &                viscAh*uFld(i,j)*gridScalingAh
80       &               -viscA4*del2u(i,j)*gridScalingA4       &               -viscA4*del2u(i,j)*gridScalingA4
81  #endif  #endif
82       &              )       &              )
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22