/[MITgcm]/MITgcm/pkg/generic_advdiff/gad_som_exchanges.F
ViewVC logotype

Contents of /MITgcm/pkg/generic_advdiff/gad_som_exchanges.F

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


Revision 1.1 - (show annotations) (download)
Tue Jan 16 04:38:34 2007 UTC (17 years, 4 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint58w_post, checkpoint58x_post, checkpoint59e, checkpoint59d, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59, checkpoint58y_post, checkpoint58v_post
2nd-Order Moment Advection Scheme (Prather, 1986): first check-in
 - enable by setting #define GAD_ALLOW_SOM_ADVECT (in GAD_OPTIONS.h)
 - used without limiter (AdvScheme=80) or with Prather limiter (AdvScheme=81)
 - still needs work (not working with some options ; efficiency to improve)
   and serious testing.

1 C $Header: /u/gcmpack/MITgcm/pkg/gmredi/gmredi_init.F,v 1.10 2005/01/04 00:20:35 jmc Exp $
2 C $Name: $
3
4 #include "GAD_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: GAD_SOM_EXCHANGES
8 C !INTERFACE:
9 SUBROUTINE GAD_SOM_EXCHANGES( myThid )
10
11 C !DESCRIPTION: \bv
12 C *==========================================================*
13 C | SUBROUTINE GAD_SOM_EXCHANGES
14 C | o Apply exchanges to update overlaps
15 C | for 2nd.Order Moment fields
16 C *==========================================================*
17 C \ev
18 C !USES:
19 IMPLICIT NONE
20
21 C === Global variables ===
22 #include "SIZE.h"
23 #include "EEPARAMS.h"
24 #include "GAD.h"
25 #include "GAD_SOM_VARS.h"
26
27 C !INPUT/OUTPUT PARAMETERS:
28 C === Routine arguments ===
29 C myThid :: my Thread Id number
30 INTEGER myThid
31 CEOP
32
33 #ifdef GAD_ALLOW_SOM_ADVECT
34 C !LOCAL VARIABLES:
35 C === Local variables ===
36 INTEGER n
37
38 C-- Apply exchanges to Temp. 2nd.O.Moments:
39 IF ( tempSOM_Advection ) THEN
40 DO n=1,nSOM
41 CALL EXCH_3D_RL( som_T(1-OLx,1-OLy,1,1,1,n), Nr, myThid )
42 ENDDO
43 ENDIF
44
45 C-- Apply exchanges to Salin. 2nd.O.Moments:
46 IF ( saltSOM_Advection ) THEN
47 DO n=1,nSOM
48 CALL EXCH_3D_RL( som_S(1-OLx,1-OLy,1,1,1,n), Nr, myThid )
49 ENDDO
50 ENDIF
51
52 #endif /* GAD_ALLOW_SOM_ADVECT */
53
54 RETURN
55 END

  ViewVC Help
Powered by ViewVC 1.1.22