model/src/ini_vertical_grid.F

C $Header: /u/gcmpack/models/MITgcmUV/model/src/ini_vertical_grid.F,v 1.7 1998/11/06 22:44:48 cnh Exp $

#include "CPP_OPTIONS.h"

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

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
      rF(1) = 0. _d 0
      DO K=1,Nr
       drF(K)     = delR(K)
       rF(K+1) = rF(K)-rkFac*delR(K)
C      drF(K)     = delZ(K)
C      rF(K+1)    = rF(K)-delZ(K)
      ENDDO
      drC(1)      = delR(1) * 0.5 _d 0
      rC(1)       = rf(1)-rkFac*delR(1) * 0.5 _d 0
C     drC(1)      = delZ(1) * 0.5 _d 0
C     rC(1)       = -delZ(1) * 0.5 _d 0
      DO K=2,Nr
       drC(K)     = 0.5 _d 0 *(delR(K-1)+delR(K))
       rC(K)      = rC(K-1) - rkFac*drC(K)
C      drC(K)     = 0.5 _d 0 *(delZ(K-1)+delZ(K))
C      rC(K)      = rC(K-1) - drC(K)
      ENDDO
      DO K=1,Nr
       saFac(K)  = 1. _d 0
       recip_drC(K)   = 1. _d 0/drC(K)
       recip_drF(K)   = 1. _d 0/drF(K)
      ENDDO
C
      RETURN
      END
1	adcroft	1.8	C $Header: /u/gcmpack/models/MITgcmUV/model/src/ini_vertical_grid.F,v 1.7 1998/11/06 22:44:48 cnh Exp $
2	adcroft	1.1
3	cnh	1.7	#include "CPP_OPTIONS.h"
4	adcroft	1.1
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	adcroft	1.8	IMPLICIT NONE
14	adcroft	1.1
15			C === Global variables ===
16			#include "SIZE.h"
17			#include "EEPARAMS.h"
18			#include "PARAMS.h"
19			#include "GRID.h"
20
21			C == Routine arguments ==
22			C myThid - Number of this instance of INI_DEPTHS
23			INTEGER myThid
24			CEndOfInterface
25
26			C == Local variables ==
27			C K
28			INTEGER K
29
30			C Calculate depths of centers and interfaces
31	cnh	1.4	rF(1) = 0. _d 0
32			DO K=1,Nr
33	cnh	1.5	drF(K) = delR(K)
34			rF(K+1) = rF(K)-rkFac*delR(K)
35			C drF(K) = delZ(K)
36			C rF(K+1) = rF(K)-delZ(K)
37	adcroft	1.1	ENDDO
38	cnh	1.5	drC(1) = delR(1) * 0.5 _d 0
39			rC(1) = rf(1)-rkFacdelR(1) 0.5 _d 0
40			C drC(1) = delZ(1) * 0.5 _d 0
41			C rC(1) = -delZ(1) * 0.5 _d 0
42	cnh	1.4	DO K=2,Nr
43	cnh	1.5	drC(K) = 0.5 _d 0 *(delR(K-1)+delR(K))
44			rC(K) = rC(K-1) - rkFac*drC(K)
45			C drC(K) = 0.5 _d 0 *(delZ(K-1)+delZ(K))
46			C rC(K) = rC(K-1) - drC(K)
47	adcroft	1.1	ENDDO
48	cnh	1.4	DO K=1,Nr
49	cnh	1.6	saFac(K) = 1. _d 0
50	cnh	1.4	recip_drC(K) = 1. _d 0/drC(K)
51			recip_drF(K) = 1. _d 0/drF(K)
52	adcroft	1.1	ENDDO
53			C
54			RETURN
55			END