/[MITgcm]/MITgcm_contrib/darwin2/pkg/darwin/darwin_diags.F
ViewVC logotype

Diff of /MITgcm_contrib/darwin2/pkg/darwin/darwin_diags.F

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

revision 1.7 by stephd, Tue Oct 23 17:55:30 2012 UTC revision 1.10 by benw, Tue May 19 14:30:15 2015 UTC
# Line 32  C     myThid -  Number of this instance Line 32  C     myThid -  Number of this instance
32        INTEGER myThid        INTEGER myThid
33  CEndOfInterface  CEndOfInterface
34    
35    #ifdef ALLOW_DARWIN
36    #ifdef ALLOW_TIMEAVE
37    
38  C     == Local variables ==  C     == Local variables ==
39        LOGICAL  DIFFERENT_MULTIPLE        LOGICAL  DIFFERENT_MULTIPLE
# Line 50  C     Functions Line 52  C     Functions
52    
53        _RL Nfixtot, PPtot, Denittot        _RL Nfixtot, PPtot, Denittot
54    
 #ifdef ALLOW_DARWIN  
   
55        IF ( writeBinaryPrec .EQ. precFloat64 ) THEN        IF ( writeBinaryPrec .EQ. precFloat64 ) THEN
56          pf(1:1) = 'D'          pf(1:1) = 'D'
57        ELSE        ELSE
58          pf(1:1) = 'R'          pf(1:1) = 'R'
59        ENDIF        ENDIF
60    
 #ifdef ALLOW_TIMEAVE  
   
