35 |
|
|
36 |
NAMELIST /STREAMICE_PARM01/ |
NAMELIST /STREAMICE_PARM01/ |
37 |
& streamice_density, streamice_density_ocean_avg, |
& streamice_density, streamice_density_ocean_avg, |
38 |
& A_glen_isothermal, n_glen, eps_glen_min, |
& A_glen_isothermal, n_glen, eps_glen_min, eps_u_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_geom_file_setup, |
55 |
|
& STREAMICE_diagnostic_only, |
56 |
! & STREAMICE_construct_matrix, |
! & STREAMICE_construct_matrix, |
57 |
& STREAMICE_lower_cg_tol, |
& STREAMICE_lower_cg_tol, |
58 |
& streamice_CFL_factor, |
& streamice_CFL_factor, |
59 |
& streamice_adjDump |
& streamice_adjDump, |
60 |
|
! & STREAMICE_hybrid_stress, |
61 |
|
& streamice_bg_surf_slope_x, streamice_bg_surf_slope_y, |
62 |
|
& streamice_kx_b_init, streamice_ky_b_init, |
63 |
|
& STREAMICEbasalTracConfig, |
64 |
|
& STREAMICEbasalTracFile, |
65 |
|
& STREAMICEvelOptimFile, |
66 |
|
& STREAMICEhmaskFile, |
67 |
|
& STREAMICEuFaceBdryFile, |
68 |
|
& STREAMICEvFaceBdryFile, |
69 |
|
& STREAMICEuDirichValsFile, |
70 |
|
& STREAMICEvDirichValsFile, |
71 |
|
& STREAMICEGlenConstFile, STREAMICEGlenConstConfig, |
72 |
|
& STREAMICE_ppm_driving_stress, |
73 |
|
& STREAMICE_h_ctrl_const_surf, |
74 |
|
& streamice_wgt_drift,streamice_wgt_surf,streamice_wgt_vel, |
75 |
|
& streamice_wgt_avthick, streamice_wgt_tikh, |
76 |
|
& streamice_addl_backstress |
77 |
|
|
78 |
|
|
79 |
NAMELIST /STREAMICE_PARM02/ |
NAMELIST /STREAMICE_PARM02/ |
80 |
& shelf_max_draft, |
& shelf_max_draft, |
106 |
& min_y_CFBC_WEST, max_y_CFBC_WEST, |
& min_y_CFBC_WEST, max_y_CFBC_WEST, |
107 |
& min_y_CFBC_EAST, max_y_CFBC_EAST, |
& min_y_CFBC_EAST, max_y_CFBC_EAST, |
108 |
& flux_bdry_val_SOUTH, flux_bdry_val_NORTH, |
& flux_bdry_val_SOUTH, flux_bdry_val_NORTH, |
109 |
& flux_bdry_val_WEST, flux_bdry_val_EAST |
& flux_bdry_val_WEST, flux_bdry_val_EAST, |
110 |
|
& STREAMICE_NS_periodic, STREAMICE_EW_periodic |
111 |
|
|
112 |
_BEGIN_MASTER(myThid) |
_BEGIN_MASTER(myThid) |
113 |
|
|
118 |
A_glen_isothermal = 9.461e-18 ! Pa (-1/3) a |
A_glen_isothermal = 9.461e-18 ! Pa (-1/3) a |
119 |
n_glen = 3. |
n_glen = 3. |
120 |
eps_glen_min = 1.0e-12 |
eps_glen_min = 1.0e-12 |
121 |
|
eps_u_min = 1.0e-6 |
122 |
C_basal_fric_const = 31.71 ! Pa (m/a)-1n |
C_basal_fric_const = 31.71 ! Pa (m/a)-1n |
123 |
n_basal_friction = 1. |
n_basal_friction = 1. |
124 |
streamice_vel_update = 169200. ! seconds |
streamice_vel_update = 169200. ! seconds |
125 |
streamice_cg_tol = 1e-6 |
streamice_cg_tol = 1e-6 |
126 |
streamice_nonlin_tol = 1e-6 |
streamice_nonlin_tol = 1e-6 |
127 |
|
streamice_nonlin_tol_fp = 1.e-14 |
128 |
streamice_max_cg_iter = 2000 |
streamice_max_cg_iter = 2000 |
129 |
streamice_max_nl_iter = 100 |
streamice_max_nl_iter = 100 |
130 |
streamice_n_sub_regularize = 4 |
streamice_n_sub_regularize = 4 |
131 |
streamice_CFL_factor = .5 |
streamice_CFL_factor = .5 |
132 |
streamice_adjDump = 0. |
streamice_adjDump = 0. |
133 |
|
streamice_bg_surf_slope_x = .0 |
134 |
|
streamice_bg_surf_slope_y = 0. |
135 |
|
streamice_kx_b_init = 1. |
136 |
|
streamice_ky_b_init = 1. |
137 |
|
streamice_wgt_drift = 0. |
138 |
|
streamice_wgt_tikh = 0. |
139 |
|
streamice_wgt_surf = 0. |
140 |
|
streamice_wgt_vel = 0. |
141 |
|
streamice_wgt_avthick = 0. |
142 |
|
streamice_addl_backstress = 0. |
143 |
|
|
144 |
STREAMICEthickInit = 'FILE' |
STREAMICEthickInit = 'FILE' |
145 |
STREAMICEthickFile = ' ' |
STREAMICEthickFile = ' ' |
146 |
STREAMICEcalveMaskFile = ' ' |
STREAMICEcalveMaskFile = ' ' |
147 |
|
STREAMICEsigcoordInit = 'UNIFORM' |
148 |
|
STREAMICEsigcoordFile = ' ' |
149 |
|
STREAMICEbasalTracConfig = 'UNIFORM' |
150 |
|
STREAMICEbasalTracFile = ' ' |
151 |
|
STREAMICEvelOptimFile = '' |
152 |
|
STREAMICEhmaskFile = '' |
153 |
|
STREAMICEuFaceBdryFile = '' |
154 |
|
STREAMICEvFaceBdryFile = '' |
155 |
|
STREAMICEuDirichValsFile = '' |
156 |
|
STREAMICEvDirichValsFile = '' |
157 |
|
STREAMICEGlenConstFile = '' |
158 |
|
STREAMICEGlenConstConfig = '' |
159 |
|
|
160 |
STREAMICEison = .TRUE. |
STREAMICEison = .TRUE. |
161 |
STREAMICE_tave_mdsio = .TRUE. |
STREAMICE_tave_mdsio = .TRUE. |
165 |
STREAMICE_GL_regularize = .FALSE. |
STREAMICE_GL_regularize = .FALSE. |
166 |
STREAMICE_move_front = .FALSE. |
STREAMICE_move_front = .FALSE. |
167 |
STREAMICE_calve_to_mask = .FALSE. |
STREAMICE_calve_to_mask = .FALSE. |
168 |
|
! STREAMICE_geom_file_setup = .FALSE. |
169 |
! STREAMICE_construct_matrix = .TRUE. |
! STREAMICE_construct_matrix = .TRUE. |
170 |
STREAMICE_lower_cg_tol = .FALSE. |
STREAMICE_lower_cg_tol = .FALSE. |
171 |
|
STREAMICE_diagnostic_only = .FALSE. |
172 |
|
STREAMICE_ppm_driving_stress = .FALSE. |
173 |
|
STREAMICE_h_ctrl_const_surf = .FALSE. |
174 |
|
! STREAMICE_hybrid_stress= .FALSE. |
175 |
|
|
176 |
min_x_noflow_NORTH = 0. |
min_x_noflow_NORTH = 0. |
177 |
max_x_noflow_NORTH = 0. |
max_x_noflow_NORTH = 0. |
214 |
flux_bdry_val_WEST = 0. |
flux_bdry_val_WEST = 0. |
215 |
flux_bdry_val_EAST = 0. |
flux_bdry_val_EAST = 0. |
216 |
|
|
217 |
|
STREAMICE_NS_periodic = .FALSE. |
218 |
|
STREAMICE_EW_periodic = .FALSE. |
219 |
|
|
220 |
WRITE(msgBuf,'(A)') 'STREAMICE_READPARMS: opening data.streamice' |
WRITE(msgBuf,'(A)') 'STREAMICE_READPARMS: opening data.streamice' |
221 |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |