pkg/icefront/ICEFRONT.h

C $Header: /u/gcmpack/MITgcm/pkg/icefront/ICEFRONT.h,v 1.5 2010/02/11 23:14:31 dimitri Exp $
C $Name:  $

#ifdef ALLOW_ICEFRONT

CBOP
C !ROUTINE: ICEFRONT.h

C !DESCRIPTION: \bv
C     /==========================================================\
C     | ICEFRONT.h                                               |
C     | o Basic header thermodnynamic shelf ice package.         |
C     |   Contains all ICEFRONT field declarations.              |
C     \==========================================================/

C-----------------------------------------------------------------------
C
C--   Constants that can be set in data.icefront
C     ICEFRONTdepthFile        - name of icefront depth file (m)
C                                2D file containing depth of the ice front
C                                at each model grid cell
C     ICEFRONTlengthFile       - name of icefront length file (m/m^2)
C                                2D file containing the ratio of the horizontal
C                                length of the ice front in each model grid cell
C                                divided by the grid cell area
C     ICEFRONTheatTransCoeff   - heat transfer coefficient that determines
C                                 heat flux into icefront (m/s)
C     ICEFRONTsaltTransCoeff   - salinity transfer coefficient that determines
C                                salt flux into icefront (m/s)
C     ICEFRONTlatentHeat       - latent heat of fusion (J/kg)
C     applyIcefrontTendT/S     -  
C
C--   Fields
C     K_icefront             - # of icefront model levels at every horizontal location (2D)
C     R_icefront             - icefront depth [m] (2D)
C     icefrontlength         - icefront horizontal length divided by grid cell area [m/m^2] (2D)
C     icefrontHeatFlux       - outward heat flux, +ve cools the ocean [W/m^2]
C     icefrontFreshWaterFlux - outward fresh water flux, +ve increases ocean salinity [m/s]
C     icefront_TendT         - temperature tendency (Kelvin/s)
C     icefront_TendS         - salinity tendency (psu/s)
C-----------------------------------------------------------------------
C \ev
CEOP

      COMMON /ICEFRONT_PARMS_I/  K_icefront
      INTEGER K_icefront (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)

      COMMON /ICEFRONT_PARMS_R/ 
     &     ICEFRONTheatTransCoeff, ICEFRONTsaltTransCoeff,
     &     rhoIcefront, ICEFRONTkappa,
     &     ICEFRONTlatentHeat, recip_ICEFRONTlatentHeat,
     &     ICEFRONTheatCapacity_Cp,
     &     ICEFRONTthetaSurface
      _RL ICEFRONTheatTransCoeff
      _RL ICEFRONTsaltTransCoeff
      _RL ICEFRONTlatentHeat
      _RL ICEFRONTheatCapacity_Cp
      _RL rhoIcefront
      _RL ICEFRONTkappa
      _RL recip_ICEFRONTlatentHeat
      _RL ICEFRONTthetaSurface

      COMMON /ICEFRONT_FIELDS_RL/ 
     &     icefront_TendT,
     &     icefront_TendS
      _RL icefront_TendT (1:sNx,1:sNy,Nr,nSx,nSy)
      _RL icefront_TendS (1:sNx,1:sNy,Nr,nSx,nSy)

      COMMON /ICEFRONT_FIELDS_RS/ 
     &     R_icefront,
     &     icefrontlength
      _RS R_icefront     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS icefrontlength (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      
      LOGICAL ICEFRONTisOn
      LOGICAL applyIcefrontTendT
      LOGICAL applyIcefrontTendS
      COMMON /ICEFRONT_PARMS_L/
     &     ICEFRONTisOn,
     &     applyIcefrontTendT,
     &     applyIcefrontTendS

      CHARACTER*(MAX_LEN_FNAM) ICEFRONTlengthFile
      CHARACTER*(MAX_LEN_FNAM) ICEFRONTdepthFile
      COMMON /ICEFRONT_PARM_C/ 
     &     ICEFRONTlengthFile, 
     &     ICEFRONTdepthFile

#endif /* ALLOW_ICEFRONT */
1	dimitri	1.6	C $Header: /u/gcmpack/MITgcm/pkg/icefront/ICEFRONT.h,v 1.5 2010/02/11 23:14:31 dimitri 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			C /==========================================================\
11			C \| ICEFRONT.h \|
12			C \| o Basic header thermodnynamic shelf ice package. \|
13			C \| Contains all ICEFRONT field declarations. \|
14			C \==========================================================/
15
16			C-----------------------------------------------------------------------
17			C
18			C-- Constants that can be set in data.icefront
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			C ICEFRONTlatentHeat - latent heat of fusion (J/kg)
31	dimitri	1.4	C applyIcefrontTendT/S -
32	dimitri	1.1	C
33			C-- Fields
34	dimitri	1.4	C K_icefront - # of icefront model levels at every horizontal location (2D)
35			C R_icefront - icefront depth [m] (2D)
36			C icefrontlength - icefront horizontal length divided by grid cell area [m/m^2] (2D)
37			C icefrontHeatFlux - outward heat flux, +ve cools the ocean [W/m^2]
38			C icefrontFreshWaterFlux - outward fresh water flux, +ve increases ocean salinity [m/s]
39			C icefront_TendT - temperature tendency (Kelvin/s)
40			C icefront_TendS - salinity tendency (psu/s)
41	dimitri	1.1	C-----------------------------------------------------------------------
42			C \ev
43			CEOP
44
45	dimitri	1.4	COMMON /ICEFRONT_PARMS_I/ K_icefront
46			INTEGER K_icefront (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
47
48	dimitri	1.1	COMMON /ICEFRONT_PARMS_R/
49			& ICEFRONTheatTransCoeff, ICEFRONTsaltTransCoeff,
50			& rhoIcefront, ICEFRONTkappa,
51			& ICEFRONTlatentHeat, recip_ICEFRONTlatentHeat,
52			& ICEFRONTheatCapacity_Cp,
53	dimitri	1.2	& ICEFRONTthetaSurface
54	dimitri	1.1	_RL ICEFRONTheatTransCoeff
55			_RL ICEFRONTsaltTransCoeff
56			_RL ICEFRONTlatentHeat
57			_RL ICEFRONTheatCapacity_Cp
58			_RL rhoIcefront
59			_RL ICEFRONTkappa
60			_RL recip_ICEFRONTlatentHeat
61			_RL ICEFRONTthetaSurface
62
63			COMMON /ICEFRONT_FIELDS_RL/
64	dimitri	1.4	& icefront_TendT,
65			& icefront_TendS
66	dimitri	1.6	_RL icefront_TendT (1:sNx,1:sNy,Nr,nSx,nSy)
67			_RL icefront_TendS (1:sNx,1:sNy,Nr,nSx,nSy)
68	dimitri	1.1
69			COMMON /ICEFRONT_FIELDS_RS/
70			& R_icefront,
71	dimitri	1.5	& icefrontlength
72	dimitri	1.6	_RS R_icefront (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
73			_RS icefrontlength (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
74	dimitri	1.1
75			LOGICAL ICEFRONTisOn
76	dimitri	1.4	LOGICAL applyIcefrontTendT
77			LOGICAL applyIcefrontTendS
78	dimitri	1.1	COMMON /ICEFRONT_PARMS_L/
79			& ICEFRONTisOn,
80	dimitri	1.4	& applyIcefrontTendT,
81			& applyIcefrontTendS
82	dimitri	1.1
83	dimitri	1.4	CHARACTER*(MAX_LEN_FNAM) ICEFRONTlengthFile
84			CHARACTER*(MAX_LEN_FNAM) ICEFRONTdepthFile
85	dimitri	1.1	COMMON /ICEFRONT_PARM_C/
86	dimitri	1.4	& ICEFRONTlengthFile,
87			& ICEFRONTdepthFile
88	dimitri	1.1
89			#endif /* ALLOW_ICEFRONT */