/[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.8 by dgoldberg, Wed Sep 26 15:16:51 2012 UTC revision 1.22 by dgoldberg, Sat Jun 8 22:15:33 2013 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         & 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
# Line 25  C     -------------------------- REAL PA Line 32  C     -------------------------- REAL PA
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
# Line 71  C     -------------------------- CHAR PA Line 84  C     -------------------------- CHAR PA
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          CHARACTER*(MAX_LEN_FNAM) STREAMICE_ADV_SCHEME
94          
95    !     THE FOLLOWING FILENAMES ARE FOR SPECIFYING IRREGULAR DOMAIN GEOMETRIES
96    !     (i.e. boundaries that do not conform with rectangular walls)
97          CHARACTER*(MAX_LEN_FNAM) STREAMICEhmaskFile
98          CHARACTER*(MAX_LEN_FNAM) STREAMICEHBCxFile
99          CHARACTER*(MAX_LEN_FNAM) STREAMICEHBCyFile
100          CHARACTER*(MAX_LEN_FNAM) STREAMICEuFaceBdryFile
101          CHARACTER*(MAX_LEN_FNAM) STREAMICEvFaceBdryFile
102          CHARACTER*(MAX_LEN_FNAM) STREAMICEuDirichValsFile
103          CHARACTER*(MAX_LEN_FNAM) STREAMICEvDirichValsFile
104    
105    #ifdef ALLOW_PETSC
106    !     CHARACTER PARAMS FOR PETSC
107          CHARACTER*(MAX_LEN_FNAM) PETSC_SOLVER_TYPE
108          CHARACTER*(MAX_LEN_FNAM) PETSC_PRECOND_TYPE
109    #endif
110        
111    #ifdef ALLOW_STREAMICE_2DTRACER
112    !     CHARACTER PARAMS FOR TRACER
113          CHARACTER*(MAX_LEN_FNAM) STREAMICETrac2DBCxFile
114          CHARACTER*(MAX_LEN_FNAM) STREAMICETrac2DBCyFile
115          CHARACTER*(MAX_LEN_FNAM) STREAMICETrac2DinitFile
116    #endif
117        
118        COMMON /STREAMICE_PARM_C/        COMMON /STREAMICE_PARM_C/
119       &     STREAMICEthickInit,       &     STREAMICEthickInit,
120       &     STREAMICEthickFile,       &     STREAMICEthickFile,
# Line 82  C     -------------------------- CHAR PA Line 124  C     -------------------------- CHAR PA
124       &     STREAMICEdelsigFile,       &     STREAMICEdelsigFile,
125       &     STREAMICEbasalTracConfig,       &     STREAMICEbasalTracConfig,
126       &     STREAMICEbasalTracFile,       &     STREAMICEbasalTracFile,
127       &     STREAMICEvelOptimFile       &     STREAMICEvelOptimFile,
128         &     STREAMICEtopogFile,
129         &     STREAMICEhmaskFile,
130         &     STREAMICEHBCxFile,
131         &     STREAMICEHBCyFile,
132         &     STREAMICEuFaceBdryFile,
133         &     STREAMICEvFaceBdryFile,
134         &     STREAMICEuDirichValsFile,
135         &     STREAMICEvDirichValsFile,
136         &     STREAMICEGlenConstFile,
137         &     STREAMICEGlenConstConfig,
138         &     STREAMICEcostMaskFile,
139         &     STREAMICE_ADV_SCHEME
140    
141    #ifdef ALLOW_PETSC
142          COMMON /PETSC_PARM_C/
143         &     PETSC_SOLVER_TYPE,
144         &     PETSC_PRECOND_TYPE
145    #endif
146    
147    #ifdef ALLOW_STREAMICE_2DTRACER
148          COMMON /STREAMICE_TRAC2D_C/
149         &     STREAMICETrac2DBCxFile,
150         &     STREAMICETrac2DBCyFile,
151         &     STREAMICETrac2DinitFile
152    #endif
153            
154  C     -------------------------- LOGICAL PARAMS ---------------------------------------------------  C     -------------------------- LOGICAL PARAMS ---------------------------------------------------
155    
# Line 96  C     -------------------------- LOGICAL Line 163  C     -------------------------- LOGICAL
163        LOGICAL STREAMICE_construct_matrix        LOGICAL STREAMICE_construct_matrix
164        LOGICAL STREAMICE_lower_cg_tol        LOGICAL STREAMICE_lower_cg_tol
165        LOGICAL STREAMICE_diagnostic_only        LOGICAL STREAMICE_diagnostic_only
166          LOGICAL STREAMICE_ppm_driving_stress
167        LOGICAL STREAMICE_h_ctrl_const_surf        LOGICAL STREAMICE_h_ctrl_const_surf
168                
169  C     The following parameters specify periodic boundary conditions.  C     The following parameters specify periodic boundary conditions.
# Line 116  C      LOGICAL STREAMICE_hybrid_stress Line 184  C      LOGICAL STREAMICE_hybrid_stress
184       & STREAMICE_construct_matrix,       & STREAMICE_construct_matrix,
185       & STREAMICE_lower_cg_tol,       & STREAMICE_lower_cg_tol,
186       & STREAMICE_NS_periodic, STREAMICE_EW_periodic,       & STREAMICE_NS_periodic, STREAMICE_EW_periodic,
187       & STREAMICE_diagnostic_only       & STREAMICE_diagnostic_only,
188         & STREAMICE_ppm_driving_stress,
189         & STREAMICE_h_ctrl_const_surf
190    
191  C     -------------------------- AND NOW ARRAYS ---------------------------------------------------  C     -------------------------- AND NOW ARRAYS ---------------------------------------------------
192    
# Line 162  C     Short arrays (e.g. masks) Line 232  C     Short arrays (e.g. masks)
232       &     STREAMICE_ufacemask_bdry,       &     STREAMICE_ufacemask_bdry,
233       &     STREAMICE_vfacemask_bdry,       &     STREAMICE_vfacemask_bdry,
234       &     STREAMICE_float_cond,       &     STREAMICE_float_cond,
235       &     STREAMICE_calve_mask       &     STREAMICE_calve_mask,
236         &     STREAMICE_ctrl_mask,
237         &     STREAMICE_cost_mask
238        _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)
239        _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)
240        _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 178  C     Short arrays (e.g. masks) Line 250  C     Short arrays (e.g. masks)
250       &  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)       &  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
251        _RS STREAMICE_calve_mask        _RS STREAMICE_calve_mask
252       & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)       & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
253          _RS STREAMICE_ctrl_mask
254         & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
255          _RS STREAMICE_cost_mask
256         & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
257            
258  C    NOTES :  C    NOTES :
259  C     STREAMICE_ufacemask_bdry, STREAMICE_vfacemask_bdry ARE CONSTANT (FIXED)  C     STREAMICE_ufacemask_bdry, STREAMICE_vfacemask_bdry ARE CONSTANT (FIXED)
# Line 201  C    REAL ARRAYS Line 277  C    REAL ARRAYS
277       &     mass_ice_streamice,       &     mass_ice_streamice,
278       &     u_flux_bdry_SI,       &     u_flux_bdry_SI,
279       &     v_flux_bdry_SI,       &     v_flux_bdry_SI,
280       &     h_bdry_values_SI,       &     h_ubdry_values_SI,
281         &     h_vbdry_values_SI,
282       &     u_bdry_values_SI,       &     u_bdry_values_SI,
283       &     v_bdry_values_SI,       &     v_bdry_values_SI,
284       &     STREAMICE_dummy_array,       &     STREAMICE_dummy_array,
285       &     C_basal_friction,       &     C_basal_friction,
286       &     A_glen,  !     &     A_glen,
287       &     BDOT_streamice,       &     B_glen,
288         &     BDOT_streamice, ADOT_streamice,  ! mass balances in meters per year
289         &     streamice_sigma_coord, streamice_delsigma,
290         &     H_streamice_prev
291    
292    #ifdef STREAMICE_HYBRID_STRESS
293          COMMON /STREAMICE_HYBRID/
294       &     streamice_taubx, streamice_tauby,       &     streamice_taubx, streamice_tauby,
295       &     streamice_u_surf, streamice_v_surf,       &     streamice_u_surf, streamice_v_surf,
296       &     visc_streamice_full, streamice_omega, streamice_basal_geom,       &     visc_streamice_full, streamice_omega, streamice_basal_geom,
297       &     streamice_vert_shear_uz, streamice_vert_shear_vz,       &     streamice_vert_shear_uz, streamice_vert_shear_vz    
298       &     streamice_sigma_coord, streamice_delsigma  #endif
299    
300    #ifdef ALLOW_STREAMICE_2DTRACER
301          COMMON /STREAMICE_TRAC2D_FIELDS_RL/
302         &     trac2d_ubdry_values_SI,
303         &     trac2d_vbdry_values_SI,
304         &     trac2d
305    #ifdef STREAMICE_TRACER_AB
306          COMMON /STREAMICE_TRAC2D_AB_RL/
307         &     GAD_trac_2d
308    #endif
309    #endif
310    
311    #ifdef USE_ALT_RLOW
312          COMMON /STREAMICE_RLOW/
313         &     R_low_si
314    #endif
315    
316    
317        _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)
318        _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)
319        _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 225  C    REAL ARRAYS Line 326  C    REAL ARRAYS
326        _RL mass_ice_streamice    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL mass_ice_streamice    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
327        _RL u_flux_bdry_SI    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL u_flux_bdry_SI    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
328        _RL v_flux_bdry_SI    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL v_flux_bdry_SI    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
329        _RL h_bdry_values_SI    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL h_ubdry_values_SI    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
330          _RL h_vbdry_values_SI    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
331        _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)
332        _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)
333        _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)
334        _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)
335                _RL B_glen    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
336  C      The following arrays are used for the hybrid stress balance              _RL streamice_sigma_coord (Nr)
337                _RL streamice_delsigma (Nr)      
338    
339    #ifdef USE_ALT_RLOW
340          _RL R_low_si    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
341    #endif
342    
343    C     The following arrays are used for the hybrid stress balance            
344    #ifdef STREAMICE_HYBRID_STRESS      
345        _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)
346        _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)
347        _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)
# Line 244  C      The following arrays are used for Line 353  C      The following arrays are used for
353       & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)       & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
354        _RL streamice_vert_shear_uz (Nr)        _RL streamice_vert_shear_uz (Nr)
355        _RL streamice_vert_shear_vz (Nr)        _RL streamice_vert_shear_vz (Nr)
356        _RL streamice_sigma_coord (Nr)  #endif      
357        _RL streamice_delsigma (Nr)      
358    #ifdef ALLOW_STREAMICE_2DTRACER
359          _RL trac2d_ubdry_values_SI
360         &   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
361          _RL trac2d_vbdry_values_SI
362         &   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
363          _RL trac2d (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
364    #ifdef STREAMICE_TRACER_AB
365          _RL GAD_trac_2d (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
366    #endif
367    #endif
368                
369          _RL ADOT_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
370  !! IMPORTANT: MELT RATE IN METERS PER YEAR  !! IMPORTANT: MELT RATE IN METERS PER YEAR
371  !! POSITIVE WHERE MELTING  !! POSITIVE WHERE MELTING
372        _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)
373                _RL H_streamice_prev (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
374        _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)
375                
376          
377              
378          COMMON /STREAMICE_COST_RL/
379         &       cost_func1_streamice
380          _RL cost_func1_streamice(nSx,nSy)
381          
382  C    NOTES :  C    NOTES :
383  C      REAL ARRAYS THAT COMPRISE "STATE":  C      REAL ARRAYS THAT COMPRISE "STATE":
384  C       H_streamice,  C       H_streamice,
# Line 266  C Line 392  C
392  C       visc & tau are now calculated based on U,V in streamice_vel_solve  C       visc & tau are now calculated based on U,V in streamice_vel_solve
393  C        but with Hybdrid stress formulation they will become part of  C        but with Hybdrid stress formulation they will become part of
394  C        velocity initial guess, so they are kept  C        velocity initial guess, so they are kept
395    
396    #ifdef ALLOW_PETSC
397          COMMON /STREAMICE_PETSC_DOFS_COMMON/
398         &      streamice_petsc_dofs_u,
399         &      streamice_petsc_dofs_v,
400         &      n_dofs_process
401          _RS streamice_petsc_dofs_u
402         & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
403          _RS streamice_petsc_dofs_v
404         & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
405          INTEGER n_dofs_process (0:nPx*nPy-1)
406    #endif
407                
408                
409  #endif /* ALLOW_STREAMICE */  #endif /* ALLOW_STREAMICE */
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22