/[MITgcm]/MITgcm/pkg/layers/layers_init_varia.F
ViewVC logotype

Diff of /MITgcm/pkg/layers/layers_init_varia.F

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

revision 1.6 by jmc, Thu Oct 18 19:51:14 2012 UTC revision 1.7 by rpa, Wed Jun 4 14:48:32 2014 UTC
# Line 25  C     myThid ::  my Thread Id number Line 25  C     myThid ::  my Thread Id number
25  #ifdef ALLOW_LAYERS  #ifdef ALLOW_LAYERS
26    
27  C     === Local variables ===  C     === Local variables ===
28        INTEGER i,j,kg,bi,bj        INTEGER i,j,kg,bi,bj,iTracer
29  c     INTEGER iLa  c     INTEGER iLa
30    
31        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
32         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
33    
34  C     Loop through 3D variables:  C     Loop through 3D layers variables:
35          DO kg=1,Nlayers          DO kg=1,Nlayers
36           DO j=1-OLy,sNy+OLy           DO j=1-OLy,sNy+OLy
37            DO i=1-OLx,sNx+OLx            DO i=1-OLx,sNx+OLx
# Line 54  C     Loop through 3D variables: Line 54  C     Loop through 3D variables:
54  #endif /* LAYERS_THICKNESS */  #endif /* LAYERS_THICKNESS */
55  #endif /* LAYERS_VFLUX */  #endif /* LAYERS_VFLUX */
56    
57    #ifdef LAYERS_THERMODYNAMICS
58               layers_TtendSurf (i,j,kg,bi,bj) = 0. _d 0
59               layers_TtendDiffh(i,j,kg,bi,bj) = 0. _d 0
60               layers_TtendDiffr(i,j,kg,bi,bj) = 0. _d 0
61               layers_StendSurf (i,j,kg,bi,bj) = 0. _d 0
62               layers_StendDiffh(i,j,kg,bi,bj) = 0. _d 0
63               layers_StendDiffr(i,j,kg,bi,bj) = 0. _d 0
64    #endif /* LAYERS_THERMODYNAMICS */
65    
66            ENDDO            ENDDO
67           ENDDO           ENDDO
68          ENDDO          ENDDO
# Line 67  C     Loop through 3D variables: Line 76  C     Loop through 3D variables:
76          ENDDO          ENDDO
77  #endif /* LAYERS_PRHO_REF */  #endif /* LAYERS_PRHO_REF */
78    
79    
80    #ifdef LAYERS_THERMODYNAMICS  
81            DO iTracer = 1,2
82             DO kg=1,Nr
83              DO j=1-OLy,sNy+OLy
84               DO i=1-OLx,sNx+OLx
85                layers_surfflux(i,j,kg,iTracer,bi,bj) = 0. _d 0
86                layers_dfx     (i,j,kg,iTracer,bi,bj) = 0. _d 0
87                layers_dfy     (i,j,kg,iTracer,bi,bj) = 0. _d 0
88                layers_dfr     (i,j,kg,iTracer,bi,bj) = 0. _d 0
89               ENDDO
90              ENDDO
91             ENDDO
92            ENDDO  
93    #endif /* LAYERS_THERMODYNAMICS */
94    
95  C- end bi,bj loops  C- end bi,bj loops
96         ENDDO         ENDDO
97        ENDDO        ENDDO
98                                      
99    
100  #endif /* ALLOW_LAYERS */  #endif /* ALLOW_LAYERS */
101    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22