| 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 | 
  | 
      & B_glen_isothermal, n_glen, eps_glen_min, eps_u_min, | 
| 11 | 
      & C_basal_fric_const, n_basal_friction, streamice_input_flux_unif, | 
      & C_basal_fric_const, n_basal_friction, streamice_input_flux_unif, | 
| 12 | 
      & streamice_vel_update, streamice_cg_tol, streamice_nonlin_tol, | 
      & streamice_vel_update, streamice_cg_tol, streamice_nonlin_tol, | 
| 13 | 
      & streamice_nonlin_tol_fp, | 
      & streamice_nonlin_tol_fp, | 
| 14 | 
      & streamice_CFL_factor, streamice_adjDump,  | 
      & streamice_CFL_factor, streamice_adjDump,  | 
| 15 | 
      & streamice_bg_surf_slope_x, streamice_bg_surf_slope_y, | 
      & streamice_bg_surf_slope_x, streamice_bg_surf_slope_y, | 
| 16 | 
      & streamice_kx_b_init, streamice_ky_b_init | 
      & streamice_kx_b_init, streamice_ky_b_init, | 
| 17 | 
  | 
      & streamice_wgt_drift, streamice_wgt_surf, | 
| 18 | 
  | 
      & streamice_wgt_avthick, streamice_wgt_vel,  | 
| 19 | 
  | 
      & streamice_wgt_tikh, | 
| 20 | 
  | 
      & streamice_addl_backstress, | 
| 21 | 
  | 
      & streamice_smooth_gl_width | 
| 22 | 
       _RL streamice_density, streamice_density_ocean_avg | 
       _RL streamice_density, streamice_density_ocean_avg | 
| 23 | 
       _RL A_glen_isothermal, n_glen, eps_glen_min | 
 !      _RL A_glen_isothermal, n_glen, eps_glen_min, eps_u_min | 
| 24 | 
  | 
       _RL B_glen_isothermal, n_glen, eps_glen_min, eps_u_min | 
| 25 | 
       _RL C_basal_fric_const | 
       _RL C_basal_fric_const | 
| 26 | 
       _RL n_basal_friction  | 
       _RL n_basal_friction  | 
| 27 | 
       _RL streamice_input_flux_unif | 
       _RL streamice_input_flux_unif | 
| 32 | 
       _RL streamice_adjDump | 
       _RL streamice_adjDump | 
| 33 | 
       _RL streamice_bg_surf_slope_x, streamice_bg_surf_slope_y | 
       _RL streamice_bg_surf_slope_x, streamice_bg_surf_slope_y | 
| 34 | 
       _RL streamice_kx_b_init, streamice_ky_b_init | 
       _RL streamice_kx_b_init, streamice_ky_b_init | 
| 35 | 
  | 
       _RL streamice_wgt_drift, streamice_wgt_surf | 
| 36 | 
  | 
       _RL streamice_wgt_avthick, streamice_wgt_vel | 
| 37 | 
  | 
       _RL streamice_wgt_tikh | 
| 38 | 
  | 
       _RL streamice_addl_backstress | 
| 39 | 
  | 
       _RL streamice_smooth_gl_width | 
| 40 | 
  | 
        | 
| 41 | 
        | 
        | 
| 42 | 
 C     parms for parameterized initial thickness | 
 C     parms for parameterized initial thickness | 
| 43 | 
 C     SHELF_MAX_DRAFT: max thickness of ice in m | 
 C     SHELF_MAX_DRAFT: max thickness of ice in m | 
| 84 | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEsigcoordFile | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEsigcoordFile | 
| 85 | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEdelsigFile | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEdelsigFile | 
| 86 | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEbasalTracConfig | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEbasalTracConfig | 
| 87 | 
  | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEGlenConstConfig | 
| 88 | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEbasalTracFile | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEbasalTracFile | 
| 89 | 
  | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEGlenConstFile | 
| 90 | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEvelOptimFile | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEvelOptimFile | 
| 91 | 
  | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEtopogFile | 
| 92 | 
  | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEcostMaskFile | 
| 93 | 
  | 
        | 
| 94 | 
  | 
 !     THE FOLLOWING FILENAMES ARE FOR SPECIFYING IRREGULAR DOMAIN GEOMETRIES  | 
| 95 | 
  | 
 !     (i.e. boundaries that do not conform with rectangular walls) | 
| 96 | 
  | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEhmaskFile | 
| 97 | 
  | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEuFaceBdryFile | 
| 98 | 
  | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEvFaceBdryFile | 
| 99 | 
  | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEuDirichValsFile | 
| 100 | 
  | 
       CHARACTER*(MAX_LEN_FNAM) STREAMICEvDirichValsFile | 
| 101 | 
  | 
        | 
| 102 | 
  | 
        | 
| 103 | 
       COMMON /STREAMICE_PARM_C/  | 
       COMMON /STREAMICE_PARM_C/  | 
| 104 | 
      &     STREAMICEthickInit, | 
      &     STREAMICEthickInit, | 
| 105 | 
      &     STREAMICEthickFile, | 
      &     STREAMICEthickFile, | 
| 109 | 
      &     STREAMICEdelsigFile, | 
      &     STREAMICEdelsigFile, | 
| 110 | 
      &     STREAMICEbasalTracConfig, | 
      &     STREAMICEbasalTracConfig, | 
| 111 | 
      &     STREAMICEbasalTracFile, | 
      &     STREAMICEbasalTracFile, | 
| 112 | 
      &     STREAMICEvelOptimFile | 
      &     STREAMICEvelOptimFile, | 
| 113 | 
  | 
      &     STREAMICEtopogFile, | 
| 114 | 
  | 
      &     STREAMICEhmaskFile, | 
| 115 | 
  | 
      &     STREAMICEuFaceBdryFile, | 
| 116 | 
  | 
      &     STREAMICEvFaceBdryFile, | 
| 117 | 
  | 
      &     STREAMICEuDirichValsFile, | 
| 118 | 
  | 
      &     STREAMICEvDirichValsFile, | 
| 119 | 
  | 
      &     STREAMICEGlenConstFile, | 
| 120 | 
  | 
      &     STREAMICEGlenConstConfig, | 
| 121 | 
  | 
      &     STREAMICEcostMaskFile | 
| 122 | 
       | 
       | 
| 123 | 
 C     -------------------------- LOGICAL PARAMS --------------------------------------------------- | 
 C     -------------------------- LOGICAL PARAMS --------------------------------------------------- | 
| 124 | 
  | 
  | 
| 201 | 
      &     STREAMICE_ufacemask_bdry, | 
      &     STREAMICE_ufacemask_bdry, | 
| 202 | 
      &     STREAMICE_vfacemask_bdry, | 
      &     STREAMICE_vfacemask_bdry, | 
| 203 | 
      &     STREAMICE_float_cond, | 
      &     STREAMICE_float_cond, | 
| 204 | 
      &     STREAMICE_calve_mask | 
      &     STREAMICE_calve_mask, | 
| 205 | 
  | 
      &     STREAMICE_ctrl_mask, | 
| 206 | 
  | 
      &     STREAMICE_cost_mask | 
| 207 | 
       _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) | 
| 208 | 
       _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) | 
| 209 | 
       _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) | 
| 219 | 
      &  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
      &  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
| 220 | 
       _RS STREAMICE_calve_mask  | 
       _RS STREAMICE_calve_mask  | 
| 221 | 
      & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
      & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
| 222 | 
  | 
       _RS STREAMICE_ctrl_mask  | 
| 223 | 
  | 
      & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) | 
| 224 | 
  | 
       _RS STREAMICE_cost_mask  | 
| 225 | 
  | 
      & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
| 226 | 
       | 
       | 
| 227 | 
 C    NOTES : | 
 C    NOTES : | 
| 228 | 
 C     STREAMICE_ufacemask_bdry, STREAMICE_vfacemask_bdry ARE CONSTANT (FIXED) | 
 C     STREAMICE_ufacemask_bdry, STREAMICE_vfacemask_bdry ARE CONSTANT (FIXED) | 
| 251 | 
      &     v_bdry_values_SI, | 
      &     v_bdry_values_SI, | 
| 252 | 
      &     STREAMICE_dummy_array, | 
      &     STREAMICE_dummy_array, | 
| 253 | 
      &     C_basal_friction, | 
      &     C_basal_friction, | 
| 254 | 
      &     A_glen, | 
 !     &     A_glen, | 
| 255 | 
      &     BDOT_streamice, | 
      &     B_glen, | 
| 256 | 
  | 
      &     BDOT_streamice, ADOT_streamice,  ! mass balances in meters per year | 
| 257 | 
  | 
      &     streamice_sigma_coord, streamice_delsigma, | 
| 258 | 
  | 
      &     H_streamice_prev | 
| 259 | 
  | 
  | 
| 260 | 
  | 
 #ifdef STREAMICE_HYBRID_STRESS | 
| 261 | 
  | 
       COMMON /STREAMICE_HYBRID/  | 
| 262 | 
      &     streamice_taubx, streamice_tauby, | 
      &     streamice_taubx, streamice_tauby, | 
| 263 | 
      &     streamice_u_surf, streamice_v_surf, | 
      &     streamice_u_surf, streamice_v_surf, | 
| 264 | 
      &     visc_streamice_full, streamice_omega, streamice_basal_geom, | 
      &     visc_streamice_full, streamice_omega, streamice_basal_geom, | 
| 265 | 
      &     streamice_vert_shear_uz, streamice_vert_shear_vz, | 
      &     streamice_vert_shear_uz, streamice_vert_shear_vz      | 
| 266 | 
  | 
 #endif | 
| 267 | 
  | 
  | 
| 268 | 
 #ifdef USE_ALT_RLOW | 
 #ifdef USE_ALT_RLOW | 
| 269 | 
      &     R_low_si, | 
       COMMON /STREAMICE_RLOW/  | 
| 270 | 
  | 
      &     R_low_si | 
| 271 | 
 #endif | 
 #endif | 
 | 
      &     streamice_sigma_coord, streamice_delsigma | 
  | 
| 272 | 
  | 
  | 
| 273 | 
  | 
  | 
| 274 | 
       _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) | 
| 287 | 
       _RL u_bdry_values_SI    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
       _RL u_bdry_values_SI    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
| 288 | 
       _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) | 
| 289 | 
       _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) | 
| 290 | 
       _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) | 
| 291 | 
  | 
       _RL B_glen    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
| 292 | 
  | 
       _RL streamice_sigma_coord (Nr) | 
| 293 | 
  | 
       _RL streamice_delsigma (Nr)       | 
| 294 | 
  | 
  | 
| 295 | 
 #ifdef USE_ALT_RLOW | 
 #ifdef USE_ALT_RLOW | 
| 296 | 
       _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) | 
| 297 | 
 #endif | 
 #endif | 
| 298 | 
  | 
  | 
| 299 | 
 C      The following arrays are used for the hybrid stress balance       | 
 C     The following arrays are used for the hybrid stress balance             | 
| 300 | 
        | 
 #ifdef STREAMICE_HYBRID_STRESS       | 
| 301 | 
       _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) | 
| 302 | 
       _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) | 
| 303 | 
       _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) | 
| 309 | 
      & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) | 
      & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) | 
| 310 | 
       _RL streamice_vert_shear_uz (Nr) | 
       _RL streamice_vert_shear_uz (Nr) | 
| 311 | 
       _RL streamice_vert_shear_vz (Nr) | 
       _RL streamice_vert_shear_vz (Nr) | 
| 312 | 
       _RL streamice_sigma_coord (Nr) | 
 #endif       | 
| 313 | 
       _RL streamice_delsigma (Nr) | 
        | 
| 314 | 
        | 
        | 
| 315 | 
  | 
       _RL ADOT_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
| 316 | 
 !! IMPORTANT: MELT RATE IN METERS PER YEAR | 
 !! IMPORTANT: MELT RATE IN METERS PER YEAR | 
| 317 | 
 !! POSITIVE WHERE MELTING | 
 !! POSITIVE WHERE MELTING | 
| 318 | 
       _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) | 
| 319 | 
        | 
       _RL H_streamice_prev (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
| 320 | 
       _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) | 
| 321 | 
        | 
        | 
| 322 | 
  | 
        | 
| 323 | 
             | 
             | 
| 324 | 
       COMMON /STREAMICE_COST_RL/ | 
       COMMON /STREAMICE_COST_RL/ | 
| 325 | 
      &       cost_func1_streamice | 
      &       cost_func1_streamice |