1 |
jmc |
1.12 |
C $Header: /u/gcmpack/MITgcm/pkg/icefront/ICEFRONT.h,v 1.11 2011/05/16 22:41:03 yunx Exp $ |
2 |
dimitri |
1.1 |
C $Name: $ |
3 |
|
|
|
4 |
|
|
#ifdef ALLOW_ICEFRONT |
5 |
|
|
|
6 |
|
|
CBOP |
7 |
|
|
C !ROUTINE: ICEFRONT.h |
8 |
|
|
|
9 |
|
|
C !DESCRIPTION: \bv |
10 |
jmc |
1.12 |
C *==========================================================* |
11 |
dimitri |
1.1 |
C | ICEFRONT.h | |
12 |
|
|
C | o Basic header thermodnynamic shelf ice package. | |
13 |
|
|
C | Contains all ICEFRONT field declarations. | |
14 |
jmc |
1.12 |
C *==========================================================* |
15 |
dimitri |
1.1 |
|
16 |
|
|
C----------------------------------------------------------------------- |
17 |
|
|
C |
18 |
|
|
C-- Constants that can be set in data.icefront |
19 |
yunx |
1.8 |
C- Namelist /ICEFRONT_PARM01/ |
20 |
dimitri |
1.4 |
C ICEFRONTdepthFile - name of icefront depth file (m) |
21 |
|
|
C 2D file containing depth of the ice front |
22 |
|
|
C at each model grid cell |
23 |
|
|
C ICEFRONTlengthFile - name of icefront length file (m/m^2) |
24 |
|
|
C 2D file containing the ratio of the horizontal |
25 |
|
|
C length of the ice front in each model grid cell |
26 |
|
|
C divided by the grid cell area |
27 |
dimitri |
1.1 |
C ICEFRONTheatTransCoeff - heat transfer coefficient that determines |
28 |
|
|
C heat flux into icefront (m/s) |
29 |
|
|
C ICEFRONTsaltTransCoeff - salinity transfer coefficient that determines |
30 |
|
|
C salt flux into icefront (m/s) |
31 |
yunx |
1.8 |
C ICEFRONTthetaSurface - surface temperature on the top of icefront (oC) |
32 |
|
|
C interior temperture of the ice changes linearly |
33 |
|
|
C from ICEFRONTthetaSurface at surface to 0 oC at |
34 |
|
|
C the bottom |
35 |
dimitri |
1.1 |
C ICEFRONTlatentHeat - latent heat of fusion (J/kg) |
36 |
jmc |
1.12 |
C applyIcefrontTendT/S - |
37 |
dimitri |
1.1 |
C |
38 |
dimitri |
1.4 |
C K_icefront - # of icefront model levels at every horizontal location (2D) |
39 |
|
|
C R_icefront - icefront depth [m] (2D) |
40 |
|
|
C icefrontlength - icefront horizontal length divided by grid cell area [m/m^2] (2D) |
41 |
|
|
C icefront_TendT - temperature tendency (Kelvin/s) |
42 |
|
|
C icefront_TendS - salinity tendency (psu/s) |
43 |
dimitri |
1.9 |
C |
44 |
|
|
C SGrunoffFile - name of subglacial runoff file (kg/s) |
45 |
|
|
C 2D file containing mass of subglacial runoff |
46 |
|
|
C added at bottom of model ocean |
47 |
|
|
C SGrunoff - subglacial runoff (kg/s) |
48 |
|
|
C Arrays *0 and *1 below are used for temporal interpolation. |
49 |
|
|
C |
50 |
dimitri |
1.1 |
C----------------------------------------------------------------------- |
51 |
|
|
C \ev |
52 |
|
|
CEOP |
53 |
|
|
|
54 |
dimitri |
1.4 |
COMMON /ICEFRONT_PARMS_I/ K_icefront |
55 |
|
|
INTEGER K_icefront (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
56 |
|
|
|
57 |
jmc |
1.12 |
COMMON /ICEFRONT_PARMS_R/ |
58 |
dimitri |
1.1 |
& rhoIcefront, ICEFRONTkappa, |
59 |
|
|
& ICEFRONTlatentHeat, recip_ICEFRONTlatentHeat, |
60 |
|
|
& ICEFRONTheatCapacity_Cp, |
61 |
dimitri |
1.2 |
& ICEFRONTthetaSurface |
62 |
dimitri |
1.1 |
_RL ICEFRONTlatentHeat |
63 |
|
|
_RL ICEFRONTheatCapacity_Cp |
64 |
|
|
_RL rhoIcefront |
65 |
|
|
_RL ICEFRONTkappa |
66 |
|
|
_RL recip_ICEFRONTlatentHeat |
67 |
|
|
_RL ICEFRONTthetaSurface |
68 |
|
|
|
69 |
jmc |
1.12 |
COMMON /ICEFRONT_FIELDS_RL/ |
70 |
dimitri |
1.4 |
& icefront_TendT, |
71 |
|
|
& icefront_TendS |
72 |
dimitri |
1.6 |
_RL icefront_TendT (1:sNx,1:sNy,Nr,nSx,nSy) |
73 |
|
|
_RL icefront_TendS (1:sNx,1:sNy,Nr,nSx,nSy) |
74 |
dimitri |
1.1 |
|
75 |
jmc |
1.12 |
COMMON /ICEFRONT_FIELDS_RS/ |
76 |
dimitri |
1.1 |
& R_icefront, |
77 |
dimitri |
1.5 |
& icefrontlength |
78 |
dimitri |
1.6 |
_RS R_icefront (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
79 |
|
|
_RS icefrontlength (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
80 |
jmc |
1.12 |
|
81 |
dimitri |
1.1 |
LOGICAL ICEFRONTisOn |
82 |
dimitri |
1.4 |
LOGICAL applyIcefrontTendT |
83 |
|
|
LOGICAL applyIcefrontTendS |
84 |
dimitri |
1.1 |
COMMON /ICEFRONT_PARMS_L/ |
85 |
|
|
& ICEFRONTisOn, |
86 |
dimitri |
1.4 |
& applyIcefrontTendT, |
87 |
|
|
& applyIcefrontTendS |
88 |
dimitri |
1.1 |
|
89 |
dimitri |
1.4 |
CHARACTER*(MAX_LEN_FNAM) ICEFRONTlengthFile |
90 |
|
|
CHARACTER*(MAX_LEN_FNAM) ICEFRONTdepthFile |
91 |
jmc |
1.12 |
CHARACTER*(MAX_LEN_FNAM) SGrunoffFile |
92 |
|
|
COMMON /ICEFRONT_PARM_C/ |
93 |
|
|
& ICEFRONTlengthFile, |
94 |
|
|
& ICEFRONTdepthFile, |
95 |
|
|
& SGrunoffFile |
96 |
dimitri |
1.1 |
|
97 |
yunx |
1.8 |
#ifdef ALLOW_SUBGLACIAL_RUNOFF |
98 |
dimitri |
1.9 |
_RL SGrunoff (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
99 |
|
|
_RL SGrunoff0(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
100 |
|
|
_RL SGrunoff1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
101 |
dimitri |
1.10 |
COMMON /SUBGLACIAL_RUNOFF/ |
102 |
|
|
& SGrunoff, SGrunoff0,SGrunoff1 |
103 |
yunx |
1.8 |
#endif /* ALLOW_SUBGLACIAL_RUNOFF */ |
104 |
|
|
|
105 |
dimitri |
1.1 |
#endif /* ALLOW_ICEFRONT */ |