/[MITgcm]/MITgcm/pkg/autodiff/adzero_adj.F
ViewVC logotype

Annotation of /MITgcm/pkg/autodiff/adzero_adj.F

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

Revision 1.7 - (hide annotations) (download)
Sun Aug 12 18:46:20 2012 UTC (11 years, 9 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint63r, checkpoint63s, checkpoint64, checkpoint65, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65o, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, HEAD
Changes since 1.6: +59 -43 lines 
- remove PARAMS.h (not needed)
- just includes CPP_EEOPTIONS.h (previously CPP_OPTIONS.h);
1 jmc 1.7 C $Header: /u/gcmpack/MITgcm/pkg/autodiff/adzero_adj.F,v 1.6 2010/08/06 18:41:03 gforget Exp $
2 jmc 1.2 C $Name: $
3 gforget 1.1
4 jmc 1.7 #include "CPP_EEOPTIONS.h"
5 gforget 1.1
6 jmc 1.7 C-- File zero_adj.F:
7     C-- Contents
8     C-- o ADZERO_ADJ_1D
9     C-- o ADZERO_ADJ_LOC
10     C-- o ADZERO_ADJ
11    
12     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
13     CBOP
14 gforget 1.6 C !ROUTINE: ADZERO_ADJ_1D
15     C !INTERFACE:
16     SUBROUTINE ADZERO_ADJ_1D(
17 jmc 1.7 I NrIn, adFieldFwd, myThid )
18 gforget 1.6 C !DESCRIPTION: \bv
19     C *==========================================================*
20     C | o zero out a 1D adjoint field
21     C | author: Gael Forget
22     C *==========================================================*
23 jmc 1.7 C \ev
24 gforget 1.6
25 jmc 1.7 C !USES:
26 gforget 1.6 IMPLICIT NONE
27 jmc 1.7 C === Global variables ===
28 gforget 1.6 #include "SIZE.h"
29     #include "EEPARAMS.h"
30    
31 jmc 1.7 C !INPUT/OUTPUT PARAMETERS:
32     C myThid :: my Thread Id number
33     INTEGER NrIn, myThid
34 gforget 1.6 _RL adFieldFwd(NrIn)
35    
36 jmc 1.7 C !LOCAL VARIABLES:
37     INTEGER k
38     CEOP
39 gforget 1.6
40     DO k=1,NrIn
41     adFieldFwd(k) = 0. _d 0
42     ENDDO
43    
44 jmc 1.7 RETURN
45     END
46 gforget 1.6
47 jmc 1.7 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
48     CBOP
49 gforget 1.3 C !ROUTINE: ADZERO_ADJ_LOC
50 gforget 1.1 C !INTERFACE:
51 gforget 1.3 SUBROUTINE ADZERO_ADJ_LOC(
52 jmc 1.7 I NrIn, adFieldFwd, myThid )
53 gforget 1.1 C !DESCRIPTION: \bv
54     C *==========================================================*
55 gforget 1.3 C | o zero out a local (no nsx,nsy dims.) adjoint field
56 gforget 1.1 C | author: Gael Forget
57     C *==========================================================*
58 jmc 1.7 C \ev
59 gforget 1.1
60 jmc 1.7 C !USES:
61 heimbach 1.4 IMPLICIT NONE
62 jmc 1.7 C === Global variables ===
63 gforget 1.1 #include "SIZE.h"
64     #include "EEPARAMS.h"
65    
66 jmc 1.7 C !INPUT/OUTPUT PARAMETERS:
67     C myThid :: my Thread Id number
68     INTEGER NrIn, myThid
69     _RL adFieldFwd(1-OLx:sNx+OLx,1-OLy:sNy+OLy,NrIn)
70    
71     C !LOCAL VARIABLES:
72     INTEGER i,j,k
73     CEOP
74 heimbach 1.4
75 gforget 1.3 DO k=1,NrIn
76 jmc 1.7 DO j=1-OLy,sNy+OLy
77     DO i=1-OLx,sNx+OLx
78 heimbach 1.4 adFieldFwd(i,j,k) = 0. _d 0
79 gforget 1.1 ENDDO
80     ENDDO
81     ENDDO
82    
83 jmc 1.7 RETURN
84     END
85 gforget 1.3
86 jmc 1.7 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
87     CBOP
88 gforget 1.3 C !ROUTINE: ADZERO_ADJ
89 gforget 1.1 C !INTERFACE:
90 gforget 1.3 SUBROUTINE ADZERO_ADJ(
91 jmc 1.7 I NrIn, adFieldFwd, myThid )
92 gforget 1.1 C !DESCRIPTION: \bv
93     C *==========================================================*
94 gforget 1.3 C | o zero out a full (incl. nsx,nsy dims.) adjoint field
95 gforget 1.1 C | author: Gael Forget
96     C *==========================================================*
97 jmc 1.7 C \ev
98 gforget 1.1
99 jmc 1.7 C !USES:
100 heimbach 1.4 IMPLICIT NONE
101 jmc 1.7 C === Global variables ===
102 gforget 1.1 #include "SIZE.h"
103     #include "EEPARAMS.h"
104    
105 jmc 1.7 C !INPUT/OUTPUT PARAMETERS:
106     C myThid :: my Thread Id number
107     INTEGER NrIn, myThid
108     _RL adFieldFwd(1-OLx:sNx+OLx,1-OLy:sNy+OLy,NrIn,nSx,nSy)
109    
110     C !LOCAL VARIABLES:
111     INTEGER bi,bj
112     INTEGER i,j,k
113     CEOP
114 heimbach 1.4
115 gforget 1.3 DO bj=myByLo(myThid),myByHi(myThid)
116     DO bi=myBxLo(myThid),myBxHi(myThid)
117     DO k=1,NrIn
118 jmc 1.7 DO j=1-OLy,sNy+OLy
119     DO i=1-OLx,sNx+OLx
120 heimbach 1.4 adFieldFwd(i,j,k,bi,bj) = 0. _d 0
121 gforget 1.3 ENDDO
122     ENDDO
123 gforget 1.1 ENDDO
124     ENDDO
125 gforget 1.3 ENDDO
126 gforget 1.1
127 jmc 1.7 RETURN
128     END
  ViewVC Help
Powered by ViewVC 1.1.22