/[MITgcm]/MITgcm/model/src/impldiff.F
ViewVC logotype

Diff of /MITgcm/model/src/impldiff.F

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

revision 1.16 by heimbach, Mon May 14 21:46:17 2001 UTC revision 1.17 by cnh, Wed Sep 26 18:09:15 2001 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6  C     |==========================================================|  CBOP
7  C     | S/R IMPLDIFF                                             |  C     !ROUTINE: IMPLDIFF
8  C     | o Solve implicit diffusion equation for vertical         |  C     !INTERFACE:
 C     |   diffusivity.                                           |  
 C     | o Recoded from 2d intermediate fields to 3d to reduce    |  
 C     |   TAMC storage                                           |  
 C     | o Fixed missing masks for fields a(), c()                |  
 C     |==========================================================|  
9        SUBROUTINE IMPLDIFF( bi, bj, iMin, iMax, jMin, jMax,        SUBROUTINE IMPLDIFF( bi, bj, iMin, iMax, jMin, jMax,
10       I                     deltaTX,KappaRX,recip_hFac,       I                     deltaTX,KappaRX,recip_hFac,
11       U                     gXNm1,       U                     gXNm1,
12       I                     myThid )       I                     myThid )
13    C     !DESCRIPTION: \bv
14    C     *==========================================================*
15    C     | S/R IMPLDIFF                                              
16    C     | o Solve implicit diffusion equation for vertical          
17    C     |   diffusivity.                                            
18    C     *==========================================================*
19    C     | o Recoded from 2d intermediate fields to 3d to reduce    
20    C     |   TAMC storage                                            
21    C     | o Fixed missing masks for fields a(), c()                
22    C     *==========================================================*
23    C     \ev
24    
25    C     !USES:
26        IMPLICIT NONE        IMPLICIT NONE
27  C     == Global data ==  C     == Global data ==
28  #include "SIZE.h"  #include "SIZE.h"
# Line 22  C     == Global data == Line 30  C     == Global data ==
30  #include "EEPARAMS.h"  #include "EEPARAMS.h"
31  #include "PARAMS.h"  #include "PARAMS.h"
32  #include "GRID.h"  #include "GRID.h"
   
33  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
34  #include "tamc_keys.h"  #include "tamc_keys.h"
35  #endif  #endif
36    
37    C     !INPUT/OUTPUT PARAMETERS:
38  C     == Routine Arguments ==  C     == Routine Arguments ==
39        INTEGER bi,bj,iMin,iMax,jMin,jMax        INTEGER bi,bj,iMin,iMax,jMin,jMax
40        _RL deltaTX        _RL deltaTX
41        _RL KappaRX(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL KappaRX(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
42        _RS recip_hFac(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr,nSx,nSy)        _RS recip_hFac(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr,nSx,nSy)
43        _RL gXnm1(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr,nSx,nSy)        _RL gXnm1(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr,nSx,nSy)
       _RL gYnm1(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr,nSx,nSy)  
44        INTEGER myThid        INTEGER myThid
45    
46    C     !LOCAL VARIABLES:
47  C     == Local variables ==  C     == Local variables ==
48        INTEGER i,j,k        INTEGER i,j,k
49          _RL gYnm1(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr,nSx,nSy)
50        _RL a(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL a(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
51        _RL b(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL b(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
52        _RL c(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL c(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
53        _RL bet(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL bet(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
54        _RL gam(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL gam(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
55    CEOP
56    
57  C--   Initialise  C--   Initialise
58        DO k=1,Nr        DO k=1,Nr
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22