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)') |
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,', |
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 |
|
|