/[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.4 - (hide annotations) (download)
Wed Jul 23 20:26:12 2008 UTC (15 years, 9 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint62, checkpoint62b, checkpoint62a, checkpoint61f, checkpoint61g, checkpoint61d, checkpoint61e, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61n, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61q, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.3: +17 -5 lines 
Clean up (no IMPLICIT NONE)
1 heimbach 1.4 C $Header: /u/gcmpack/MITgcm/pkg/autodiff/adzero_adj.F,v 1.3 2008/06/11 18:29:11 gforget 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 heimbach 1.4 IMPLICIT NONE
17    
18     C == Global arguments ==
19 gforget 1.1 #include "SIZE.h"
20     #include "EEPARAMS.h"
21     #include "PARAMS.h"
22    
23     C == Routine arguments ==
24     C myThid - Thread number for this instance of the routine.
25 heimbach 1.4 integer NrIn,myThid
26 gforget 1.3 _RL adFieldFwd(1-Olx:sNx+Olx,1-Oly:sNy+Oly,NrIn)
27 gforget 1.1
28 heimbach 1.4 C == Local arguments ==
29     integer bi,bj,i,j,k
30    
31 gforget 1.3 DO k=1,NrIn
32 gforget 1.1 DO j=1-Oly,sNy+Oly
33     DO i=1-Olx,sNx+Olx
34 heimbach 1.4 adFieldFwd(i,j,k) = 0. _d 0
35 gforget 1.1 ENDDO
36     ENDDO
37     ENDDO
38    
39     end
40    
41 gforget 1.3
42     C !ROUTINE: ADZERO_ADJ
43 gforget 1.1 C !INTERFACE:
44 gforget 1.3 SUBROUTINE ADZERO_ADJ(
45     I NrIn,myThid,adFieldFwd)
46 gforget 1.1 C !DESCRIPTION: \bv
47     C *==========================================================*
48 gforget 1.3 C | o zero out a full (incl. nsx,nsy dims.) adjoint field
49 gforget 1.1 C | author: Gael Forget
50     C *==========================================================*
51    
52 heimbach 1.4 IMPLICIT NONE
53    
54     C == Global arguments ==
55 gforget 1.1 #include "SIZE.h"
56     #include "EEPARAMS.h"
57     #include "PARAMS.h"
58    
59     C == Routine arguments ==
60     C myThid - Thread number for this instance of the routine.
61 heimbach 1.4 integer NrIn,myThid
62 gforget 1.3 _RL adFieldFwd(1-Olx:sNx+Olx,1-Oly:sNy+Oly,NrIn,nSx,nSy)
63 gforget 1.1
64 heimbach 1.4 C == Local arguments ==
65     integer bi,bj,i,j,k
66    
67 gforget 1.3 DO bj=myByLo(myThid),myByHi(myThid)
68     DO bi=myBxLo(myThid),myBxHi(myThid)
69     DO k=1,NrIn
70     DO j=1-Oly,sNy+Oly
71     DO i=1-Olx,sNx+Olx
72 heimbach 1.4 adFieldFwd(i,j,k,bi,bj) = 0. _d 0
73 gforget 1.3 ENDDO
74     ENDDO
75 gforget 1.1 ENDDO
76     ENDDO
77 gforget 1.3 ENDDO
78 gforget 1.1
79     end
80    
  ViewVC Help
Powered by ViewVC 1.1.22