6 |
|
|
7 |
COMMON /STREAMICE_PARMS_R/ |
COMMON /STREAMICE_PARMS_R/ |
8 |
& streamice_density, streamice_density_ocean_avg, |
& streamice_density, streamice_density_ocean_avg, |
9 |
& A_glen_isothermal, n_glen, eps_glen_min, |
& A_glen_isothermal, n_glen, eps_glen_min, eps_u_min, |
10 |
& C_basal_fric_const, n_basal_friction, streamice_input_flux_unif, |
& C_basal_fric_const, n_basal_friction, streamice_input_flux_unif, |
11 |
& streamice_vel_update, streamice_cg_tol, streamice_nonlin_tol, |
& streamice_vel_update, streamice_cg_tol, streamice_nonlin_tol, |
12 |
& streamice_nonlin_tol_fp, |
& streamice_nonlin_tol_fp, |
13 |
& streamice_CFL_factor, streamice_adjDump, |
& streamice_CFL_factor, streamice_adjDump, |
14 |
& streamice_bg_surf_slope_x, streamice_bg_surf_slope_y, |
& streamice_bg_surf_slope_x, streamice_bg_surf_slope_y, |
15 |
& streamice_kx_b_init, streamice_ky_b_init |
& streamice_kx_b_init, streamice_ky_b_init, |
16 |
|
& streamice_wgt_drift, streamice_wgt_surf, |
17 |
|
& streamice_wgt_avthick, streamice_wgt_vel, |
18 |
|
& streamice_addl_backstress |
19 |
_RL streamice_density, streamice_density_ocean_avg |
_RL streamice_density, streamice_density_ocean_avg |
20 |
_RL A_glen_isothermal, n_glen, eps_glen_min |
_RL A_glen_isothermal, n_glen, eps_glen_min, eps_u_min |
21 |
_RL C_basal_fric_const |
_RL C_basal_fric_const |
22 |
_RL n_basal_friction |
_RL n_basal_friction |
23 |
_RL streamice_input_flux_unif |
_RL streamice_input_flux_unif |
28 |
_RL streamice_adjDump |
_RL streamice_adjDump |
29 |
_RL streamice_bg_surf_slope_x, streamice_bg_surf_slope_y |
_RL streamice_bg_surf_slope_x, streamice_bg_surf_slope_y |
30 |
_RL streamice_kx_b_init, streamice_ky_b_init |
_RL streamice_kx_b_init, streamice_ky_b_init |
31 |
|
_RL streamice_wgt_drift, streamice_wgt_surf |
32 |
|
_RL streamice_wgt_avthick, streamice_wgt_vel |
33 |
|
_RL streamice_addl_backstress |
34 |
|
|
35 |
|
|
36 |
C parms for parameterized initial thickness |
C parms for parameterized initial thickness |
37 |
C SHELF_MAX_DRAFT: max thickness of ice in m |
C SHELF_MAX_DRAFT: max thickness of ice in m |
172 |
& STREAMICE_ufacemask_bdry, |
& STREAMICE_ufacemask_bdry, |
173 |
& STREAMICE_vfacemask_bdry, |
& STREAMICE_vfacemask_bdry, |
174 |
& STREAMICE_float_cond, |
& STREAMICE_float_cond, |
175 |
& STREAMICE_calve_mask |
& STREAMICE_calve_mask, |
176 |
|
& STREAMICE_ctrl_mask |
177 |
_RS STREAMICE_hmask (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS STREAMICE_hmask (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
178 |
_RS STREAMICE_umask (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS STREAMICE_umask (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
179 |
_RS STREAMICE_vmask (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS STREAMICE_vmask (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
189 |
& (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
& (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
190 |
_RS STREAMICE_calve_mask |
_RS STREAMICE_calve_mask |
191 |
& (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
& (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
192 |
|
_RS STREAMICE_ctrl_mask |
193 |
|
& (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy,Nr) |
194 |
|
|
195 |
C NOTES : |
C NOTES : |
196 |
C STREAMICE_ufacemask_bdry, STREAMICE_vfacemask_bdry ARE CONSTANT (FIXED) |
C STREAMICE_ufacemask_bdry, STREAMICE_vfacemask_bdry ARE CONSTANT (FIXED) |
221 |
& C_basal_friction, |
& C_basal_friction, |
222 |
& A_glen, |
& A_glen, |
223 |
& BDOT_streamice, |
& BDOT_streamice, |
224 |
|
& streamice_sigma_coord, streamice_delsigma, |
225 |
|
& H_streamice_prev |
226 |
|
|
227 |
|
#ifdef STREAMICE_HYBRID_STRESS |
228 |
|
COMMON /STREAMICE_HYBRID/ |
229 |
& streamice_taubx, streamice_tauby, |
& streamice_taubx, streamice_tauby, |
230 |
& streamice_u_surf, streamice_v_surf, |
& streamice_u_surf, streamice_v_surf, |
231 |
& visc_streamice_full, streamice_omega, streamice_basal_geom, |
& visc_streamice_full, streamice_omega, streamice_basal_geom, |
232 |
& streamice_vert_shear_uz, streamice_vert_shear_vz, |
& streamice_vert_shear_uz, streamice_vert_shear_vz |
233 |
|
#endif |
234 |
|
|
235 |
#ifdef USE_ALT_RLOW |
#ifdef USE_ALT_RLOW |
236 |
& R_low_si, |
COMMON /STREAMICE_RLOW/ |
237 |
|
& R_low_si |
238 |
#endif |
#endif |
|
& streamice_sigma_coord, streamice_delsigma |
|
239 |
|
|
240 |
|
|
241 |
_RL H_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL H_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
255 |
_RL v_bdry_values_SI (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL v_bdry_values_SI (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
256 |
_RL C_basal_friction (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL C_basal_friction (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
257 |
_RL A_glen (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL A_glen (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
258 |
|
_RL streamice_sigma_coord (Nr) |
259 |
|
_RL streamice_delsigma (Nr) |
260 |
|
|
261 |
#ifdef USE_ALT_RLOW |
#ifdef USE_ALT_RLOW |
262 |
_RL R_low_si (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL R_low_si (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
263 |
#endif |
#endif |
264 |
|
|
265 |
C The following arrays are used for the hybrid stress balance |
C The following arrays are used for the hybrid stress balance |
266 |
|
#ifdef STREAMICE_HYBRID_STRESS |
267 |
_RL streamice_taubx (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL streamice_taubx (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
268 |
_RL streamice_tauby (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL streamice_tauby (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
269 |
_RL streamice_u_surf (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL streamice_u_surf (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
275 |
& (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
& (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
276 |
_RL streamice_vert_shear_uz (Nr) |
_RL streamice_vert_shear_uz (Nr) |
277 |
_RL streamice_vert_shear_vz (Nr) |
_RL streamice_vert_shear_vz (Nr) |
278 |
_RL streamice_sigma_coord (Nr) |
#endif |
279 |
_RL streamice_delsigma (Nr) |
|
280 |
|
|
281 |
!! IMPORTANT: MELT RATE IN METERS PER YEAR |
!! IMPORTANT: MELT RATE IN METERS PER YEAR |
282 |
!! POSITIVE WHERE MELTING |
!! POSITIVE WHERE MELTING |
283 |
_RL BDOT_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL BDOT_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
284 |
|
_RL H_streamice_prev (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
285 |
_RL STREAMICE_dummy_array (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL STREAMICE_dummy_array (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
286 |
|
|
287 |
|
|
288 |
|
|
289 |
COMMON /STREAMICE_COST_RL/ |
COMMON /STREAMICE_COST_RL/ |
290 |
& cost_func1_streamice |
& cost_func1_streamice |