** Warning **

Issuing rollback() due to DESTROY without explicit disconnect() of DBD::mysql::db handle dbname=MITgcm at /usr/local/share/lxr/lib/LXR/Common.pm line 1224.

Last-Modified: Wed, 14 May 2024 05:11:28 GMT Content-Type: text/html; charset=utf-8 MITgcm/MITgcm/pkg/mom_fluxform/mom_u_del2u.F
Back to home page

MITgcm
 
 

    

File indexing completed on 2018-03-02 18:42:13 UTC

view on githubraw file Latest commit add29e06 on 2018-01-31 20:35:05 UTC
6d54cf9ca1 Ed H*0001 #include "MOM_FLUXFORM_OPTIONS.h"
aea29c8517 Alis*0002 
71207ba845 Alis*0003 CBOP
                0004 C !ROUTINE: MOM_U_DEL2U
                0005 
                0006 C !INTERFACE: ==========================================================
aea29c8517 Alis*0007       SUBROUTINE MOM_U_DEL2U(
fd362e9f7c Jean*0008      I        bi, bj, k,
                0009      I        uFld, hFacZ, h0FacZ,
aea29c8517 Alis*0010      O        del2u,
fd362e9f7c Jean*0011      I        myThid )
aea29c8517 Alis*0012 
71207ba845 Alis*0013 C !DESCRIPTION:
                0014 C Calculates the Laplacian of zonal flow
                0015 
                0016 C !USES: ===============================================================
                0017       IMPLICIT NONE
aea29c8517 Alis*0018 #include "SIZE.h"
                0019 #include "EEPARAMS.h"
                0020 #include "PARAMS.h"
                0021 #include "GRID.h"
bf5496678a Jean*0022 #include "SURFACE.h"
aea29c8517 Alis*0023 
71207ba845 Alis*0024 C !INPUT PARAMETERS: ===================================================
                0025 C  bi,bj                :: tile indices
                0026 C  k                    :: vertical level
                0027 C  uFld                 :: zonal flow
                0028 C  hFacZ                :: fractional thickness at vorticity points
