/[MITgcm]/MITgcm_contrib/dgoldberg/streamice/STREAMICE.h
ViewVC logotype

Diff of /MITgcm_contrib/dgoldberg/streamice/STREAMICE.h

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by dgoldberg, Mon May 14 16:47:10 2012 UTC revision 1.14 by dgoldberg, Wed Oct 10 02:26:19 2012 UTC
# Line 1  Line 1 
1  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1--+-+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
2    
3  #ifdef ALLOW_STREAMICE  #ifdef ALLOW_STREAMICE
4    
# Line 9  C     -------------------------- REAL PA Line 9  C     -------------------------- REAL PA
9       & A_glen_isothermal, n_glen, eps_glen_min,       & A_glen_isothermal, n_glen, eps_glen_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_CFL_factor       & streamice_nonlin_tol_fp,
13         & streamice_CFL_factor, streamice_adjDump,
14         & streamice_bg_surf_slope_x, streamice_bg_surf_slope_y,
15         & streamice_kx_b_init, streamice_ky_b_init,
16         & streamice_wgt_drift, streamice_wgt_surf    
17        _RL streamice_density, streamice_density_ocean_avg        _RL streamice_density, streamice_density_ocean_avg
18        _RL A_glen_isothermal, n_glen, eps_glen_min        _RL A_glen_isothermal, n_glen, eps_glen_min
19        _RL C_basal_fric_const        _RL C_basal_fric_const
# Line 17  C     -------------------------- REAL PA Line 21  C     -------------------------- REAL PA
21        _RL streamice_input_flux_unif        _RL streamice_input_flux_unif
22        _RL streamice_vel_update        _RL streamice_vel_update
23        _RL streamice_cg_tol, streamice_nonlin_tol        _RL streamice_cg_tol, streamice_nonlin_tol
24          _RL streamice_nonlin_tol_fp
25        _RL streamice_CFL_factor        _RL streamice_CFL_factor
26          _RL streamice_adjDump
27          _RL streamice_bg_surf_slope_x, streamice_bg_surf_slope_y
28          _RL streamice_kx_b_init, streamice_ky_b_init
29          _RL streamice_wgt_drift, streamice_wgt_surf    
30          
31                
32  C     parms for parameterized initial thickness  C     parms for parameterized initial thickness
33  C     SHELF_MAX_DRAFT: max thickness of ice in m  C     SHELF_MAX_DRAFT: max thickness of ice in m
# Line 44  C     FLOW_DIR: 1.0=west, 2.0=east, 3.0= Line 54  C     FLOW_DIR: 1.0=west, 2.0=east, 3.0=
54    
55  C     -------------------------- INT PARAMS ---------------------------------------------------  C     -------------------------- INT PARAMS ---------------------------------------------------
56    
57          INTEGER streamice_max_nl
58          PARAMETER ( streamice_max_nl = 100 )
59    
60        COMMON /STREAMICE_PARMS_I/        COMMON /STREAMICE_PARMS_I/
61       &     streamice_max_cg_iter, streamice_max_nl_iter,       &     streamice_max_cg_iter, streamice_max_nl_iter,
62       &     streamice_vel_upd_counter, streamice_nstep_velocity,       &     streamice_vel_upd_counter, streamice_nstep_velocity,
# Line 57  C     -------------------------- CHAR PA Line 70  C     -------------------------- CHAR PA
70        CHARACTER*(MAX_LEN_FNAM) STREAMICEthickFile        CHARACTER*(MAX_LEN_FNAM) STREAMICEthickFile
71        CHARACTER*(MAX_LEN_FNAM) STREAMICEthickInit        CHARACTER*(MAX_LEN_FNAM) STREAMICEthickInit
72        CHARACTER*(MAX_LEN_FNAM) STREAMICEcalveMaskFile        CHARACTER*(MAX_LEN_FNAM) STREAMICEcalveMaskFile
73          CHARACTER*(MAX_LEN_FNAM) STREAMICEsigcoordInit
74          CHARACTER*(MAX_LEN_FNAM) STREAMICEsigcoordFile
75          CHARACTER*(MAX_LEN_FNAM) STREAMICEdelsigFile
76          CHARACTER*(MAX_LEN_FNAM) STREAMICEbasalTracConfig
77          CHARACTER*(MAX_LEN_FNAM) STREAMICEbasalTracFile
78          CHARACTER*(MAX_LEN_FNAM) STREAMICEvelOptimFile
79        COMMON /STREAMICE_PARM_C/        COMMON /STREAMICE_PARM_C/
80       &     STREAMICEthickInit,       &     STREAMICEthickInit,
81       &     STREAMICEthickFile,       &     STREAMICEthickFile,
82       &     STREAMICEcalveMaskFile       &     STREAMICEcalveMaskFile,
83         &     STREAMICEsigcoordInit,
84         &     STREAMICEsigcoordFile,
85         &     STREAMICEdelsigFile,
86         &     STREAMICEbasalTracConfig,
87         &     STREAMICEbasalTracFile,
88         &     STREAMICEvelOptimFile
89            
90  C     -------------------------- LOGICAL PARAMS ---------------------------------------------------  C     -------------------------- LOGICAL PARAMS ---------------------------------------------------
91    
# Line 71  C     -------------------------- LOGICAL Line 96  C     -------------------------- LOGICAL
96        LOGICAL STREAMICE_tave_mnc        LOGICAL STREAMICE_tave_mnc
97        LOGICAL STREAMICE_GL_regularize, STREAMICE_move_front        LOGICAL STREAMICE_GL_regularize, STREAMICE_move_front
98        LOGICAL STREAMICE_calve_to_mask        LOGICAL STREAMICE_calve_to_mask
99  !       LOGICAL STREAMICE_construct_matrix        LOGICAL STREAMICE_construct_matrix
100        LOGICAL STREAMICE_lower_cg_tol        LOGICAL STREAMICE_lower_cg_tol
101          LOGICAL STREAMICE_diagnostic_only
102          LOGICAL STREAMICE_ppm_driving_stress
103          LOGICAL STREAMICE_h_ctrl_const_surf
104          
105    C     The following parameters specify periodic boundary conditions.
106    C     For now this will completely override all other boundary conditions
107    C     and apply to the entire boundary
108    
109          LOGICAL STREAMICE_NS_periodic
110          LOGICAL STREAMICE_EW_periodic
111          
112    C      LOGICAL STREAMICE_hybrid_stress
113          
114        COMMON /STREAMICE_PARM_L/        COMMON /STREAMICE_PARM_L/
115       & STREAMICEison,       & STREAMICEison,
116       & STREAMICE_dump_mdsio, STREAMICE_tave_mdsio,       & STREAMICE_dump_mdsio, STREAMICE_tave_mdsio,
117       & STREAMICE_dump_mnc, STREAMICE_tave_mnc,       & STREAMICE_dump_mnc, STREAMICE_tave_mnc,
118       & STREAMICE_GL_regularize, STREAMICE_move_front,       & STREAMICE_GL_regularize, STREAMICE_move_front,
119       & STREAMICE_calve_to_mask,       & STREAMICE_calve_to_mask,
120  !     & STREAMICE_construct_matrix,       & STREAMICE_construct_matrix,
121       & STREAMICE_lower_cg_tol       & STREAMICE_lower_cg_tol,
122         & STREAMICE_NS_periodic, STREAMICE_EW_periodic,
123         & STREAMICE_diagnostic_only,
124         & STREAMICE_ppm_driving_stress,
125         & STREAMICE_h_ctrl_const_surf
126    
127  C     -------------------------- AND NOW ARRAYS ---------------------------------------------------  C     -------------------------- AND NOW ARRAYS ---------------------------------------------------
128    
# Line 169  C    REAL ARRAYS Line 211  C    REAL ARRAYS
211       &     u_bdry_values_SI,       &     u_bdry_values_SI,
212       &     v_bdry_values_SI,       &     v_bdry_values_SI,
213       &     STREAMICE_dummy_array,       &     STREAMICE_dummy_array,
214       &     C_basal_friction       &     C_basal_friction,
215         &     A_glen,
216         &     BDOT_streamice,
217         &     streamice_taubx, streamice_tauby,
218         &     streamice_u_surf, streamice_v_surf,
219         &     visc_streamice_full, streamice_omega, streamice_basal_geom,
220         &     streamice_vert_shear_uz, streamice_vert_shear_vz,
221         &     streamice_sigma_coord, streamice_delsigma,
222         &     H_streamice_prev
223    
224    #ifdef USE_ALT_RLOW
225          COMMON /STREAMICE_RLOW/
226         &     R_low_si,
227    #endif
228    
229    
230        _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)
231        _RL U_streamice           (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL U_streamice           (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
232        _RL V_streamice           (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL V_streamice           (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
# Line 186  C    REAL ARRAYS Line 243  C    REAL ARRAYS
243        _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)
244        _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)
245        _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)
246          _RL A_glen    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
247    #ifdef USE_ALT_RLOW
248          _RL R_low_si    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
249    #endif
250    
251    C      The following arrays are used for the hybrid stress balance      
252          
253          _RL streamice_taubx (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
254          _RL streamice_tauby (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
255          _RL streamice_u_surf (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
256          _RL streamice_v_surf (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
257          _RL streamice_omega (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
258          _RL streamice_basal_geom
259         & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
260          _RL visc_streamice_full
261         & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
262          _RL streamice_vert_shear_uz (Nr)
263          _RL streamice_vert_shear_vz (Nr)
264          _RL streamice_sigma_coord (Nr)
265          _RL streamice_delsigma (Nr)
266                
267    !! IMPORTANT: MELT RATE IN METERS PER YEAR
268    !! POSITIVE WHERE MELTING
269          _RL BDOT_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
270          _RL H_streamice_prev (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
271        _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)
272                
273          
274              
275          COMMON /STREAMICE_COST_RL/
276         &       cost_func1_streamice
277          _RL cost_func1_streamice(nSx,nSy)
278          
279  C    NOTES :  C    NOTES :
280  C      REAL ARRAYS THAT COMPRISE "STATE":  C      REAL ARRAYS THAT COMPRISE "STATE":
281  C       H_streamice,  C       H_streamice,
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22