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

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

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

revision 1.16 by mlosch, Wed Aug 7 16:55:52 2002 UTC revision 1.17 by adcroft, Thu Aug 15 17:25:31 2002 UTC
# Line 207  C     abbreviations Line 207  C     abbreviations
207              t3 = t2*t              t3 = t2*t
208              t4 = t3*t              t4 = t3*t
209    
210              if ( sFld(i,j,k,bi,bj) .lt. 0. _d 0 ) then              s  = sFld(i,j,k,bi,bj)
211                if ( s .lt. 0. _d 0 ) then
212  C     issue a warning  C     issue a warning
213                 write(*,'(a,i3,a,i3,a,i3,a,e13.5)')                 write(*,'(a,i3,a,i3,a,i3,a,e13.5)')
214       &           ' FIND_RHOP0:   WARNING, s(',i,',',j,',',k,') = ', s       &           ' FIND_RHOP0:   WARNING, s(',i,',',j,',',k,') = ', s
215                 s = 0. _d 0                 s = 0. _d 0
             else  
                s  = sFld(i,j,k,bi,bj)  
216              end if              end if
217              s3o2 = s*sqrt(s)              s3o2 = s*sqrt(s)
218                            
# Line 300  C     abbreviations Line 299  C     abbreviations
299              t3 = t2*t              t3 = t2*t
300              t4 = t3*t              t4 = t3*t
301    
302              if ( sFld(i,j,k,bi,bj) .lt. 0. _d 0 ) then              s  = sFld(i,j,k,bi,bj)
303                if ( s .lt. 0. _d 0 ) then
304  C     issue a warning  C     issue a warning
305                 write(*,'(a,i3,a,i3,a,i3,a,e13.5)')                 write(*,'(a,i3,a,i3,a,i3,a,e13.5)')
306       &           ' FIND_BULKMOD: WARNING, s(',i,',',j,',',k,') = ', s       &           ' FIND_BULKMOD: WARNING, s(',i,',',j,',',k,') = ', s
307                 s = 0. _d 0                 s = 0. _d 0
             else  
                s  = sFld(i,j,k,bi,bj)  
308              end if              end if
309              s3o2 = s*sqrt(s)              s3o2 = s*sqrt(s)
310  C  C
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22