pkg/autodiff/adzero_adj.F

C $Header: /u/gcmpack/MITgcm/pkg/autodiff/adzero_adj.F,v 1.2 2007/10/08 23:50:53 jmc Exp $
C $Name:  $

#include "CPP_OPTIONS.h"

C     !ROUTINE: ADZERO_ADJ_LOC
C     !INTERFACE:
      SUBROUTINE ADZERO_ADJ_LOC(
     I             NrIn,myThid,adFieldFwd)
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | o zero out a local (no nsx,nsy dims.) adjoint field
C     | author: Gael Forget
C     *==========================================================*

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

C     == Routine arguments ==
C     myThid - Thread number for this instance of the routine.
      integer myThid,bi,bj,i,j,k
      _RL adFieldFwd(1-Olx:sNx+Olx,1-Oly:sNy+Oly,NrIn)

      DO k=1,NrIn
       DO j=1-Oly,sNy+Oly
        DO i=1-Olx,sNx+Olx
         adFieldFwd(i,j,k) = 0.
        ENDDO
       ENDDO
      ENDDO

      end


C     !ROUTINE: ADZERO_ADJ
C     !INTERFACE:
      SUBROUTINE ADZERO_ADJ(
     I             NrIn,myThid,adFieldFwd)
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | o zero out a full (incl. nsx,nsy dims.) adjoint field
C     | author: Gael Forget
C     *==========================================================*

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

C     == Routine arguments ==
C     myThid - Thread number for this instance of the routine.
      integer myThid,bi,bj,i,j,k
      _RL adFieldFwd(1-Olx:sNx+Olx,1-Oly:sNy+Oly,NrIn,nSx,nSy)

      DO bj=myByLo(myThid),myByHi(myThid)
       DO bi=myBxLo(myThid),myBxHi(myThid)
        DO k=1,NrIn
         DO j=1-Oly,sNy+Oly
          DO i=1-Olx,sNx+Olx
           adFieldFwd(i,j,k,bi,bj) = 0.
          ENDDO
         ENDDO
        ENDDO
       ENDDO
      ENDDO

      end

1	gforget	1.3	C $Header: /u/gcmpack/MITgcm/pkg/autodiff/adzero_adj.F,v 1.2 2007/10/08 23:50:53 jmc Exp $
2	jmc	1.2	C $Name: $
3	gforget	1.1
4			#include "CPP_OPTIONS.h"
5
6	gforget	1.3	C !ROUTINE: ADZERO_ADJ_LOC
7	gforget	1.1	C !INTERFACE:
8	gforget	1.3	SUBROUTINE ADZERO_ADJ_LOC(
9			I NrIn,myThid,adFieldFwd)
10	gforget	1.1	C !DESCRIPTION: \bv
11			C ==========================================================
12	gforget	1.3	C \| o zero out a local (no nsx,nsy dims.) adjoint field
13	gforget	1.1	C \| author: Gael Forget
14			C ==========================================================
15
16			#include "SIZE.h"
17			#include "EEPARAMS.h"
18			#include "PARAMS.h"
19
20			C == Routine arguments ==
21			C myThid - Thread number for this instance of the routine.
22			integer myThid,bi,bj,i,j,k
23	gforget	1.3	_RL adFieldFwd(1-Olx:sNx+Olx,1-Oly:sNy+Oly,NrIn)
24	gforget	1.1
25	gforget	1.3	DO k=1,NrIn
26	gforget	1.1	DO j=1-Oly,sNy+Oly
27			DO i=1-Olx,sNx+Olx
28			adFieldFwd(i,j,k) = 0.
29			ENDDO
30			ENDDO
31			ENDDO
32
33			end
34
35	gforget	1.3
36			C !ROUTINE: ADZERO_ADJ
37	gforget	1.1	C !INTERFACE:
38	gforget	1.3	SUBROUTINE ADZERO_ADJ(
39			I NrIn,myThid,adFieldFwd)
40	gforget	1.1	C !DESCRIPTION: \bv
41			C ==========================================================
42	gforget	1.3	C \| o zero out a full (incl. nsx,nsy dims.) adjoint field
43	gforget	1.1	C \| author: Gael Forget
44			C ==========================================================
45
46			#include "SIZE.h"
47			#include "EEPARAMS.h"
48			#include "PARAMS.h"
49
50			C == Routine arguments ==
51			C myThid - Thread number for this instance of the routine.
52	gforget	1.3	integer myThid,bi,bj,i,j,k
53			_RL adFieldFwd(1-Olx:sNx+Olx,1-Oly:sNy+Oly,NrIn,nSx,nSy)
54	gforget	1.1
55	gforget	1.3	DO bj=myByLo(myThid),myByHi(myThid)
56			DO bi=myBxLo(myThid),myBxHi(myThid)
57			DO k=1,NrIn
58			DO j=1-Oly,sNy+Oly
59			DO i=1-Olx,sNx+Olx
60			adFieldFwd(i,j,k,bi,bj) = 0.
61			ENDDO
62			ENDDO
63	gforget	1.1	ENDDO
64			ENDDO
65	gforget	1.3	ENDDO
66	gforget	1.1
67			end
68