/[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.9 by mlosch, Tue May 30 22:48:59 2006 UTC revision 1.10 by mlosch, Mon Jun 5 22:32:53 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 .OR. Tice(2).GT.0. _d 0) THEN        IF (Tice(1).GT.0. _d 0 ) THEN
213            WRITE (6,*) 'BBerr Tice(1) > 0 = ',Tice(1)         WRITE (standardMessageUnit,*) 'BBerr: Tice(1) > 0'
214            WRITE (6,*) 'BBerr Tice(2) > 0 = ',Tice(2)         WRITE (standardMessageUnit,'(A,4I5,2F11.4)')
215         &      'BBerr: i,j,bi,bj,Tice = ',i,j,bi,bj,Tice
216          ENDIF
217          IF ( Tice(2).GT.0. _d 0) THEN
218           WRITE (standardMessageUnit,*) 'BBerr: Tice(2) > 0'
219           WRITE (standardMessageUnit,'(A,4I5,2F11.4)')
220         &      'BBerr: i,j,bi,bj,Tice = ',i,j,bi,bj,Tice
221        ENDIF        ENDIF
222  #ifdef ALLOW_DBUG_THSICE  #ifdef ALLOW_DBUG_THSICE
223        IF ( dBug(i,j,bi,bj) ) WRITE(6,1010)        IF ( dBug(i,j,bi,bj) ) WRITE(6,1010)
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22