1 |
C $Header$ |
C $Header$ |
2 |
|
|
3 |
#include "CPP_EEOPTIONS.h" |
#include "CPP_OPTIONS.h" |
4 |
|
|
5 |
CStartOfInterface |
CStartOfInterface |
6 |
SUBROUTINE INI_CARTESIAN_GRID( myThid ) |
SUBROUTINE INI_CARTESIAN_GRID( myThid ) |
33 |
C | and Y are in metres. Disktance in Z are in m or Pa | |
C | and Y are in metres. Disktance in Z are in m or Pa | |
34 |
C | depending on the vertical gridding mode. | |
C | depending on the vertical gridding mode. | |
35 |
C \==========================================================/ |
C \==========================================================/ |
36 |
|
IMPLICIT NONE |
37 |
|
|
38 |
C === Global variables === |
C === Global variables === |
39 |
#include "SIZE.h" |
#include "SIZE.h" |
48 |
|
|
49 |
C == Local variables == |
C == Local variables == |
50 |
C xG, yG - Global coordinate location. |
C xG, yG - Global coordinate location. |
|
C zG |
|
51 |
C xBase - South-west corner location for process. |
C xBase - South-west corner location for process. |
52 |
C yBase |
C yBase |
|
C zUpper - Work arrays for upper and lower |
|
|
C zLower cell-face heights. |
|
|
C phi - Temporary scalar |
|
53 |
C xBase - Temporaries for lower corner coordinate |
C xBase - Temporaries for lower corner coordinate |
54 |
C yBase |
C yBase |
55 |
C iG, jG - Global coordinate index. Usually used to hold |
C iG, jG - Global coordinate index. Usually used to hold |
58 |
C zUpper - Temporary arrays holding z coordinates of |
C zUpper - Temporary arrays holding z coordinates of |
59 |
C zLower upper and lower faces. |
C zLower upper and lower faces. |
60 |
C I,J,K |
C I,J,K |
61 |
_RL xG, yG, zG |
_RL xG, yG |
|
_RL phi |
|
|
_RL zUpper(Nr), zLower(Nr) |
|
62 |
_RL xBase, yBase |
_RL xBase, yBase |
63 |
INTEGER iG, jG |
INTEGER iG, jG |
64 |
INTEGER bi, bj |
INTEGER bi, bj |
65 |
INTEGER I, J, K |
INTEGER I, J |
66 |
|
|
67 |
C-- Simple example of inialisation on cartesian grid |
C-- Simple example of inialisation on cartesian grid |
68 |
C-- First set coordinates of cell centers |
C-- First set coordinates of cell centers |
151 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
152 |
DO J=1,sNy |
DO J=1,sNy |
153 |
DO I=1,sNx |
DO I=1,sNx |
154 |
rA(I,J,bi,bj) = dxF(I,J,bi,bj)*dyF(I,J,bi,bj) |
rA (I,J,bi,bj) = dxF(I,J,bi,bj)*dyF(I,J,bi,bj) |
155 |
|
rAw(I,J,bi,bj) = dxC(I,J,bi,bj)*dyG(I,J,bi,bj) |
156 |
|
rAs(I,J,bi,bj) = dxG(I,J,bi,bj)*dyC(I,J,bi,bj) |
157 |
tanPhiAtU(I,J,bi,bj) = 0. _d 0 |
tanPhiAtU(I,J,bi,bj) = 0. _d 0 |
158 |
tanPhiAtV(I,J,bi,bj) = 0. _d 0 |
tanPhiAtV(I,J,bi,bj) = 0. _d 0 |
159 |
ENDDO |
ENDDO |
161 |
ENDDO |
ENDDO |
162 |
ENDDO |
ENDDO |
163 |
_EXCH_XY_R4 (rA , myThid ) |
_EXCH_XY_R4 (rA , myThid ) |
164 |
|
_EXCH_XY_R4 (rAw , myThid ) |
165 |
|
_EXCH_XY_R4 (rAs , myThid ) |
166 |
_EXCH_XY_R4 (tanPhiAtU , myThid ) |
_EXCH_XY_R4 (tanPhiAtU , myThid ) |
167 |
_EXCH_XY_R4 (tanPhiAtV , myThid ) |
_EXCH_XY_R4 (tanPhiAtV , myThid ) |
168 |
|
|