/[MITgcm]/MITgcm/model/src/initialise_fixed.F
ViewVC logotype

Contents of /MITgcm/model/src/initialise_fixed.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.24 - (show annotations) (download)
Thu Feb 26 22:31:58 2004 UTC (20 years, 3 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint52l_pre, hrcube5, checkpoint52l_post
Changes since 1.23: +1 -46 lines
 o fix the calling sequence for MNC
 o add coordinates to the variables within the "state" file
 o reads still unfinished

1 C $Header: /u/u3/gcmpack/MITgcm/model/src/initialise_fixed.F,v 1.23 2004/02/05 19:38:55 edhill Exp $
2 C $Name: $
3
4 #include "PACKAGES_CONFIG.h"
5 #include "CPP_OPTIONS.h"
6
7 CBOP
8 C !ROUTINE: INITIALISE_FIXED
9 C !INTERFACE:
10 SUBROUTINE INITIALISE_FIXED(myThid)
11
12 C !DESCRIPTION: \bv
13 C *==========================================================*
14 C | SUBROUTINE INITIALISE_FIXED
15 C | o Routine for setting fixed model arrays such as
16 C | topography, grid, solver matrices, etc.
17 C *==========================================================*
18 C | INITIALISE_FIXED is invoked at the start of the model to
19 C | set fixed model arrays. It reads data from an input file
20 C | and from various binary files.
21 C | Each thread invokes an instance of this routine as does
22 C | each process in a multi-process parallel environment like
23 C | MPI.
24 C *==========================================================*
25 C \ev
26
27 C !CALLING SEQUENCE:
28 C INITIALISE_FIXED
29 C |
30 C |-- INI_PARMS
31 C |
32 C |-- MON_INIT
33 C |
34 C |-- INI_GRID
35 C |
36 C |-- INI_EOS
37 C |
38 C |-- INI_DEPTHS
39 C |
40 C |-- INI_MASKS_ETC
41 C |
42 C |-- PACKAGES_BOOT
43 C |
44 C |-- PACKAGES_READPARMS
45 C |
46 C |-- PACKAGES_INIT_FIXED
47 C |
48 C |-- PACKAGES_CHECK
49 C |
50 C |-- INI_LINEAR_PHSURF
51 C |
52 C |-- INI_CORI
53 C |
54 C |-- INI_CG2D
55 C |
56 C |-- INI_CG3D
57 C |
58 C |-- CONFIG_CHECK
59 C |
60 C |-- CONFIG_SUMMARY
61 C |
62 C |-- MNC_INIT
63
64 C !USES:
65 IMPLICIT NONE
66 C == Global variables ==
67 #include "SIZE.h"
68 #include "EEPARAMS.h"
69 #include "PARAMS.h"
70
71 C !INPUT/OUTPUT PARAMETERS:
72 C == Routine arguments ==
73 INTEGER myThid
74 CEOP
75
76 C-- Set model parameters.
77 C Parameters are set to defaults and then updates are read from
78 C an input file called data.
79 CALL INI_PARMS( myThid )
80 _BARRIER
81
82 #ifdef ALLOW_MONITOR
83 C-- Initialise MONITOR I/O streams so we can report config. info
84 CALL MON_INIT( myThid )
85 _BARRIER
86 #endif
87
88
89 C-- Set model grid.
90 C Variables defining model grid spacing are defined.
91 CALL INI_GRID( myThid )
92 _BARRIER
93
94 C-- Set equation of state parameters.
95 CALL INI_EOS( myThid )
96 _BARRIER
97
98 C-- Initialise map of depths
99 CALL INI_DEPTHS( myThid )
100 _BARRIER
101
102 C-- Derive masks, lopping factors and recipricols of quantities.
103 C Volumes and areas are set according to grid and depth map.
104 CALL INI_MASKS_ETC( myThid )
105 _BARRIER
106
107 C-- Configure packages
108 CALL PACKAGES_BOOT( myThid )
109
110 C-- Read configuration parameters for packages
111 CALL PACKAGES_READPARMS( myThid )
112
113 C-- Call fixed data initialization phase of packages
114 CALL PACKAGES_INIT_FIXED( myThid )
115
116 C-- Check dependances between packages
117 CALL PACKAGES_CHECK( myThid )
118
119 C-- Set Bo_surf => define the Linear Relation: Phi_surf(eta)
120 CALL INI_LINEAR_PHISURF( myThid )
121
122 C-- Set coriolis operators
123 CALL INI_CORI( myThid )
124
125 C-- Set laplace operators for use in 2D conjugate gradient solver.
126 CALL INI_CG2D( myThid )
127
128 #ifdef ALLOW_NONHYDROSTATIC
129 C-- Set laplace operators for use in 3D conjugate gradient solver.
130 ceh3 needs an IF ( useNONHYDROSTATIC ) THEN
131 CALL INI_CG3D( myThid )
132 #endif
133
134 C-- Check parameters and model cofiguration
135 CALL CONFIG_CHECK( myThid )
136
137 C-- Finally summarise the model cofiguration
138 CALL CONFIG_SUMMARY( myThid )
139
140 #ifdef COMPONENT_MODULE
141 C-- Post component-model configuration information to coupler
142 C and get config. info for other component(s).
143 IF ( useCoupler ) CALL CPL_EXCH_CONFIGS( myThid )
144 #endif
145
146 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
147
148 RETURN
149 END

  ViewVC Help
Powered by ViewVC 1.1.22