pkg/autodiff/autodiff_check.F

C $Header: /u/gcmpack/MITgcm/pkg/autodiff/autodiff_check.F,v 1.3 2005/08/10 03:34:48 heimbach Exp $
C $Name:  $

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

#include "AUTODIFF_OPTIONS.h"

CBOP
C     !ROUTINE: AUTODIFF_CHECK
C     !INTERFACE:
      SUBROUTINE AUTODIFF_CHECK( myThid )

C     !DESCRIPTION: \bv
C     \ev

      IMPLICIT NONE
#include "SIZE.h"
#include "GRID.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#ifdef ALLOW_AUTODIFF
#include "tamc.h"
#endif

C     !INPUT/OUTPUT PARAMETERS:
C     myThid -  Number of this instances
      INTEGER myThid
CEOP

#ifdef ALLOW_AUTODIFF_TAMC

#if (defined (AUTODIFF_2_LEVEL_CHECKPOINT))
      if (nchklev_1*nchklev_2 .lt. nTimeSteps) then
        print*, ' the_main_loop: TAMC checkpointing parameters'
        print*, '     nchklev_1*nchklev_2 = ',
     &       nchklev_1*nchklev_2
        print*, '            are not consistent with nTimeSteps = ',
     &       nTimeSteps
        stop    ' ... stopped in autodiff_check'
      endif
#elif (defined (AUTODIFF_4_LEVEL_CHECKPOINT))
      if (nchklev_1*nchklev_2*nchklev_3*nchklev_4 .lt. nTimeSteps) then
        print*, ' the_main_loop: TAMC checkpointing parameters'
        print*, '     nchklev_1*nchklev_2*nchklev_3*nchklev_4 = ',
     &       nchklev_1*nchklev_2*nchklev_3*nchklev_4
        print*, '            are not consistent with nTimeSteps = ',
     &       nTimeSteps
        stop    ' ... stopped in autodiff_check'
      endif
#else
c--   Check the choice of the checkpointing parameters in relation
c--   to nTimeSteps: (nchklev_1*nchklev_2*nchklev_3 .ge. nTimeSteps)
      if (nchklev_1*nchklev_2*nchklev_3 .lt. nTimeSteps) then
        print*, ' the_main_loop: TAMC checkpointing parameters'
        print*, '     nchklev_1*nchklev_2*nchklev_3 = ',
     &       nchklev_1*nchklev_2*nchklev_3
        print*, '            are not consistent with nTimeSteps = ',
     &       nTimeSteps
        stop    ' ... stopped in autodiff_check'
      endif
#endif

      IF ( select_rStar.NE.0 ) THEN
        STOP 'in autodiff_check: adjoint not implemented for rStar'
      ENDIF

#endif

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

      RETURN
      END
1	heimbach	1.4	C $Header: /u/gcmpack/MITgcm/pkg/autodiff/autodiff_check.F,v 1.3 2005/08/10 03:34:48 heimbach Exp $
2	heimbach	1.1	C $Name: $
3
4			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
5
6	jmc	1.2	#include "AUTODIFF_OPTIONS.h"
7	heimbach	1.1
8			CBOP
9			C !ROUTINE: AUTODIFF_CHECK
10			C !INTERFACE:
11			SUBROUTINE AUTODIFF_CHECK( myThid )
12
13			C !DESCRIPTION: \bv
14			C \ev
15
16			IMPLICIT NONE
17			#include "SIZE.h"
18			#include "GRID.h"
19			#include "EEPARAMS.h"
20			#include "PARAMS.h"
21			#ifdef ALLOW_AUTODIFF
22			#include "tamc.h"
23			#endif
24
25			C !INPUT/OUTPUT PARAMETERS:
26			C myThid - Number of this instances
27			INTEGER myThid
28			CEOP
29
30			#ifdef ALLOW_AUTODIFF_TAMC
31
32	heimbach	1.3	#if (defined (AUTODIFF_2_LEVEL_CHECKPOINT))
33	heimbach	1.1	if (nchklev_1*nchklev_2 .lt. nTimeSteps) then
34			print*, ' the_main_loop: TAMC checkpointing parameters'
35	heimbach	1.3	print, ' nchklev_1nchklev_2 = ',
36	heimbach	1.1	& nchklev_1*nchklev_2
37			print*, ' are not consistent with nTimeSteps = ',
38			& nTimeSteps
39			stop ' ... stopped in autodiff_check'
40			endif
41	heimbach	1.3	#elif (defined (AUTODIFF_4_LEVEL_CHECKPOINT))
42			if (nchklev_1nchklev_2nchklev_3*nchklev_4 .lt. nTimeSteps) then
43			print*, ' the_main_loop: TAMC checkpointing parameters'
44			print, ' nchklev_1nchklev_2nchklev_3nchklev_4 = ',
45			& nchklev_1nchklev_2nchklev_3*nchklev_4
46			print*, ' are not consistent with nTimeSteps = ',
47			& nTimeSteps
48			stop ' ... stopped in autodiff_check'
49			endif
50	heimbach	1.1	#else
51			c-- Check the choice of the checkpointing parameters in relation
52			c-- to nTimeSteps: (nchklev_1nchklev_2nchklev_3 .ge. nTimeSteps)
53			if (nchklev_1nchklev_2nchklev_3 .lt. nTimeSteps) then
54			print*, ' the_main_loop: TAMC checkpointing parameters'
55	heimbach	1.3	print, ' nchklev_1nchklev_2*nchklev_3 = ',
56	heimbach	1.1	& nchklev_1nchklev_2nchklev_3
57			print*, ' are not consistent with nTimeSteps = ',
58			& nTimeSteps
59			stop ' ... stopped in autodiff_check'
60			endif
61			#endif
62
63	heimbach	1.4	IF ( select_rStar.NE.0 ) THEN
64			STOP 'in autodiff_check: adjoint not implemented for rStar'
65			ENDIF
66
67	heimbach	1.1	#endif
68
69			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
70
71			RETURN
72			END