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

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

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

Revision 1.9 - (show annotations) (download)
Wed Sep 26 18:09:16 2001 UTC (22 years, 7 months ago) by cnh
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint51k_post, checkpoint47e_post, checkpoint57m_post, checkpoint52l_pre, checkpoint44e_post, hrcube4, hrcube5, checkpoint46l_post, checkpoint57g_pre, checkpoint46g_pre, checkpoint47c_post, release1_p13_pre, checkpoint50c_post, checkpoint57s_post, checkpoint58b_post, checkpoint57b_post, checkpoint46f_post, checkpoint52d_pre, checkpoint57g_post, checkpoint48e_post, checkpoint56b_post, checkpoint50g_post, checkpoint57y_post, checkpoint46b_post, checkpoint52j_pre, checkpoint43a-release1mods, checkpoint51o_pre, checkpoint44g_post, checkpoint54d_post, checkpoint48c_post, checkpoint54e_post, ecco_c50_e32, ecco_c50_e33, ecco_c50_e30, ecco_c50_e31, release1_p13, checkpoint62c, checkpoint51l_post, checkpoint48i_post, checkpoint57r_post, checkpoint46l_pre, checkpoint57d_post, checkpoint57i_post, checkpoint50d_pre, checkpoint52k_post, chkpt44d_post, checkpoint59, checkpoint58, checkpoint55, checkpoint54, checkpoint57, checkpoint56, checkpoint51, checkpoint53, checkpoint52, release1_p8, release1_p9, checkpoint50d_post, release1_p1, release1_p2, release1_p3, release1_p4, release1_p5, release1_p6, release1_p7, checkpoint58f_post, checkpoint52f_post, checkpoint57n_post, checkpoint58d_post, checkpoint58a_post, checkpoint50b_pre, checkpoint44e_pre, checkpoint62p, checkpoint57z_post, checkpoint54f_post, checkpoint51f_post, checkpoint62a, release1_b1, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, ecco_c51_e34d, ecco_c51_e34e, ecco_c51_e34f, ecco_c51_e34g, ecco_c51_e34a, ecco_c51_e34b, ecco_c51_e34c, checkpoint58y_post, checkpoint48b_post, checkpoint43, checkpoint51d_post, checkpoint48c_pre, checkpoint55a_post, checkpoint51t_post, checkpoint58t_post, checkpoint51n_post, release1_chkpt44d_post, checkpoint55i_post, checkpoint57l_post, checkpoint52i_pre, checkpoint57h_post, hrcube_2, hrcube_3, checkpoint51s_post, checkpoint57t_post, checkpoint55c_post, checkpoint48d_pre, checkpoint51j_post, checkpoint47i_post, checkpoint52e_pre, checkpoint57v_post, release1_p11, checkpoint57f_post, checkpoint52e_post, checkpoint51n_pre, checkpoint47d_post, icebear5, icebear4, icebear3, icebear2, checkpoint53d_post, checkpoint46d_pre, checkpoint60, checkpoint61, checkpoint62, checkpoint57a_post, checkpoint48d_post, release1-branch_tutorials, checkpoint57h_pre, checkpoint57x_post, checkpoint48f_post, checkpoint45d_post, checkpoint52b_pre, checkpoint54b_post, checkpoint46j_pre, checkpoint58w_post, ecco_c50_e28, checkpoint51l_pre, checkpoint52m_post, checkpoint47d_pre, checkpoint57y_pre, chkpt44a_post, checkpoint55g_post, checkpoint44h_pre, checkpoint48h_post, checkpoint51q_post, ecco_c50_e29, checkpoint51b_pre, checkpoint46a_post, checkpoint47g_post, checkpoint52b_post, checkpoint52c_post, checkpoint46j_post, checkpoint51h_pre, checkpoint46k_post, checkpoint46b_pre, chkpt44c_pre, checkpoint58o_post, checkpoint52h_pre, checkpoint45a_post, checkpoint57c_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, hrcube_1, checkpoint58p_post, checkpoint58q_post, checkpoint51m_post, checkpoint55d_post, ecco_c44_e19, ecco_c44_e18, ecco_c44_e17, ecco_c44_e16, release1_p12, checkpoint58e_post, release1_p10, mitgcm_mapl_00, release1_p16, release1_p17, release1_p14, release1_p15, checkpoint58m_post, checkpoint47a_post, ecco_c50_e33a, checkpoint53c_post, checkpoint55d_pre, checkpoint57c_pre, checkpoint58r_post, checkpoint55j_post, branchpoint-genmake2, checkpoint54a_post, checkpoint46e_pre, checkpoint55h_post, checkpoint58n_post, checkpoint51r_post, checkpoint45b_post, checkpoint51i_post, checkpoint57e_post, release1-branch-end, release1_final_v1, checkpoint55b_post, checkpoint51b_post, release1_p12_pre, checkpoint46c_pre, checkpoint53a_post, checkpoint44f_post, checkpoint47b_post, checkpoint44b_post, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint55f_post, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, ecco_c51_e34, checkpoint46h_pre, checkpoint52d_post, checkpoint53g_post, checkpoint46m_post, checkpoint57p_post, checkpint57u_post, checkpoint46a_pre, checkpoint50c_pre, checkpoint45c_post, checkpoint57q_post, ecco_ice2, ecco_ice1, checkpoint44h_post, eckpoint57e_pre, checkpoint46g_post, checkpoint51c_post, checkpoint58k_post, checkpoint52a_pre, checkpoint62b, ecco_c44_e22, checkpoint58v_post, checkpoint50h_post, checkpoint52i_post, checkpoint50e_pre, checkpoint50i_post, ecco_c44_e25, checkpoint51i_pre, checkpoint48a_post, checkpoint56a_post, checkpoint58l_post, checkpoint53f_post, checkpoint47j_post, checkpoint54a_pre, checkpoint53b_pre, branch-exfmods-tag, checkpoint57h_done, checkpoint52j_post, checkpoint47f_post, checkpoint50e_post, chkpt44a_pre, checkpoint46i_post, ecco_c44_e23, ecco_c44_e20, ecco_c44_e21, ecco_c44_e26, ecco_c44_e27, ecco_c44_e24, checkpoint57j_post, checkpoint57f_pre, checkpoint61f, checkpoint46c_post, checkpoint58g_post, ecco-branch-mod1, ecco-branch-mod2, ecco-branch-mod3, ecco-branch-mod4, ecco-branch-mod5, branch-netcdf, checkpoint52l_post, checkpoint58x_post, checkpoint61n, checkpoint52n_post, checkpoint46e_post, checkpoint59j, release1_beta1, checkpoint58h_post, checkpoint56c_post, checkpoint58j_post, checkpoint51e_post, checkpoint44b_pre, checkpoint42, checkpoint57a_pre, checkpoint41, checkpoint46, checkpoint47, checkpoint44, checkpoint45, checkpoint48, checkpoint49, checkpoint57o_post, checkpoint46h_post, checkpoint51o_post, checkpoint61q, checkpoint50, checkpoint57k_post, checkpoint51f_pre, chkpt44c_post, checkpoint61z, checkpoint53b_post, checkpoint47h_post, checkpoint52a_post, checkpoint57w_post, checkpoint61e, checkpoint44f_pre, checkpoint58i_post, checkpoint51g_post, ecco_c52_e35, checkpoint46d_post, checkpoint50b_post, checkpoint58c_post, checkpoint58u_post, release1-branch_branchpoint, checkpoint52f_pre, checkpoint53d_pre, checkpoint58s_post, checkpoint55e_post, checkpoint61g, checkpoint61d, checkpoint54c_post, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint51a_post, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint51p_post, checkpoint48g_post, checkpoint61x, checkpoint61y, checkpoint51u_post
Branch point for: c24_e25_ice, branch-exfmods-curt, release1_coupled, release1_final, release1-branch, branch-genmake2, release1, branch-nonh, tg2-branch, ecco-branch, release1_50yr, netcdf-sm0, icebear, checkpoint51n_branch
Changes since 1.8: +24 -14 lines 
Bringing comments up to data and formatting for document extraction.
1 C $Header: /u/gcmpack/models/MITgcmUV/model/src/modeldata_example.F,v 1.8 2001/02/04 14:38:48 cnh Exp $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: MODELDATA_EXAMPLE
8 C !INTERFACE:
9 SUBROUTINE MODELDATA_EXAMPLE( myThid )
10
11 C !DESCRIPTION: \bv
12 C *==========================================================*
13 C | S/R MODELDATA_EXAMPLE
14 C | o Write example data file
15 C *==========================================================*
16 C | Notes
17 C | =====
18 C | Some systems require & as the namelist terminator.
19 C | Other systems use a / character.
20 C *==========================================================*
21 C \ev
22
23 C !USES:
24 IMPLICIT NONE
25 #include "SIZE.h"
26 #include "EEPARAMS.h"
27 #include "PARAMS.h"
28
29 C !INPUT/OUTPUT PARAMETERS:
30 C == Routine arguments ==
31 INTEGER myThid
32
33 C !LOCAL VARIABLES:
34 C == Local variables ==
35 CHARACTER*(MAX_LEN_MBUF) msgBuf
36 CEOP
37
38 WRITE(msgBuf,'(A)') '// Shown below is an example "data" file.'
39 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
40 WRITE(msgBuf,'(A)') '// To use this example copy and paste the '
41 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
42 WRITE(msgBuf,'(A)') '// ">" lines. Then remove the text up to'
43 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
44 WRITE(msgBuf,'(A)') '// and including the ">".'
45 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
46 WRITE(msgBuf,'(A)') '># Example "data" file'
47 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
48 WRITE(msgBuf,'(A)') '># Lines beginning "#" are comments'
49 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
50
51 WRITE(msgBuf,'(A)') '># '
52 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
53 WRITE(msgBuf,'(A)') '># o Continuous equation parameters'
54 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
55 WRITE(msgBuf,'(A)') '># gravity - Accel due to gravity (m.s^2)'
56 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
57 WRITE(msgBuf,'(A)') '># rhonil - Reference density (kg/m^3)'
58 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
59 WRITE(msgBuf,'(A)')
60 & '># tAlpha - Thermal expansion coefficient (1/oC)'
61 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
62 WRITE(msgBuf,'(A)')
63 & '># sBeta - Haline contraction coefficient (1/ppt)'
64 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
65 WRITE(msgBuf,'(A)')
66 & '># f0 - Reference coriolis parameter ( 1/s ).'
67 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
68 WRITE(msgBuf,'(A)')
69 & '># ( South edge f on beta plane.)'
70 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
71 WRITE(msgBuf,'(A)') '># beta - df/dy ( s^-1.m^-1 ).'
72 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
73 WRITE(msgBuf,'(A)')
74 & '># viscAh - Horizontal eddy viscosity coefficient '
75 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
76 WRITE(msgBuf,'(A)') '># ( m^2/s ).'
77 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
78 WRITE(msgBuf,'(A)')
79 & '># viscAz - Vertical eddy viscosity coefficient '
80 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
81 WRITE(msgBuf,'(A)') '># ( m^2/s ).'
82 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
83 WRITE(msgBuf,'(A)')
84 & '># viscA4 - Biharmonic eddy viscosity coefficient '
85 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
86 WRITE(msgBuf,'(A)') '># ( m^4/s ).'
87 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
88 WRITE(msgBuf,'(A)')
89 & '># diffKhT - Horizontal temperature diffusivity '
90 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
91 WRITE(msgBuf,'(A)') '># ( m^2/s ).'
92 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
93 WRITE(msgBuf,'(A)')
94 & '># diffKzT - Vertical temperature diffusivity '
95 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
96 WRITE(msgBuf,'(A)') '># ( m^2/s ).'
97 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
98 WRITE(msgBuf,'(A)')
99 & '># diffK4T - Biharmonic temperature diffusivity '
100 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
101 WRITE(msgBuf,'(A)') '># ( m^4/s ).'
102 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
103 WRITE(msgBuf,'(A)') '># diffKhS - Horizontal salt diffusivity '
104 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
105 WRITE(msgBuf,'(A)') '># ( m^2/s ).'
106 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
107 WRITE(msgBuf,'(A)') '># diffKzS - Vertical salt diffusivity '
108 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
109 WRITE(msgBuf,'(A)') '># ( m^2/s ).'
110 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
111 WRITE(msgBuf,'(A)') '># diffK4S - Biharmonic salt diffusivity '
112 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
113 WRITE(msgBuf,'(A)') '># ( m^4/s ).'
114 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
115 WRITE(msgBuf,'(A)') '># momStepping - On/Off flag for momentum'
116 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
117 WRITE(msgBuf,'(A)') '># equation. '
118 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
119 WRITE(msgBuf,'(A)') '># momViscosity - On/Off flag for momentum'
120 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
121 WRITE(msgBuf,'(A)') '># mixing. '
122 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
123 WRITE(msgBuf,'(A)') '># momAdvection - On/Off flag for momentum'
124 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
125 WRITE(msgBuf,'(A)') '># self transport. '
126 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
127 WRITE(msgBuf,'(A)')
128 & '># momPressureForcing - On/Off flag for momentum'
129 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
130 WRITE(msgBuf,'(A)')
131 & '># pressure terms. '
132 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
133 WRITE(msgBuf,'(A)') '># useCoriolis - On/Off flag for momentum'
134 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
135 WRITE(msgBuf,'(A)') '># equation coriolis term. '
136 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
137 WRITE(msgBuf,'(A)')
138 & '># tempStepping - On/Off flag for temperature'
139 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
140 WRITE(msgBuf,'(A)') '># eqaution. '
141 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
142 WRITE(msgBuf,'(A)')
143 & '># tempDiffusion- On/Off flag for temperature'
144 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
145 WRITE(msgBuf,'(A)') '># mixing. '
146 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
147 WRITE(msgBuf,'(A)')
148 & '># tempAdvection- On/Off flag for temperature'
149 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
150 WRITE(msgBuf,'(A)') '># transport. '
151 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
152 WRITE(msgBuf,'(A)')
153 & '># tempForcing - On/Off flag for temperature'
154 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
155 WRITE(msgBuf,'(A)') '># forcing.'
156 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
157 WRITE(msgBuf,'(A)') '># saltDiffusion- On/Off flag for salt'
158 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
159 WRITE(msgBuf,'(A)') '># mixing. '
160 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
161 WRITE(msgBuf,'(A)') '># saltAdvection- On/Off flag for salt'
162 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
163 WRITE(msgBuf,'(A)') '># transport. '
164 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
165 WRITE(msgBuf,'(A)') '># saltForcing - On/Off flag for salt'
166 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
167 WRITE(msgBuf,'(A)') '># forcing.'
168 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
169 WRITE(msgBuf,'(A)')
170 & '># tRef - Reference vertical pot. temp'
171 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
172 WRITE(msgBuf,'(A)')
173 & '># sRef - Reference vertical salinity'
174 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
175 WRITE(msgBuf,'(A)') '>&PARM01'
176 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
177 WRITE(msgBuf,'(A)') '> gravity=9.81,'
178 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
179 WRITE(msgBuf,'(A)') '> rhonil=999.8,'
180 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
181 WRITE(msgBuf,'(A)') '> tAlpha=2.e-4,'
182 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
183 WRITE(msgBuf,'(A)') '> sBeta=7e-4'
184 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
185 WRITE(msgBuf,'(A)') '> f0=1.e-4'
186 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
187 WRITE(msgBuf,'(A)') '> viscAh=1.e3'
188 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
189 WRITE(msgBuf,'(A)') '> viscAz=1.e-5'
190 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
191 WRITE(msgBuf,'(A)') '> viscA4=0.'
192 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
193 WRITE(msgBuf,'(A)') '> diffKhT=1.e3'
194 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
195 WRITE(msgBuf,'(A)') '> diffKzT=1.e-5'
196 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
197 WRITE(msgBuf,'(A)') '> diffK4T=0.'
198 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
199 WRITE(msgBuf,'(A)') '> diffKhS=1.e3'
200 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
201 WRITE(msgBuf,'(A)') '> diffKzS=1.e-5'
202 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
203 WRITE(msgBuf,'(A)') '> diffK4S=0.'
204 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
205 WRITE(msgBuf,'(A)') '> momStepping=.TRUE.'
206 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
207 WRITE(msgBuf,'(A)') '> momViscosity=.TRUE.'
208 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
209 WRITE(msgBuf,'(A)') '> momAdvection=.TRUE.'
210 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
211 WRITE(msgBuf,'(A)') '> momPressureForcing=.TRUE.'
212 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
213 WRITE(msgBuf,'(A)') '> momForcing=.TRUE.'
214 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
215 WRITE(msgBuf,'(A)') '> useCoriolis=.TRUE.'
216 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
217 WRITE(msgBuf,'(A)') '> tempStepping=.TRUE.'
218 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
219 WRITE(msgBuf,'(A)') '> tempDiffusion=.TRUE.'
220 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
221 WRITE(msgBuf,'(A)') '> tempAdvection=.TRUE.'
222 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
223 WRITE(msgBuf,'(A)') '> tempForcing=.TRUE.'
224 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
225 WRITE(msgBuf,'(A)') '> saltDiffusion=.TRUE.'
226 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
227 WRITE(msgBuf,'(A)') '> saltAdvection=.TRUE.'
228 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
229 WRITE(msgBuf,'(A)') '> saltForcing=.TRUE.'
230 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
231 WRITE(msgBuf,'(A)') '> tRef=20.,'
232 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
233 WRITE(msgBuf,'(A)') '> sRef=35.,'
234 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
235 WRITE(msgBuf,'(A)') '> implicitFreeSurface=.TRUE.,'
236 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
237 WRITE(msgBuf,'(A)') '> rigidLid=.FALSE.,'
238 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
239 WRITE(msgBuf,'(A)') '> GMmaxSlope=1.d-2,'
240 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
241 WRITE(msgBuf,'(A)') '> GMlength=200.d3,'
242 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
243 WRITE(msgBuf,'(A)') '> GMalpha=200.d3,'
244 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
245 WRITE(msgBuf,'(A)') '> GMdepth=1000.,'
246 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
247 WRITE(msgBuf,'(A)') '> GMkBackground=0.,'
248 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
249 WRITE(msgBuf,'(A)') '>/ '
250 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
251
252 WRITE(msgBuf,'(A)') '># '
253 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
254 WRITE(msgBuf,'(A)') '># o Elliptic solver parameters'
255 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
256 WRITE(msgBuf,'(A)') '># cg2dMaxIters - Maximum number of 2d '
257 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
258 WRITE(msgBuf,'(A)') '># solver iterations. '
259 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
260 WRITE(msgBuf,'(A)')
261 & '># cg2dChkReqFreq - Frequency solver tests '
262 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
263 WRITE(msgBuf,'(A)') '># convergence. '
264 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
265 WRITE(msgBuf,'(A)') '># cg2dTargetResidual - Solver target'
266 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
267 WRITE(msgBuf,'(A)') '># residual. '
268 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
269 WRITE(msgBuf,'(A)') '>&PARM02'
270 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
271 WRITE(msgBuf,'(A)') '> cg2dMaxIters=200,'
272 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
273 WRITE(msgBuf,'(A)') '> cg2dChkResFreq=5,'
274 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
275 WRITE(msgBuf,'(A)') '> cg2dTargetResidual=1.e-7,'
276 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
277 WRITE(msgBuf,'(A)') '>/ '
278 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
279
280 WRITE(msgBuf,'(A)') '># '
281 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
282 WRITE(msgBuf,'(A)') '># o Timestepping parameters'
283 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
284 WRITE(msgBuf,'(A)') '># nIter0 - Start timestep index'
285 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
286 WRITE(msgBuf,'(A)')
287 & '># nTimeSteps - Number of timesteps in run.'
288 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
289 WRITE(msgBuf,'(A)') '># delT - Timestep ( s ).'
290 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
291 WRITE(msgBuf,'(A)') '># deltaTtracer - Tracer timestep ( s ).'
292 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
293 WRITE(msgBuf,'(A)') '># abEps - Adams-Bashforth stabilising '
294 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
295 WRITE(msgBuf,'(A)') '># factor. '
296 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
297 WRITE(msgBuf,'(A)')
298 & '># tauCD - CD scheme coupling timescale (s)'
299 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
300 WRITE(msgBuf,'(A)')
301 & '># startTime - Integration starting time (s)'
302 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
303 WRITE(msgBuf,'(A)') '># endTime - Integration ending time (s)'
304 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
305 WRITE(msgBuf,'(A)')
306 & '># cAdjFreq - Convective adjustment period (s)'
307 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
308 WRITE(msgBuf,'(A)') '># chkPtFreq - Frequency at which check '
309 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
310 WRITE(msgBuf,'(A)') '># pointing is done ( s ). '
311 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
312 WRITE(msgBuf,'(A)') '># dumpFreq - Frequency at which model '
313 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
314 WRITE(msgBuf,'(A)') '># state is stored ( s ). '
315 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
316 WRITE(msgBuf,'(A)') '>&PARM03'
317 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
318 WRITE(msgBuf,'(A)') '> nIter0=0'
319 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
320 WRITE(msgBuf,'(A)') '> nTimeSteps=5000'
321 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
322 WRITE(msgBuf,'(A)') '> delT=3600.'
323 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
324 WRITE(msgBuf,'(A)') '> deltaTtracer=3600.'
325 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
326 WRITE(msgBuf,'(A)') '> abEps=0.1'
327 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
328 WRITE(msgBuf,'(A)') '> tauCD=345600.'
329 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
330 WRITE(msgBuf,'(A)') '> startTime=0.,'
331 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
332 WRITE(msgBuf,'(A)') '> endTime=31104000.,'
333 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
334 WRITE(msgBuf,'(A)') '> chkPtFreq=864000.,'
335 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
336 WRITE(msgBuf,'(A)') '> dumpFreq=2592000.,'
337 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
338 WRITE(msgBuf,'(A)') '> cAdjFreq=86400.,'
339 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
340 WRITE(msgBuf,'(A)') '>/ '
341 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
342
343 WRITE(msgBuf,'(A)') '># '
344 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
345 WRITE(msgBuf,'(A)') '># o Gridding parameters'
346 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
347 WRITE(msgBuf,'(A)') '># l - Global domain grid-points in X'
348 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
349 WRITE(msgBuf,'(A)') '># m - Global domain grid-points in Y'
350 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
351 WRITE(msgBuf,'(A)') '># n - Grid-points in Z'
352 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
353 WRITE(msgBuf,'(A)')
354 & '># usingSphericalPolarGrid - On/Off flag for'
355 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
356 WRITE(msgBuf,'(A)')
357 & '># selecting spherical polar coordinates'
358 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
359 WRITE(msgBuf,'(A)') '># usingCartesianGrid - On/Off flag for'
360 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
361 WRITE(msgBuf,'(A)') '># selecting cartesian coordinates'
362 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
363 WRITE(msgBuf,'(A)') '># delX - Zonal grid spacing. Degrees'
364 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
365 WRITE(msgBuf,'(A)') '># for spherical polar and m for'
366 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
367 WRITE(msgBuf,'(A)')
368 & '># cartesian. A value for each point'
369 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
370 WRITE(msgBuf,'(A)') '># in X can be specified.'
371 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
372 WRITE(msgBuf,'(A)') '># delY - Meridional grid spacing. Degrees'
373 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
374 WRITE(msgBuf,'(A)') '># for spherical polar and m for'
375 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
376 WRITE(msgBuf,'(A)')
377 & '># cartesian. A value for each point'
378 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
379 WRITE(msgBuf,'(A)') '># in Y can be specified.'
380 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
381 WRITE(msgBuf,'(A)') '># delZ - Vertical grid spacing (m).'
382 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
383 WRITE(msgBuf,'(A)') '># delP - Vertical grid spacing (Pa).'
384 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
385 WRITE(msgBuf,'(A)') '># phiMin - Southern boundary latitude'
386 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
387 WRITE(msgBuf,'(A)') '># (spherical polar grid). '
388 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
389 WRITE(msgBuf,'(A)') '># rSphere- Radius of globe '
390 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
391 WRITE(msgBuf,'(A)') '># (spherical polar grid). '
392 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
393 WRITE(msgBuf,'(A)') '>&PARM04'
394 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
395 WRITE(msgBuf,'(A)') '> n=20,'
396 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
397 WRITE(msgBuf,'(A)') '> l=122,'
398 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
399 WRITE(msgBuf,'(A)') '> m=86,'
400 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
401 WRITE(msgBuf,'(A)') '> usingSphericalPolarGrid=.TRUE.,'
402 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
403 WRITE(msgBuf,'(A)') '> usingCartesianGrid=.FALSE.,'
404 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
405 WRITE(msgBuf,'(A)') '> delx=1.,'
406 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
407 WRITE(msgBuf,'(A)') '> dely=1.,'
408 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
409 WRITE(msgBuf,'(A)') '> delz= 100., 100., 100., 100., 100.,'
410 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
411 WRITE(msgBuf,'(A)') '> 100., 100., 100., 100., 100.,'
412 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
413 WRITE(msgBuf,'(A)') '> 100., 100., 100., 100., 100.,'
414 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
415 WRITE(msgBuf,'(A)') '> 100., 100., 100., 100., 100.,'
416 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
417 WRITE(msgBuf,'(A)') '> phiMin=-80.,'
418 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
419 WRITE(msgBuf,'(A)') '> rSphere=6430.E3'
420 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
421 WRITE(msgBuf,'(A)') '>/ '
422 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
423
424 WRITE(msgBuf,'(A)') '># Note: Some systems use & as the '
425 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
426 WRITE(msgBuf,'(A)') '># namelist terminator. Other systems'
427 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
428 WRITE(msgBuf,'(A)') '># use a / character (as shown here).'
429 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
430 WRITE(msgBuf,'(A)') ' '
431 CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
432 C
433 RETURN
434 END
435
  ViewVC Help
Powered by ViewVC 1.1.22