339 |
#endif |
#endif |
340 |
IF ( useBlkFlx .AND. k.EQ.nitMaxTsf |
IF ( useBlkFlx .AND. k.EQ.nitMaxTsf |
341 |
& .AND. ABS(dTsf).GE.Terrmax ) THEN |
& .AND. ABS(dTsf).GE.Terrmax ) THEN |
342 |
WRITE (6,*) 'BB: thermw conv err ',i,j,bi,bj,dTsf |
WRITE (6,'(A,4I4,I12,F15.9)') |
343 |
WRITE (6,*) 'BB: thermw conv err, iceheight ', hi |
& ' BB: thermw conv err: i,j,it,hi=',i,j,bi,bj, |
344 |
WRITE (6,*) 'BB: thermw conv err: Tsf, flx0', Tsf,flx0 |
& myIter,hi |
345 |
|
WRITE (6,*) 'BB: thermw conv err: Tsf, dTsf=', Tsf,dTsf |
346 |
|
WRITE (6,*) 'BB: thermw conv err: flx0,dfdT=',flx0,df0dT |
347 |
IF (Tsf.LT.-70. _d 0) STOP |
IF (Tsf.LT.-70. _d 0) STOP |
348 |
ENDIF |
ENDIF |
349 |
|
|
|
100 continue ! surface temperature iteration |
|
350 |
ENDIF |
ENDIF |
351 |
ENDDO |
ENDDO |
|
150 continue |
|
352 |
C ------ end iteration ------------ |
C ------ end iteration ------------ |
353 |
|
|
354 |
C Compute new bottom layer temperature. |
C Compute new bottom layer temperature. |