| 15 |
C by each thread and initializes only the region of the domain it is |
C by each thread and initializes only the region of the domain it is |
| 16 |
C "responsible" for. |
C "responsible" for. |
| 17 |
|
|
| 18 |
|
C !CALLING SEQUENCE: |
| 19 |
|
C INI_GRID |
| 20 |
|
C | -- LOAD_GRID_SPACING |
| 21 |
|
C | -- INI_VERTICAL_GRID |
| 22 |
|
C | / INI_CARTESIAN_GRID |
| 23 |
|
C | / INI_SPHERICAL_POLAR_GRID |
| 24 |
|
C | \ INI_CURVILINEAR_GRID |
| 25 |
|
C | \ INI_CYLINDER_GRID |
| 26 |
|
|
| 27 |
C !USES: |
C !USES: |
| 28 |
IMPLICIT NONE |
IMPLICIT NONE |
| 29 |
#include "SIZE.h" |
#include "SIZE.h" |
| 56 |
C load grid spacing (vector) from files |
C load grid spacing (vector) from files |
| 57 |
CALL LOAD_GRID_SPACING( myThid ) |
CALL LOAD_GRID_SPACING( myThid ) |
| 58 |
|
|
|
C Set up reference vertical profile (vector) for state variables |
|
|
C jmc: this call will not stay here but will move to S/R INITIALISE_FIXED |
|
|
C (a better place) once anelatic initialisation is cleaned-up. |
|
|
CALL INI_REFERENCE_STATE( myThid ) |
|
|
|
|
| 59 |
C Set up vertical grid and coordinate system |
C Set up vertical grid and coordinate system |
| 60 |
CALL INI_VERTICAL_GRID( myThid ) |
CALL INI_VERTICAL_GRID( myThid ) |
| 61 |
|
|
| 93 |
_END_MASTER(myThid) |
_END_MASTER(myThid) |
| 94 |
ENDIF |
ENDIF |
| 95 |
|
|
| 96 |
|
#ifdef ALLOW_OBCS |
| 97 |
|
IF ( useOBCS ) THEN |
| 98 |
|
C Fix overlap regions to avoid discontinuities in dxc, dyc, etc. |
| 99 |
|
CALL OBCS_CHECK_GRID( myThid ) |
| 100 |
|
ENDIF |
| 101 |
|
#endif /* ALLOW_OBCS */ |
| 102 |
|
|
| 103 |
#ifdef ALLOW_MONITOR |
#ifdef ALLOW_MONITOR |
| 104 |
IF ( MASTER_CPU_IO(myThid) ) THEN |
IF ( MASTER_CPU_IO(myThid) ) THEN |
| 105 |
C-- only the master thread is allowed to switch On/Off mon_write_stdout |
C-- only the master thread is allowed to switch On/Off mon_write_stdout |