model/src/calc_buoyancy.F

C $Header: /u/gcmpack/models/MITgcmUV/model/src/calc_buoyancy.F,v 1.1 1998/08/23 15:35:10 cnh Exp $

#include "CPP_EEOPTIONS.h"

      SUBROUTINE CALC_BUOYANCY(
     I      bi, bj, iMin, iMax, jMin, jMax,  k, rhoLoc,
     O      buoy,
     I      myThid )

C     /==========================================================\
C     | o SUBROUTINE BUOYANCY                                    |
C     |   Calculates buoyancy for an XY slice of a tile.         |
C     |==========================================================|
C     |                                                          |
C     | k - is the density level                                 |
C     | kRef - determines pressure reference level               |
C     |        (not used in 'LINEAR' mode)                       |
C     |                                                          |
C     \==========================================================/
      IMPLICIT NONE

C--   == Global data ==
#include "SIZE.h"
#include "DYNVARS.h"
#include "EEPARAMS.h"
#include "PARAMS.h"

C--   == Routine arguments ==
      INTEGER bi,bj,iMin,iMax,jMin,jMax
      INTEGER k                 
      INTEGER kRef      
      _RL rhoLoc(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
      _RL buoy  (1-Olx:sNx+Olx,1-Oly:sNy+Oly)
      INTEGER myThid

C--   == Local variables ==
      INTEGER i,j
      _RL refTemp,refSalt,sigRef,tP,sP,deltaSig
! ------------------------------------------------------------------------------

      IF    ( buoyancyRelation .eq. 'ATMOSPHERIC'  ) THEN

       DO j=jMin,jMax
        DO i=iMin,iMax
         buoy(i,j)=-specVol_S(K)/theta_S(K)*(theta(i,j,k,bi,bj)-tRef(K))
        ENDDO
       ENDDO
       
      ELSEIF ( buoyancyRelation .eq. 'OCEANIC' ) THEN

       DO j=jMin,jMax
        DO i=iMin,iMax
         buoy(i,j)=-Gravity*recip_rhoConst*rholoc(i,j)
        ENDDO
       ENDDO

      ELSE

       STOP 'CALC_BUOANCY: variable "buoyancyRelation" has an illegal value'

      ENDIF


      RETURN
      END
1	C $Header: /u/gcmpack/models/MITgcmUV/model/src/calc_buoyancy.F,v 1.1 1998/08/23 15:35:10 cnh Exp $
2
3	#include "CPP_EEOPTIONS.h"
4
5	SUBROUTINE CALC_BUOYANCY(
6	I bi, bj, iMin, iMax, jMin, jMax, k, rhoLoc,
7	O buoy,
8	I myThid )
9
10	C /==========================================================\
11	C \| o SUBROUTINE BUOYANCY \|
12	C \| Calculates buoyancy for an XY slice of a tile. \|
13	C \|==========================================================\|
14	C \| \|
15	C \| k - is the density level \|
16	C \| kRef - determines pressure reference level \|
17	C \| (not used in 'LINEAR' mode) \|
18	C \| \|
19	C \==========================================================/
20	IMPLICIT NONE
21
22	C-- == Global data ==
23	#include "SIZE.h"
24	#include "DYNVARS.h"
25	#include "EEPARAMS.h"
26	#include "PARAMS.h"
27
28	C-- == Routine arguments ==
29	INTEGER bi,bj,iMin,iMax,jMin,jMax
30	INTEGER k
31	INTEGER kRef
32	_RL rhoLoc(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
33	_RL buoy (1-Olx:sNx+Olx,1-Oly:sNy+Oly)
34	INTEGER myThid
35
36	C-- == Local variables ==
37	INTEGER i,j
38	_RL refTemp,refSalt,sigRef,tP,sP,deltaSig
39	! ------------------------------------------------------------------------------
40
41	IF ( buoyancyRelation .eq. 'ATMOSPHERIC' ) THEN
42
43	DO j=jMin,jMax
44	DO i=iMin,iMax
45	buoy(i,j)=-specVol_S(K)/theta_S(K)*(theta(i,j,k,bi,bj)-tRef(K))
46	ENDDO
47	ENDDO
48
49	ELSEIF ( buoyancyRelation .eq. 'OCEANIC' ) THEN
50
51	DO j=jMin,jMax
52	DO i=iMin,iMax
53	buoy(i,j)=-Gravityrecip_rhoConstrholoc(i,j)
54	ENDDO
55	ENDDO
56
57	ELSE
58
59	STOP 'CALC_BUOANCY: variable "buoyancyRelation" has an illegal value'
60
61	ENDIF
62
63
64	RETURN
65	END