/[MITgcm]/MITgcm/pkg/thsice/thsice_solve4temp.F
ViewVC logotype

Diff of /MITgcm/pkg/thsice/thsice_solve4temp.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.11 by jmc, Tue Jun 6 22:10:52 2006 UTC revision 1.12 by jmc, Tue Jun 20 21:33:51 2006 UTC
# Line 209  C compute ice temperatures Line 209  C compute ice temperatures
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)
# Line 340  C If no convergence, then repeat. Line 342  C If no convergence, then repeat.
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    
# Line 400  C (This should not happen for reasonable Line 402  C (This should not happen for reasonable
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-|--+----|
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22