forked from altMITgcm/MITgcm
-
Notifications
You must be signed in to change notification settings - Fork 0
/
autodiff_inadmode_set_ad.F
120 lines (107 loc) · 4.24 KB
/
autodiff_inadmode_set_ad.F
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
C $Header: /u/gcmpack/MITgcm/pkg/autodiff/autodiff_inadmode_set_ad.F,v 1.14 2016/01/22 15:33:21 mlosch Exp $
C $Name: $
#include "AUTODIFF_OPTIONS.h"
#ifdef ALLOW_CTRL
# include "CTRL_OPTIONS.h"
#endif
SUBROUTINE ADAUTODIFF_INADMODE_SET( myThid )
C *==========================================================*
C | SUBROUTINE ADAUTODIFF_INADMODE_SET
C *==========================================================*
IMPLICIT NONE
C == Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "AUTODIFF_PARAMS.h"
#ifdef ALLOW_SEAICE
#include "SEAICE_SIZE.h"
#include "SEAICE_PARAMS.h"
#endif
#ifdef ALLOW_CTRL
#include "ctrl.h"
#endif
C == Routine arguments ==
C myThid :: my Thread Id number
INTEGER myThid
CHARACTER*(MAX_LEN_MBUF) msgBuf
inAdMode = inAdTrue
useKPP = useKPPinAdMode
useGMRedi = useGMRediInAdMode
useSEAICE = useSEAICEinAdMode
useGGL90 = useGGL90inAdMode
useSALT_PLUME = useSALT_PLUMEinAdMode
#if (defined (ALLOW_CTRL) && defined (ECCO_CTRL_DEPRECATED))
ctrlSmoothCorrel2D = useSmoothCorrel2DinAdMode
#endif
#ifdef ALLOW_SEAICE
IF ( SEAICEuseFREEDRIFTswitchInAd ) THEN
SEAICEuseFREEDRIFT = .NOT.SEAICEuseFREEDRIFTinFwdMode
SEAICEuseLSR = .NOT.SEAICEuseFREEDRIFT
ENDIF
IF ( SEAICEuseDYNAMICSswitchInAd ) THEN
SEAICEuseDYNAMICS = .NOT.SEAICEuseDYNAMICSinFwdMode
ENDIF
SEAICEadjMODE=SEAICEapproxLevInAd
#endif /* ALLOW_SEAICE */
viscFacAdj=viscFacInAd
IF ( debugLevel.GE.debLevC ) THEN
WRITE(msgBuf,'(A,L5)')
& 'S/R ADAUTODIFF_INADMODE_SET: useKPP =', useKPP
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
& SQUEEZE_RIGHT, myThid )
WRITE(msgBuf,'(A,L5)')
& 'S/R ADAUTODIFF_INADMODE_SET: useGMRedi =', useGMRedi
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
& SQUEEZE_RIGHT, myThid )
WRITE(msgBuf,'(A,L5)')
& 'S/R ADAUTODIFF_INADMODE_SET: useSEAICE =', useSEAICE
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
& SQUEEZE_RIGHT, myThid )
WRITE(msgBuf,'(A,L5)')
& 'S/R ADAUTODIFF_INADMODE_SET: useGGL90 =', useGGL90
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
& SQUEEZE_RIGHT, myThid )
WRITE(msgBuf,'(A,L5)')
& 'S/R ADAUTODIFF_INADMODE_SET: useSALT_PLUME =',
& useSALT_PLUME
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
& SQUEEZE_RIGHT, myThid )
WRITE(msgBuf,'(A,1PE21.14)')
& 'S/R ADAUTODIFF_INADMODE_SET: viscFacAdj =',viscFacAdj
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
& SQUEEZE_RIGHT, myThid )
#if (defined (ALLOW_CTRL) && defined (ECCO_CTRL_DEPRECATED))
WRITE(msgBuf,'(A,L5)')
& 'S/R ADAUTODIFF_INADMODE_SET: ctrlSmoothCorrel2D =',
& ctrlSmoothCorrel2D
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
& SQUEEZE_RIGHT, myThid )
#endif
#ifdef ALLOW_SEAICE
IF ( SEAICEuseFREEDRIFTswitchInAd ) THEN
WRITE(msgBuf,'(2A,L5)') 'S/R ADAUTODIFF_INADMODE_SET: ',
& 'SEAICEuseFREEDRIFT =', SEAICEuseFREEDRIFT
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
& SQUEEZE_RIGHT, myThid )
WRITE(msgBuf,'(2A,L5)') 'S/R ADAUTODIFF_INADMODE_SET: ',
& 'SEAICEuseLSR =', SEAICEuseLSR
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
& SQUEEZE_RIGHT, myThid )
ENDIF
IF ( SEAICEuseDYNAMICSswitchInAd ) THEN
WRITE(msgBuf,'(2A,L5)') 'S/R ADAUTODIFF_INADMODE_SET: ',
& 'SEAICEuseDYNAMICS =', SEAICEuseDYNAMICS
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
& SQUEEZE_RIGHT, myThid )
ENDIF
IF ( SEAICEapproxLevInAd.NE.0 ) THEN
WRITE(msgBuf,'(2A,I2)') 'S/R ADAUTODIFF_INADMODE_SET: ',
& 'SEAICEadjMODE =', SEAICEadjMODE
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
& SQUEEZE_RIGHT, myThid )
ENDIF
#endif /* ALLOW_SEAICE */
ENDIF
RETURN
END