1 |
cnh |
1.3 |
C $Header: find_rho.F,v 1.2 1998/04/22 19:25:35 cnh Exp $ |
2 |
cnh |
1.1 |
|
3 |
|
|
#include "CPP_EEOPTIONS.h" |
4 |
|
|
|
5 |
|
|
! ============================================================================== |
6 |
|
|
subroutine FIND_RHO( bi, bj, iMin, iMax, jMin, jMax, myThid ) |
7 |
|
|
implicit none |
8 |
|
|
! Arguments |
9 |
|
|
integer bi,bj,iMin,iMax,jMin,jMax |
10 |
|
|
integer myThid |
11 |
|
|
! Common |
12 |
|
|
#include "SIZE.h" |
13 |
|
|
#include "DYNVARS.h" |
14 |
|
|
#include "PARAMS.h" |
15 |
|
|
! Local |
16 |
cnh |
1.2 |
INTEGER i,j,k |
17 |
cnh |
1.1 |
! ------------------------------------------------------------------------------ |
18 |
|
|
|
19 |
|
|
do k=1,Nz |
20 |
|
|
do j=jMin,jMax |
21 |
|
|
do i=iMin,iMax |
22 |
|
|
rho(i,j,k,bi,bj)=rhonil*( |
23 |
|
|
& sBeta*( salt(i,j,k,bi,bj)-sRef(k)) |
24 |
|
|
& -tAlpha*(theta(i,j,k,bi,bj)-tRef(k)) |
25 |
|
|
& ) |
26 |
|
|
enddo |
27 |
|
|
enddo |
28 |
|
|
enddo |
29 |
|
|
|
30 |
|
|
! ------------------------------------------------------------------------------ |
31 |
|
|
return |
32 |
|
|
end |
33 |
|
|
! ============================================================================== |