49 |
INTEGER i, j |
INTEGER i, j |
50 |
#ifndef SEAICE_EXTERNAL_FLUXES |
#ifndef SEAICE_EXTERNAL_FLUXES |
51 |
INTEGER ITER |
INTEGER ITER |
52 |
_RL QS1, C1, C2, C3, C4, C5, TB, D1, D1W, D1I, D3 |
_RL QS1, TB, D1, D1W, D3, TMELT |
|
_RL TMELT, TMELTP, XKI, XKS, HCUT, ASNOW, XIO |
|
53 |
C effective conductivity of combined ice and snow |
C effective conductivity of combined ice and snow |
54 |
_RL effConduct |
_RL effConduct |
55 |
C specific humidity at ice surface |
C specific humidity at ice surface |
72 |
C NOW DEFINE ASSORTED CONSTANTS |
C NOW DEFINE ASSORTED CONSTANTS |
73 |
C SATURATION VAPOR PRESSURE CONSTANT |
C SATURATION VAPOR PRESSURE CONSTANT |
74 |
QS1=0.622 _d +00/1013.0 _d +00 |
QS1=0.622 _d +00/1013.0 _d +00 |
|
C MAYKUTS CONSTANTS FOR SAT. VAP. PRESSURE TEMP. POLYNOMIAL |
|
|
C1=2.7798202 _d -06 |
|
|
C2=-2.6913393 _d -03 |
|
|
C3=0.97920849 _d +00 |
|
|
C4=-158.63779 _d +00 |
|
|
C5=9653.1925 _d +00 |
|
75 |
C FREEZING TEMPERATURE OF SEAWATER |
C FREEZING TEMPERATURE OF SEAWATER |
76 |
TB=271.2 _d +00 |
TB=271.2 _d +00 |
77 |
C SENSIBLE HEAT CONSTANT |
C SENSIBLE HEAT CONSTANT |
78 |
D1=SEAICE_sensHeat |
D1=SEAICE_sensHeat |
79 |
C WATER LATENT HEAT CONSTANT |
C WATER LATENT HEAT CONSTANT |
80 |
D1W=SEAICE_latentWater |
D1W=SEAICE_latentWater |
|
C ICE LATENT HEAT CONSTANT |
|
|
D1I=SEAICE_latentIce |
|
81 |
C STEFAN BOLTZMAN CONSTANT TIMES 0.97 EMISSIVITY |
C STEFAN BOLTZMAN CONSTANT TIMES 0.97 EMISSIVITY |
82 |
D3=SEAICE_emissivity |
D3=SEAICE_emissivity |
83 |
C MELTING TEMPERATURE OF ICE |
C MELTING TEMPERATURE OF ICE |
84 |
TMELT=273.16 _d +00 |
TMELT=273.16 _d +00 |
|
TMELTP=273.159 _d +00 |
|
|
C ICE CONDUCTIVITY |
|
|
XKI=SEAICE_iceConduct |
|
|
C SNOW CONDUCTIVITY |
|
|
XKS=SEAICE_snowConduct |
|
|
C CUTOFF SNOW THICKNESS |
|
|
HCUT=SEAICE_snowThick |
|
|
C PENETRATION SHORTWAVE RADIATION FACTOR |
|
|
XIO=SEAICE_shortwave |
|
85 |
|
|
86 |
DO J=1,sNy |
DO J=1,sNy |
87 |
DO I=1,sNx |
DO I=1,sNx |