/[MITgcm]/MITgcm/model/src/ini_eos.F
ViewVC logotype

Diff of /MITgcm/model/src/ini_eos.F

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

revision 1.5 by mlosch, Wed Sep 18 16:38:02 2002 UTC revision 1.6 by mlosch, Wed Sep 25 19:36:50 2002 UTC
# Line 112  C Line 112  C
112                    do K=1,Nr                    do K=1,Nr
113                       do J=1-Oly,sNy+Oly                       do J=1-Oly,sNy+Oly
114                          do I=1-Olx,sNx+Olx                          do I=1-Olx,sNx+Olx
115                             pressure(i,j,k,bi,bj) = rhonil * (                             pressure(i,j,k,bi,bj) = rhoConst * (
116       &                          - gravity*rC(k)       &                          - gravity*rC(k)
117       &                          )       &                          )
118                          end do                          end do
# Line 430  C     check nonlinear EOS Line 430  C     check nonlinear EOS
430       &           SQUEEZE_RIGHT , 1)       &           SQUEEZE_RIGHT , 1)
431              write(msgBuf,'(a14,a22,f10.5,1x,f11.5)')              write(msgBuf,'(a14,a22,f10.5,1x,f11.5)')
432       &           'rho(find_rho) ',       &           'rho(find_rho) ',
433       &           ' = ', rholoc(i,j)+rhoNil, bulkMod(i,j)       &           ' = ', rholoc(i,j)+rhoConst, bulkMod(i,j)
434              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
435       &           SQUEEZE_RIGHT , 1)       &           SQUEEZE_RIGHT , 1)
436    
# Line 439  C     check nonlinear EOS Line 439  C     check nonlinear EOS
439              bulkMod(i,j) = 0. _d 0              bulkMod(i,j) = 0. _d 0
440              write(msgBuf,'(a21,a15,f10.5,1x,f11.5)')              write(msgBuf,'(a21,a15,f10.5,1x,f11.5)')
441       &           'rho(find_rho_scalar) ',       &           'rho(find_rho_scalar) ',
442       &           ' = ', rholoc(i,j)+rhoNil, bulkMod(i,j)       &           ' = ', rholoc(i,j)+rhoConst, bulkMod(i,j)
443              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
444       &           SQUEEZE_RIGHT , 1)       &           SQUEEZE_RIGHT , 1)
445                            
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22