--- MITgcm_contrib/dgoldberg/streamice/STREAMICE.h 2012/10/03 14:54:31 1.12 +++ MITgcm_contrib/dgoldberg/streamice/STREAMICE.h 2012/10/10 02:26:19 1.14 @@ -12,7 +12,8 @@ & streamice_nonlin_tol_fp, & streamice_CFL_factor, streamice_adjDump, & streamice_bg_surf_slope_x, streamice_bg_surf_slope_y, - & streamice_kx_b_init, streamice_ky_b_init + & streamice_kx_b_init, streamice_ky_b_init, + & streamice_wgt_drift, streamice_wgt_surf _RL streamice_density, streamice_density_ocean_avg _RL A_glen_isothermal, n_glen, eps_glen_min _RL C_basal_fric_const @@ -25,6 +26,8 @@ _RL streamice_adjDump _RL streamice_bg_surf_slope_x, streamice_bg_surf_slope_y _RL streamice_kx_b_init, streamice_ky_b_init + _RL streamice_wgt_drift, streamice_wgt_surf + C parms for parameterized initial thickness C SHELF_MAX_DRAFT: max thickness of ice in m @@ -215,10 +218,13 @@ & streamice_u_surf, streamice_v_surf, & visc_streamice_full, streamice_omega, streamice_basal_geom, & streamice_vert_shear_uz, streamice_vert_shear_vz, + & streamice_sigma_coord, streamice_delsigma, + & H_streamice_prev + #ifdef USE_ALT_RLOW + COMMON /STREAMICE_RLOW/ & R_low_si, #endif - & streamice_sigma_coord, streamice_delsigma _RL H_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) @@ -261,9 +267,10 @@ !! IMPORTANT: MELT RATE IN METERS PER YEAR !! POSITIVE WHERE MELTING _RL BDOT_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) - + _RL H_streamice_prev (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL STREAMICE_dummy_array (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) + COMMON /STREAMICE_COST_RL/ & cost_func1_streamice