105 |
#include "CG3D.h" |
#include "CG3D.h" |
106 |
#endif |
#endif |
107 |
|
|
108 |
|
#ifdef ALLOW_AUTODIFF_TAMC |
109 |
|
#include "tamc.h" |
110 |
|
#include "cost.h" |
111 |
|
#include "adcost.h" |
112 |
|
#endif |
113 |
|
|
114 |
C == Routine arguments == |
C == Routine arguments == |
115 |
C myThid - Thread number for this instance of the routine. |
C myThid - Thread number for this instance of the routine. |
116 |
INTEGER myThid |
INTEGER myThid |
135 |
CALL INITIALISE_FIXED( myThid ) |
CALL INITIALISE_FIXED( myThid ) |
136 |
CALL TIMER_STOP ('INITIALISE_FIXED [THE_MODEL_MAIN]',myThid) |
CALL TIMER_STOP ('INITIALISE_FIXED [THE_MODEL_MAIN]',myThid) |
137 |
|
|
|
|
|
138 |
myCurrentTime = startTime |
myCurrentTime = startTime |
139 |
myCurrentIter = nIter0 |
myCurrentIter = nIter0 |
140 |
|
|
141 |
|
#ifdef ALLOW_MIT_ADJOINT_RUN |
142 |
|
|
143 |
|
_BEGIN_MASTER( mythid ) |
144 |
|
IF (myProcId .eq. 0) THEN |
145 |
|
CALL CTRL_UNPACK( mycurrentiter, mycurrenttime, mythid ) |
146 |
|
ENDIF |
147 |
|
_END_MASTER( mythid ) |
148 |
|
_BARRIER |
149 |
|
|
150 |
|
fc = 0.0 |
151 |
|
adfc = 1.0 |
152 |
|
CALL TIMER_START('ADTHE_MAIN_SENS [ADJOINT RUN]', mythid) |
153 |
|
call ADTHE_MAIN_SENS ( mythid ) |
154 |
|
call TIMER_STOP ('ADTHE_MAIN_SENS [ADJOINT RUN]', mythid) |
155 |
|
_BARRIER |
156 |
|
|
157 |
|
_BEGIN_MASTER( mythid ) |
158 |
|
IF (myProcId .eq. 0) THEN |
159 |
|
call CTRL_PACK( mycurrentiter, mycurrenttime, mythid ) |
160 |
|
ENDIF |
161 |
|
_END_MASTER( mythid ) |
162 |
|
_BARRIER |
163 |
|
|
164 |
|
#else /* ALLOW_MIT_ADJOINT_RUN undef */ |
165 |
|
|
166 |
C-- Call time stepping loop of full model |
C-- Call time stepping loop of full model |
167 |
CALL TIMER_START('THE_MAIN_LOOP [THE_MODEL_MAIN]',myThid) |
CALL TIMER_START('THE_MAIN_LOOP [THE_MODEL_MAIN]',myThid) |
168 |
CALL THE_MAIN_LOOP(myCurrentTime, myCurrentIter, myThid ) |
CALL THE_MAIN_SENS(myCurrentTime, myCurrentIter, myThid ) |
169 |
CALL TIMER_STOP ('THE_MAIN_LOOP [THE_MODEL_MAIN]',myThid) |
CALL TIMER_STOP ('THE_MAIN_LOOP [THE_MODEL_MAIN]',myThid) |
170 |
|
|
171 |
|
#endif /* ALLOW_MIT_ADJOINT_RUN */ |
172 |
|
|
173 |
CALL TIMER_STOP ('ALL [THE_MODEL_MAIN]',myThid) |
CALL TIMER_STOP ('ALL [THE_MODEL_MAIN]',myThid) |
174 |
|
|