/[MITgcm]/MITgcm/pkg/my82/my82_calc.F
ViewVC logotype

Diff of /MITgcm/pkg/my82/my82_calc.F

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

revision 1.6 by jmc, Thu Oct 8 20:07:53 2009 UTC revision 1.7 by jmc, Wed Jun 27 22:39:09 2012 UTC
# Line 8  C !ROUTINE: MY82_CALC Line 8  C !ROUTINE: MY82_CALC
8    
9  C !INTERFACE: ======================================================  C !INTERFACE: ======================================================
10        subroutine MY82_CALC(        subroutine MY82_CALC(
11       I     bi, bj, myTime, myThid )       I                bi, bj, sigmaR, myTime, myIter, myThid )
12    
13  C !DESCRIPTION: \bv  C !DESCRIPTION: \bv
14  C     *==========================================================*  C     *==========================================================*
# Line 31  C !USES: =============================== Line 31  C !USES: ===============================
31  #include "SIZE.h"  #include "SIZE.h"
32  #include "EEPARAMS.h"  #include "EEPARAMS.h"
33  #include "PARAMS.h"  #include "PARAMS.h"
 c#include "DYNVARS.h"  
34  #include "MY82.h"  #include "MY82.h"
 c#include "FFIELDS.h"  
35  #include "GRID.h"  #include "GRID.h"
36    
37  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
38  c Routine arguments  C Routine arguments
39  c     bi, bj - array indices on which to apply calculations  C     bi, bj :: Current tile indices
40  c     myTime - Current time in simulation  C     sigmaR :: Vertical gradient of iso-neutral density
41    C     myTime :: Current time in simulation
42    C     myIter :: Current time-step number
43    C     myThid :: My Thread Id number
44        INTEGER bi, bj        INTEGER bi, bj
45          _RL     sigmaR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
46        _RL     myTime        _RL     myTime
47          INTEGER myIter
48        INTEGER myThid        INTEGER myThid
49    
50  #ifdef ALLOW_MY82  #ifdef ALLOW_MY82
   
51  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
52  c Local constants  c Local constants
53  C     imin, imax, jmin, jmax  - array computation indices  C     imin, imax, jmin, jmax  - array computation indices
# Line 72  CEOP Line 73  CEOP
73        jMax = sNy+OLy-1        jMax = sNy+OLy-1
74    
75  C     Initialize local fields  C     Initialize local fields
76        DO J=1-Oly,sNy+Oly        DO J=1-OLy,sNy+OLy
77         DO I=1-Olx,sNx+Olx         DO I=1-OLx,sNx+OLx
78          GH(I,J) = 0. _d 0          GH(I,J) = 0. _d 0
79          GM(I,J) = 0. _d 0          GM(I,J) = 0. _d 0
80         ENDDO         ENDDO
81        ENDDO        ENDDO
82        DO K = 1, Nr        DO K = 1, Nr
83         DO J=1-Oly,sNy+Oly         DO J=1-OLy,sNy+OLy
84          DO I=1-Olx,sNx+Olx          DO I=1-OLx,sNx+OLx
85           SH(I,J,K)  = 0. _d 0           SH(I,J,K)  = 0. _d 0
86           SM(I,J,K)  = 0. _d 0           SM(I,J,K)  = 0. _d 0
87           tke(I,J,K) = 0. _d 0           tke(I,J,K) = 0. _d 0
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22