OpenAD/code_ad_moc/openad_checkpointInit.F

C $Header: /u/gcmpack/MITgcm_contrib/heimbach/OpenAD/code_ad_moc/openad_checkpointInit.F,v 1.1 2008/09/15 22:16:01 utke Exp $

      SUBROUTINE openad_checkpointinit(uCheckLev1, 
     +     uCheckLev2, 
     +     uCheckLev3,
     +     uCheckLev4 )

      IMPLICIT NONE

#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"

      integer iloop
      logical haveRemainder, allStepsCovered
      integer  uCheckLev1, uCheckLev2, uCheckLev3,uCheckLev4
      integer  nCheckLev1, nCheckLev2, nCheckLev3,nCheckLev4

      uCheckLev1=0
      uCheckLev2=0
      uCheckLev3=0
      uCheckLev4=0

      open(2,action='read',file='ncheckLev.conf')
      read(2,'(I5,/,I5,/I5,/I5)') nCheckLev1, nCheckLev2, nCheckLev3, 
     +nCheckLev4
      close(2)

      write(*,'(A,I4,A,I4,A,I4,A,I8)') 
     +'JU: the_main_loop checkpoints allowed 1st:' ,nCheckLev1,
     +', 2nd:', nCheckLev2,
     +', 3rd:', nCheckLev3,
     +', 4th:', nCheckLev4

cju in this current scheme 
cju there is no point in  specifying 1 as the number of checkpoints
cju because for now don't avoid the incidental store of 
cju the same checkpoint across multiple levels. 
      if (nCheckLev4 .eq. 1 ) then 
         uCheckLev4=0
         write(*,'(A)') 'JU: NOTE: resetting level 4 to 0'
      else
         uCheckLev4=nCheckLev4
      end if
      if (nCheckLev3 .eq. 1 ) then 
         uCheckLev3=0
         write(*,'(A)') 'JU: NOTE: resetting level 3 to 0'
      else
         uCheckLev3=nCheckLev3
      end if
      if (nCheckLev2 .eq. 1 ) then 
         uCheckLev2=0
         write(*,'(A)') 'JU: NOTE: resetting level 2 to 0'
      else
         uCheckLev2=nCheckLev2
      end if
      if (nCheckLev1 .eq. 1 ) then 
         uCheckLev1=0
         write(*,'(A)') 'JU: NOTE: resetting level 1 to 0'
      else
         uCheckLev1=nCheckLev1
      end if
cju check sanity
      if (uCheckLev4 .gt. 0 ) then 
        if ( uCheckLev3*uCheckLev2*uCheckLev1 .eq. 0) then 
           print *, 'JU: ERROR: uCheckLev4 > 0 implies ',
     +'all lower levels must be nonzero too.'
        end if
      end if
      if (uCheckLev3 .gt. 0 ) then 
        if ( uCheckLev2*uCheckLev1 .eq. 0) then 
           print *, 'JU: ERROR: uCheckLev3 > 0 implies ',
     +'all lower levels must be nonzero too.'
        end if
      end if
      if (uCheckLev2 .gt. 0 ) then 
        if ( uCheckLev1 .eq. 0) then 
           print *, 'JU: ERROR: uCheckLev2 > 0 implies ',
     +'all lower levels must be nonzero too.'
        end if
      end if
      if (uCheckLev1 .eq. 0 ) then 
         print *, 'JU: ERROR: uCheckLev1 must be nonzero '
      end if
      
cju uCheckLev1
      allStepsCovered = .false.
      if (nTimeSteps .lt. uCheckLev1 ) then 
         uCheckLev1=nTimeSteps
         allStepsCovered =.true.
      end if

cju uCheckLev2
      if (.NOT. allStepsCovered )  then 
         if (MOD(nTimeSteps,uCheckLev1).gt.0)  then 
            if((nTimeSteps/uCheckLev1)+1 .le. uCheckLev2 ) then 
               uCheckLev2=(nTimeSteps/uCheckLev1)+1
               allStepsCovered =.true.
            end if
         else
            if((nTimeSteps/uCheckLev1) .le. uCheckLev2 ) then 
               uCheckLev2=(nTimeSteps/uCheckLev1)
               allStepsCovered =.true.
            end if
         end if
      else
         uCheckLev2 = 0 
      end if

