model/src/ini_vertical_grid.F

C $Header: /u/gcmpack/models/MITgcmUV/model/src/ini_vertical_grid.F,v 1.11 2001/02/04 14:38:48 cnh Exp $
C $Name:  $

#include "CPP_OPTIONS.h"

CBOP
C     !ROUTINE: INI_VERTICAL_GRID
C     !INTERFACE:
      SUBROUTINE INI_VERTICAL_GRID( myThid )
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE INI_VERTICAL_GRID                              
C     | o Initialise vertical gridding arrays                     
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE
C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine arguments ==
C     myThid -  Number of this instance of INI_DEPTHS
      INTEGER myThid

C     !LOCAL VARIABLES:
C     == Local variables ==
C     K
      INTEGER K
CEOP

C     Calculate depths of centers and interfaces
      rF(1) = Ro_SeaLevel
      DO K=1,Nr
       drF(K)     = delR(K)
       rF(K+1) = rF(K)-rkFac*delR(K)
      ENDDO
      drC(1)      = delR(1) * 0.5 _d 0
      rC(1)       = rf(1)-rkFac*delR(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)
      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	cnh	1.12	C $Header: /u/gcmpack/models/MITgcmUV/model/src/ini_vertical_grid.F,v 1.11 2001/02/04 14:38:48 cnh Exp $
2			C $Name: $
3	adcroft	1.1
4	cnh	1.7	#include "CPP_OPTIONS.h"
5	adcroft	1.1
6	cnh	1.12	CBOP
7			C !ROUTINE: INI_VERTICAL_GRID
8			C !INTERFACE:
9	adcroft	1.1	SUBROUTINE INI_VERTICAL_GRID( myThid )
10	cnh	1.12	C !DESCRIPTION: \bv
11			C ==========================================================
12			C \| SUBROUTINE INI_VERTICAL_GRID
13			C \| o Initialise vertical gridding arrays
14			C ==========================================================
15			C \ev
16
17			C !USES:
18	adcroft	1.8	IMPLICIT NONE
19	adcroft	1.1	C === Global variables ===
20			#include "SIZE.h"
21			#include "EEPARAMS.h"
22			#include "PARAMS.h"
23			#include "GRID.h"
24
25	cnh	1.12	C !INPUT/OUTPUT PARAMETERS:
26	adcroft	1.1	C == Routine arguments ==
27			C myThid - Number of this instance of INI_DEPTHS
28			INTEGER myThid
29
30	cnh	1.12	C !LOCAL VARIABLES:
31	adcroft	1.1	C == Local variables ==
32			C K
33			INTEGER K
34	cnh	1.12	CEOP
35	adcroft	1.1
36			C Calculate depths of centers and interfaces
37	adcroft	1.10	rF(1) = Ro_SeaLevel
38	cnh	1.4	DO K=1,Nr
39	cnh	1.5	drF(K) = delR(K)
40			rF(K+1) = rF(K)-rkFac*delR(K)
41	adcroft	1.1	ENDDO
42	cnh	1.5	drC(1) = delR(1) * 0.5 _d 0
43			rC(1) = rf(1)-rkFacdelR(1) 0.5 _d 0
44	cnh	1.4	DO K=2,Nr
45	cnh	1.5	drC(K) = 0.5 _d 0 *(delR(K-1)+delR(K))
46			rC(K) = rC(K-1) - rkFac*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