model/src/ini_vertical_grid.F

C $Header: /u/gcmpack/models/MITgcmUV/model/src/ini_vertical_grid.F,v 1.2 1998/07/29 18:33:48 adcroft Exp $

#include "CPP_EEOPTIONS.h"

CStartOfInterface
      SUBROUTINE INI_VERTICAL_GRID( myThid )
C     /==========================================================\
C     | SUBROUTINE INI_VERTICAL_GRID                             |
C     | o Initialise vertical gridding arrays                    |
C     |==========================================================|
C     |                                                          |
C     \==========================================================/

C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"

C     == Routine arguments ==
C     myThid -  Number of this instance of INI_DEPTHS
      INTEGER myThid
CEndOfInterface

C     == Local variables ==
C     K
      INTEGER K

C     Calculate depths of centers and interfaces
      zFace(1) = 0. _d 0
      DO K=1,Nz
       dzF(K)     = delZ(K)
       zFace(K+1) = zFace(K)-delZ(K)
      ENDDO
      dzC(1)      = delZ(1) * 0.5 _d 0
      zC(1)       = -delZ(1) * 0.5 _d 0
      DO K=2,Nz
       dzC(K)     = 0.5 _d 0 *(delZ(K-1)+delZ(K))
       zC(K)      = zC(K-1) - dzC(K)
      ENDDO
      DO K=1,Nz
       saFac(K)  = 1. D0
       rDzC(K)   = 1. _d 0/dzC(K)
       rDzF(K)   = 1. _d 0/dzF(K)
      ENDDO
C
      RETURN
      END
1	adcroft	1.3	C $Header: /u/gcmpack/models/MITgcmUV/model/src/ini_vertical_grid.F,v 1.2 1998/07/29 18:33:48 adcroft Exp $
2	adcroft	1.1
3			#include "CPP_EEOPTIONS.h"
4
5			CStartOfInterface
6			SUBROUTINE INI_VERTICAL_GRID( myThid )
7			C /==========================================================\
8			C \| SUBROUTINE INI_VERTICAL_GRID \|
9			C \| o Initialise vertical gridding arrays \|
10			C \|==========================================================\|
11			C \| \|
12			C \==========================================================/
13
14			C === Global variables ===
15			#include "SIZE.h"
16			#include "EEPARAMS.h"
17			#include "PARAMS.h"
18			#include "GRID.h"
19
20			C == Routine arguments ==
21			C myThid - Number of this instance of INI_DEPTHS
22			INTEGER myThid
23			CEndOfInterface
24
25			C == Local variables ==
26			C K
27			INTEGER K
28
29			C Calculate depths of centers and interfaces
30			zFace(1) = 0. _d 0
31			DO K=1,Nz
32			dzF(K) = delZ(K)
33	adcroft	1.2	zFace(K+1) = zFace(K)-delZ(K)
34	adcroft	1.1	ENDDO
35			dzC(1) = delZ(1) * 0.5 _d 0
36	adcroft	1.3	zC(1) = -delZ(1) * 0.5 _d 0
37	adcroft	1.1	DO K=2,Nz
38			dzC(K) = 0.5 _d 0 *(delZ(K-1)+delZ(K))
39	adcroft	1.3	zC(K) = zC(K-1) - dzC(K)
40	adcroft	1.1	ENDDO
41			DO K=1,Nz
42			saFac(K) = 1. D0
43			rDzC(K) = 1. _d 0/dzC(K)
44			rDzF(K) = 1. _d 0/dzF(K)
45			ENDDO
46			C
47			RETURN
48			END