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

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

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


Revision 1.4 - (hide annotations) (download)
Fri Jul 27 21:07:13 2012 UTC (12 years, 11 months ago) by heimbach
Branch: MAIN
Changes since 1.3: +2 -1 lines
File MIME type: text/plain
Enable transient I/O

1 heimbach 1.1 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
2    
3     #ifdef ALLOW_STREAMICE
4    
5     C -------------------------- REAL PARAMS ---------------------------------------------------
6    
7     COMMON /STREAMICE_PARMS_R/
8     & streamice_density, streamice_density_ocean_avg,
9     & A_glen_isothermal, n_glen, eps_glen_min,
10     & C_basal_fric_const, n_basal_friction, streamice_input_flux_unif,
11     & streamice_vel_update, streamice_cg_tol, streamice_nonlin_tol,
12 heimbach 1.4 & streamice_CFL_factor, streamice_adjDump
13 heimbach 1.1 _RL streamice_density, streamice_density_ocean_avg
14     _RL A_glen_isothermal, n_glen, eps_glen_min
15     _RL C_basal_fric_const
16     _RL n_basal_friction
17     _RL streamice_input_flux_unif
18     _RL streamice_vel_update
19     _RL streamice_cg_tol, streamice_nonlin_tol
20     _RL streamice_CFL_factor
21 heimbach 1.4 _RL streamice_adjDump
22 heimbach 1.1
23     C parms for parameterized initial thickness
24     C SHELF_MAX_DRAFT: max thickness of ice in m
25     C SHELF_MIN_DRAFT: min thickness of ice in m
26     C SHELF_EDGE_POS: extent of ice shelf in (km?)
27     C SHELF_SLOPE_SCALE: dist over which shelf slopes (km?)
28     C SHELF_FLAT_WIDTH: width of flat shelf (km?)
29     C also must be aware of units (m for cartesian, deg for curvilinear, m/deg for CYLINDRICAL POLAR)
30     C FLOW_DIR: 1.0=west, 2.0=east, 3.0=south, 4.0=north
31    
32     COMMON /STREAMICE_H_INIT_R/
33     & shelf_max_draft,
34     & shelf_min_draft,
35     & shelf_edge_pos,
36     & shelf_slope_scale,
37     & shelf_flat_width,
38     & flow_dir
39     _RL shelf_max_draft
40     _RL shelf_min_draft
41     _RL shelf_edge_pos
42     _RL shelf_slope_scale
43     _RL shelf_flat_width
44     _RL flow_dir
45    
46     C -------------------------- INT PARAMS ---------------------------------------------------
47    
48 heimbach 1.3 INTEGER streamice_max_nl
49     PARAMETER ( streamice_max_nl = 100 )
50    
51 heimbach 1.1 COMMON /STREAMICE_PARMS_I/
52     & streamice_max_cg_iter, streamice_max_nl_iter,
53     & streamice_vel_upd_counter, streamice_nstep_velocity,
54     & streamice_n_sub_regularize
55     INTEGER streamice_max_cg_iter, streamice_max_nl_iter
56     INTEGER streamice_vel_upd_counter, streamice_nstep_velocity
57     INTEGER streamice_n_sub_regularize
58    
59     C -------------------------- CHAR PARAMS ---------------------------------------------------
60    
61     CHARACTER*(MAX_LEN_FNAM) STREAMICEthickFile
62     CHARACTER*(MAX_LEN_FNAM) STREAMICEthickInit
63     CHARACTER*(MAX_LEN_FNAM) STREAMICEcalveMaskFile
64     COMMON /STREAMICE_PARM_C/
65     & STREAMICEthickInit,
66     & STREAMICEthickFile,
67     & STREAMICEcalveMaskFile
68    
69     C -------------------------- LOGICAL PARAMS ---------------------------------------------------
70    
71     LOGICAL STREAMICEison
72     LOGICAL STREAMICE_dump_mdsio
73     LOGICAL STREAMICE_tave_mdsio
74     LOGICAL STREAMICE_dump_mnc
75     LOGICAL STREAMICE_tave_mnc
76     LOGICAL STREAMICE_GL_regularize, STREAMICE_move_front
77     LOGICAL STREAMICE_calve_to_mask
78 heimbach 1.3 LOGICAL STREAMICE_construct_matrix
79 heimbach 1.1 LOGICAL STREAMICE_lower_cg_tol
80     COMMON /STREAMICE_PARM_L/
81     & STREAMICEison,
82     & STREAMICE_dump_mdsio, STREAMICE_tave_mdsio,
83     & STREAMICE_dump_mnc, STREAMICE_tave_mnc,
84     & STREAMICE_GL_regularize, STREAMICE_move_front,
85     & STREAMICE_calve_to_mask,
86 heimbach 1.3 & STREAMICE_construct_matrix,
87 heimbach 1.1 & STREAMICE_lower_cg_tol
88    
89     C -------------------------- AND NOW ARRAYS ---------------------------------------------------
90    
91     C EXPLANATION OF MASKS
92    
93     C STREAMICE_hmask VALUES 1=ice-covered cell
94     C 2=partially ice-covered cell (no dynamics)
95     C 0=ice-free cell (for now, it means the cell
96     C is treated as an open-ocean cell
97     C that ice shelf can flow into)
98     C -1=outside computational domain; will not change
99     C
100     C STREAMICE_umask VALUES 1=degree of freedom;
101     C 0=homogeneous dirich condition
102     C 3=inhomogeneous dirich condition
103    
104     C STREAMICE_vmask similar to umask
105    
106     C STREAMICE_ufacemask VALUES -1=unset,
107     C 0=no-flow boundary,
108     C 1=no-stress bdry
109     C 2=stress bdry condition,
110     C 3=inhomogeneous dirichlet boundary,
111     C 4=flux boundary: at these faces a flux will be specified by u_flux_bdry_SI
112     C
113     C STREAMICE_vfacemask similar to ufacemask
114     C STREAMICE_ufacemask_bdry field initialized at the beginning of simulation
115     C specified all ufacemask values except for calving front
116     C CONSTANT FOR A SIMULATION (ie not changes after streamice_init_fixed)
117     C STREAMICE_vfacemask_bdry CONSTANT FOR A SIMULATION
118     C STREAMICE_calve_mask specified allowed extent of ice shelf
119     C (should be integer, but don't know howvi sre to read ints from file)
120     C not necessarily used, but CONSTANT FOR A SIMULATION
121     C STREAMICE_float_cond will only be used if partial floatation is implemented
122    
123     C Short arrays (e.g. masks)
124     COMMON /STREAMICE_FIELDS_RS/
125     & STREAMICE_hmask,
126     & STREAMICE_umask,
127     & STREAMICE_vmask,
128     & STREAMICE_ufacemask,
129     & STREAMICE_vfacemask,
130     & STREAMICE_ufacemask_bdry,
131     & STREAMICE_vfacemask_bdry,
132     & STREAMICE_float_cond,
133     & STREAMICE_calve_mask
134     _RS STREAMICE_hmask (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
135     _RS STREAMICE_umask (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
136     _RS STREAMICE_vmask (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
137     _RS STREAMICE_ufacemask
138     & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
139     _RS STREAMICE_vfacemask
140     & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
141     _RS STREAMICE_ufacemask_bdry
142     & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
143     _RS STREAMICE_vfacemask_bdry
144     & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
145     _RS STREAMICE_float_cond
146     & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
147     _RS STREAMICE_calve_mask
148     & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
149    
150     C NOTES :
151     C STREAMICE_ufacemask_bdry, STREAMICE_vfacemask_bdry ARE CONSTANT (FIXED)
152     C STREAMICE_hmask IS PART OF **STATE**
153     C All other masks are updated within a timestep BASED ON STREAMICE_hmask
154    
155     C Number of quadrature points are hardcoded.. could turn into a CPP macro
156    
157     C REAL ARRAYS
158    
159     COMMON /STREAMICE_FIELDS_RL/
160     & H_streamice,
161     & U_streamice,
162     & V_streamice,
163     & visc_streamice,
164     & tau_beta_eff_streamice,
165     & float_frac_streamice,
166     & base_el_streamice,
167     & surf_el_streamice,
168     & area_shelf_streamice,
169     & mass_ice_streamice,
170     & u_flux_bdry_SI,
171     & v_flux_bdry_SI,
172     & h_bdry_values_SI,
173     & u_bdry_values_SI,
174     & v_bdry_values_SI,
175     & STREAMICE_dummy_array,
176     & C_basal_friction
177     _RL H_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
178     _RL U_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
179     _RL V_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
180     _RL visc_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
181     _RL tau_beta_eff_streamice(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
182     _RL float_frac_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
183     _RL surf_el_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
184     _RL base_el_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
185     _RL area_shelf_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
186     _RL mass_ice_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
187     _RL u_flux_bdry_SI (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
188     _RL v_flux_bdry_SI (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
189     _RL h_bdry_values_SI (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
190     _RL u_bdry_values_SI (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
191     _RL v_bdry_values_SI (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
192     _RL C_basal_friction (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
193    
194     _RL STREAMICE_dummy_array (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
195    
196     C NOTES :
197     C REAL ARRAYS THAT COMPRISE "STATE":
198     C H_streamice,
199     C U_streamice,
200     C V_streamice,
201     C visc_streamice,
202     C tau_beta_eff_streamice,
203     C area_shelf_streamice
204     C (and don't forget STREAMICE_hmask)
205     C
206     C visc & tau are now calculated based on U,V in streamice_vel_solve
207     C but with Hybdrid stress formulation they will become part of
208     C velocity initial guess, so they are kept
209    
210    
211     #endif /* ALLOW_STREAMICE */
212    
213     CEH3 ;;; Local Variables: ***
214     CEH3 ;;; mode:fortran ***
215     CEH3 ;;; End: ***

  ViewVC Help
Powered by ViewVC 1.1.22