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 |
|
|
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 |
|
|