1 |
jmc |
1.15 |
C $Header: /u/gcmpack/MITgcm/pkg/icefront/ICEFRONT.h,v 1.14 2016/09/15 00:16:03 jmc 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 |
jmc |
1.14 |
C \ev |
16 |
dimitri |
1.1 |
|
17 |
|
|
C----------------------------------------------------------------------- |
18 |
jmc |
1.14 |
C-- Constants that can be set in data.icefront, namelist /ICEFRONT_PARM01/ |
19 |
dimitri |
1.4 |
C ICEFRONTdepthFile - name of icefront depth file (m) |
20 |
|
|
C 2D file containing depth of the ice front |
21 |
|
|
C at each model grid cell |
22 |
|
|
C ICEFRONTlengthFile - name of icefront length file (m/m^2) |
23 |
|
|
C 2D file containing the ratio of the horizontal |
24 |
|
|
C length of the ice front in each model grid cell |
25 |
|
|
C divided by the grid cell area |
26 |
dimitri |
1.1 |
C ICEFRONTheatTransCoeff - heat transfer coefficient that determines |
27 |
|
|
C heat flux into icefront (m/s) |
28 |
|
|
C ICEFRONTsaltTransCoeff - salinity transfer coefficient that determines |
29 |
|
|
C salt flux into icefront (m/s) |
30 |
yunx |
1.8 |
C ICEFRONTthetaSurface - surface temperature on the top of icefront (oC) |
31 |
|
|
C interior temperture of the ice changes linearly |
32 |
|
|
C from ICEFRONTthetaSurface at surface to 0 oC at |
33 |
|
|
C the bottom |
34 |
dimitri |
1.1 |
C ICEFRONTlatentHeat - latent heat of fusion (J/kg) |
35 |
jmc |
1.12 |
C applyIcefrontTendT/S - |
36 |
dimitri |
1.1 |
C |
37 |
jmc |
1.14 |
C K_icefront - # of icefront model levels at every horizontal location (2D) |
38 |
|
|
C R_icefront - icefront depth [m] (2D) |
39 |
|
|
C icefrontlength - icefront horizontal length divided |
40 |
|
|
C 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 |
jmc |
1.14 |
C-- Constants that can be set in data.icefront, namelist /ICEFRONT_EXF_PARM02/ |
45 |
|
|
C SGRunOffFile - name of subglacial runoff file (kg/s) |
46 |
dimitri |
1.9 |
C 2D file containing mass of subglacial runoff |
47 |
|
|
C added at bottom of model ocean |
48 |
jmc |
1.14 |
C SGRunOff - subglacial runoff (kg/s) |
49 |
dimitri |
1.9 |
C Arrays *0 and *1 below are used for temporal interpolation. |
50 |
dimitri |
1.1 |
C----------------------------------------------------------------------- |
51 |
|
|
CEOP |
52 |
|
|
|
53 |
dimitri |
1.4 |
COMMON /ICEFRONT_PARMS_I/ K_icefront |
54 |
|
|
INTEGER K_icefront (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
55 |
|
|
|
56 |
jmc |
1.12 |
COMMON /ICEFRONT_PARMS_R/ |
57 |
dimitri |
1.1 |
& rhoIcefront, ICEFRONTkappa, |
58 |
|
|
& ICEFRONTlatentHeat, recip_ICEFRONTlatentHeat, |
59 |
|
|
& ICEFRONTheatCapacity_Cp, |
60 |
dimitri |
1.2 |
& ICEFRONTthetaSurface |
61 |
dimitri |
1.1 |
_RL ICEFRONTlatentHeat |
62 |
|
|
_RL ICEFRONTheatCapacity_Cp |
63 |
|
|
_RL rhoIcefront |
64 |
|
|
_RL ICEFRONTkappa |
65 |
|
|
_RL recip_ICEFRONTlatentHeat |
66 |
|
|
_RL ICEFRONTthetaSurface |
67 |
|
|
|
68 |
jmc |
1.12 |
COMMON /ICEFRONT_FIELDS_RL/ |
69 |
dimitri |
1.4 |
& icefront_TendT, |
70 |
|
|
& icefront_TendS |
71 |
dimitri |
1.6 |
_RL icefront_TendT (1:sNx,1:sNy,Nr,nSx,nSy) |
72 |
|
|
_RL icefront_TendS (1:sNx,1:sNy,Nr,nSx,nSy) |
73 |
dimitri |
1.1 |
|
74 |
jmc |
1.12 |
COMMON /ICEFRONT_FIELDS_RS/ |
75 |
dimitri |
1.1 |
& R_icefront, |
76 |
dimitri |
1.5 |
& icefrontlength |
77 |
dimitri |
1.6 |
_RS R_icefront (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
78 |
|
|
_RS icefrontlength (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
79 |
jmc |
1.12 |
|
80 |
dimitri |
1.1 |
LOGICAL ICEFRONTisOn |
81 |
dimitri |
1.4 |
LOGICAL applyIcefrontTendT |
82 |
|
|
LOGICAL applyIcefrontTendS |
83 |
dimitri |
1.1 |
COMMON /ICEFRONT_PARMS_L/ |
84 |
|
|
& ICEFRONTisOn, |
85 |
dimitri |
1.4 |
& applyIcefrontTendT, |
86 |
|
|
& applyIcefrontTendS |
87 |
dimitri |
1.1 |
|
88 |
dimitri |
1.4 |
CHARACTER*(MAX_LEN_FNAM) ICEFRONTlengthFile |
89 |
|
|
CHARACTER*(MAX_LEN_FNAM) ICEFRONTdepthFile |
90 |
jmc |
1.12 |
COMMON /ICEFRONT_PARM_C/ |
91 |
|
|
& ICEFRONTlengthFile, |
92 |
jmc |
1.14 |
& ICEFRONTdepthFile |
93 |
|
|
|
94 |
|
|
#ifdef ALLOW_EXF |
95 |
|
|
C the following variables are used in conjunction |
96 |
|
|
C with pkg/exf to specify sub-glacial runoff |
97 |
|
|
INTEGER SGRunOffstartdate1 |
98 |
|
|
INTEGER SGRunOffstartdate2 |
99 |
jmc |
1.15 |
_RL SGRunOffStartTime |
100 |
jmc |
1.14 |
_RL SGRunOffperiod |
101 |
|
|
_RL SGRunOffconst |
102 |
|
|
_RL SGRunOff_inscal |
103 |
|
|
_RL SGRunOff_remov_intercept |
104 |
|
|
_RL SGRunOff_remov_slope |
105 |
|
|
CHARACTER*(MAX_LEN_FNAM) SGRunOffFile |
106 |
|
|
CHARACTER*1 SGRunOffmask |
107 |
|
|
|
108 |
|
|
COMMON /ICEFRONT_EXF_PAR_I/ |
109 |
|
|
& SGRunOffstartdate1, SGRunOffstartdate2 |
110 |
|
|
COMMON /ICEFRONT_EXF_PAR_R/ |
111 |
jmc |
1.15 |
& SGRunOffStartTime, SGRunOffperiod, |
112 |
jmc |
1.14 |
& SGRunOffconst, SGRunOff_inscal, |
113 |
|
|
& SGRunOff_remov_intercept, SGRunOff_remov_slope |
114 |
|
|
COMMON /ICEFRONT_EXF_PAR_C/ |
115 |
|
|
& SGRunOffFile, SGRunOffmask |
116 |
|
|
#endif /* ALLOW_EXF */ |
117 |
dimitri |
1.1 |
|
118 |
|
|
#endif /* ALLOW_ICEFRONT */ |