60 |
& fracEnMelt, fracEnFreez, hThinIce, hThickIce, hNewIceMax |
& fracEnMelt, fracEnFreez, hThinIce, hThickIce, hNewIceMax |
61 |
|
|
62 |
NAMELIST /THSICE_PARM01/ |
NAMELIST /THSICE_PARM01/ |
63 |
& startIceModel, stepFwd_oceMxL, |
& startIceModel, stepFwd_oceMxL, thSIce_calc_albNIR, |
64 |
& thSIce_deltaT, thSIce_dtTemp, |
& thSIce_deltaT, thSIce_dtTemp, |
65 |
& ocean_deltaT, tauRelax_MxL, |
& ocean_deltaT, tauRelax_MxL, |
66 |
& hMxL_default, sMxL_default, vMxL_default, |
& hMxL_default, sMxL_default, vMxL_default, |
127 |
|
|
128 |
C-- Default values (parameters) |
C-- Default values (parameters) |
129 |
stepFwd_oceMxL = .FALSE. |
stepFwd_oceMxL = .FALSE. |
130 |
|
thSIce_calc_albNIR = .FALSE. |
131 |
startIceModel = 0 |
startIceModel = 0 |
132 |
thSIce_deltaT = dTtracerLev(1) |
thSIce_deltaT = dTtracerLev(1) |
133 |
thSIce_dtTemp = UNSET_RL |
thSIce_dtTemp = UNSET_RL |
187 |
|
|
188 |
C- Define other constants (from previous ones): |
C- Define other constants (from previous ones): |
189 |
Tmlt1=-mu_Tf*S_winton |
Tmlt1=-mu_Tf*S_winton |
190 |
rhoiw = rhosw - rhoi |
floodFac = (rhosw - rhoi)/rhos |
191 |
|
|
192 |
C Set I/O parameters |
C Set I/O parameters |
193 |
thSIce_tave_mdsio = .TRUE. |
thSIce_tave_mdsio = .TRUE. |
251 |
WRITE(iUnit,*) 'ThSI: rhoi =',rhoi |
WRITE(iUnit,*) 'ThSI: rhoi =',rhoi |
252 |
WRITE(iUnit,*) 'ThSI: rhosw =',rhosw |
WRITE(iUnit,*) 'ThSI: rhosw =',rhosw |
253 |
WRITE(iUnit,*) 'ThSI: rhofw =',rhofw |
WRITE(iUnit,*) 'ThSI: rhofw =',rhofw |
254 |
WRITE(iUnit,*) 'ThSI: rhoiw =',rhoiw |
WRITE(iUnit,*) 'ThSI: floodFac =',floodFac |
255 |
WRITE(iUnit,*) 'ThSI: cpIce =',cpIce |
WRITE(iUnit,*) 'ThSI: cpIce =',cpIce |
256 |
WRITE(iUnit,*) 'ThSI: cpWater =',cpWater |
WRITE(iUnit,*) 'ThSI: cpWater =',cpWater |
257 |
WRITE(iUnit,*) 'ThSI: kIce =',kIce |
WRITE(iUnit,*) 'ThSI: kIce =',kIce |