/[MITgcm]/MITgcm/model/src/the_model_main.F

Diff of /MITgcm/model/src/the_model_main.F

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

-revision 1.6 by cnh,
Fri Apr 24 03:45:12 1998 UTC
+revision 1.39 by heimbach,
Fri Jul 13 14:26:57 2001 UTC
 Line 1
  C $Header$
+ C $Name$
- #include "CPP_EEOPTIONS.h"
+ #include "CPP_OPTIONS.h"
        SUBROUTINE THE_MODEL_MAIN(myThid)
  C     /==========================================================\
-Line 22 
 C     | =====
+Line 23 
 C     | =====
  C     | C*P* comments indicating place holders for which code is |
  C     |      presently being developed.                          |
  C     \==========================================================/
+       IMPLICIT NONE
+ C
  C     == Global variables ===
  #include "SIZE.h"
  #include "EEPARAMS.h"
  #include "PARAMS.h"
- #include "CG2D.h"
  #include "DYNVARS.h"
+ #ifdef ALLOW_NONHYDROSTATIC
+ #include "CG3D.h"
+ #endif
+ #ifdef ALLOW_AUTODIFF_TAMC
+ #include "tamc.h"
+ #include "cost.h"
+ #include "adcost.h"
+ #endif
  C     == Routine arguments ==
  C     myThid - Thread number for this instance of the routine.
        INTEGER myThid
  C     == Local variables ==
- C     Note: Under the multi-threaded model myCurrentIter and myCurrentTime
+ C     Note: Under the multi-threaded model myCurrentIter and
- C           are local variables passed around as routine arguments. Although
+ C           myCurrentTime are local variables passed around as routine
- C           this is fiddly it saves the need to impose additional synchronisation
+ C           arguments. Although this is fiddly it saves the need to
- C           points when they are updated.
+ C           impose additional synchronisation points when they are
- C     myCurrentIter - Iteration counter for this thread
+ C           updated.
  C     myCurrentTime - Time counter for this thread
- C     I             - Loop counter
+ C     myCurrentIter - Iteration counter for this thread
-       INTEGER I, myCurrentIter
+       INTEGER myCurrentIter
-       REAL    myCurrentTime
+       _RL     myCurrentTime
- C--   Set model initial conditions
-       CALL INITIALISE( myThid )
-       myCurrentTime = startTime
-       myCurrentIter = nIter0
- C--   Begin time stepping loop
+ C--   This timer encompasses the whole code
-       DO I=1, nTimeSteps
+       CALL TIMER_START('ALL                [THE_MODEL_MAIN]',myThid)
- C--    Do IO if needed.
-        CALL DO_THE_MODEL_IO( myCurrentIter, myThid )
- C--    Load offline tracer fields
+ C--   Set model configuration (fixed arrays)
-        IF ( MOD(myCurrentIter,numStepsPerPickup) .EQ. 1 ) THEN
+       CALL TIMER_START('INITIALISE_FIXED   [THE_MODEL_MAIN]',myThid)
- C       I/O
+       CALL INITIALISE_FIXED( myThid )
- C       o Writes to arrays are performed by their own thread (to ensure
+       CALL TIMER_STOP ('INITIALISE_FIXED   [THE_MODEL_MAIN]',myThid)
- C         physical memory allocation will pair with thread).
- C       o Thread 1 reads into a buffer.
- C       CALL LOAD_OFFLINE_FIELDS( myCurrentTime, myCurrentIter, myThid )
-        ENDIF
- C--    Step forward fields and calculate time tendency terms
-        CALL DYNAMICS( myThid )
- C--    Solve elliptic equation(s).
- C      Two-dimensional only for conventional hydrostatic or three-dimensional
- C      for non-hydrostatic and/or IGW scheme.
-        CALL SOLVE_FOR_PRESSURE( myThid )
- C--    Do "blocking" sends and receives for tendency "overlap" terms
-        CALL DO_GTERM_BLOCKING_EXCHANGES( myThid )
-        myCurrentIter = myCurrentIter + 1
+       myCurrentTime = startTime
-        myCurrentTime = myCurrentTime  + deltaTtracer
+       myCurrentIter = nIter0
- C--    Save state for restarts
+ #ifdef ALLOW_ADJOINT_RUN
- C*P*   CALL CHECKPOINT( myThid )
-       ENDDO
+       _BEGIN_MASTER( mythid )
+       IF (myProcId .eq. 0) THEN
+          CALL CTRL_UNPACK( mycurrentiter, mycurrenttime, mythid )
+       ENDIF
+       _END_MASTER( mythid )
+       _BARRIER
+       fc   = 0.0
+       adfc = 1.0
+       CALL TIMER_START('ADTHE_MAIN_LOOP          [ADJOINT RUN]', mythid)
+       CALL ADTHE_MAIN_LOOP ( mythid )
+       CALL TIMER_STOP ('ADTHE_MAIN_LOOP          [ADJOINT RUN]', mythid)
+       _BARRIER
+       _BEGIN_MASTER( mythid )
+       IF (myProcId .eq. 0) THEN
+          call CTRL_PACK( mycurrentiter, mycurrenttime, mythid )
+       ENDIF
+       _END_MASTER( mythid )
+       _BARRIER
+ #ifdef ALLOW_GRADIENT_CHECK
+       if (useGrdchk) then
+          CALL TIMER_START('GRDCHK_MAIN          [ADJOINT RUN]', mythid)
+          CALL GRDCHK_MAIN( mythid )
+          CALL TIMER_STOP ('GRDCHK_MAIN          [ADJOINT RUN]', mythid)
+          _BARRIER
+       end if
+ #endif
+ #else /* ALLOW_ADJOINT_RUN undef */
+ C--   Call time stepping loop of full model
+       CALL TIMER_START('THE_MAIN_LOOP      [THE_MODEL_MAIN]',myThid)
+       CALL THE_MAIN_LOOP(myCurrentTime, myCurrentIter, myThid )
+       CALL TIMER_STOP ('THE_MAIN_LOOP      [THE_MODEL_MAIN]',myThid)
+ #endif /* ALLOW_ADJOINT_RUN */
+ C--   Final checkpoint (incase the in-loop checkpoint was missed)
+       CALL TIMER_START('WRITE_CHECKPOINT   [THE_MODEL_MAIN]',myThid)
+       CALL WRITE_CHECKPOINT(
+      &        .TRUE., myCurrentTime, myCurrentIter, myThid )
+       CALL TIMER_STOP ('WRITE_CHECKPOINT   [THE_MODEL_MAIN]',myThid)
+       CALL TIMER_STOP ('ALL                [THE_MODEL_MAIN]',myThid)
+ C--   Write timer statistics
+       IF ( myThid .EQ. 1 ) THEN
+        CALL TIMER_PRINTALL( myThid )
+        CALL COMM_STATS
+       ENDIF
        RETURN
        END

 Legend:



Removed from v.1.6
 


changed lines


 
Added in v.1.39
 Legend:



Removed from v.1.6
 


changed lines


 
Added in v.1.39
-Removed from v.1.6
+Added in v.1.39

	ViewVC Help
Powered by ViewVC 1.1.22