/[MITgcm]/MITgcm/pkg/seaice/seaice_fake.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/seaice_fake.F

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

revision 1.2 by jmc, Fri Aug 9 22:59:38 2013 UTC revision 1.3 by gforget, Thu Jan 9 19:36:27 2014 UTC
# Line 38  C     EndOfInterface(global-font-lock-mo Line 38  C     EndOfInterface(global-font-lock-mo
38  C     === Local variables ===  C     === Local variables ===
39  C     i,j,bi,bj - Loop counters  C     i,j,bi,bj - Loop counters
40    
41        INTEGER i, j, k, bi, bj        INTEGER i, j, bi, bj
42        _RL fac, tempFrz        _RL fac, tempFrz
43          CHARACTER*(MAX_LEN_MBUF) msgBuf
44    
45          WRITE(msgBuf,'(2A)') 'SEAICE_FAKE:',
46         &  ' forward code is not meant to be used (adj only)'
47          CALL PRINT_ERROR( msgBuf, myThid )
48          STOP 'ABNORMAL END: S/R SEAICE_FAKE'
49    
50        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
51         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
# Line 65  c relaxation to freezing point Line 71  c relaxation to freezing point
71              fac=MIN(1. _d 0, MAX(0. _d 0 , area(i,j,bi,bj)))              fac=MIN(1. _d 0, MAX(0. _d 0 , area(i,j,bi,bj)))
72              tempFrz = SEAICE_tempFrz0 +              tempFrz = SEAICE_tempFrz0 +
73       &                SEAICE_dTempFrz_dS *salt(i,j,1,bi,bj)       &                SEAICE_dTempFrz_dS *salt(i,j,1,bi,bj)
74              theta(i,j,k,bi,bj)=theta(i,j,k,bi,bj) + fac *              theta(i,j,1,bi,bj)=theta(i,j,1,bi,bj) + fac *
75       &        ( tempFrz-theta(i,j,k,bi,bj) ) *       &        ( tempFrz-theta(i,j,1,bi,bj) ) *
76       &        SEAICE_mcPheePiston/drF(1)*SEAICE_deltaTtherm       &        SEAICE_mcPheePiston/drF(1)*SEAICE_deltaTtherm
77           ENDDO           ENDDO
78          ENDDO          ENDDO
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22