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

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

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


Revision 1.3 - (show annotations) (download)
Wed Apr 6 18:20:22 2005 UTC (19 years, 2 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint57t_post, checkpoint57m_post, checkpoint57v_post, checkpoint57g_pre, checkpoint57f_post, checkpoint57s_post, checkpoint57j_post, checkpoint58b_post, checkpoint57g_post, checkpoint57h_pre, checkpoint57y_post, checkpoint57x_post, checkpoint57h_post, checkpoint57y_pre, checkpoint57o_post, checkpoint57r_post, checkpoint57k_post, checkpoint57i_post, checkpoint58, checkpoint57h_done, checkpoint57n_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint58a_post, checkpoint57q_post, checkpoint57z_post, checkpoint57l_post
Changes since 1.2: +2 -1 lines
include "GRID.h"

1 C $Header: /u/gcmpack/MITgcm/model/src/quasihydrostaticterms.F,v 1.2 2002/11/14 22:36:38 heimbach Exp $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5
6 SUBROUTINE QUASIHYDROSTATICTERMS(
7 I bi,bj,k,
8 O effectiveBuoy,
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 "EEPARAMS.h"
24 #include "PARAMS.h"
25 #include "GRID.h"
26 #include "DYNVARS.h"
27
28 C-- == Routine arguments ==
29 INTEGER bi,bj,k
30 _RL effectiveBuoy(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
31 INTEGER myThid
32
33 C-- == Local variables ==
34 INTEGER i,j
35 _RL scalingFactor
36
37 scalingFactor=rhoConst/gravity
38 scalingFactor=1./gravity
39
40 DO j=1-Oly,sNy+Oly
41 DO i=1-Olx,sNx+Olx-1
42 effectiveBuoy(i,j)=effectiveBuoy(i,j)
43 & -scalingFactor*
44 & fCoriCos(i,j,bi,bj)*
45 & 0.5*( uVel(i,j,k,bi,bj)+uVel(i+1,j,k,bi,bj) )
46 ENDDO
47 ENDDO
48
49 IF (useNHMTerms) THEN
50 DO j=1-Oly,sNy+Oly-1
51 DO i=1-Olx,sNx+Olx-1
52 effectiveBuoy(i,j)=effectiveBuoy(i,j)
53 & -scalingFactor*
54 & (
55 & 0.5*( uVel(i,j,k,bi,bj)**2+uVel(i+1,j,k,bi,bj)**2 )
56 & +0.5*( vVel(i,j,k,bi,bj)**2+vVel(i,j+1,k,bi,bj)**2 )
57 & )*recip_RSphere
58 ENDDO
59 ENDDO
60 ENDIF
61
62 RETURN
63 END

  ViewVC Help
Powered by ViewVC 1.1.22