pkg/layers/layers_init_varia.F

108a00eab9 Ryan*0001 #include "LAYERS_OPTIONS.h"
                0002
                0003 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0004
                0005       SUBROUTINE LAYERS_INIT_VARIA( myThid )
                0006
                0007 C ===================================================================
                0008 C     Initialize LAYERS variables.
                0009 C ===================================================================
                0010
                0011       IMPLICIT NONE
                0012 #include "EEPARAMS.h"
                0013 #include "SIZE.h"
                0014 #include "PARAMS.h"
                0015 #include "LAYERS_SIZE.h"
                0016 #include "LAYERS.h"
                0017
                0018 C  INPUT/OUTPUT PARAMETERS:
                0019 C     myThid ::  my Thread Id number
                0020       INTEGER myThid
                0021
                0022 #ifdef ALLOW_LAYERS
                0023
                0024 C     === Local variables ===
62b1f3e57f Jean*0025       INTEGER bi, bj
d21fbd4209 Jean*0026 #ifdef LAYERS_THERMODYNAMICS
62b1f3e57f Jean*0027       INTEGER i,j,kg
d21fbd4209 Jean*0028       INTEGER iTracer
                0029 #endif
9e68be709e Jean*0030
108a00eab9 Ryan*0031       DO bj = myByLo(myThid), myByHi(myThid)
                0032        DO bi = myBxLo(myThid), myBxHi(myThid)
                0033
cf336ab6c5 Ryan*0034 C     Loop through 3D layers variables:
62b1f3e57f Jean*0035 #ifdef LAYERS_THERMODYNAMICS
108a00eab9 Ryan*0036         DO kg=1,Nlayers
9e68be709e Jean*0037          DO j=1-OLy,sNy+OLy
                0038           DO i=1-OLx,sNx+OLx
50d8304171 Ryan*0039            layers_Hc (i,j,kg,bi,bj) = 0. _d 0
                0040            layers_PIc(i,j,kg,bi,bj) = 0. _d 0
                0041           ENDDO
                0042          ENDDO
                0043         ENDDO
                0044         DO kg=1,Nlayers-1
                0045          DO j=1-OLy,sNy+OLy
                0046           DO i=1-OLx,sNx+OLx
                0047            layers_Hc        (i,j,kg,bi,bj) = 0. _d 0
d21fbd4209 Jean*0048            layers_TtendSurf (i,j,kg,bi,bj) = 0. _d 0
                0049            layers_TtendDiffh(i,j,kg,bi,bj) = 0. _d 0
                0050            layers_TtendDiffr(i,j,kg,bi,bj) = 0. _d 0
50d8304171 Ryan*0051            layers_TtendAdvh (i,j,kg,bi,bj) = 0. _d 0
                0052            layers_TtendAdvr (i,j,kg,bi,bj) = 0. _d 0
d21fbd4209 Jean*0053            layers_StendSurf (i,j,kg,bi,bj) = 0. _d 0
                0054            layers_StendDiffh(i,j,kg,bi,bj) = 0. _d 0
                0055            layers_StendDiffr(i,j,kg,bi,bj) = 0. _d 0
50d8304171 Ryan*0056            layers_StendAdvh(i,j,kg,bi,bj) = 0. _d 0
                0057            layers_StendAdvr(i,j,kg,bi,bj) = 0. _d 0
108a00eab9 Ryan*0058           ENDDO
                0059          ENDDO
                0060         ENDDO
1f346f78e1 Jean*0061 #ifdef LAYERS_PRHO_REF
                0062         DO kg=1,Nr
                0063          DO j=1-OLy,sNy+OLy
                0064           DO i=1-OLx,sNx+OLx
50d8304171 Ryan*0065            layers_alpha(i,j,kg,bi,bj) = 0. _d 0
                0066            layers_beta(i,j,kg,bi,bj) = 0. _d 0
1f346f78e1 Jean*0067           ENDDO
                0068          ENDDO
                0069         ENDDO
                0070 #endif /* LAYERS_PRHO_REF */
108a00eab9 Ryan*0071
cf336ab6c5 Ryan*0072         DO iTracer = 1,2
62b1f3e57f Jean*0073          DO j=1-OLy,sNy+OLy
50d8304171 Ryan*0074           DO i=1-OLx,sNx+OLx
                0075            layers_surfflux(i,j,1,iTracer,bi,bj) = 0. _d 0
                0076           ENDDO
                0077          ENDDO
cf336ab6c5 Ryan*0078          DO kg=1,Nr
                0079           DO j=1-OLy,sNy+OLy
                0080            DO i=1-OLx,sNx+OLx
50d8304171 Ryan*0081             layers_tottend (i,j,kg,iTracer,bi,bj) = 0. _d 0
d21fbd4209 Jean*0082             layers_dfx     (i,j,kg,iTracer,bi,bj) = 0. _d 0
cf336ab6c5 Ryan*0083             layers_dfy     (i,j,kg,iTracer,bi,bj) = 0. _d 0
                0084             layers_dfr     (i,j,kg,iTracer,bi,bj) = 0. _d 0
50d8304171 Ryan*0085             layers_afx     (i,j,kg,iTracer,bi,bj) = 0. _d 0
                0086             layers_afy     (i,j,kg,iTracer,bi,bj) = 0. _d 0
                0087             layers_afr     (i,j,kg,iTracer,bi,bj) = 0. _d 0
                0088 #ifdef SHORTWAVE_HEATING
                0089             layers_sw      (i,j,kg,1      ,bi,bj) = 0. _d 0
                0090 #endif /* SHORTWAVE_HEATING */
cf336ab6c5 Ryan*0091            ENDDO
                0092           ENDDO
                0093          ENDDO
d21fbd4209 Jean*0094         ENDDO
cf336ab6c5 Ryan*0095 #endif /* LAYERS_THERMODYNAMICS */
                0096
108a00eab9 Ryan*0097 C- end bi,bj loops
                0098        ENDDO
                0099       ENDDO
                0100
                0101 #endif /* ALLOW_LAYERS */
                0102
                0103       RETURN
                0104       END