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

Annotation of /MITgcm_contrib/natl_12/code/mom_u_yviscflux.F

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


Revision 1.1 - (hide annotations) (download)
Tue Aug 5 21:22:44 2003 UTC (21 years, 11 months ago) by cnh
Branch: MAIN
Adding set of files for 1/12 Atlantic configuration

1 cnh 1.1 C $Header: /u/u0/gcmpack/MITgcm/pkg/mom_fluxform/mom_u_yviscflux.F,v 1.3 2001/09/26 19:05:21 adcroft Exp $
2     C $Name: $
3    
4     #include "CPP_OPTIONS.h"
5    
6     CBOP
7     C !ROUTINE: MOM_U_YVISCFLUX
8    
9     C !INTERFACE: ==========================================================
10     SUBROUTINE MOM_U_YVISCFLUX(
11     I bi,bj,k,
12     I uFld, del2u, hFacZ,
13     O yViscFluxU,
14     I myThid)
15    
16     C !DESCRIPTION:
17     C Calculates the area integrated meridional viscous fluxes of U:
18     C \begin{equation*}
19     C F^y = - \frac{ \Delta y_v \Delta r_f h_z }{\Delta y_u}
20     C ( A_h \delta_j u - A_4 \delta_j \nabla^2 u )
21     C \end{equation*}
22    
23     C !USES: ===============================================================
24     IMPLICIT NONE
25     #include "SIZE.h"
26     #include "EEPARAMS.h"
27     #include "PARAMS.h"
28     #include "GRID.h"
29    
30     C !INPUT PARAMETERS: ===================================================
31     C bi,bj :: tile indices
32     C k :: vertical level
33     C uFld :: zonal flow
34     C del2u :: Laplacian of zonal flow
35     C myThid :: thread number
36     INTEGER bi,bj,k
37     _RL uFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
38     _RL del2u(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
39     _RS hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
40     INTEGER myThid
41    
42     C !OUTPUT PARAMETERS: ==================================================
43     C yViscFluxU :: viscous fluxes
44     _RL yViscFluxU(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
45    
46     C !LOCAL VARIABLES: ====================================================
47     C i,j :: loop indices
48     INTEGER I,J
49     _RS gridScalingA4
50     CEOP
51    
52     C - Laplacian and bi-harmonic terms
53     DO j=1-Oly+1,sNy+Oly
54     DO i=1-Olx,sNx+Olx
55     gridScalingA4 = (_dyU(i,j,bi,bj)**4)/(5000.**4)
56     yViscFluxU(i,j) =
57     & _dxV(i,j,bi,bj)*drF(k)*hFacZ(i,j)
58     & *(
59     & -viscAh*(uFld(i,j)-uFld(i,j-1))
60     #ifdef ISOTROPIC_COS_SCALING
61     & *cosFacV(J,bi,bj)
62     #endif
63     & +viscA4*(del2u(i,j) -del2u(i,j-1) )
64     & *gridScalingA4
65     #ifdef ISOTROPIC_COS_SCALING
66     #ifdef COSINEMETH_III
67     & *sqcosFacV(J,bi,bj)
68     #else
69     & *cosFacV(J,bi,bj)
70     #endif
71     #endif
72     & )*_recip_dyU(i,j,bi,bj)
73     ENDDO
74     ENDDO
75    
76     RETURN
77     END

  ViewVC Help
Powered by ViewVC 1.1.22