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

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

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

revision 1.7 by adcroft, Thu Aug 15 18:48:47 2002 UTC revision 1.8 by mlosch, Mon Aug 19 14:21:30 2002 UTC
# Line 119  C     abbreviations Line 119  C     abbreviations
119                 t2 = t*t                 t2 = t*t
120                 t3 = t2*t                 t3 = t2*t
121                                
122                 if ( salt(i,j,k,bi,bj) .lt. 0. _d 0 ) then                 s  = salt(i,j,k,bi,bj)
123                   if ( s .lt. 0. _d 0 ) then
124  C     issue a warning  C     issue a warning
125                    write(*,'(a,i3,a,i3,a,i3,a,e13.5)')                    write(*,'(a,i3,a,i3,a,i3,a,e13.5)')
126       &             ' FIND_BETA: WARNING, s(',i,',',j,',',k,') = ', s       &             ' FIND_ALPHA: WARNING, salinity(',
127         &                 i,',',j,',',k,') = ', s
128                    s = 0. _d 0                    s = 0. _d 0
                else  
                   s  = salt(i,j,k,bi,bj)  
129                 end if                 end if
130                 s3o2 = sqrt(s*s*s)                 s3o2 = sqrt(s*s*s)
131                                
# Line 310  C     abbreviations Line 310  C     abbreviations
310                 t2 = t*t                 t2 = t*t
311                 t3 = t2*t                 t3 = t2*t
312                                
313                 if ( salt(i,j,k,bi,bj) .lt. 0. _d 0 ) then                 s  = salt(i,j,k,bi,bj)
314                   if ( s .lt. 0. _d 0 ) then
315  C     issue a warning  C     issue a warning
316                    write(*,'(a,i3,a,i3,a,i3,a,e13.5)')                    write(*,'(a,i3,a,i3,a,i3,a,e13.5)')
317       &             ' FIND_BETA: WARNING, s(',i,',',j,',',k,') = ', s       &             ' FIND_BETA: WARNING, salinity(',
318         &                 i,',',j,',',k,') = ', s
319                    s = 0. _d 0                    s = 0. _d 0
                else  
                   s  = salt(i,j,k,bi,bj)  
320                 end if                 end if
321                 s3o2 = 1.5*sqrt(s)                 s3o2 = 1.5*sqrt(s)
322    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22