cju uCheckLev3
      if (.NOT. allStepsCovered )  then 
         if (MOD(nTimeSteps,uCheckLev1*uCheckLev2).gt.0)  then 
            if((nTimeSteps/(uCheckLev1*uCheckLev2))+1 
     +           .le. uCheckLev3 ) then 
               uCheckLev3=(nTimeSteps/(uCheckLev1*uCheckLev2))+1
               allStepsCovered =.true.
            end if
         else
           if((nTimeSteps/(uCheckLev1*uCheckLev2)) .le. uCheckLev3 )then
               uCheckLev3=(nTimeSteps/(uCheckLev1*uCheckLev2))
               allStepsCovered =.true.
            end if
         end if
      else
         uCheckLev3 = 0 
      end if

cju uCheckLev4
      if (.NOT. allStepsCovered )  then 
         if (MOD(nTimeSteps,uCheckLev1*uCheckLev2*uCheckLev3)
     +        .gt.0)  then 
            if((nTimeSteps/(uCheckLev1*uCheckLev2*uCheckLev3))+1 
     +           .le. uCheckLev4 ) then 
               uCheckLev4=
     +              (nTimeSteps/(uCheckLev1*uCheckLev2*uCheckLev3))+1
               allStepsCovered =.true.
            end if
         else
            if((nTimeSteps/(uCheckLev1*uCheckLev2*uCheckLev3)) 
     +           .le. uCheckLev4 ) then 
               uCheckLev4=
     +              (nTimeSteps/(uCheckLev1*uCheckLev2*uCheckLev3))
               allStepsCovered =.true.
            end if
         end if
      else
         uCheckLev4 = 0 
      end if

      if (.NOT. allStepsCovered) then 
         print *, 'JU: ERROR: not enough checkpoints allowed'
      end if

      write(*,'(A,I6)') 'JU: timesteps required:',ntimesteps
      write(*,'(A,I4,A,I4,A,I4,A,I4)') 
     +'JU: the_main_loop checkpoints used    1st:',uCheckLev1,
     +', 2nd:', uCheckLev2,
     +', 3rd:', uCheckLev3,
     +', 4th:', uCheckLev4

      END

