209 |
Tice(1) = 0.5 _d 0 *(-b1 - SQRT(b1*b1-4. _d 0*a1*c1))/a1 |
Tice(1) = 0.5 _d 0 *(-b1 - SQRT(b1*b1-4. _d 0*a1*c1))/a1 |
210 |
Tice(2) = (Lfresh-qicen(2)) / cpice |
Tice(2) = (Lfresh-qicen(2)) / cpice |
211 |
|
|
212 |
IF (Tice(1).GT.0. _d 0 ) THEN |
IF (Tice(1).GT.0. _d 0 ) THEN |
213 |
WRITE (standardMessageUnit,*) 'BBerr: Tice(1) > 0' |
WRITE (standardMessageUnit,'(A,I12,1PE14.6)') |
214 |
|
& ' BBerr: Tice(1) > 0 ; it=', myIter, qicen(1) |
215 |
WRITE (standardMessageUnit,'(A,4I5,2F11.4)') |
WRITE (standardMessageUnit,'(A,4I5,2F11.4)') |
216 |
& 'BBerr: i,j,bi,bj,Tice = ',i,j,bi,bj,Tice |
& ' BBerr: i,j,bi,bj,Tice = ',i,j,bi,bj,Tice |
217 |
ENDIF |
ENDIF |
218 |
IF ( Tice(2).GT.0. _d 0) THEN |
IF ( Tice(2).GT.0. _d 0) THEN |
219 |
WRITE (standardMessageUnit,*) 'BBerr: Tice(2) > 0' |
WRITE (standardMessageUnit,'(A,I12,1PE14.6)') |
220 |
|
& ' BBerr: Tice(2) > 0 ; it=', myIter, qicen(2) |
221 |
WRITE (standardMessageUnit,'(A,4I5,2F11.4)') |
WRITE (standardMessageUnit,'(A,4I5,2F11.4)') |
222 |
& 'BBerr: i,j,bi,bj,Tice = ',i,j,bi,bj,Tice |
& ' BBerr: i,j,bi,bj,Tice = ',i,j,bi,bj,Tice |
223 |
ENDIF |
ENDIF |
224 |
#ifdef ALLOW_DBUG_THSICE |
#ifdef ALLOW_DBUG_THSICE |
225 |
IF ( dBug(i,j,bi,bj) ) WRITE(6,1010) |
IF ( dBug(i,j,bi,bj) ) WRITE(6,1010) |
342 |
IF ( useBlkFlx .AND. k.EQ.nitMaxTsf |
IF ( useBlkFlx .AND. k.EQ.nitMaxTsf |
343 |
& .AND. ABS(dTsf).GE.Terrmax ) THEN |
& .AND. ABS(dTsf).GE.Terrmax ) THEN |
344 |
WRITE (6,'(A,4I4,I12,F15.9)') |
WRITE (6,'(A,4I4,I12,F15.9)') |
345 |
& ' BB: thermw conv err: i,j,it,hi=',i,j,bi,bj, |
& ' BB: not converge: i,j,it,hi=',i,j,bi,bj, |
346 |
& myIter,hi |
& myIter,hi |
347 |
WRITE (6,*) 'BB: thermw conv err: Tsf, dTsf=', Tsf,dTsf |
WRITE (6,*) 'BB: not converge: Tsf, dTsf=', Tsf,dTsf |
348 |
WRITE (6,*) 'BB: thermw conv err: flx0,dfdT=',flx0,df0dT |
WRITE (6,*) 'BB: not converge: flx0,dfdT=',flx0,df0dT |
349 |
IF (Tsf.LT.-70. _d 0) STOP |
IF (Tsf.LT.-70. _d 0) STOP |
350 |
ENDIF |
ENDIF |
351 |
|
|
402 |
|
|
403 |
IF (Tice(1) .GE. Tmlt1) THEN |
IF (Tice(1) .GE. Tmlt1) THEN |
404 |
WRITE (6,'(A,2I4,2I3,1P2E14.6)') |
WRITE (6,'(A,2I4,2I3,1P2E14.6)') |
405 |
& 'BBerr - Bug: IceT(1) > Tmlt',i,j,bi,bj,Tice(1),Tmlt1 |
& ' BBerr - Bug: IceT(1) > Tmlt',i,j,bi,bj,Tice(1),Tmlt1 |
406 |
ENDIF |
ENDIF |
407 |
IF (Tice(2) .GE. 0. _d 0) THEN |
IF (Tice(2) .GE. 0. _d 0) THEN |
408 |
WRITE (6,'(A,2I4,2I3,1P2E14.6)') |
WRITE (6,'(A,2I4,2I3,1P2E14.6)') |
409 |
& 'BBerr - Bug: IceT(2) > 0',i,j,bi,bj,Tice(2) |
& ' BBerr - Bug: IceT(2) > 0',i,j,bi,bj,Tice(2) |
410 |
ENDIF |
ENDIF |
411 |
|
|
412 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |