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

Contents of /MITgcm/model/src/calc_buoyancy.F

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

Revision 1.7 - (show annotations) (download)
Sun Feb 4 14:38:45 2001 UTC (23 years, 4 months ago) by cnh
Branch: MAIN
CVS Tags: checkpoint44e_post, checkpoint46g_pre, release1_p13_pre, checkpoint46f_post, checkpoint46b_post, checkpoint43a-release1mods, checkpoint44g_post, release1_p13, checkpoint40pre3, checkpoint40pre1, checkpoint40pre7, checkpoint40pre6, checkpoint40pre9, checkpoint40pre8, chkpt44d_post, release1_p8, release1_p9, release1_p1, release1_p2, release1_p3, release1_p4, release1_p5, release1_p6, release1_p7, checkpoint44e_pre, release1_b1, checkpoint43, checkpoint38, release1_chkpt44d_post, release1_p11, icebear5, icebear4, icebear3, icebear2, checkpoint46d_pre, checkpoint40pre2, release1-branch_tutorials, checkpoint45d_post, chkpt44a_post, checkpoint44h_pre, checkpoint40pre4, checkpoint46a_post, checkpoint46b_pre, chkpt44c_pre, checkpoint45a_post, ecco_c44_e19, ecco_c44_e18, ecco_c44_e17, ecco_c44_e16, release1_p12, release1_p10, release1_p16, release1_p17, release1_p14, release1_p15, pre38tag1, checkpoint46e_pre, checkpoint45b_post, release1-branch-end, c37_adj, release1_final_v1, release1_p12_pre, checkpoint46c_pre, checkpoint46, checkpoint44f_post, checkpoint44b_post, checkpoint46h_pre, checkpoint46a_pre, checkpoint45c_post, ecco_ice2, ecco_ice1, checkpoint44h_post, pre38-close, checkpoint46g_post, checkpoint39, checkpoint37, checkpoint36, checkpoint35, ecco_c44_e22, ecco_c44_e25, checkpoint40pre5, chkpt44a_pre, ecco_c44_e23, ecco_c44_e20, ecco_c44_e21, ecco_c44_e26, ecco_c44_e27, ecco_c44_e24, checkpoint46c_post, ecco-branch-mod1, ecco-branch-mod2, ecco-branch-mod3, ecco-branch-mod4, ecco-branch-mod5, checkpoint46e_post, release1_beta1, checkpoint44b_pre, checkpoint42, checkpoint40, checkpoint41, checkpoint44, checkpoint45, chkpt44c_post, checkpoint44f_pre, checkpoint46d_post, release1-branch_branchpoint
Branch point for: c24_e25_ice, release1_coupled, release1_final, release1-branch, release1, ecco-branch, release1_50yr, icebear, pre38
Changes since 1.6: +2 -1 lines 
Made sure each .F and .h file had
the CVS keywords Header and Name at its start.
Most had header but very few currently have Name, so
lots of changes!
1 C $Header: /u/gcmpack/models/MITgcmUV/model/src/calc_buoyancy.F,v 1.6 2000/03/27 22:25:43 adcroft Exp $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5
6 SUBROUTINE CALC_BUOYANCY(
7 I bi, bj, iMin, iMax, jMin, jMax, k, rhoLoc,
8 O buoy,
9 I myThid )
10
11 C /==========================================================\
12 C | o SUBROUTINE BUOYANCY |
13 C | Calculates buoyancy for an XY slice of a tile. |
14 C |==========================================================|
15 C | |
16 C | k - is the density level |
17 C | |
18 C \==========================================================/
19 IMPLICIT NONE
20
21 C-- == Global data ==
22 #include "SIZE.h"
23 #include "DYNVARS.h"
24 #include "EEPARAMS.h"
25 #include "PARAMS.h"
26
27 C-- == Routine arguments ==
28 INTEGER bi,bj,iMin,iMax,jMin,jMax
29 INTEGER k
30 _RL rhoLoc(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
31 _RL buoy (1-Olx:sNx+Olx,1-Oly:sNy+Oly)
32 INTEGER myThid
33
34 C-- == Local variables ==
35 INTEGER i,j
36
37 IF ( buoyancyRelation .eq. 'ATMOSPHERIC' ) THEN
38
39 DO j=jMin,jMax
40 DO i=iMin,iMax
41 buoy(i,j)=(theta(i,j,k,bi,bj)-tRef(K))/tRef(K)
42 ENDDO
43 ENDDO
44
45 ELSEIF ( buoyancyRelation .eq. 'OCEANIC' ) THEN
46
47 DO j=jMin,jMax
48 DO i=iMin,iMax
49 buoy(i,j)=-Gravity*recip_rhoConst*rholoc(i,j)
50 ENDDO
51 ENDDO
52
53 ELSE
54
55 STOP
56 &'CALC_BUOANCY: variable "buoyancyRelation" has an illegal value'
57
58 ENDIF
59
60
61 RETURN
62 END
  ViewVC Help
Powered by ViewVC 1.1.22