1	C $Header: /u/gcmpack/MITgcm_contrib/heimbach/OpenAD/code_ad_moc/openad_checkpointInit.F,v 1.1 2008/09/15 22:16:01 utke Exp $
2
3	SUBROUTINE openad_checkpointinit(uCheckLev1,
4	+ uCheckLev2,
5	+ uCheckLev3,
6	+ uCheckLev4 )
7
8	IMPLICIT NONE
9
10	#include "SIZE.h"
11	#include "EEPARAMS.h"
12	#include "PARAMS.h"
13
14	integer iloop
15	logical haveRemainder, allStepsCovered
16	integer uCheckLev1, uCheckLev2, uCheckLev3,uCheckLev4
17	integer nCheckLev1, nCheckLev2, nCheckLev3,nCheckLev4
18
19	uCheckLev1=0
20	uCheckLev2=0
21	uCheckLev3=0
22	uCheckLev4=0
23
24	open(2,action='read',file='ncheckLev.conf')
25	read(2,'(I5,/,I5,/I5,/I5)') nCheckLev1, nCheckLev2, nCheckLev3,
26	+nCheckLev4
27	close(2)
28
29	write(*,'(A,I4,A,I4,A,I4,A,I8)')
30	+'JU: the_main_loop checkpoints allowed 1st:' ,nCheckLev1,
31	+', 2nd:', nCheckLev2,
32	+', 3rd:', nCheckLev3,
33	+', 4th:', nCheckLev4
34
35	cju in this current scheme
36	cju there is no point in specifying 1 as the number of checkpoints
37	cju because for now don't avoid the incidental store of
38	cju the same checkpoint across multiple levels.
39	if (nCheckLev4 .eq. 1 ) then
40	uCheckLev4=0
41	write(*,'(A)') 'JU: NOTE: resetting level 4 to 0'
42	else
43	uCheckLev4=nCheckLev4
44	end if
45	if (nCheckLev3 .eq. 1 ) then
46	uCheckLev3=0
47	write(*,'(A)') 'JU: NOTE: resetting level 3 to 0'
48	else
49	uCheckLev3=nCheckLev3
50	end if
51	if (nCheckLev2 .eq. 1 ) then
52	uCheckLev2=0
53	write(*,'(A)') 'JU: NOTE: resetting level 2 to 0'
54	else
55	uCheckLev2=nCheckLev2
56	end if
57	if (nCheckLev1 .eq. 1 ) then
58	uCheckLev1=0
59	write(*,'(A)') 'JU: NOTE: resetting level 1 to 0'
60	else
61	uCheckLev1=nCheckLev1
62	end if
63	cju check sanity
64	if (uCheckLev4 .gt. 0 ) then
65	if ( uCheckLev3uCheckLev2uCheckLev1 .eq. 0) then
66	print *, 'JU: ERROR: uCheckLev4 > 0 implies ',
67	+'all lower levels must be nonzero too.'
68	end if
69	end if
70	if (uCheckLev3 .gt. 0 ) then
71	if ( uCheckLev2*uCheckLev1 .eq. 0) then
72	print *, 'JU: ERROR: uCheckLev3 > 0 implies ',
73	+'all lower levels must be nonzero too.'
74	end if
75	end if
76	if (uCheckLev2 .gt. 0 ) then
77	if ( uCheckLev1 .eq. 0) then
78	print *, 'JU: ERROR: uCheckLev2 > 0 implies ',
79	+'all lower levels must be nonzero too.'
80	end if
81	end if
82	if (uCheckLev1 .eq. 0 ) then
83	print *, 'JU: ERROR: uCheckLev1 must be nonzero '
84	end if
85
86	cju uCheckLev1
87	allStepsCovered = .false.
88	if (nTimeSteps .lt. uCheckLev1 ) then
89	uCheckLev1=nTimeSteps
90	allStepsCovered =.true.
91	end if
92
93	cju uCheckLev2
94	if (.NOT. allStepsCovered ) then
95	if (MOD(nTimeSteps,uCheckLev1).gt.0) then
96	if((nTimeSteps/uCheckLev1)+1 .le. uCheckLev2 ) then
97	uCheckLev2=(nTimeSteps/uCheckLev1)+1
98	allStepsCovered =.true.
99	end if
100	else
101	if((nTimeSteps/uCheckLev1) .le. uCheckLev2 ) then
102	uCheckLev2=(nTimeSteps/uCheckLev1)
103	allStepsCovered =.true.
104	end if
105	end if
106	else
107	uCheckLev2 = 0
108	end if
109
110	cju uCheckLev3
111	if (.NOT. allStepsCovered ) then
112	if (MOD(nTimeSteps,uCheckLev1*uCheckLev2).gt.0) then
113	if((nTimeSteps/(uCheckLev1*uCheckLev2))+1
114	+ .le. uCheckLev3 ) then
115	uCheckLev3=(nTimeSteps/(uCheckLev1*uCheckLev2))+1
116	allStepsCovered =.true.
117	end if
118	else
119	if((nTimeSteps/(uCheckLev1*uCheckLev2)) .le. uCheckLev3 )then
120	uCheckLev3=(nTimeSteps/(uCheckLev1*uCheckLev2))
121	allStepsCovered =.true.
122	end if
123	end if
124	else
125	uCheckLev3 = 0
126	end if
127
128	cju uCheckLev4
129	if (.NOT. allStepsCovered ) then
130	if (MOD(nTimeSteps,uCheckLev1uCheckLev2uCheckLev3)
131	+ .gt.0) then
132	if((nTimeSteps/(uCheckLev1uCheckLev2uCheckLev3))+1
133	+ .le. uCheckLev4 ) then
134	uCheckLev4=
135	+ (nTimeSteps/(uCheckLev1uCheckLev2uCheckLev3))+1
136	allStepsCovered =.true.
137	end if
138	else
139	if((nTimeSteps/(uCheckLev1uCheckLev2uCheckLev3))
140	+ .le. uCheckLev4 ) then
141	uCheckLev4=
142	+ (nTimeSteps/(uCheckLev1uCheckLev2uCheckLev3))
143	allStepsCovered =.true.
144	end if
145	end if
146	else
147	uCheckLev4 = 0
148	end if
149
150	if (.NOT. allStepsCovered) then
151	print *, 'JU: ERROR: not enough checkpoints allowed'
152	end if
153
154	write(*,'(A,I6)') 'JU: timesteps required:',ntimesteps
155	write(*,'(A,I4,A,I4,A,I4,A,I4)')
156	+'JU: the_main_loop checkpoints used 1st:',uCheckLev1,
157	+', 2nd:', uCheckLev2,
158	+', 3rd:', uCheckLev3,
159	+', 4th:', uCheckLev4
160
161	END
162