/[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.1 by heimbach, Thu Mar 29 15:59:20 2012 UTC revision 1.17 by dgoldberg, Mon Dec 10 02:34:45 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 6  C     -------------------------- REAL PA Line 6  C     -------------------------- REAL PA
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_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         & 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
24        _RL streamice_vel_update        _RL streamice_vel_update
25        _RL streamice_cg_tol, streamice_nonlin_tol        _RL streamice_cg_tol, streamice_nonlin_tol
26          _RL streamice_nonlin_tol_fp
27        _RL streamice_CFL_factor        _RL streamice_CFL_factor
28          _RL streamice_adjDump
29          _RL streamice_bg_surf_slope_x, streamice_bg_surf_slope_y
30          _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
# Line 44  C     FLOW_DIR: 1.0=west, 2.0=east, 3.0= Line 58  C     FLOW_DIR: 1.0=west, 2.0=east, 3.0=
58    
59  C     -------------------------- INT PARAMS ---------------------------------------------------  C     -------------------------- INT PARAMS ---------------------------------------------------
60    
61          INTEGER streamice_max_nl
62          PARAMETER ( streamice_max_nl = 100 )
63    
64        COMMON /STREAMICE_PARMS_I/        COMMON /STREAMICE_PARMS_I/
65       &     streamice_max_cg_iter, streamice_max_nl_iter,       &     streamice_max_cg_iter, streamice_max_nl_iter,
66       &     streamice_vel_upd_counter, streamice_nstep_velocity,       &     streamice_vel_upd_counter, streamice_nstep_velocity,
# Line 57  C     -------------------------- CHAR PA Line 74  C     -------------------------- CHAR PA
74        CHARACTER*(MAX_LEN_FNAM) STREAMICEthickFile        CHARACTER*(MAX_LEN_FNAM) STREAMICEthickFile
75        CHARACTER*(MAX_LEN_FNAM) STREAMICEthickInit        CHARACTER*(MAX_LEN_FNAM) STREAMICEthickInit
76        CHARACTER*(MAX_LEN_FNAM) STREAMICEcalveMaskFile        CHARACTER*(MAX_LEN_FNAM) STREAMICEcalveMaskFile
77          CHARACTER*(MAX_LEN_FNAM) STREAMICEsigcoordInit
78          CHARACTER*(MAX_LEN_FNAM) STREAMICEsigcoordFile
79          CHARACTER*(MAX_LEN_FNAM) STREAMICEdelsigFile
80          CHARACTER*(MAX_LEN_FNAM) STREAMICEbasalTracConfig
81          CHARACTER*(MAX_LEN_FNAM) STREAMICEbasalTracFile
82          CHARACTER*(MAX_LEN_FNAM) STREAMICEvelOptimFile
83        COMMON /STREAMICE_PARM_C/        COMMON /STREAMICE_PARM_C/
84       &     STREAMICEthickInit,       &     STREAMICEthickInit,
85       &     STREAMICEthickFile,       &     STREAMICEthickFile,
86       &     STREAMICEcalveMaskFile       &     STREAMICEcalveMaskFile,
87         &     STREAMICEsigcoordInit,
88         &     STREAMICEsigcoordFile,
89         &     STREAMICEdelsigFile,
90         &     STREAMICEbasalTracConfig,
91         &     STREAMICEbasalTracFile,
92         &     STREAMICEvelOptimFile
93            
94  C     -------------------------- LOGICAL PARAMS ---------------------------------------------------  C     -------------------------- LOGICAL PARAMS ---------------------------------------------------
95    
# Line 73  C     -------------------------- LOGICAL Line 102  C     -------------------------- LOGICAL
102        LOGICAL STREAMICE_calve_to_mask        LOGICAL STREAMICE_calve_to_mask
103        LOGICAL STREAMICE_construct_matrix        LOGICAL STREAMICE_construct_matrix
104        LOGICAL STREAMICE_lower_cg_tol        LOGICAL STREAMICE_lower_cg_tol
105          LOGICAL STREAMICE_diagnostic_only
106          LOGICAL STREAMICE_ppm_driving_stress
107          LOGICAL STREAMICE_h_ctrl_const_surf
108          
109    C     The following parameters specify periodic boundary conditions.
110    C     For now this will completely override all other boundary conditions
111    C     and apply to the entire boundary
112    
113          LOGICAL STREAMICE_NS_periodic
114          LOGICAL STREAMICE_EW_periodic
115          
116    C      LOGICAL STREAMICE_hybrid_stress
117          
118        COMMON /STREAMICE_PARM_L/        COMMON /STREAMICE_PARM_L/
119       & STREAMICEison,       & STREAMICEison,
120       & STREAMICE_dump_mdsio, STREAMICE_tave_mdsio,       & STREAMICE_dump_mdsio, STREAMICE_tave_mdsio,
# Line 80  C     -------------------------- LOGICAL Line 122  C     -------------------------- LOGICAL
122       & STREAMICE_GL_regularize, STREAMICE_move_front,       & STREAMICE_GL_regularize, STREAMICE_move_front,
123       & STREAMICE_calve_to_mask,       & STREAMICE_calve_to_mask,
124       & STREAMICE_construct_matrix,       & STREAMICE_construct_matrix,
125       & STREAMICE_lower_cg_tol       & STREAMICE_lower_cg_tol,
126         & STREAMICE_NS_periodic, STREAMICE_EW_periodic,
127         & STREAMICE_diagnostic_only,
128         & STREAMICE_ppm_driving_stress,
129         & STREAMICE_h_ctrl_const_surf
130    
131  C     -------------------------- AND NOW ARRAYS ---------------------------------------------------  C     -------------------------- AND NOW ARRAYS ---------------------------------------------------
132    
# Line 126  C     Short arrays (e.g. masks) Line 172  C     Short arrays (e.g. masks)
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)
# Line 142  C     Short arrays (e.g. masks) Line 189  C     Short arrays (e.g. masks)
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)
# Line 169  C    REAL ARRAYS Line 218  C    REAL ARRAYS
218       &     u_bdry_values_SI,       &     u_bdry_values_SI,
219       &     v_bdry_values_SI,       &     v_bdry_values_SI,
220       &     STREAMICE_dummy_array,       &     STREAMICE_dummy_array,
221       &     C_basal_friction       &     C_basal_friction,
222         &     A_glen,
223         &     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,
230         &     streamice_u_surf, streamice_v_surf,
231         &     visc_streamice_full, streamice_omega, streamice_basal_geom,
232         &     streamice_vert_shear_uz, streamice_vert_shear_vz    
233    #endif
234    
235    #ifdef USE_ALT_RLOW
236          COMMON /STREAMICE_RLOW/
237         &     R_low_si
238    #endif
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)
242        _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)
243        _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 254  C    REAL ARRAYS
254        _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)
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)
258          _RL streamice_sigma_coord (Nr)
259          _RL streamice_delsigma (Nr)      
260    
261    #ifdef USE_ALT_RLOW
262          _RL R_low_si    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
263    #endif
264    
265    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)
268          _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)
270          _RL streamice_v_surf (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
271          _RL streamice_omega (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
272          _RL streamice_basal_geom
273         & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
274          _RL visc_streamice_full
275         & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
276          _RL streamice_vert_shear_uz (Nr)
277          _RL streamice_vert_shear_vz (Nr)
278    #endif      
279                
280          
281    !! IMPORTANT: MELT RATE IN METERS PER YEAR
282    !! POSITIVE WHERE MELTING
283          _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/
290         &       cost_func1_streamice
291          _RL cost_func1_streamice(nSx,nSy)
292          
293  C    NOTES :  C    NOTES :
294  C      REAL ARRAYS THAT COMPRISE "STATE":  C      REAL ARRAYS THAT COMPRISE "STATE":
295  C       H_streamice,  C       H_streamice,
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22