/[MITgcm]/MITgcm/pkg/ptracers/ptracers_init_fixed.F
ViewVC logotype

Diff of /MITgcm/pkg/ptracers/ptracers_init_fixed.F

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

revision 1.2 by jmc, Mon Nov 5 18:48:04 2007 UTC revision 1.3 by jmc, Mon Dec 17 22:03:15 2007 UTC
# Line 20  C     !USES: Line 20  C     !USES:
20  #include "PARAMS.h"  #include "PARAMS.h"
21  #include "PTRACERS_SIZE.h"  #include "PTRACERS_SIZE.h"
22  #include "PTRACERS_PARAMS.h"  #include "PTRACERS_PARAMS.h"
23    #include "PTRACERS_RESTART.h"
24  #include "GAD.h"  #include "GAD.h"
25    
26  C     !INPUT PARAMETERS:  C     !INPUT PARAMETERS:
# Line 35  C     msgBuf      - Informational/error Line 36  C     msgBuf      - Informational/error
36        INTEGER minOlSize        INTEGER minOlSize
37        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
38    
39          _BEGIN_MASTER( myThid )
40    
41    C     Initialise internal parameter in common block:
42          DO iTracer = 1, PTRACERS_num
43            PTRACERS_MultiDimAdv(iTracer) = multiDimAdvection
44            PTRACERS_AdamsBashGtr(iTracer) = .FALSE.
45            PTRACERS_startAB(iTracer) = nIter0 - PTRACERS_Iter0
46          ENDDO
47    
48  C     Loop over tracers  C     Loop over tracers
49        DO iTracer = 1, PTRACERS_numInUse        DO iTracer = 1, PTRACERS_numInUse
50    
51          IF (          IF (
52       &       multiDimAdvection .AND.       &       PTRACERS_advScheme(iTracer).EQ.ENUM_CENTERED_2ND .OR.
53       &       PTRACERS_advScheme(iTracer).NE.ENUM_CENTERED_2ND .AND.       &       PTRACERS_advScheme(iTracer).EQ.ENUM_UPWIND_3RD .OR.
54       &       PTRACERS_advScheme(iTracer).NE.ENUM_UPWIND_3RD .AND.       &       PTRACERS_advScheme(iTracer).EQ.ENUM_CENTERED_4TH
55       &       PTRACERS_advScheme(iTracer).NE.ENUM_CENTERED_4TH ) THEN       &     ) PTRACERS_MultiDimAdv(iTracer) = .FALSE.
56  c         PTRACERS_multiDimAdv(iTracer) = .TRUE.          useMultiDimAdvec = useMultiDimAdvec
57            useMultiDimAdvec = .TRUE.       &                .OR. PTRACERS_MultiDimAdv(iTracer)
58          ENDIF          PTRACERS_AdamsBashGtr(iTracer) =
59         &       PTRACERS_advScheme(iTracer).EQ.ENUM_CENTERED_2ND .OR.
60         &       PTRACERS_advScheme(iTracer).EQ.ENUM_UPWIND_3RD .OR.
61         &       PTRACERS_advScheme(iTracer).EQ.ENUM_CENTERED_4TH
62    
63  C       end of Tracer loop  C       end of Tracer loop
64        ENDDO        ENDDO
65          _END_MASTER( myThid )
66          _BARRIER
67    
68  C--   Check size of the overlap :  C--   Check size of the overlap :
69        IF ( useCubedSphereExchange .AND. useMultiDimAdvec ) THEN        IF ( useCubedSphereExchange .AND. useMultiDimAdvec ) THEN
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22