109 |
|
|
110 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
111 |
|
|
112 |
IF ( startIceModel.LE.0 .AND. nIter0.NE.0 ) THEN |
IF ( startIceModel.LE.0 .AND. |
113 |
|
& ( nIter0.NE.0 .OR. pickupSuff.NE.' ' ) |
114 |
|
& ) THEN |
115 |
C-- Read ice pickup fields |
C-- Read ice pickup fields |
|
_BARRIER |
|
116 |
prec=precFloat64 |
prec=precFloat64 |
117 |
CALL THSICE_READ_PICKUP( prec, nIter0, myThid ) |
CALL THSICE_READ_PICKUP( prec, nIter0, myThid ) |
118 |
|
|
119 |
ELSE |
ELSE |
120 |
C-- Read initial conditions: |
C-- Read initial conditions: |
|
_BARRIER |
|
121 |
IF ( thSIceFract_InitFile .NE. ' ' ) THEN |
IF ( thSIceFract_InitFile .NE. ' ' ) THEN |
122 |
CALL READ_REC_XY_RL(thSIceFract_InitFile,iceMask,1,0,myThid) |
CALL READ_REC_XY_RL(thSIceFract_InitFile,iceMask,1,0,myThid) |
123 |
ENDIF |
ENDIF |
138 |
CALL READ_REC_XY_RS(thSIceTsurf_InitFile,Tsrf,1,0,myThid) |
CALL READ_REC_XY_RS(thSIceTsurf_InitFile,Tsrf,1,0,myThid) |
139 |
ENDIF |
ENDIF |
140 |
IF ( thSIceEnthp_InitFile .EQ. ' ' ) THEN |
IF ( thSIceEnthp_InitFile .EQ. ' ' ) THEN |
|
_BARRIER |
|
141 |
C- enthalpy of new ice in J/kg, taken from thsice_extend.F with Tf beeing |
C- enthalpy of new ice in J/kg, taken from thsice_extend.F with Tf beeing |
142 |
C the freezing Temp of seawater computed from a fixed salinity (31.5 psu) |
C the freezing Temp of seawater computed from a fixed salinity (31.5 psu) |
143 |
C Tf = -mu_Tf*salinity = -1.70 deg C; Qice1 ~ 3.2e5; Qice2 ~ 3.4e5 |
C Tf = -mu_Tf*salinity = -1.70 deg C; Qice1 ~ 3.2e5; Qice2 ~ 3.4e5 |