/[MITgcm]/MITgcm_contrib/dgoldberg/streamice/streamice_readparms.F
ViewVC logotype

Contents of /MITgcm_contrib/dgoldberg/streamice/streamice_readparms.F

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


Revision 1.15 - (show annotations) (download)
Sat Jun 8 22:15:34 2013 UTC (12 years, 1 month ago) by dgoldberg
Branch: MAIN
Changes since 1.14: +83 -9 lines
new advected scalar; new advection scheme for thickness update; corresponding TAF directives

1 C $Header: /u/gcmpack/MITgcm_contrib/dgoldberg/streamice/streamice_readparms.F,v 1.13 2013/04/06 17:43:41 dgoldberg Exp $
2 C $Name: $
3
4
5 C this needs changes
6
7 #include "STREAMICE_OPTIONS.h"
8
9 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
10 CBOP 0
11 SUBROUTINE STREAMICE_READPARMS( myThid )
12
13 C !DESCRIPTION:
14 C Initialize STREAMICE variables and constants.
15
16 C !USES:
17 IMPLICIT NONE
18 #include "SIZE.h"
19 #include "EEPARAMS.h"
20 #include "PARAMS.h"
21 #include "STREAMICE.h"
22 #include "STREAMICE_BDRY.h"
23 #ifdef ALLOW_STREAMICE_FLUX_CONTROL
24 #include "STREAMICE_CTRL_FLUX.h"
25 #endif
26
27 C !INPUT PARAMETERS:
28 INTEGER myThid
29 CEOP
30
31 #ifdef ALLOW_STREAMICE
32
33 C !LOCAL VARIABLES:
34 C msgBuf :: Informational/error message buffer
35 C iUnit :: Work variable for IO unit number
36 CHARACTER*(MAX_LEN_MBUF) msgBuf
37 INTEGER iUnit, iarr, tarr
38
39 NAMELIST /STREAMICE_PARM01/
40 & streamice_density, streamice_density_ocean_avg,
41 & B_glen_isothermal, n_glen, eps_glen_min, eps_u_min,
42 & C_basal_fric_const, n_basal_friction,
43 & streamice_vel_update,streamice_cg_tol,streamice_nonlin_tol,
44 & streamice_nonlin_tol_fp,
45 & streamice_max_cg_iter, streamice_max_nl_iter,
46 & STREAMICE_GL_regularize,
47 & STREAMICEthickInit,
48 & STREAMICEsigcoordInit,
49 & STREAMICEsigcoordFile,
50 & STREAMICEthickFile,
51 & STREAMICEcalveMaskFile,
52 & STREAMICEcostMaskFile,
53 & STREAMICEison,
54 & STREAMICE_dump_mdsio, STREAMICE_tave_mdsio,
55 & STREAMICE_dump_mnc, STREAMICE_tave_mnc,
56 & STREAMICE_GL_regularize, STREAMICE_move_front,
57 & STREAMICE_calve_to_mask,
58 & STREAMICE_diagnostic_only,
59 & STREAMICE_lower_cg_tol,
60 & streamice_CFL_factor,
61 & streamice_adjDump,
62 & streamice_bg_surf_slope_x, streamice_bg_surf_slope_y,
63 & streamice_kx_b_init, streamice_ky_b_init,
64 & STREAMICEbasalTracConfig,
65 & STREAMICEbasalTracFile,
66 & STREAMICEvelOptimFile,
67 & STREAMICEtopogFile,
68 & STREAMICEhmaskFile,
69 & STREAMICEHBCyFile,
70 & STREAMICEHBCxFile,
71 & STREAMICEuFaceBdryFile,
72 & STREAMICEvFaceBdryFile,
73 & STREAMICEuDirichValsFile,
74 & STREAMICEvDirichValsFile,
75 & STREAMICEGlenConstFile, STREAMICEGlenConstConfig,
76 & STREAMICE_ppm_driving_stress,
77 & STREAMICE_h_ctrl_const_surf,
78 & streamice_wgt_drift,streamice_wgt_surf,streamice_wgt_vel,
79 & streamice_wgt_avthick, streamice_wgt_tikh,
80 & streamice_addl_backstress,
81 & streamice_smooth_gl_width,
82 & STREAMICE_ADV_SCHEME
83
84 #ifdef ALLOW_STREAMICE_2DTRACER
85 NAMELIST /STREAMICE_PARMTRACER/
86 & STREAMICETrac2DBCxFile,
87 & STREAMICETrac2DBCyFile,
88 & STREAMICETrac2DINITFile
89 #endif
90
91 #ifdef ALLOW_PETSC
92 NAMELIST /STREAMICE_PARMPETSC/
93 & PETSC_PRECOND_TYPE, PETSC_SOLVER_TYPE
94 #endif
95
96 #ifdef ALLOW_STREAMICE_FLUX_CONTROL
97 NAMELIST /STREAMICE_PARMFLUXCTRL/
98 & n_fluxes, n_epochs,
99 & streamice_ctrl_flux_id,
100 & streamice_ctrl_flux_scaleVel
101 #endif
102
103
104 NAMELIST /STREAMICE_PARM02/
105 & shelf_max_draft,
106 & shelf_min_draft,
107 & shelf_edge_pos,
108 & shelf_slope_scale,
109 & shelf_flat_width,
110 & flow_dir
111
112 NAMELIST /STREAMICE_PARM03/
113 & min_x_noflow_NORTH, max_x_noflow_NORTH,
114 & min_x_noflow_SOUTH, max_x_noflow_SOUTH,
115 & min_y_noflow_WEST, max_y_noflow_WEST,
116 & min_y_noflow_EAST, max_y_noflow_EAST,
117 & min_x_noStress_NORTH, max_x_noStress_NORTH,
118 & min_x_noStress_SOUTH, max_x_noStress_SOUTH,
119 & min_y_noStress_WEST, max_y_noStress_WEST,
120 & min_y_noStress_EAST, max_y_noStress_EAST,
121 & min_x_FluxBdry_NORTH, max_x_FluxBdry_NORTH,
122 & min_x_FluxBdry_SOUTH, max_x_FluxBdry_SOUTH,
123 & min_y_FluxBdry_WEST, max_y_FluxBdry_WEST,
124 & min_y_FluxBdry_EAST, max_y_FluxBdry_EAST,
125 & min_x_Dirich_NORTH, max_x_Dirich_NORTH,
126 & min_x_Dirich_SOUTH, max_x_Dirich_SOUTH,
127 & min_y_Dirich_WEST, max_y_Dirich_WEST,
128 & min_y_Dirich_EAST, max_y_Dirich_EAST,
129 & min_x_CFBC_NORTH, max_x_CFBC_NORTH,
130 & min_x_CFBC_SOUTH, max_x_CFBC_SOUTH,
131 & min_y_CFBC_WEST, max_y_CFBC_WEST,
132 & min_y_CFBC_EAST, max_y_CFBC_EAST,
133 & flux_bdry_val_SOUTH, flux_bdry_val_NORTH,
134 & flux_bdry_val_WEST, flux_bdry_val_EAST,
135 & STREAMICE_NS_periodic, STREAMICE_EW_periodic
136
137 _BEGIN_MASTER(myThid)
138
139 C-- Default values for STREAMICE
140
141 streamice_density = 917.
142 streamice_density_ocean_avg = 1024.
143 B_glen_isothermal = 9.461e-18 ! Pa (-1/3) a
144 n_glen = 3.
145 eps_glen_min = 1.0e-12
146 eps_u_min = 1.0e-6
147 C_basal_fric_const = 31.71 ! Pa (m/a)-1n
148 n_basal_friction = 1.
149 streamice_vel_update = 169200. ! seconds
150 streamice_cg_tol = 1e-6
151 streamice_nonlin_tol = 1e-6
152 streamice_nonlin_tol_fp = 1.e-14
153 streamice_max_cg_iter = 2000
154 streamice_max_nl_iter = 100
155 streamice_n_sub_regularize = 4
156 streamice_CFL_factor = .5
157 streamice_adjDump = 0.
158 streamice_bg_surf_slope_x = .0
159 streamice_bg_surf_slope_y = 0.
160 streamice_kx_b_init = 1.
161 streamice_ky_b_init = 1.
162 streamice_wgt_drift = 0.
163 streamice_wgt_tikh = 0.
164 streamice_wgt_surf = 0.
165 streamice_wgt_vel = 0.
166 streamice_wgt_avthick = 0.
167 streamice_addl_backstress = 0.0
168 streamice_smooth_gl_width = 0.0
169
170 STREAMICEthickInit = 'FILE'
171 STREAMICEthickFile = ' '
172 STREAMICEcalveMaskFile = ' '
173 STREAMICEsigcoordInit = 'UNIFORM'
174 STREAMICEsigcoordFile = ' '
175 STREAMICEbasalTracConfig = 'UNIFORM'
176 STREAMICEbasalTracFile = ' '
177 STREAMICEvelOptimFile = ''
178 STREAMICEtopogFile = ''
179 STREAMICEhmaskFile = ''
180 STREAMICEHBCyFile = ''
181 STREAMICEHBCxFile = ''
182 #ifdef ALLOW_STREAMICE_2DTRACER
183 STREAMICETrac2DBCxFile = ''
184 STREAMICETrac2DBCyFile = ''
185 STREAMICETrac2DInitFile = ''
186 #endif
187 STREAMICEuFaceBdryFile = ''
188 STREAMICEvFaceBdryFile = ''
189 STREAMICEuDirichValsFile = ''
190 STREAMICEvDirichValsFile = ''
191 STREAMICEGlenConstFile = ''
192 STREAMICEcostMaskFile = ''
193 STREAMICEGlenConstConfig = 'UNIFORM'
194 #ifdef ALLOW_PETSC
195 PETSC_PRECOND_TYPE = 'PCBJACOBI'
196 PETSC_SOLVER_TYPE = 'KSPCG'
197 #endif
198 STREAMICE_ADV_SCHEME = ''
199
200 #ifdef ALLOW_STREAMICE_FLUX_CONTROL
201 n_fluxes = 0
202 n_epochs = 0
203 DO iarr=1,n_fluxes_max
204 streamice_ctrl_flux_id(iarr) = 0
205 DO tarr=1,n_epochs_max
206 streamice_ctrl_flux_scaleVel(iarr,tarr) = 0. _d 0
207 ENDDO
208 ENDDO
209 #endif
210
211 STREAMICEison = .TRUE.
212 STREAMICE_tave_mdsio = .TRUE.
213 STREAMICE_dump_mdsio = .TRUE.
214 STREAMICE_dump_mnc = .FALSE.
215 STREAMICE_tave_mnc = .FALSE.
216 STREAMICE_GL_regularize = .FALSE.
217 STREAMICE_move_front = .FALSE.
218 STREAMICE_calve_to_mask = .FALSE.
219 ! STREAMICE_geom_file_setup = .FALSE.
220 ! STREAMICE_construct_matrix = .TRUE.
221 STREAMICE_lower_cg_tol = .FALSE.
222 STREAMICE_diagnostic_only = .FALSE.
223 STREAMICE_ppm_driving_stress = .FALSE.
224 STREAMICE_h_ctrl_const_surf = .FALSE.
225 ! STREAMICE_hybrid_stress= .FALSE.
226
227 min_x_noflow_NORTH = 0.
228 max_x_noflow_NORTH = 0.
229 min_x_noflow_SOUTH = 0.
230 max_x_noflow_SOUTH = 0.
231 min_y_noflow_WEST = 0.
232 max_y_noflow_WEST = 0.
233 min_y_noflow_EAST = 0.
234 max_y_noflow_EAST = 0.
235 min_x_noStress_NORTH = 0.
236 max_x_noStress_NORTH = 0.
237 min_x_noStress_SOUTH = 0.
238 max_x_noStress_SOUTH = 0.
239 min_y_noStress_WEST = 0.
240 max_y_noStress_WEST = 0.
241 min_y_noStress_EAST = 0.
242 max_y_noStress_EAST = 0.
243 min_x_FluxBdry_NORTH = 0.
244 max_x_FluxBdry_NORTH = 0.
245 min_x_FluxBdry_SOUTH = 0.
246 max_x_FluxBdry_SOUTH = 0.
247 min_y_FluxBdry_WEST = 0.
248 max_y_FluxBdry_WEST = 0.
249 min_y_FluxBdry_EAST = 0.
250 max_y_FluxBdry_EAST = 0.
251 min_x_Dirich_NORTH = 0.
252 max_x_Dirich_NORTH = 0.
253 min_x_Dirich_SOUTH = 0.
254 max_x_Dirich_SOUTH = 0.
255 min_y_Dirich_WEST = 0.
256 max_y_Dirich_WEST = 0.
257 min_y_Dirich_EAST = 0.
258 max_y_Dirich_EAST = 0.
259 min_y_CFBC_WEST = 0.
260 max_y_CFBC_WEST = 0.
261 min_y_CFBC_EAST = 0.
262 max_y_CFBC_EAST = 0.
263 flux_bdry_val_SOUTH = 0.
264 flux_bdry_val_NORTH = 0.
265 flux_bdry_val_WEST = 0.
266 flux_bdry_val_EAST = 0.
267
268 STREAMICE_NS_periodic = .FALSE.
269 STREAMICE_EW_periodic = .FALSE.
270
271 WRITE(msgBuf,'(A)') 'STREAMICE_READPARMS: opening data.streamice'
272 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
273 & SQUEEZE_RIGHT , 1)
274 CALL OPEN_COPY_DATA_FILE(
275 I 'data.streamice', 'STREAMICE_READPARMS',
276 O iUnit,
277 I myThid )
278
279 C Read parameters from open data file
280 READ(UNIT=iUnit,NML=STREAMICE_PARM01)
281 WRITE(msgBuf,'(A)')
282 & 'STREAMICE_READPARMS: read first param block'
283 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
284 & SQUEEZE_RIGHT , 1)
285
286 IF (TRIM(STREAMICEthickInit) .eq. "PARAM") THEN
287 READ(UNIT=iUnit,NML=STREAMICE_PARM02)
288 WRITE(msgBuf,'(A)')
289 & 'STREAMICE_READPARMS: read second param block'
290 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
291 & SQUEEZE_RIGHT , 1)
292 ENDIF
293
294 #ifdef ALLOW_STREAMICE_2DTRACER
295 READ(UNIT=iUnit,NML=STREAMICE_PARMTRACER)
296 WRITE(msgBuf,'(A)')
297 & 'STREAMICE_READPARMS: read tracer param block'
298 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
299 & SQUEEZE_RIGHT , 1)
300 #endif
301
302 #ifdef ALLOW_PETSC
303 READ(UNIT=iUnit,NML=STREAMICE_PARMPETSC)
304 WRITE(msgBuf,'(A)')
305 & 'STREAMICE_READPARMS: read petsc param block'
306 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
307 & SQUEEZE_RIGHT , 1)
308 #endif
309
310
311 READ(UNIT=iUnit,NML=STREAMICE_PARM03)
312 WRITE(msgBuf,'(A)')
313 & 'STREAMICE_READPARMS: read third param block'
314 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
315 & SQUEEZE_RIGHT , 1)
316 CLOSE(iUnit)
317
318
319 #ifdef ALLOW_STREAMICE_FLUX_CONTROL
320
321 CALL OPEN_COPY_DATA_FILE(
322 I 'data.strmctrlflux', 'STREAMICE_READPARMS',
323 O iUnit,
324 I myThid )
325
326
327 READ(UNIT=iUnit,NML=STREAMICE_PARMFLUXCTRL)
328 WRITE(msgBuf,'(A)')
329 & 'STREAMICE_READPARMS: read flux_ctrl param block'
330 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
331 & SQUEEZE_RIGHT , 1)
332 CLOSE(iUnit)
333 #endif
334
335
336
337 streamice_nstep_velocity = NINT (streamice_vel_update / deltaT)
338
339 C- Set Output type flags :
340
341 #ifdef ALLOW_MNC
342 IF (useMNC) THEN
343 IF ( .NOT.outputTypesInclusive
344 & .AND. STREAMICE_tave_mnc ) STREAMICE_tave_mdsio = .FALSE.
345 IF ( .NOT.outputTypesInclusive
346 & .AND. STREAMICE_dump_mnc ) STREAMICE_dump_mdsio = .FALSE.
347 ENDIF
348 #endif
349
350 _END_MASTER(myThid)
351
352 C-- Everyone else must wait for the parameters to be loaded
353 _BARRIER
354
355 #endif /* ALLOW_STREAMICE */
356
357 RETURN
358 END

  ViewVC Help
Powered by ViewVC 1.1.22