1 |
C |
C |
2 |
C Invocation from WRAPPER level... |
C Invocation from WRAPPER level... |
3 |
C : |
C |
|
C : |
|
4 |
C | |
C | |
5 |
C |-THE_MODEL_MAIN :: Primary driver for the MITgcm algorithm |
C |-THE_MODEL_MAIN :: Primary driver for the MITgcm algorithm |
6 |
C | :: Called from WRAPPER level numerical |
C | :: Called from WRAPPER level numerical |
7 |
C | :: code innvocation routine. On entry |
C | :: code invocation routine. On entry |
8 |
C | :: to THE_MODEL_MAIN separate thread and |
C | :: to THE_MODEL_MAIN separate thread and |
9 |
C | :: separate processes will have been established. |
C | :: separate processes will have been established. |
10 |
C | :: Each thread and process will have a unique ID |
C | :: Each thread and process will have a unique ID |
11 |
C | :: but as yet it will not be associated with a |
C | :: but as yet it will not be associated with a |
12 |
C | :: specific region in decomposed discrete space. |
C | :: specific region in decomposed discrete space. |
13 |
C | |
C | |
14 |
C |-INITIALISE_FIXED :: Set fixed model arrays such as topography, |
C |-INITIALISE_FIXED :: Set fixed model arrays such as topography, |
15 |
C | | :: grid, solver matrices etc.. |
C | | :: grid, solver matrices etc.. |
16 |
C | | |
C | | |
17 |
C | |-INI_PARMS :: Routine to set kernel model parameters. |
C | |-INI_PARMS :: Routine to set kernel model parameters. |
18 |
C | | :: By default kernel parameters are read from file |
C | | :: Kernel parameters are read from file "data" |
19 |
C | | :: "data" in directory in which code executes. |
C | | :: in directory in which code executes. |
20 |
C | | |
C | | |
21 |
C | |-MON_INIT :: Initialises monitor pacakge ( see pkg/monitor ) |
C | |-PACKAGES_BOOT :: Start up the optional package environment. |
22 |
|
C | | :: Runtime selection of active packages. |
23 |
|
C | |-PACKAGES_READPARMS :: read all packages input parameter file |
24 |
|
C | | |- ${PKG}_READPARMS |
25 |
C | | |
C | | |
26 |
C | |-INI_GRID :: Control grid array (vert. and hori.) initialisation. |
C | |-SET_PARMS :: Finalise model parameter setting (if fct of pkg usage) |
27 |
C | | | :: Grid arrays are held and described in GRID.h. |
C | | |
28 |
C | | | |
C | |-INI_MODEL_IO :: Initialise Input/Ouput setting |
29 |
C | | |-INI_VERTICAL_GRID :: Initialise vertical grid arrays. |
C | | |-MNC_INIT |
30 |
C | | | |
C | | |-MNC_CW_INIT |
31 |
|
C | | | |
32 |
|
C | | |-MON_INIT :: Initialises monitor pacakge ( see pkg/monitor ) |
33 |
|
C | | |
34 |
|
C | |-INI_GRID :: Control grid array (vert. and hori.) initialisation. |
35 |
|
C | | | :: Grid arrays are held and described in GRID.h. |
36 |
|
C | | |-LOAD_GRID_SPACING :: Load grid spacing (vector) from files |
37 |
|
C | | |-INI_VERTICAL_GRID :: Set up vertical grid and coordinate |
38 |
C | | |-INI_CARTESIAN_GRID :: Cartesian horiz. grid initialisation |
C | | |-INI_CARTESIAN_GRID :: Cartesian horiz. grid initialisation |
39 |
C | | | :: (calculate grid from kernel parameters). |
C | | | :: (calculate grid from kernel parameters). |
40 |
C | | | |
C | | |-INI_SPHERICAL_POLAR_GRID :: Spherical polar horiz. grid setting |
41 |
C | | |-INI_SPHERICAL_POLAR_GRID :: Spherical polar horiz. grid |
C | | | :: (calculate grid from kernel parameters). |
|
C | | | :: initialisation (calculate grid from |
|
|
C | | | :: kernel parameters). |
|
|
C | | | |
|
42 |
C | | |-INI_CURVILINEAR_GRID :: General orthogonal, structured horiz. |
C | | |-INI_CURVILINEAR_GRID :: General orthogonal, structured horiz. |
43 |
C | | :: grid initialisations. ( input from raw |
C | | | :: grid initialisations. ( input from raw |
44 |
C | | :: grid files, LONC.bin, DXF.bin etc... ) |
C | | | :: grid files, LONC.bin, DXF.bin etc... ) |
45 |
C | | |
C | | |-INI_CYLINDER_GRID :: Cylindrical horiz. grid setting |
46 |
C | |-INI_DEPTHS :: Read (from "bathyFile") or set bathymetry/orgography. |
C | | |
47 |
|
C | |-LOAD_REF_FILES :: |
48 |
|
C | |-INI_EOS :: |
49 |
|
C | |-SET_REF_STATE :: |
50 |
|
C | |-SET_GRID_FACTORS :: |
51 |
C | | |
C | | |
52 |
|
C | |-INI_DEPTHS :: Read (from "bathyFile") or set bathymetry/orography. |
53 |
C | |-INI_MASKS_ETC :: Derive horizontal and vertical cell fractions and |
C | |-INI_MASKS_ETC :: Derive horizontal and vertical cell fractions and |
54 |
C | | :: land masking for solid-fluid boundaries. |
C | | :: land masking for solid-fluid boundaries. |
55 |
C | | |
C | | |
56 |
|
C | |-PACKAGES_INIT_FIXED :: do all packages fixed-initialisation setting |
57 |
|
C | | |- ${PKG}_INIT_FIXED |
58 |
|
C | | |
59 |
|
C | |-INI_GLOBAL_DOMAIN |
60 |
C | |-INI_LINEAR_PHSURF :: Set ref. surface Bo_surf |
C | |-INI_LINEAR_PHSURF :: Set ref. surface Bo_surf |
61 |
C | | |
C | | |
62 |
C | |-INI_CORI :: Set coriolis term. zero, f-plane, beta-plane, |
C | |-INI_CORI :: Set coriolis term. zero, f-plane, beta-plane, |
63 |
C | | :: sphere optins are coded. |
C | | :: sphere optins are coded. |
64 |
|
C | |-INI_CG2D :: 2d con. grad solver initialisation. |
65 |
|
C | |-INI_CG3D :: 3d con. grad solver initialisation. |
66 |
C | | |
C | | |
67 |
C | |-PACAKGES_BOOT :: Start up the optional package environment. |
C | |-CONFIG_SUMMARY :: Provide synopsis of kernel setup. Includes |
68 |
C | | :: Runtime selection of active packages. |
C | | :: annotated table of kernel parameter settings. |
69 |
C | | |
C | | |
70 |
C | |-PACKAGES_READPARMS :: Call active package internal parameter load. |
C | |-PACKAGES_CHECK :: call each package configuration checking S/R |
71 |
C | | | |
C | | |- ${PKG}_CHECK |
72 |
C | | |-GMREDI_READPARMS :: GM Package. see pkg/gmredi |
C | | |
73 |
C | | |-KPP_READPARMS :: KPP Package. see pkg/kpp |
C | |-CONFIG_CHECK :: Check config and parameter consistency. |
74 |
C | | |-SHAP_FILT_READPARMS :: Shapiro filter package. see pkg/shap_filt |
C | | |
75 |
C | | |-OBCS_READPARMS :: Open bndy package. see pkg/obcs |
C | |-WRITE_GRID |
76 |
C | | |-AIM_READPARMS :: Intermediate Atmos. pacakage. see pkg/aim |
C | |-CPL_EXCH_CONFIGS |
|
C | | |-COST_READPARMS :: Cost function package. see pkg/cost |
|
|
C | | |-CTRL_INIT :: Control vector support package. see pkg/ctrl |
|
|
C | | |-OPTIM_READPARMS :: Optimisation support package. see pkg/ctrl |
|
|
C | | |-GRDCHK_READPARMS :: Gradient check package. see pkg/grdchk |
|
|
C | | |-ECCO_READPARMS :: ECCO Support Package. see pkg/ecco |
|
|
C | | |
|
|
C | |-PACKAGES_CHECK |
|
|
C | | | |
|
|
C | | |-KPP_CHECK :: KPP Package. pkg/kpp |
|
|
C | | |-OBCS_CHECK :: Open bndy Pacakge. pkg/obcs |
|
|
C | | |-GMREDI_CHECK :: GM Package. pkg/gmredi |
|
|
C | | |
|
|
C | |-PACKAGES_INIT_FIXED |
|
|
C | | |-OBCS_INIT_FIXED :: Open bndy Package. see pkg/obcs |
|
|
C | | |-FLT_INIT :: Floats Package. see pkg/flt |
|
|
C | | |
|
|
C | |-ZONAL_FILT_INIT :: FFT filter Package. see pkg/zonal_filt |
|
|
C | | |
|
|
C | |-INI_CG2D :: 2d con. grad solver initialisation. |
|
|
C | | |
|
|
C | |-INI_CG3D :: 3d con. grad solver initialisation. |
|
|
C | | |
|
|
C | |-CONFIG_SUMMARY :: Provide synopsis of kernel setup. |
|
|
C | :: Includes annotated table of kernel |
|
|
C | :: parameter settings. |
|
77 |
C | |
C | |
78 |
C |-CTRL_UNPACK :: Control vector support package. see pkg/ctrl |
C |-CTRL_UNPACK :: Control vector support package. see pkg/ctrl |
79 |
|
C |-COST_DEPENDENT_INIT :: |
80 |
C | |
C | |
81 |
C |-ADTHE_MAIN_LOOP :: Derivative evaluating form of main time stepping loop |
C |-ADTHE_MAIN_LOOP :: Derivative evaluating form of main time stepping loop |
82 |
C ! :: Auotmatically gerenrated by TAMC/TAF. |
C ! :: Automatically generated by TAMC/TAF. |
|
C | |
|
|
C |-CTRL_PACK :: Control vector support package. see pkg/ctrl |
|
83 |
C | |
C | |
84 |
C |-GRDCHK_MAIN :: Gradient check package. see pkg/grdchk |
C |-THE_MAIN_LOOP :: Main timestepping loop routine. |
|
C | |
|
|
C |-THE_MAIN_LOOP :: Main timestepping loop routine. |
|
85 |
C | | |
C | | |
86 |
C | |-INITIALISE_VARIA :: Set the initial conditions for time evolving |
C | |-INITIALISE_VARIA :: Set the initial conditions for time evolving |
|
C | | | :: variables |
|
|
C | | | |
|
|
C | | |-INI_LINEAR_PHISURF :: Set ref. surface Bo_surf |
|
87 |
C | | | |
C | | | |
88 |
C | | |-INI_CORI :: Set coriolis term. zero, f-plane, beta-plane, |
C #ifdef ALLOW_AUTODIFF |
89 |
C | | | :: sphere optins are coded. |
C | | |-INI_DEPTHS \ |
90 |
C | | | |
C | | |-CTRL_DEPTH_INI \ |
91 |
C | | |-INI_CG2D :: 2d con. grad solver initialisation. |
C | | |-UPDATE_MASKS_ETC } ALLOW_DEPTH_CONTROL case |
92 |
C | | |-INI_CG3D :: 3d con. grad solver initialisation. |
C | | |-UPDATE_CG2D / |
93 |
C | | |-INI_MIXING :: Initialise diapycnal diffusivity. |
C #endif |
94 |
C | | |-INI_DYNVARS :: Initialise to zero all DYNVARS.h arrays (dynamical |
C | | |-INI_NLFS_VARS :: Initialise all Non-Lin Free-Surf arrays (SURFACE.h) |
95 |
C | | | :: fields). |
C | | |-INI_DYNVARS :: Initialise to zero all DYNVARS.h arrays |
96 |
|
C | | |-INI_NH_VARS :: Initialise to zero all NH_VARS.h arrays |
97 |
C | | | |
C | | | |
98 |
C | | |-INI_FIELDS :: Control initialising model fields to non-zero |
C | | |-INI_FIELDS :: Control initialising model fields to non-zero |
99 |
C | | | |-INI_VEL :: Initialize 3D flow field. |
C | | | |-INI_VEL :: Initialize 3D flow field. |
100 |
C | | | |-INI_THETA :: Set model initial temperature field. |
C | | | |-INI_THETA :: Set model initial temperature field. |
101 |
C | | | |-INI_SALT :: Set model initial salinity field. |
C | | | |-INI_SALT :: Set model initial salinity field. |
102 |
C | | | |-INI_PSURF :: Set model initial free-surface height/pressure. |
C | | | |-INI_PSURF :: Set model initial free-surface height/pressure. |
103 |
|
C | | | |-READ_PICKUP |
104 |
C | | | |
C | | | |
105 |
C | | |-INI_TR1 :: Set initial tracer 1 distribution. |
C | | |-INI_MIXING :: Initialise diapycnal diffusivity. |
106 |
|
C | | |-TAUEDDY_INIT_VARIA |
107 |
C | | | |
C | | | |
108 |
C | | |-THE_CORRECTION_STEP :: Step forward to next time step. |
C | | |-INI_FORCING :: initialise forcing fields |
109 |
C | | | | :: Here applied to move restart conditions |
C | | |-INI_FORCING :: Set model initial forcing fields. |
110 |
C | | | | :: (saved in mid timestep) to correct level in |
C | | | | :: Either set in-line or from file as shown. |
|
C | | | | :: time (only used for pre-c35). |
|
|
C | | | | |
|
|
C | | | |-CALC_GRAD_PHI_SURF :: Return DDx and DDy of surface pressure |
|
|
C | | | |-CORRECTION_STEP :: Pressure correction to momentum |
|
|
C | | | |-CYCLE_TRACER :: Move tracers forward in time. |
|
|
C | | | |-OBCS_APPLY :: Open bndy package. see pkg/obcs |
|
|
C | | | |-SHAP_FILT_APPLY :: Shapiro filter package. see pkg/shap_filt |
|
|
C | | | |-ZONAL_FILT_APPLY :: FFT filter package. see pkg/zonal_filt |
|
|
C | | | |-CONVECTIVE_ADJUSTMENT :: Control static instability mixing. |
|
|
C | | | | |-FIND_RHO :: Find adjacent densities. |
|
|
C | | | | |-CONVECT :: Mix static instability. |
|
|
C | | | | |-TIMEAVE_CUMULATE :: Update convection statistics. |
|
|
C | | | | |
|
|
C | | | |-CALC_EXACT_ETA :: Change SSH to flow divergence. |
|
|
C | | | |
|
|
C | | |-CONVECTIVE_ADJUSTMENT_INI :: Control static instability mixing |
|
|
C | | | | :: Extra time history interactions. |
|
|
C | | | | |
|
|
C | | | |-FIND_RHO :: Find adjacent densities. |
|
|
C | | | |-CONVECT :: Mix static instability. |
|
|
C | | | |-TIMEAVE_CUMULATE :: Update convection statistics. |
|
|
C | | | |
|
|
C | | |-PACKAGES_INIT_VARIABLES :: Does initialisation of time evolving |
|
|
C | | | | :: package data. |
|
|
C | | | | |
|
|
C | | | |-GMREDI_INIT :: GM package. ( see pkg/gmredi ) |
|
|
C | | | |-KPP_INIT :: KPP package. ( see pkg/kpp ) |
|
|
C | | | |-KPP_OPEN_DIAGS |
|
|
C | | | |-OBCS_INIT_VARIABLES :: Open bndy. package. ( see pkg/obcs ) |
|
|
C | | | |-AIM_INIT :: Interm. atmos package. ( see pkg/aim ) |
|
|
C | | | |-CTRL_MAP_INI :: Control vector package.( see pkg/ctrl ) |
|
|
C | | | |-COST_INIT :: Cost function package. ( see pkg/cost ) |
|
|
C | | | |-ECCO_INIT :: ECCO support package. ( see pkg/ecco ) |
|
|
C | | | |-INI_FORCING :: Set model initial forcing fields. |
|
|
C | | | | :: Either set in-line or from file as shown. |
|
111 |
C | | | |-READ_FLD_XY_RS(zonalWindFile) |
C | | | |-READ_FLD_XY_RS(zonalWindFile) |
112 |
C | | | |-READ_FLD_XY_RS(meridWindFile) |
C | | | |-READ_FLD_XY_RS(meridWindFile) |
113 |
C | | | |-READ_FLD_XY_RS(surfQFile) |
C | | | |-READ_FLD_XY_RS(surfQnetFile) |
114 |
C | | | |-READ_FLD_XY_RS(EmPmRfile) |
C | | | |-READ_FLD_XY_RS(EmPmRfile) |
115 |
C | | | |-READ_FLD_XY_RS(thetaClimFile) |
C | | | |-READ_FLD_XY_RS(thetaClimFile) |
116 |
C | | | |-READ_FLD_XY_RS(saltClimFile) |
C | | | |-READ_FLD_XY_RS(saltClimFile) |
117 |
C | | | |-READ_FLD_XY_RS(surfQswFile) |
C | | | |-READ_FLD_XY_RS(surfQswFile) |
118 |
C | | | |
C | | | |
119 |
|
C | | |-AUTODIFF_INIT_VARIA |
120 |
|
C | | | |
121 |
|
C | | |-PACKAGES_INIT_VARIABLES :: Does initialisation of time evolving |
122 |
|
C | | | | ${PKG}_INIT_VARIA :: package data. |
123 |
|
C | | | |
124 |
|
C | | |-COST_INIT_VARIA |
125 |
|
C | | |-CONVECTIVE_ADJUSTMENT_INI |
126 |
|
C | | | |
127 |
|
C | | |-CALC_R_STAR :: Calculate the new level thickness factor (r* coord) |
128 |
|
C | | |-UPDATE_R_STAR :: Update the level thickness fraction (r* coord). |
129 |
|
C | | |-UPDATE_SIGMA :: Update the level thickness fraction (sigma-coord). |
130 |
C | | |-CALC_SURF_DR :: Calculate the new surface level thickness. |
C | | |-CALC_SURF_DR :: Calculate the new surface level thickness. |
131 |
C | | |-UPDATE_SURF_DR :: Update the surface-level thickness fraction. |
C | | |-UPDATE_SURF_DR :: Update the surface-level thickness fraction. |
132 |
|
C | | | |
133 |
C | | |-UPDATE_CG2D :: Update 2d conjugate grad. for Free-Surf. |
C | | |-UPDATE_CG2D :: Update 2d conjugate grad. for Free-Surf. |
134 |
C | | |-STATE_SUMMARY :: Summarize model prognostic variables. |
C | | | |
135 |
C | | |-TIMEAVE_STATVARS :: Time averaging package ( see pkg/timeave ). |
C | | |-INTEGR_CONTINUITY :: Integrate the continuity Eqiuation |
136 |
C | | |
C | | | |-INTEGRATE_FOR_W :: Integrate for vertical velocity |
137 |
C | |-WRITE_STATE :: Controlling routine for IO to dump model state. |
C | | | |-OBCS_APPLY_W :: Open boundary package (see pkg/obcs). |
138 |
C | | |-WRITE_REC_XYZ_RL :: Single file I/O |
C | | | |-UPDATE_ETAH :: Update Surface height/pressure |
139 |
C | | |-WRITE_FLD_XYZ_RL :: Multi-file I/O |
C | | | |
140 |
C | | |
C | | |-CALC_R_STAR :: Calculate the new level thickness factor (r* coord) |
141 |
C | |-MONITOR :: Monitor state ( see pkg/monitor ) |
C | | |-CALC_SURF_DR :: Calculate the new surface level thickness. |
142 |
C | |-CTRL_MAP_FORCING :: Control vector support package. ( see pkg/ctrl ) |
C | | | |
143 |
C====|>| |
C | | |-STATE_SUMMARY :: Summarise model prognostic variables. |
144 |
|
C | | | |
145 |
|
C | | |-MONITOR :: Monitor state (see pkg/monitor) |
146 |
|
C | | | |
147 |
|
C | | |-DO_STATEVARS_TAVE :: Time averaging package ( see pkg/timeave ). |
148 |
|
C | | | |-TIMEAVE_STATVARS :: |
149 |
|
C | | | |-PTRACERS_TIMEAVE :: |
150 |
|
C | | | |
151 |
|
C | | |-DO_THE_MODEL_IO :: Controlling routine for IO |
152 |
|
C | | | |-WRITE_STATE :: Write model state variables. |
153 |
|
C | | | |-TIMEAVE_STATV_WRITE :: Write Time averaged output (see pkg/timeave) |
154 |
|
C | | | |-FIZHI_WRITE_STATE |
155 |
|
C | | | |-AIM_WRITE_TAVE |
156 |
|
C | | | |-LAND_OUTPUT |
157 |
|
C | | | |-OBCS_OUTPUT |
158 |
|
C | | | |-GMREDI_OUTPUT |
159 |
|
C | | | |-KPP_OUTPUT |
160 |
|
C | | | |-PP81_OUTPUT |
161 |
|
C | | | |-KL10_OUTPUT |
162 |
|
C | | | |-MY82_OUTPUT |
163 |
|
C | | | |-OPPS_OUTPUT |
164 |
|
C | | | |-GGL90_OUTPUT |
165 |
|
C | | | |-SBO_CALC |
166 |
|
C | | | |-SBO_OUTPUT |
167 |
|
C | | | |-SEAICE_OUTPUT |
168 |
|
C | | | |-SHELFICE_OUTPUT |
169 |
|
C | | | |-BULKF_OUTPUT |
170 |
|
C | | | |-THSICE_OUTPUT |
171 |
|
C | | | |-PTRACERS_OUTPUT |
172 |
|
C | | | |-MATRIX_OUTPUT |
173 |
|
C | | | |-GCHEM_OUTPUT |
174 |
|
C | | | |-CPL_OUTPUT |
175 |
|
C | | | |-LAYERS_CALC |
176 |
|
C | | | |-LAYERS_OUTPUT |
177 |
|
C | | | |-DIAGNOSTICS_WRITE :: Write pkg/diagnostics output |
178 |
|
C | | | |
179 |
C====|>| **************************** |
C====|>| **************************** |
180 |
C====|>| BEGIN MAIN TIMESTEPPING LOOP |
C====|>| BEGIN MAIN TIMESTEPPING LOOP |
181 |
C====|>| **************************** |
C====|>| **************************** |
182 |
C====|>| |
C | |-COST_AVERAGESFIELDS :: time-averaged Cost function terms (see pkg/cost) |
183 |
C/\ | |-FORWARD_STEP :: Step forward a time-step ( AT LAST !!! ) |
C | |-PROFILES_INLOOP :: |
184 |
|
C | / |
185 |
|
C | |-MAIN_DO_LOOP :: Open-AD case: Main timestepping loop routine |
186 |
|
C | \ otherwise: just call FORWARD_STEP |
187 |
|
C | | |
188 |
|
C/\ | |-FORWARD_STEP :: Step forward a time-step ( AT LAST !!! ) |
189 |
C/\ | | | |
C/\ | | | |
190 |
C/\ | | |-DUMMY_IN_STEPPING :: autodiff package ( pkg/autoduff ). |
C/\ | | |-LOAD_FIELDS_DRIVER :: control loading of input fields from files |
|
C/\ | | |-CALC_EXACT_ETA :: Change SSH to flow divergence. |
|
|
C/\ | | |-CALC_SURF_DR :: Calculate the new surface level thickness. |
|
|
C/\ | | |-EXF_GETFORCING :: External forcing package. ( pkg/exf ) |
|
|
C/\ | | |-EXTERNAL_FIELDS_LOAD :: Control loading time dep. external data. |
|
|
C/\ | | | | :: Simple interpolcation between end-points |
|
|
C/\ | | | | :: for forcing datasets. |
|
|
C/\ | | | | |
|
|
C/\ | | | |-EXCH :: Sync forcing. in overlap regions. |
|
191 |
C/\ | | | |
C/\ | | | |
192 |
C/\ | | |-THERMODYNAMICS :: theta, salt + tracer equations driver. |
C/\ | | |-CTRL_MAP_FORCING :: Control vector support package. (see pkg/ctrl) |
193 |
C/\ | | | | |
C/\ | | | |
194 |
C/\ | | | |-INTEGRATE_FOR_W :: Integrate for vertical velocity. |
C/\ | | |-CPL_EXPORT_MY_DATA :: Send coupling fields to coupler |
195 |
C/\ | | | |-OBCS_APPLY_W :: Open bndy. package ( see pkg/obcs ). |
C/\ | | |-CPL_IMPORT_EXTERNAL_DATA :: Receive coupling fields from coupler |
196 |
C/\ | | | |-FIND_RHO :: Calculates [rho(S,T,z)-Rhonil] of a slice |
C/\ | | | |
197 |
C/\ | | | |-GRAD_SIGMA :: Calculate isoneutral gradients |
C/\ | | |-DO_ATMOSPHERIC_PHYS :: Atmospheric physics computation |
198 |
C/\ | | | |-CALC_IVDC :: Set Implicit Vertical Diffusivity for Convection |
C/\ | | | |
199 |
C/\ | | | | |
C/\ | | |-DO_OCEANIC_PHYS :: Oceanic (& seaice) physics computation |
200 |
C/\ | | | |-OBCS_CALC :: Open bndy. package ( see pkg/obcs ). |
C/\ | | | |-OBCS_CALC :: Open boundary. package (see pkg/obcs). |
201 |
C/\ | | | |-EXTERNAL_FORCING_SURF:: Accumulates appropriately dimensioned |
C/\ | | | |-EXTERNAL_FORCING_SURF:: Accumulates appropriately dimensioned |
202 |
C/\ | | | | :: forcing terms. |
C/\ | | | | :: forcing terms. |
203 |
|
C/\ | | | |
204 |
|
C/\ | | |-GCHEM_CALC_TENDENCY :: geochemistry driver routine (see pkg/gchem) |
205 |
|
C/\ | | | |
206 |
|
C/\ | | |-THERMODYNAMICS :: (synchronous time-stepping) |
207 |
|
C/\ | | | theta, salt + tracer equations driver. |
208 |
|
C/\ | | | |-GMREDI_RESIDUAL_FLOW :: Get the flow field used to advect tracer |
209 |
|
C/\ | | | | |
210 |
|
C/\ | | | |-TEMP_INTEGRATE :: Step forward Prognostic Eq for Temperature. |
211 |
|
C/\ | | | | |-GAD_ADVECTION :: Generalised advection driver (multi-dim |
212 |
|
C/\ | | | | | advection case) (see pkg/gad). |
213 |
|
C/\ | | | | |-GAD_CALC_RHS :: Calculate Advection-Diffusion tendency terms |
214 |
|
C/\ | | | | |-TIMESTEP_TRACER :: Step tracer field forward in time |
215 |
|
C/\ | | | | |-IMPLDIFF :: Solve vertical implicit diffusion equation. |
216 |
|
C/\ | | | | |-CYCLE_AB_TRACER :: Cycle time-stepping arrays for tracer field |
217 |
|
C/\ | | | | |
218 |
|
C/\ | | | |-SALT_INTEGRATE :: Step forward Prognostic Eq for Salinity. |
219 |
|
C/\ | | | | | same sequence of calls as in TEMP_INTEGRATE |
220 |
|
C/\ | | | | |
221 |
|
C/\ | | | |-PTRACERS_INTEGRATE :: Integrate other tracer(s) (see pkg/ptracers). |
222 |
|
C/\ | | | | | same sequence of calls as in TEMP_INTEGRATE |
223 |
|
C/\ | | | | |-OBCS_APPLY_PTRACER :: Open boundary package for pTracers |
224 |
C/\ | | | | |
C/\ | | | | |
225 |
C/\ | | | |-GMREDI_CALC_TENSOR :: GM package ( see pkg/gmredi ). |
C/\ | | | |-OBCS_APPLY_TS :: Open boundary package (see pkg/obcs ). |
|
C/\ | | | |-GMREDI_CALC_TENSOR_DUMMY :: GM package ( see pkg/gmredi ). |
|
|
C/\ | | | |-KPP_CALC :: KPP package ( see pkg/kpp ). |
|
|
C/\ | | | |-KPP_CALC_DUMMY :: KPP package ( see pkg/kpp ). |
|
|
C/\ | | | |-AIM_DO_ATMOS_PHYSICS :: Intermed. atmos package ( see pkg/aim ). |
|
|
C/\ | | | |-GAD_ADVECTION :: Generalised advection driver (multi-dim |
|
|
C/\ | | | | advection case) (see pkg/gad). |
|
|
C/\ | | | |-CALC_COMMON_FACTORS :: Calculate common data (such as volume flux) |
|
|
C/\ | | | |-CALC_DIFFUSIVITY :: Calculate net vertical diffusivity |
|
|
C/\ | | | | | |
|
|
C/\ | | | | |-GMREDI_CALC_DIFF :: GM package ( see pkg/gmredi ). |
|
|
C/\ | | | | |-KPP_CALC_DIFF :: KPP package ( see pkg/kpp ). |
|
|
C/\ | | | | |
|
|
C/\ | | | |-CALC_GT :: Calculate the temperature tendency terms |
|
|
C/\ | | | | | |
|
|
C/\ | | | | |-GAD_CALC_RHS :: Generalised advection package |
|
|
C/\ | | | | | :: ( see pkg/gad ) |
|
|
C/\ | | | | |-EXTERNAL_FORCING_T :: Problem specific forcing for temperature. |
|
|
C/\ | | | | |-ADAMS_BASHFORTH2 :: Extrapolate tendencies forward in time. |
|
|
C/\ | | | | |-FREESURF_RESCALE_G :: Re-scale Gt for free-surface height. |
|
|
C/\ | | | | |
|
|
C/\ | | | |-TIMESTEP_TRACER :: Step tracer field forward in time |
|
|
C/\ | | | | |
|
|
C/\ | | | |-CALC_GS :: Calculate the salinity tendency terms |
|
|
C/\ | | | | | |
|
|
C/\ | | | | |-GAD_CALC_RHS :: Generalised advection package |
|
|
C/\ | | | | | :: ( see pkg/gad ) |
|
|
C/\ | | | | |-EXTERNAL_FORCING_S :: Problem specific forcing for salt. |
|
|
C/\ | | | | |-ADAMS_BASHFORTH2 :: Extrapolate tendencies forward in time. |
|
|
C/\ | | | | |-FREESURF_RESCALE_G :: Re-scale Gs for free-surface height. |
|
|
C/\ | | | | |
|
|
C/\ | | | |-TIMESTEP_TRACER :: Step tracer field forward in time |
|
|
C/\ | | | | |
|
|
C/\ | | | |-CALC_GTR1 :: Calculate other tracer(s) tendency terms |
|
|
C/\ | | | | | |
|
|
C/\ | | | | |-GAD_CALC_RHS :: Generalised advection package |
|
|
C/\ | | | | | :: ( see pkg/gad ) |
|
|
C/\ | | | | |-EXTERNAL_FORCING_TR:: Problem specific forcing for tracer. |
|
|
C/\ | | | | |-ADAMS_BASHFORTH2 :: Extrapolate tendencies forward in time. |
|
|
C/\ | | | | |-FREESURF_RESCALE_G :: Re-scale Gs for free-surface height. |
|
|
C/\ | | | | |
|
|
C/\ | | | |-TIMESTEP_TRACER :: Step tracer field forward in time |
|
|
C/\ | | | |-OBCS_APPLY_TS :: Open bndy. package (see pkg/obcs ). |
|
|
C/\ | | | |-FREEZE :: Limit range of temperature. |
|
|
C/\ | | | | |
|
|
C/\ | | | |-IMPLDIFF :: Solve vertical implicit diffusion equation. |
|
|
C/\ | | | |-OBCS_APPLY_TS :: Open bndy. package (see pkg/obcs ). |
|
|
C/\ | | | | |
|
|
C/\ | | | |-AIM_AIM2DYN_EXCHANGES :: Inetermed. atmos (see pkg/aim). |
|
|
C/\ | | | |-EXCH :: Update overlaps |
|
226 |
C/\ | | | |
C/\ | | | |
227 |
C/\ | | |-DYNAMICS :: Momentum equations driver. |
C/\ | | |-DYNAMICS :: Momentum equations driver. |
228 |
C/\ | | | | |
C/\ | | | | |
229 |
C/\ | | | |-CALC_GRAD_PHI_SURF :: Calculate the gradient of the surface |
C/\ | | | |-CALC_GRAD_PHI_SURF :: Calculate the gradient of the surface |
230 |
C/\ | | | | Potential anomaly. |
C/\ | | | | Potential anomaly. |
231 |
C/\ | | | |-CALC_VISCOSITY :: Calculate net vertical viscosity |
C/\ | | | |-CALC_VISCOSITY :: Calculate net vertical viscosity |
232 |
C/\ | | | | |-KPP_CALC_VISC :: KPP package ( see pkg/kpp ). |
C/\ | | | |-CALC_PHI_HYD :: Integrate the hydrostatic relation. |
233 |
C/\ | | | | |
C/\ | | | |-MOM_FLUXFORM :: Flux form mom eqn. package (pkg/mom_fluxform) |
234 |
C/\ | | | |-CALC_PHI_HYD :: Integrate the hydrostatic relation. |
C/\ | | | |-MOM_VECINV :: Vector invariant form mom eqn (pkg/mom_vecinv) |
235 |
C/\ | | | |-MOM_FLUXFORM :: Flux form mom eqn. package ( see |
C/\ | | | |-TIMESTEP :: Step momentum fields forward in time |
236 |
C/\ | | | | pkg/mom_fluxform ). |
C/\ | | | |-OBCS_APPLY_UV :: Open boundary package (see pkg/obcs). |
237 |
C/\ | | | |-MOM_VECINV :: Vector invariant form mom eqn. package ( see |
C/\ | | | |-IMPLDIFF :: Solve vertical implicit diffusion equation. |
238 |
C/\ | | | | pkg/mom_vecinv ). |
C/\ | | | |-CALC_GW :: vert. momentum tendency terms (Non-Hydrostatic) |
|
C/\ | | | |-TIMESTEP :: Step momentum fields forward in time |
|
|
C/\ | | | |-OBCS_APPLY_UV :: Open bndy. package (see pkg/obcs ). |
|
|
C/\ | | | | |
|
|
C/\ | | | |-IMPLDIFF :: Solve vertical implicit diffusion equation. |
|
|
C/\ | | | |-OBCS_APPLY_UV :: Open bndy. package (see pkg/obcs ). |
|
|
C/\ | | | | |
|
|
C/\ | | | |-TIMEAVE_CUMUL_1T :: Time averaging package ( see pkg/timeave ). |
|
|
C/\ | | | |-TIMEAVE_CUMUATE :: Time averaging package ( see pkg/timeave ). |
|
|
C/\ | | | |-DEBUG_STATS_RL :: Quick debug package ( see pkg/debug ). |
|
|
C/\ | | | |
|
|
C/\ | | |-CALC_GW :: vert. momentum tendency terms ( NH, QH only ). |
|
239 |
C/\ | | | |
C/\ | | | |
240 |
C/\ | | |-UPDATE_SURF_DR :: Update the surface-level thickness fraction. |
C/\ | | |-UPDATE_SURF_DR :: Update the surface-level thickness fraction. |
241 |
C/\ | | | |
C/\ | | |-UPDATE_R_STAR :: Update the level thickness fraction. |
242 |
C/\ | | |-UPDATE_CG2D :: Update 2d conjugate grad. for Free-Surf. |
C/\ | | |-UPDATE_CG2D :: Update 2d conjugate grad. for Free-Surf. |
243 |
C/\ | | | |
C/\ | | | |
244 |
C/\ | | |-SOLVE_FOR_PRESSURE :: Find surface pressure. |
C/\ | | |-SOLVE_FOR_PRESSURE :: Find surface pressure. |
|
C/\ | | | |-CALC_DIV_GHAT :: Form the RHS of the surface pressure eqn. |
|
245 |
C/\ | | | |-CG2D :: Two-dim pre-con. conjugate-gradient. |
C/\ | | | |-CG2D :: Two-dim pre-con. conjugate-gradient. |
246 |
C/\ | | | |-CG3D :: Three-dim pre-con. conjugate-gradient solver. |
C/\ | | | |-CG3D :: Three-dim pre-con. conjugate-gradient solver. |
247 |
C/\ | | | |
C/\ | | | |
248 |
C/\ | | |-THE_CORRECTION_STEP :: Step forward to next time step. |
C/\ | | |-MOMENTUM_CORRECTION_STEP :: Finalise momentum stepping |
|
C/\ | | | | |
|
249 |
C/\ | | | |-CALC_GRAD_PHI_SURF :: Return DDx and DDy of surface pressure |
C/\ | | | |-CALC_GRAD_PHI_SURF :: Return DDx and DDy of surface pressure |
250 |
C/\ | | | |-CORRECTION_STEP :: Pressure correction to momentum |
C/\ | | | |-CORRECTION_STEP :: Pressure correction to momentum |
251 |
C/\ | | | |-CYCLE_TRACER :: Move tracers forward in time. |
C/\ | | | |-OBCS_APPLY_UV :: Open boundary package (see pkg/obcs). |
|
C/\ | | | |-OBCS_APPLY :: Open bndy package. see pkg/obcs |
|
|
C/\ | | | |-SHAP_FILT_APPLY :: Shapiro filter package. see pkg/shap_filt |
|
|
C/\ | | | |-ZONAL_FILT_APPLY :: FFT filter package. see pkg/zonal_filt |
|
|
C/\ | | | |-CONVECTIVE_ADJUSTMENT :: Control static instability mixing. |
|
|
C/\ | | | | |-FIND_RHO :: Find adjacent densities. |
|
|
C/\ | | | | |-CONVECT :: Mix static instability. |
|
|
C/\ | | | | |-TIMEAVE_CUMULATE :: Update convection statistics. |
|
|
C/\ | | | | |
|
|
C/\ | | | |-CALC_EXACT_ETA :: Change SSH to flow divergence. |
|
252 |
C/\ | | | |
C/\ | | | |
253 |
C/\ | | |-DO_FIELDS_BLOCKING_EXCHANGES :: Sync up overlap regions. |
C/\ | | |-INTEGR_CONTINUITY :: Integrate continuity equation |
|
C/\ | | | |-EXCH |
|
254 |
C/\ | | | |
C/\ | | | |
255 |
C/\ | | |-FLT_MAIN :: Float package ( pkg/flt ). |
C/\ | | |-THERMODYNAMICS :: (staggered time-stepping) |
256 |
|
C/\ | | | theta, salt + tracer equations driver. |
257 |
C/\ | | | |
C/\ | | | |
258 |
C/\ | | |-MONITOR :: Monitor package ( pkg/monitor ). |
C/\ | | |-TRACERS_CORRECTION_STEP :: Finalise tracer steppin: |
259 |
|
C/\ | | | apply filter, conv.adjustment |
260 |
C/\ | | | |
C/\ | | | |
261 |
|
C/\ | | |-GCHEM_FORCING_SEP :: tracer forcing for gchem pkg (if tracer |
262 |
|
C/\ | | | dependent tendencies calculated separately) |
263 |
|
C/\ | | | |
264 |
|
C/\ | | |-DO_FIELDS_BLOCKING_EXCHANGES :: Sync up overlap regions. |
265 |
|
C/\ | | | |
266 |
|
C/\ | | |-MONITOR :: Monitor package (pkg/monitor). |
267 |
C/\ | | |-DO_THE_MODEL_IO :: Standard diagnostic I/O. |
C/\ | | |-DO_THE_MODEL_IO :: Standard diagnostic I/O. |
268 |
C/\ | | | |-WRITE_STATE :: Core state I/O |
C/\ | | | |
269 |
C/\ | | | |-TIMEAVE_STATV_WRITE :: Time averages. see pkg/timeave |
C/\ | | |-DO_WRITE_PICKUP :: Write restart files. |
270 |
C/\ | | | |-AIM_WRITE_DIAGS :: Intermed. atmos diags. see pkg/aim |
C | | |
|
C/\ | | | |-GMREDI_DIAGS :: GM diags. see pkg/gmredi |
|
|
C/\ | | | |-KPP_DO_DIAGS :: KPP diags. see pkg/kpp |
|
|
C/\ | | | |
|
|
C/\ | | |-WRITE_CHECKPOINT :: Do I/O for restart files. |
|
|
C/\ | | |
|
|
C/\ | |-COST_TILE :: Cost function package. ( see pkg/cost ) |
|
|
C<===|=| |
|
271 |
C<===|=| ************************** |
C<===|=| ************************** |
272 |
C<===|=| END MAIN TIMESTEPPING LOOP |
C<===|=| END MAIN TIMESTEPPING LOOP |
273 |
C<===|=| ************************** |
C<===|=| ************************** |
274 |
C<===|=| |
C | | |
275 |
C | |-COST_FINAL :: Cost function package. ( see pkg/cost ) |
C | |-COST_AVERAGESFIELDS :: time-averaged Cost function terms (see pkg/cost) |
276 |
|
C | |-PROFILES_INLOOP :: |
277 |
|
C | |-COST_FINAL :: Cost function package. (see pkg/cost) |
278 |
|
C | |
279 |
|
C |-CTRL_PACK :: Control vector support package. see pkg/ctrl |
280 |
C | |
C | |
281 |
C |-WRITE_CHECKPOINT :: Final state storage, for restart. |
C |-GRDCHK_MAIN :: Gradient check package. see pkg/grdchk |
282 |
C | |
C | |
283 |
C |-TIMER_PRINTALL :: Computational timing summary |
C |-TIMER_PRINTALL :: Computational timing summary |
284 |
C | |
C | |
285 |
C |-COMM_STATS :: Summarise inter-proc and inter-thread communication |
C |-COMM_STATS :: Summarise inter-proc and inter-thread communication |
286 |
C :: events. |
C :: events. |