37 |
|
|
38 |
_BEGIN_MASTER(myThid) |
_BEGIN_MASTER(myThid) |
39 |
|
|
40 |
|
C-- Set factors required for mixing pressure and meters as vertical coordinate. |
41 |
|
C rkSign is a "sign" parameter which is used where the orientation of the vertical |
42 |
|
C coordinate (pressure or meters) relative to the vertical index (K) is important. |
43 |
|
C rkSign = -1 applies when K and the coordinate are in the opposite sense. |
44 |
|
C rkSign = 1 applies when K and the coordinate are in the same sense. |
45 |
|
rkSign = -1. _d 0 |
46 |
|
gravitySign = -1. _d 0 |
47 |
|
IF ( usingPCoords ) THEN |
48 |
|
gravitySign = 1. _d 0 |
49 |
|
ENDIF |
50 |
|
|
51 |
IF (setCenterDr) THEN |
IF (setCenterDr) THEN |
52 |
C-- Interface at middle between 2 centers : |
C-- Interface at middle between 2 centers : |
53 |
|
|
113 |
|
|
114 |
C- Calculate reciprol vertical grid spacing : |
C- Calculate reciprol vertical grid spacing : |
115 |
DO K=1,Nr |
DO K=1,Nr |
|
saFac(K) = 1. _d 0 |
|
116 |
recip_drC(K) = 1. _d 0/drC(K) |
recip_drC(K) = 1. _d 0/drC(K) |
117 |
recip_drF(K) = 1. _d 0/drF(K) |
recip_drF(K) = 1. _d 0/drF(K) |
118 |
ENDDO |
ENDDO |