fd362e9f7c Jean*0029 C  h0FacZ               :: fixed fractional thickness at vorticity points
                0030 C  myThid               :: my Thread Id number
                0031       INTEGER bi, bj, k
                0032       _RL uFld  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
                0033       _RS hFacZ (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
                0034       _RS h0FacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
aea29c8517 Alis*0035       INTEGER myThid
                0036 
71207ba845 Alis*0037 C !OUTPUT PARAMETERS: ==================================================
                0038 C  del2u                :: Laplacian
                0039       _RL del2u(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
                0040 
                0041 C !LOCAL VARIABLES: ====================================================
                0042 C  i,j                  :: loop indices
aea29c8517 Alis*0043       INTEGER I,J
                0044       _RL fZon(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
                0045       _RL fMer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
fd362e9f7c Jean*0046       _RS hFacZClosedS, hFacZClosedN
71207ba845 Alis*0047 CEOP
aea29c8517 Alis*0048 
                0049 C     Zonal flux d/dx U
fd362e9f7c Jean*0050       DO j=1-OLy+1,sNy+OLy-1
                0051        DO i=1-OLx,sNx+OLx-1
aea29c8517 Alis*0052         fZon(i,j) = drF(k)*_hFacC(i,j,k,bi,bj)
                0053      &   *_dyF(i,j,bi,bj)
                0054      &   *_recip_dxF(i,j,bi,bj)
                0055      &   *(uFld(i+1,j)-uFld(i,j))
                0056 #ifdef COSINEMETH_III
eaba2fd266 Jean*0057      &   *sqCosFacU(J,bi,bj)
aea29c8517 Alis*0058 #endif
c19923126b Jean*0059 #ifdef ALLOW_OBCS
                0060      &   *maskInC(i,j,bi,bj)
                0061 #endif
eaba2fd266 Jean*0062 c    &   *deepFacC(k)        ! dyF scaling factor
                0063 c    &   *recip_deepFacC(k)  ! recip_dxF scaling factor
aea29c8517 Alis*0064        ENDDO
                0065       ENDDO
                0066 
                0067 C     Meridional flux d/dy U
fd362e9f7c Jean*0068       DO j=1-OLy+1,sNy+OLy
                0069        DO i=1-OLx+1,sNx+OLx-1
aea29c8517 Alis*0070         fMer(i,j) = drF(k)*hFacZ(i,j)
                0071      &   *_dxV(i,j,bi,bj)
                0072      &   *_recip_dyU(i,j,bi,bj)
                0073      &   *(uFld(i,j)-uFld(i,j-1))
5bc29fe512 Mart*0074 #if (defined (ISOTROPIC_COS_SCALING) && defined (COSINEMETH_III))
eaba2fd266 Jean*0075      &   *sqCosFacV(J,bi,bj)
5bc29fe512 Mart*0076 #endif
eaba2fd266 Jean*0077 c    &   *deepFacC(k)        ! dxV scaling factor
                0078 c    &   *recip_deepFacC(k)  ! recip_dyU scaling factor
aea29c8517 Alis*0079        ENDDO
                0080       ENDDO
                0081 
                0082 C     del^2 U
fd362e9f7c Jean*0083       DO j=1-OLy+1,sNy+OLy-1
                0084        DO i=1-OLx+1,sNx+OLx-1
aea29c8517 Alis*0085         del2u(i,j) =
eaba2fd266 Jean*0086      &   recip_drF(k)*_recip_hFacW(i,j,k,bi,bj)
                0087      &  *recip_rAw(i,j,bi,bj)*recip_deepFac2C(k)
aea29c8517 Alis*0088      &  *( fZon(i,j  )    - fZon(i-1,j)
                0089      &    +fMer(i,j+1)    - fMer(i  ,j)
                0090      &   )*_maskW(i,j,k,bi,bj)
c19923126b Jean*0091 #ifdef ALLOW_OBCS
                0092      &    *maskInW(i,j,bi,bj)
                0093 #endif
aea29c8517 Alis*0094        ENDDO
                0095       ENDDO
                0096 
                0097       IF (no_slip_sides) THEN
                0098 C-- No-slip BCs impose a drag at walls...
fd362e9f7c Jean*0099        DO j=1-OLy+1,sNy+OLy-1
                0100         DO i=1-OLx+1,sNx+OLx-1
bf5496678a Jean*0101 #ifdef NONLIN_FRSURF
fd362e9f7c Jean*0102          hFacZClosedS = h0FacW(i,j,k,bi,bj) - h0FacZ(i,j)
                0103          hFacZClosedN = h0FacW(i,j,k,bi,bj) - h0FacZ(i,j+1)
bf5496678a Jean*0104 #else
fd362e9f7c Jean*0105          hFacZClosedS = _hFacW(i,j,k,bi,bj) - h0FacZ(i,j)
                0106          hFacZClosedN = _hFacW(i,j,k,bi,bj) - h0FacZ(i,j+1)
bf5496678a Jean*0107 #endif
                0108          del2u(i,j) = del2u(i,j)
                0109      &    -_recip_hFacW(i,j,k,bi,bj)
                0110      &       *recip_rAw(i,j,bi,bj)*recip_deepFac2C(k)
                0111      &       *( hFacZClosedS*dxV(i, j ,bi,bj)
                0112      &                      *_recip_dyU(i, j ,bi,bj)
                0113      &         +hFacZClosedN*dxV(i,j+1,bi,bj)
                0114      &                      *_recip_dyU(i,j+1,bi,bj)
                0115      &        )*uFld(i,j)*sideDragFactor
                0116      &         *_maskW(i,j,k,bi,bj)
                0117         ENDDO
aea29c8517 Alis*0118        ENDDO
                0119       ENDIF
                0120 
                0121       RETURN
                0122       END
[ Source navigation ] [ Diff markup ] [ Identifier search ] [ general search ]

This page was automatically generated from https://github.com/MITgcm/MITgcm by the 2.2.1-MITgcm-0.1 LXR engine.

The LXR team
Valid CSS 2.1! Valid HTML 4.01!