38 |
& A_glen_isothermal, n_glen, eps_glen_min, |
& A_glen_isothermal, n_glen, eps_glen_min, |
39 |
& C_basal_fric_const, n_basal_friction, |
& C_basal_fric_const, n_basal_friction, |
40 |
& streamice_vel_update,streamice_cg_tol,streamice_nonlin_tol, |
& streamice_vel_update,streamice_cg_tol,streamice_nonlin_tol, |
41 |
|
& streamice_nonlin_tol_fp, |
42 |
& streamice_max_cg_iter, streamice_max_nl_iter, |
& streamice_max_cg_iter, streamice_max_nl_iter, |
43 |
& STREAMICE_GL_regularize, |
& STREAMICE_GL_regularize, |
44 |
& STREAMICEthickInit, |
& STREAMICEthickInit, |
45 |
|
& STREAMICEsigcoordInit, |
46 |
|
& STREAMICEsigcoordFile, |
47 |
& STREAMICEthickFile, |
& STREAMICEthickFile, |
48 |
& STREAMICEcalveMaskFile, |
& STREAMICEcalveMaskFile, |
49 |
& STREAMICEison, |
& STREAMICEison, |
51 |
& STREAMICE_dump_mnc, STREAMICE_tave_mnc, |
& STREAMICE_dump_mnc, STREAMICE_tave_mnc, |
52 |
& STREAMICE_GL_regularize, STREAMICE_move_front, |
& STREAMICE_GL_regularize, STREAMICE_move_front, |
53 |
& STREAMICE_calve_to_mask, |
& STREAMICE_calve_to_mask, |
54 |
|
& STREAMICE_diagnostic_only, |
55 |
! & STREAMICE_construct_matrix, |
! & STREAMICE_construct_matrix, |
56 |
& STREAMICE_lower_cg_tol, |
& STREAMICE_lower_cg_tol, |
57 |
& streamice_CFL_factor |
& streamice_CFL_factor, |
58 |
|
& streamice_adjDump, |
59 |
|
! & STREAMICE_hybrid_stress, |
60 |
|
& streamice_bg_surf_slope_x, streamice_bg_surf_slope_y, |
61 |
|
& streamice_kx_b_init, streamice_ky_b_init, |
62 |
|
& STREAMICEbasalTracConfig, |
63 |
|
& STREAMICEbasalTracFile, |
64 |
|
& STREAMICEvelOptimFile, |
65 |
|
& STREAMICE_ppm_driving_stress, |
66 |
|
& STREAMICE_h_ctrl_const_surf, |
67 |
|
& streamice_wgt_drift,streamice_wgt_surf,streamice_wgt_vel, |
68 |
|
& streamice_wgt_avthick |
69 |
|
|
70 |
|
|
71 |
NAMELIST /STREAMICE_PARM02/ |
NAMELIST /STREAMICE_PARM02/ |
72 |
& shelf_max_draft, |
& shelf_max_draft, |
98 |
& min_y_CFBC_WEST, max_y_CFBC_WEST, |
& min_y_CFBC_WEST, max_y_CFBC_WEST, |
99 |
& min_y_CFBC_EAST, max_y_CFBC_EAST, |
& min_y_CFBC_EAST, max_y_CFBC_EAST, |
100 |
& flux_bdry_val_SOUTH, flux_bdry_val_NORTH, |
& flux_bdry_val_SOUTH, flux_bdry_val_NORTH, |
101 |
& flux_bdry_val_WEST, flux_bdry_val_EAST |
& flux_bdry_val_WEST, flux_bdry_val_EAST, |
102 |
|
& STREAMICE_NS_periodic, STREAMICE_EW_periodic |
103 |
|
|
104 |
_BEGIN_MASTER(myThid) |
_BEGIN_MASTER(myThid) |
105 |
|
|
115 |
streamice_vel_update = 169200. ! seconds |
streamice_vel_update = 169200. ! seconds |
116 |
streamice_cg_tol = 1e-6 |
streamice_cg_tol = 1e-6 |
117 |
streamice_nonlin_tol = 1e-6 |
streamice_nonlin_tol = 1e-6 |
118 |
|
streamice_nonlin_tol_fp = 1.e-14 |
119 |
streamice_max_cg_iter = 2000 |
streamice_max_cg_iter = 2000 |
120 |
streamice_max_nl_iter = 100 |
streamice_max_nl_iter = 100 |
121 |
streamice_n_sub_regularize = 4 |
streamice_n_sub_regularize = 4 |
122 |
streamice_CFL_factor = .5 |
streamice_CFL_factor = .5 |
123 |
|
streamice_adjDump = 0. |
124 |
|
streamice_bg_surf_slope_x = .0 |
125 |
|
streamice_bg_surf_slope_y = 0. |
126 |
|
streamice_kx_b_init = 1. |
127 |
|
streamice_ky_b_init = 1. |
128 |
|
streamice_wgt_drift = 0. |
129 |
|
streamice_wgt_surf = 0. |
130 |
|
streamice_wgt_vel = 0. |
131 |
|
streamice_wgt_avthick = 0. |
132 |
|
|
133 |
STREAMICEthickInit = 'FILE' |
STREAMICEthickInit = 'FILE' |
134 |
STREAMICEthickFile = ' ' |
STREAMICEthickFile = ' ' |
135 |
STREAMICEcalveMaskFile = ' ' |
STREAMICEcalveMaskFile = ' ' |
136 |
|
STREAMICEsigcoordInit = 'UNIFORM' |
137 |
|
STREAMICEsigcoordFile = ' ' |
138 |
|
STREAMICEbasalTracConfig = 'UNIFORM' |
139 |
|
STREAMICEbasalTracFile = ' ' |
140 |
|
STREAMICEvelOptimFile = '' |
141 |
|
|
142 |
STREAMICEison = .TRUE. |
STREAMICEison = .TRUE. |
143 |
STREAMICE_tave_mdsio = .TRUE. |
STREAMICE_tave_mdsio = .TRUE. |
149 |
STREAMICE_calve_to_mask = .FALSE. |
STREAMICE_calve_to_mask = .FALSE. |
150 |
! STREAMICE_construct_matrix = .TRUE. |
! STREAMICE_construct_matrix = .TRUE. |
151 |
STREAMICE_lower_cg_tol = .FALSE. |
STREAMICE_lower_cg_tol = .FALSE. |
152 |
|
STREAMICE_diagnostic_only = .FALSE. |
153 |
|
STREAMICE_ppm_driving_stress = .FALSE. |
154 |
|
STREAMICE_h_ctrl_const_surf = .FALSE. |
155 |
|
! STREAMICE_hybrid_stress= .FALSE. |
156 |
|
|
157 |
min_x_noflow_NORTH = 0. |
min_x_noflow_NORTH = 0. |
158 |
max_x_noflow_NORTH = 0. |
max_x_noflow_NORTH = 0. |
195 |
flux_bdry_val_WEST = 0. |
flux_bdry_val_WEST = 0. |
196 |
flux_bdry_val_EAST = 0. |
flux_bdry_val_EAST = 0. |
197 |
|
|
198 |
|
STREAMICE_NS_periodic = .FALSE. |
199 |
|
STREAMICE_EW_periodic = .FALSE. |
200 |
|
|
201 |
WRITE(msgBuf,'(A)') 'STREAMICE_READPARMS: opening data.streamice' |
WRITE(msgBuf,'(A)') 'STREAMICE_READPARMS: opening data.streamice' |
202 |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |
212 |
& 'STREAMICE_READPARMS: read first param block' |
& 'STREAMICE_READPARMS: read first param block' |
213 |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |
214 |
& SQUEEZE_RIGHT , 1) |
& SQUEEZE_RIGHT , 1) |
|
READ(UNIT=iUnit,NML=STREAMICE_PARM02) |
|
215 |
|
|
216 |
IF (TRIM(STREAMICEthickInit) .eq. "PARAM") THEN |
IF (TRIM(STREAMICEthickInit) .eq. "PARAM") THEN |
217 |
|
READ(UNIT=iUnit,NML=STREAMICE_PARM02) |
218 |
WRITE(msgBuf,'(A)') |
WRITE(msgBuf,'(A)') |
219 |
& 'STREAMICE_READPARMS: read second param block' |
& 'STREAMICE_READPARMS: read second param block' |
220 |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |
221 |
& SQUEEZE_RIGHT , 1) |
& SQUEEZE_RIGHT , 1) |
|
READ(UNIT=iUnit,NML=STREAMICE_PARM03) |
|
222 |
ENDIF |
ENDIF |
223 |
|
|
224 |
|
READ(UNIT=iUnit,NML=STREAMICE_PARM03) |
225 |
WRITE(msgBuf,'(A)') |
WRITE(msgBuf,'(A)') |
226 |
& 'STREAMICE_READPARMS: read third param block' |
& 'STREAMICE_READPARMS: read third param block' |
227 |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |