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

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

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

revision 1.2 by adcroft, Wed Jul 29 18:33:48 1998 UTC revision 1.12 by cnh, Wed Sep 26 18:09:16 2001 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "CPP_EEOPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6  CStartOfInterface  CBOP
7    C     !ROUTINE: INI_VERTICAL_GRID
8    C     !INTERFACE:
9        SUBROUTINE INI_VERTICAL_GRID( myThid )        SUBROUTINE INI_VERTICAL_GRID( myThid )
10  C     /==========================================================\  C     !DESCRIPTION: \bv
11  C     | SUBROUTINE INI_VERTICAL_GRID                             |  C     *==========================================================*
12  C     | o Initialise vertical gridding arrays                    |  C     | SUBROUTINE INI_VERTICAL_GRID                              
13  C     |==========================================================|  C     | o Initialise vertical gridding arrays                    
14  C     |                                                          |  C     *==========================================================*
15  C     \==========================================================/  C     \ev
16    
17    C     !USES:
18          IMPLICIT NONE
19  C     === Global variables ===  C     === Global variables ===
20  #include "SIZE.h"  #include "SIZE.h"
21  #include "EEPARAMS.h"  #include "EEPARAMS.h"
22  #include "PARAMS.h"  #include "PARAMS.h"
23  #include "GRID.h"  #include "GRID.h"
24    
25    C     !INPUT/OUTPUT PARAMETERS:
26  C     == Routine arguments ==  C     == Routine arguments ==
27  C     myThid -  Number of this instance of INI_DEPTHS  C     myThid -  Number of this instance of INI_DEPTHS
28        INTEGER myThid        INTEGER myThid
 CEndOfInterface  
29    
30    C     !LOCAL VARIABLES:
31  C     == Local variables ==  C     == Local variables ==
32  C     K  C     K
33        INTEGER K        INTEGER K
34    CEOP
35    
36  C     Calculate depths of centers and interfaces  C     Calculate depths of centers and interfaces
37        zFace(1) = 0. _d 0        rF(1) = Ro_SeaLevel
38        DO K=1,Nz        DO K=1,Nr
39         dzF(K)     = delZ(K)         drF(K)     = delR(K)
40         zFace(K+1) = zFace(K)-delZ(K)         rF(K+1) = rF(K)-rkFac*delR(K)
41        ENDDO        ENDDO
42        dzC(1)      = delZ(1) * 0.5 _d 0        drC(1)      = delR(1) * 0.5 _d 0
43        zC(1)       = delZ(1) * 0.5 _d 0        rC(1)       = rf(1)-rkFac*delR(1) * 0.5 _d 0
44        DO K=2,Nz        DO K=2,Nr
45         dzC(K)     = 0.5 _d 0 *(delZ(K-1)+delZ(K))         drC(K)     = 0.5 _d 0 *(delR(K-1)+delR(K))
46         zC(K)      = zC(K) - dzC(K)         rC(K)      = rC(K-1) - rkFac*drC(K)
47        ENDDO        ENDDO
48        DO K=1,Nz        DO K=1,Nr
49         saFac(K)  = 1. D0         saFac(K)  = 1. _d 0
50         rDzC(K)   = 1. _d 0/dzC(K)         recip_drC(K)   = 1. _d 0/drC(K)
51         rDzF(K)   = 1. _d 0/dzF(K)         recip_drF(K)   = 1. _d 0/drF(K)
52        ENDDO        ENDDO
53  C  C
54        RETURN        RETURN
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22