61  #ifdef ALLOW_QUOTA  #ifdef ALLOW_QUOTA
62        iPhy = ibiomass        iPhy = ibiomass
63  #endif  #endif
# Line 67  C     Functions Line 65  C     Functions
65  C     Dump files and restart average computation if needed  C     Dump files and restart average computation if needed
66        IF ( DIFFERENT_MULTIPLE(taveFreq,myTime,deltaTClock).        IF ( DIFFERENT_MULTIPLE(taveFreq,myTime,deltaTClock).
67       &                               and.myIter.ne.nIter0 ) THEN       &                               and.myIter.ne.nIter0 ) THEN
68            
69  C      Normalize by integrated time  C      Normalize by integrated time
70         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
71          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
72           CALL TIMEAVE_NORMALIZ(PARave, DAR_timeave, nR ,           CALL TIMEAVE_NORMALIZE(PARave, DAR_timeave, nR ,
73       &                 bi,bj,myThid)       &                 bi,bj,myThid)
74           CALL TIMEAVE_NORMALIZ(PPave, DAR_timeave, nR ,           CALL TIMEAVE_NORMALIZE(PPave, DAR_timeave, nR ,
75       &                 bi,bj,myThid)       &                 bi,bj,myThid)
76  #ifdef QUOTA_DIAG_LIMIT  #ifdef QUOTA_DIAG_LIMIT
77           do np=1,npmax           do np=1,npmax
78             CALL TIMEAVE_NORMALIZ(Nlimave(1-Olx,1-Oly,1,1,1,np),  ! Carbon
79       &                          DAR_timeave, nR ,bi,bj,myThid)             CALL TIMEAVE_NORMALIZE(AP_C_ave(1-Olx,1-Oly,1,1,1,np),
80             CALL TIMEAVE_NORMALIZ(Flimave(1-Olx,1-Oly,1,1,1,np),       &                           DAR_timeave, nR ,bi,bj,myThid)
81       &                          DAR_timeave, nR ,bi,bj,myThid)             CALL TIMEAVE_NORMALIZE(HP_C_ave(1-Olx,1-Oly,1,1,1,np),
82             CALL TIMEAVE_NORMALIZ(Ilimave(1-Olx,1-Oly,1,1,1,np),       &                           DAR_timeave, nR ,bi,bj,myThid)
83       &                          DAR_timeave, nR ,bi,bj,myThid)  ! Nitrogen
84               CALL TIMEAVE_NORMALIZE(AP_N_ave(1-Olx,1-Oly,1,1,1,np),
85         &                           DAR_timeave, nR ,bi,bj,myThid)
86               CALL TIMEAVE_NORMALIZE(HP_N_ave(1-Olx,1-Oly,1,1,1,np),
87         &                           DAR_timeave, nR ,bi,bj,myThid)
88               CALL TIMEAVE_NORMALIZE(Nlimave(1-Olx,1-Oly,1,1,1,np),
89         &                           DAR_timeave, nR ,bi,bj,myThid)
90    ! Phosphorus
91    #ifdef PQUOTA
92               CALL TIMEAVE_NORMALIZE(AP_P_ave(1-Olx,1-Oly,1,1,1,np),
93         &                           DAR_timeave, nR ,bi,bj,myThid)
94               CALL TIMEAVE_NORMALIZE(HP_P_ave(1-Olx,1-Oly,1,1,1,np),
95         &                           DAR_timeave, nR ,bi,bj,myThid)
96               CALL TIMEAVE_NORMALIZE(Plimave(1-Olx,1-Oly,1,1,1,np),
97         &                           DAR_timeave, nR ,bi,bj,myThid)
98    #endif
99    ! Iron
100    #ifdef FQUOTA
101               CALL TIMEAVE_NORMALIZE(AP_F_ave(1-Olx,1-Oly,1,1,1,np),
102         &                           DAR_timeave, nR ,bi,bj,myThid)
103               CALL TIMEAVE_NORMALIZE(HP_F_ave(1-Olx,1-Oly,1,1,1,np),
104         &                           DAR_timeave, nR ,bi,bj,myThid)
105               CALL TIMEAVE_NORMALIZE(Flimave(1-Olx,1-Oly,1,1,1,np),
106         &                           DAR_timeave, nR ,bi,bj,myThid)
107    #endif
108    ! Light
109               CALL TIMEAVE_NORMALIZE(Ilimave(1-Olx,1-Oly,1,1,1,np),
110         &                           DAR_timeave, nR ,bi,bj,myThid)
111           enddo           enddo
112           CALL TIMEAVE_NORMALIZ(Tlimave,DAR_timeave, nR ,  ! Temperature
113             CALL TIMEAVE_NORMALIZE(Tlimave,DAR_timeave, nR ,
114       &                                           bi,bj,myThid)       &                                           bi,bj,myThid)
115  #endif  #endif
116  #ifndef ALLOW_QUOTA  #ifndef ALLOW_QUOTA
117           CALL TIMEAVE_NORMALIZ(Chlave, DAR_timeave, nR ,           CALL TIMEAVE_NORMALIZE(Chlave, DAR_timeave, nR ,
118       &                 bi,bj,myThid)       &                 bi,bj,myThid)
119  #endif  #endif
120           CALL TIMEAVE_NORMALIZ(Nfixave, DAR_timeave, nR ,           CALL TIMEAVE_NORMALIZE(Nfixave, DAR_timeave, nR ,
121       &                 bi,bj,myThid)       &                 bi,bj,myThid)
122           CALL TIMEAVE_NORMALIZ(Denitave, DAR_timeave, nR ,           CALL TIMEAVE_NORMALIZE(Denitave, DAR_timeave, nR ,
123       &                 bi,bj,myThid)       &                 bi,bj,myThid)
124    
125  #ifdef DAR_DIAG_PARW  #ifdef DAR_DIAG_PARW
126           do ii = 1,tlam           do ii = 1,tlam
127            CALL TIMEAVE_NORMALIZ(PARwave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(PARwave(1-Olx,1-Oly,1,1,1,ii),
128       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
129           enddo           enddo
130           do np = 1,npmax           do np = 1,npmax
131            CALL TIMEAVE_NORMALIZ(chl2cave(1-Olx,1-Oly,1,1,1,np),            CALL TIMEAVE_NORMALIZE(chl2cave(1-Olx,1-Oly,1,1,1,np),
132       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
133           enddo           enddo
134  #endif  #endif
135  #ifdef DAR_DIAG_ACDOM  #ifdef DAR_DIAG_ACDOM
136           CALL TIMEAVE_NORMALIZ(aCDOMave, DAR_timeave, nR ,           CALL TIMEAVE_NORMALIZE(aCDOMave, DAR_timeave, nR ,
137       &                 bi,bj,myThid)       &                 bi,bj,myThid)
138  #endif  #endif
139  #ifdef DAR_DIAG_IRR  #ifdef DAR_DIAG_IRR
140           do ii = 1,tlam           do ii = 1,tlam
141            CALL TIMEAVE_NORMALIZ(Edave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(Edave(1-Olx,1-Oly,1,1,1,ii),
142       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
143            CALL TIMEAVE_NORMALIZ(Esave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(Esave(1-Olx,1-Oly,1,1,1,ii),
144       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
145            CALL TIMEAVE_NORMALIZ(Euave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(Euave(1-Olx,1-Oly,1,1,1,ii),
146       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
147            CALL TIMEAVE_NORMALIZ(Estave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(Estave(1-Olx,1-Oly,1,1,1,ii),
148       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
149            CALL TIMEAVE_NORMALIZ(Eutave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(Eutave(1-Olx,1-Oly,1,1,1,ii),
150       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
151           enddo           enddo
152  #endif  #endif
153  #ifdef DAR_DIAG_IRR_AMPS  #ifdef DAR_DIAG_IRR_AMPS
154           do ii = 1,tlam           do ii = 1,tlam
155            CALL TIMEAVE_NORMALIZ(amp1ave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(amp1ave(1-Olx,1-Oly,1,1,1,ii),
156       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
157            CALL TIMEAVE_NORMALIZ(amp2ave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(amp2ave(1-Olx,1-Oly,1,1,1,ii),
158       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
159           enddo           enddo
160  #endif  #endif
161  #ifdef DAR_DIAG_ABSORP  #ifdef DAR_DIAG_ABSORP
162           do ii = 1,tlam           do ii = 1,tlam
163            CALL TIMEAVE_NORMALIZ(aave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(aave(1-Olx,1-Oly,1,1,1,ii),
164       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
165           enddo           enddo
166  #endif  #endif
167  #ifdef DAR_DIAG_SCATTER  #ifdef DAR_DIAG_SCATTER
168           do ii = 1,tlam           do ii = 1,tlam
169            CALL TIMEAVE_NORMALIZ(btave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(btave(1-Olx,1-Oly,1,1,1,ii),
170       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
171            CALL TIMEAVE_NORMALIZ(bbave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(bbave(1-Olx,1-Oly,1,1,1,ii),
172       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
173           enddo           enddo
174  #endif  #endif
175  #ifdef DAR_DIAG_PART_SCATTER  #ifdef DAR_DIAG_PART_SCATTER
176           do ii = 1,tlam           do ii = 1,tlam
177            CALL TIMEAVE_NORMALIZ(apartave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(apartave(1-Olx,1-Oly,1,1,1,ii),
178       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
179            CALL TIMEAVE_NORMALIZ(btpartave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(btpartave(1-Olx,1-Oly,1,1,1,ii),
180       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
181            CALL TIMEAVE_NORMALIZ(bbpartave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(bbpartave(1-Olx,1-Oly,1,1,1,ii),
182       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
183           enddo           enddo
184  #endif  #endif
185  #ifdef DAR_RADTRANS  #ifdef DAR_RADTRANS
186           CALL TIMEAVE_NORMALIZ(rmudave(1-Olx,1-Oly,1,1),           CALL TIMEAVE_NORMALIZE(rmudave(1-Olx,1-Oly,1,1),
187       &                         DAR_timeave,1,bi,bj,myThid)       &                         DAR_timeave,1,bi,bj,myThid)
188  #endif  #endif
189  #ifdef DAR_DIAG_EK  #ifdef DAR_DIAG_EK
190           do np = 1,npmax           do np = 1,npmax
191            CALL TIMEAVE_NORMALIZ(Ekave(1-Olx,1-Oly,1,1,1,np),            CALL TIMEAVE_NORMALIZE(Ekave(1-Olx,1-Oly,1,1,1,np),
192       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
193            CALL TIMEAVE_NORMALIZ(EkoverEave(1-Olx,1-Oly,1,1,1,np),            CALL TIMEAVE_NORMALIZE(EkoverEave(1-Olx,1-Oly,1,1,1,np),
194       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
195            CALL TIMEAVE_NORMALIZ(acclimave(1-Olx,1-Oly,1,1,1,np),            CALL TIMEAVE_NORMALIZE(acclimave(1-Olx,1-Oly,1,1,1,np),
196       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
197            do ii = 1,tlam            do ii = 1,tlam
198              CALL TIMEAVE_NORMALIZ(Ek_nlave(1-Olx,1-Oly,1,1,1,np,ii),              CALL TIMEAVE_NORMALIZE(Ek_nlave(1-Olx,1-Oly,1,1,1,np,ii),
199       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
200              CALL TIMEAVE_NORMALIZ(              CALL TIMEAVE_NORMALIZE(
201       &                       EkoverE_nlave(1-Olx,1-Oly,1,1,1,np,ii),       &                       EkoverE_nlave(1-Olx,1-Oly,1,1,1,np,ii),
202       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
203            enddo            enddo
# Line 181  C      Normalize by integrated time Line 207  C      Normalize by integrated time
207  c ANNA_TAVE  c ANNA_TAVE
208  #ifdef WAVES_DIAG_PCHL  #ifdef WAVES_DIAG_PCHL
209           do np=1,npmax           do np=1,npmax
210             CALL TIMEAVE_NORMALIZ(Pchlave(1-Olx,1-Oly,1,1,1,np),             CALL TIMEAVE_NORMALIZE(Pchlave(1-Olx,1-Oly,1,1,1,np),
211       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
212       &                 bi,bj,myThid)       &                 bi,bj,myThid)
213  c           CALL TIMEAVE_NORMALIZ(Pchlave(1-Olx,1-Oly,1,1,1,np),  c           CALL TIMEAVE_NORMALIZE(Pchlave(1-Olx,1-Oly,1,1,1,np),
214  c     &                          DAR_timeave, nR ,  c     &                          DAR_timeave, nR ,
215  c     &                 bi,bj,myThid)  c     &                 bi,bj,myThid)
216           enddo           enddo
# Line 193  c ANNA end TAVE Line 219  c ANNA end TAVE
219    
220  #ifdef DAR_DIAG_RSTAR  #ifdef DAR_DIAG_RSTAR
221           do np=1,npmax           do np=1,npmax
222             CALL TIMEAVE_NORMALIZ(Rstarave(1-Olx,1-Oly,1,1,1,np),             CALL TIMEAVE_NORMALIZE(Rstarave(1-Olx,1-Oly,1,1,1,np),
223       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
224       &                 bi,bj,myThid)       &                 bi,bj,myThid)
225             CALL TIMEAVE_NORMALIZ(RNstarave(1-Olx,1-Oly,1,1,1,np),             CALL TIMEAVE_NORMALIZE(RNstarave(1-Olx,1-Oly,1,1,1,np),
226       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
227       &                 bi,bj,myThid)       &                 bi,bj,myThid)
228           enddo           enddo
229  #endif  #endif
230  #ifdef DAR_DIAG_DIVER  #ifdef DAR_DIAG_DIVER
231            CALL TIMEAVE_NORMALIZ(Diver1ave(1-Olx,1-Oly,1,1,1),            CALL TIMEAVE_NORMALIZE(Diver1ave(1-Olx,1-Oly,1,1,1),
232       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
233       &                 bi,bj,myThid)       &                 bi,bj,myThid)
234            CALL TIMEAVE_NORMALIZ(Diver2ave(1-Olx,1-Oly,1,1,1),            CALL TIMEAVE_NORMALIZE(Diver2ave(1-Olx,1-Oly,1,1,1),
235       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
236       &                 bi,bj,myThid)       &                 bi,bj,myThid)
237            CALL TIMEAVE_NORMALIZ(Diver3ave(1-Olx,1-Oly,1,1,1),            CALL TIMEAVE_NORMALIZE(Diver3ave(1-Olx,1-Oly,1,1,1),
238       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
239       &                 bi,bj,myThid)       &                 bi,bj,myThid)
240            CALL TIMEAVE_NORMALIZ(Diver4ave(1-Olx,1-Oly,1,1,1),            CALL TIMEAVE_NORMALIZE(Diver4ave(1-Olx,1-Oly,1,1,1),
241       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
242       &                 bi,bj,myThid)       &                 bi,bj,myThid)
243  #endif  #endif
244  c  c
245  #ifdef DAR_DIAG_GROW  #ifdef DAR_DIAG_GROW
246           do np=1,npmax           do np=1,npmax
247             CALL TIMEAVE_NORMALIZ(Growave(1-Olx,1-Oly,1,1,1,np),             CALL TIMEAVE_NORMALIZE(Growave(1-Olx,1-Oly,1,1,1,np),
248       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
249       &                 bi,bj,myThid)       &                 bi,bj,myThid)
250             CALL TIMEAVE_NORMALIZ(Growsqave(1-Olx,1-Oly,1,1,1,np),             CALL TIMEAVE_NORMALIZE(Growsqave(1-Olx,1-Oly,1,1,1,np),
251       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
252       &                 bi,bj,myThid)       &                 bi,bj,myThid)
253           enddo           enddo
# Line 230  c Line 256  c
256  c  c
257  #ifdef DAR_DIAG_NFIXP  #ifdef DAR_DIAG_NFIXP
258           do np=1,npmax           do np=1,npmax
259             CALL TIMEAVE_NORMALIZ(NfixPave(1-Olx,1-Oly,1,1,1,np),             CALL TIMEAVE_NORMALIZE(NfixPave(1-Olx,1-Oly,1,1,1,np),
260       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
261       &                 bi,bj,myThid)       &                 bi,bj,myThid)
262           enddo           enddo
# Line 252  c Line 278  c
278       &             Denitave, nR, Denittot, myThid )       &             Denitave, nR, Denittot, myThid )
279         Denittot=Denittot*14. _d 0 *1. _d -3* 31536000. _d 0         Denittot=Denittot*14. _d 0 *1. _d -3* 31536000. _d 0
280       &                      *1 _d -12       &                      *1 _d -12
281           _BEGIN_MASTER( myThid )
282         write(dar_ratetot_unit,*) myIter, PPtot, Nfixtot, Denittot         write(dar_ratetot_unit,*) myIter, PPtot, Nfixtot, Denittot
283           _END_MASTER( myThid )
284  #endif  #endif
285    
286          
287         IF ( PTRACERS_timeave_mdsio )  THEN         IF ( PTRACERS_timeave_mdsio )  THEN
288           WRITE(suff,'(I10.10)') myIter           WRITE(suff,'(I10.10)') myIter
289           CALL WRITE_FLD_XYZ_RL('DAR_PARtave.',suff,PARave,                   CALL WRITE_FLD_XYZ_RL('DAR_PARtave.',suff,PARave,
290       &        myIter,myThid)       &        myIter,myThid)
291           CALL WRITE_FLD_XYZ_RL('DAR_PPtave.',suff,PPave,           CALL WRITE_FLD_XYZ_RL('DAR_PPtave.',suff,PPave,
292       &        myIter,myThid)       &        myIter,myThid)
# Line 495  c Line 523  c
523       &        pf,'dar_tave',0,0,'dar_PP_ave',PPave,myThid)       &        pf,'dar_tave',0,0,'dar_PP_ave',PPave,myThid)
524  #ifdef QUOTA_DIAG_LIMIT  #ifdef QUOTA_DIAG_LIMIT
525           do np=1,npmax           do np=1,npmax
526              WRITE(name,'(a,i3.3)') 'dar_Nlim_',np  ! Carbon
527              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,             WRITE(name,'(a,i3.3)') 'dar_AP_C_',np
528       &           Nlimave(1-OLx,1-OLy,1,1,1,np),myThid)             CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
529              WRITE(name,'(a,i3.3)') 'dar_Flim_',np       &          AP_C_ave(1-OLx,1-OLy,1,1,1,np),myThid)
530              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,             WRITE(name,'(a,i3.3)') 'dar_HP_C_',np
531       &           Flimave(1-OLx,1-OLy,1,1,1,np),myThid)             CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
532              WRITE(name,'(a,i3.3)') 'dar_Ilim_',np       &          HP_C_ave(1-OLx,1-OLy,1,1,1,np),myThid)
533              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,  ! Nitrogen
534       &           Ilimave(1-OLx,1-OLy,1,1,1,np),myThid)             WRITE(name,'(a,i3.3)') 'dar_Nlim_',np
535               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
536         &          Nlimave(1-OLx,1-OLy,1,1,1,np),myThid)
537               WRITE(name,'(a,i3.3)') 'dar_AP_N_',np
538               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
539         &          AP_N_ave(1-OLx,1-OLy,1,1,1,np),myThid)
540               WRITE(name,'(a,i3.3)') 'dar_HP_N_',np
541               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
542         &          HP_N_ave(1-OLx,1-OLy,1,1,1,np),myThid)
543    ! Phosphorus
544    #ifdef PQUOTA
545               WRITE(name,'(a,i3.3)') 'dar_Plim_',np
546               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
547         &          Plimave(1-OLx,1-OLy,1,1,1,np),myThid)
548               WRITE(name,'(a,i3.3)') 'dar_AP_P_',np
549               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
550         &          AP_P_ave(1-OLx,1-OLy,1,1,1,np),myThid)
551               WRITE(name,'(a,i3.3)') 'dar_HP_P_',np
552               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
553         &          HP_P_ave(1-OLx,1-OLy,1,1,1,np),myThid)
554    #endif
555    ! Iron
556    #ifdef FQUOTA
557               WRITE(name,'(a,i3.3)') 'dar_Flim_',np
558               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
559         &          Flimave(1-OLx,1-OLy,1,1,1,np),myThid)
560               WRITE(name,'(a,i3.3)') 'dar_AP_F_',np
561               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
562         &          AP_F_ave(1-OLx,1-OLy,1,1,1,np),myThid)
563               WRITE(name,'(a,i3.3)') 'dar_HP_F_',np
564               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
565         &          HP_F_ave(1-OLx,1-OLy,1,1,1,np),myThid)
566    #endif
567    ! Light
568              WRITE(name,'(a,i3.3)') 'dar_Ilim_',np
569              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
570         &         Ilimave(1-OLx,1-OLy,1,1,1,np),myThid)
571           enddo           enddo
572    ! Temperature
573           CALL MNC_CW_RL_W(           CALL MNC_CW_RL_W(
574       &        pf,'dar_tave',0,0,'dar_Tlim_ave',Tlimave,myThid)       &        pf,'dar_tave',0,0,'dar_Tlim_ave',Tlimave,myThid)
575  #endif  #endif
# Line 737  C      Reset averages to zero Line 802  C      Reset averages to zero
802           CALL TIMEAVE_RESET(PARave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(PARave,Nr,bi,bj,myThid)
803           CALL TIMEAVE_RESET(PPave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(PPave,Nr,bi,bj,myThid)
804  #ifdef QUOTA_DIAG_LIMIT  #ifdef QUOTA_DIAG_LIMIT
805           CALL TIMEAVE_RESET(Nlimave,Nr,bi,bj,myThid)  ! Carbon
806           CALL TIMEAVE_RESET(Flimave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(AP_C_ave,Nr,bi,bj,myThid)
807             CALL TIMEAVE_RESET(HP_C_ave,Nr,bi,bj,myThid)
808    ! Nitrogen
809             CALL TIMEAVE_RESET(AP_N_ave,Nr,bi,bj,myThid)
810             CALL TIMEAVE_RESET(HP_N_ave,Nr,bi,bj,myThid)
811             CALL TIMEAVE_RESET(Nlimave ,Nr,bi,bj,myThid)
812    ! Phosphorus
813    #ifdef PQUOTA
814             CALL TIMEAVE_RESET(AP_P_ave,Nr,bi,bj,myThid)
815             CALL TIMEAVE_RESET(HP_P_ave,Nr,bi,bj,myThid)
816             CALL TIMEAVE_RESET(Plimave ,Nr,bi,bj,myThid)
817    #endif
818    ! Iron
819    #ifdef FQUOTA
820             CALL TIMEAVE_RESET(AP_F_ave,Nr,bi,bj,myThid)
821             CALL TIMEAVE_RESET(HP_F_ave,Nr,bi,bj,myThid)
822             CALL TIMEAVE_RESET(Flimave ,Nr,bi,bj,myThid)
823    #endif
824           CALL TIMEAVE_RESET(Ilimave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(Ilimave,Nr,bi,bj,myThid)
825           CALL TIMEAVE_RESET(Tlimave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(Tlimave,Nr,bi,bj,myThid)
826  #endif  #endif
# Line 826  C      Reset averages to zero Line 908  C      Reset averages to zero
908             enddo             enddo
909           enddo           enddo
910  #endif  #endif
911  c ANNA_TAVE  c ANNA_TAVE
912  #ifdef WAVES_DIAG_PCHL  #ifdef WAVES_DIAG_PCHL
913           do np=1,npmax           do np=1,npmax
914             CALL TIMEAVE_RESET(Pchlave(1-OLx,1-OLy,1,1,1,np),             CALL TIMEAVE_RESET(Pchlave(1-OLx,1-OLy,1,1,1,np),
# Line 870  c Line 952  c
952           enddo           enddo
953  #endif  #endif
954  #endif  #endif
955           do k=1,Nr           DAR_TimeAve(bi,bj)=0.d0
            DAR_TimeAve(bi,bj,k)=0.d0  
          enddo  
956          ENDDO          ENDDO
957         ENDDO         ENDDO
958            
959        ENDIF        ENDIF
960    
961  #endif /* ALLOW_TIMEAVE */  #endif /* ALLOW_TIMEAVE */
         
962  #endif /* ALLOW_DARWIN */  #endif /* ALLOW_DARWIN */
963          
964        RETURN        RETURN
965        END        END
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22