/[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.4 by mlosch, Thu Sep 5 20:49:33 2002 UTC revision 1.5 by mlosch, Wed Sep 18 16:38:02 2002 UTC
# Line 202  C     5. secant bulk modulus K of sea wa Line 202  C     5. secant bulk modulus K of sea wa
202              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
203       &           SQUEEZE_RIGHT , 1)       &           SQUEEZE_RIGHT , 1)
204              write(msgBuf,'(a)')              write(msgBuf,'(a)')
205       &           'WARNING: (see Jackett and McDougall 1995, JPO)'       &           'WARNING: (see Jackett and McDougall 1995, JAOT)'
206              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
207       &           SQUEEZE_RIGHT , 1)       &           SQUEEZE_RIGHT , 1)
208              write(msgBuf,'(a)')              write(msgBuf,'(a)')
# Line 414  C     check nonlinear EOS Line 414  C     check nonlinear EOS
414                    
415              call find_rho(              call find_rho(
416       &           bi, bj, iMin, iMax, jMin, jMax,  k, k,       &           bi, bj, iMin, iMax, jMin, jMax,  k, k,
      &           equationOfState,  
417       &           tFld, sFld, rholoc, myThid )       &           tFld, sFld, rholoc, myThid )
418    
419              call find_bulkmod(              call find_bulkmod(
420       &           bi, bj, imin, imax, jmin, jmax, k, k,       &           bi, bj, imin, imax, jmin, jmax, k, k,
421       &           tFld, sFld, bulkMod, myThid )       &           tFld, sFld, bulkMod, myThid )
422                
         
423              write(msgBuf,              write(msgBuf,
424       &           '(a4,f4.1,a5,f4.1,a6,f5.0,a5,a3,f10.5,1x,f11.5)')       &           '(a4,f4.1,a5,f4.1,a6,f5.0,a5,a3,f10.5,1x,f11.5)')
425       &           'rho(', sFld(i,j,k,bi,bj), ' PSU,',       &           'rho(', sFld(i,j,k,bi,bj), ' PSU,',
# Line 430  C     check nonlinear EOS Line 428  C     check nonlinear EOS
428       &           rloc(kcheck), bloc(kcheck)       &           rloc(kcheck), bloc(kcheck)
429              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
430       &           SQUEEZE_RIGHT , 1)       &           SQUEEZE_RIGHT , 1)
431              write(msgBuf,'(a4,a32,f10.5,1x,f11.5)')              write(msgBuf,'(a14,a22,f10.5,1x,f11.5)')
432       &           'rho ', ' = ', rholoc(i,j)+rhoNil, bulkMod(i,j)       &           'rho(find_rho) ',
433         &           ' = ', rholoc(i,j)+rhoNil, bulkMod(i,j)
434                CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
435         &           SQUEEZE_RIGHT , 1)
436    
437                call find_rho_scalar( tFld(i,j,k,bi,bj), sLoc(kcheck),
438         &           pLoc(kcheck), rhoLoc(i,j), myThid )
439                bulkMod(i,j) = 0. _d 0
440                write(msgBuf,'(a21,a15,f10.5,1x,f11.5)')
441         &           'rho(find_rho_scalar) ',
442         &           ' = ', rholoc(i,j)+rhoNil, 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.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22