111 |
|
|
112 |
IF ( startIceModel.LE.0 .AND. nIter0.NE.0 ) THEN |
IF ( startIceModel.LE.0 .AND. nIter0.NE.0 ) THEN |
113 |
C-- Read ice pickup fields |
C-- Read ice pickup fields |
114 |
_BARRIER |
_BARRIER |
115 |
prec=precFloat64 |
prec=precFloat64 |
116 |
CALL THSICE_READ_PICKUP( prec, nIter0, myThid ) |
CALL THSICE_READ_PICKUP( prec, nIter0, myThid ) |
117 |
|
|
118 |
ELSE |
ELSE |
119 |
C-- Read initial conditions: |
C-- Read initial conditions: |
120 |
_BEGIN_MASTER( myThid ) |
_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 |
137 |
IF ( thSIceTsurf_InitFile .NE. ' ' ) THEN |
IF ( thSIceTsurf_InitFile .NE. ' ' ) THEN |
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 |
|
_END_MASTER(myThid) |
|
140 |
IF ( thSIceEnthp_InitFile .EQ. ' ' ) THEN |
IF ( thSIceEnthp_InitFile .EQ. ' ' ) THEN |
141 |
|
_BARRIER |
142 |
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 |
143 |
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) |
144 |
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 |