/[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.2 by stephd, Wed May 11 18:11:56 2011 UTC revision 1.9 by jahn, Fri Dec 27 17:29:00 2013 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
77             do np=1,npmax
78               CALL TIMEAVE_NORMALIZE(Nlimave(1-Olx,1-Oly,1,1,1,np),
79         &                          DAR_timeave, nR ,bi,bj,myThid)
80               CALL TIMEAVE_NORMALIZE(Flimave(1-Olx,1-Oly,1,1,1,np),
81         &                          DAR_timeave, nR ,bi,bj,myThid)
82               CALL TIMEAVE_NORMALIZE(Ilimave(1-Olx,1-Oly,1,1,1,np),
83         &                          DAR_timeave, nR ,bi,bj,myThid)
84             enddo
85             CALL TIMEAVE_NORMALIZE(Tlimave,DAR_timeave, nR ,
86         &                                           bi,bj,myThid)
87    #endif
88  #ifndef ALLOW_QUOTA  #ifndef ALLOW_QUOTA
89           CALL TIMEAVE_NORMALIZ(Chlave, DAR_timeave, nR ,           CALL TIMEAVE_NORMALIZE(Chlave, DAR_timeave, nR ,
90       &                 bi,bj,myThid)       &                 bi,bj,myThid)
91  #endif  #endif
92           CALL TIMEAVE_NORMALIZ(Nfixave, DAR_timeave, nR ,           CALL TIMEAVE_NORMALIZE(Nfixave, DAR_timeave, nR ,
93       &                 bi,bj,myThid)       &                 bi,bj,myThid)
94           CALL TIMEAVE_NORMALIZ(Denitave, DAR_timeave, nR ,           CALL TIMEAVE_NORMALIZE(Denitave, DAR_timeave, nR ,
95       &                 bi,bj,myThid)       &                 bi,bj,myThid)
96    
97    #ifdef DAR_DIAG_PARW
98             do ii = 1,tlam
99              CALL TIMEAVE_NORMALIZE(PARwave(1-Olx,1-Oly,1,1,1,ii),
100         &                          DAR_timeave,nR,bi,bj,myThid)
101             enddo
102             do np = 1,npmax
103              CALL TIMEAVE_NORMALIZE(chl2cave(1-Olx,1-Oly,1,1,1,np),
104         &                          DAR_timeave,nR,bi,bj,myThid)
105             enddo
106    #endif
107  #ifdef DAR_DIAG_ACDOM  #ifdef DAR_DIAG_ACDOM
108           CALL TIMEAVE_NORMALIZ(aCDOMave, DAR_timeave, nR ,           CALL TIMEAVE_NORMALIZE(aCDOMave, DAR_timeave, nR ,
109       &                 bi,bj,myThid)       &                 bi,bj,myThid)
110  #endif  #endif
111  #ifdef DAR_DIAG_IRR  #ifdef DAR_DIAG_IRR
112           do ii = 1,tlam           do ii = 1,tlam
113            CALL TIMEAVE_NORMALIZ(Edave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(Edave(1-Olx,1-Oly,1,1,1,ii),
114       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
115            CALL TIMEAVE_NORMALIZ(Esave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(Esave(1-Olx,1-Oly,1,1,1,ii),
116       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
117            CALL TIMEAVE_NORMALIZ(Euave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(Euave(1-Olx,1-Oly,1,1,1,ii),
118       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
119            CALL TIMEAVE_NORMALIZ(Eutave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(Estave(1-Olx,1-Oly,1,1,1,ii),
120         &                          DAR_timeave,nR,bi,bj,myThid)
121              CALL TIMEAVE_NORMALIZE(Eutave(1-Olx,1-Oly,1,1,1,ii),
122         &                          DAR_timeave,nR,bi,bj,myThid)
123             enddo
124    #endif
125    #ifdef DAR_DIAG_IRR_AMPS
126             do ii = 1,tlam
127              CALL TIMEAVE_NORMALIZE(amp1ave(1-Olx,1-Oly,1,1,1,ii),
128         &                          DAR_timeave,nR,bi,bj,myThid)
129              CALL TIMEAVE_NORMALIZE(amp2ave(1-Olx,1-Oly,1,1,1,ii),
130       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
131           enddo           enddo
132  #endif  #endif
133  #ifdef DAR_DIAG_ABSORP  #ifdef DAR_DIAG_ABSORP
134           do ii = 1,tlam           do ii = 1,tlam
135            CALL TIMEAVE_NORMALIZ(aave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(aave(1-Olx,1-Oly,1,1,1,ii),
136       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
137           enddo           enddo
138  #endif  #endif
139  #ifdef DAR_DIAG_SCATTER  #ifdef DAR_DIAG_SCATTER
140           do ii = 1,tlam           do ii = 1,tlam
141            CALL TIMEAVE_NORMALIZ(btave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(btave(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(bbave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(bbave(1-Olx,1-Oly,1,1,1,ii),
144       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
145           enddo           enddo
146  #endif  #endif
147  #ifdef DAR_DIAG_PART_SCATTER  #ifdef DAR_DIAG_PART_SCATTER
148           do ii = 1,tlam           do ii = 1,tlam
149            CALL TIMEAVE_NORMALIZ(apartave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(apartave(1-Olx,1-Oly,1,1,1,ii),
150       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
151            CALL TIMEAVE_NORMALIZ(btpartave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(btpartave(1-Olx,1-Oly,1,1,1,ii),
152       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
153            CALL TIMEAVE_NORMALIZ(bbpartave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZE(bbpartave(1-Olx,1-Oly,1,1,1,ii),
154       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
155           enddo           enddo
156  #endif  #endif
157    #ifdef DAR_RADTRANS
158             CALL TIMEAVE_NORMALIZE(rmudave(1-Olx,1-Oly,1,1),
159         &                         DAR_timeave,1,bi,bj,myThid)
160    #endif
161    #ifdef DAR_DIAG_EK
162             do np = 1,npmax
163              CALL TIMEAVE_NORMALIZE(Ekave(1-Olx,1-Oly,1,1,1,np),
164         &                          DAR_timeave,nR,bi,bj,myThid)
165              CALL TIMEAVE_NORMALIZE(EkoverEave(1-Olx,1-Oly,1,1,1,np),
166         &                          DAR_timeave,nR,bi,bj,myThid)
167              CALL TIMEAVE_NORMALIZE(acclimave(1-Olx,1-Oly,1,1,1,np),
168         &                          DAR_timeave,nR,bi,bj,myThid)
169              do ii = 1,tlam
170                CALL TIMEAVE_NORMALIZE(Ek_nlave(1-Olx,1-Oly,1,1,1,np,ii),
171         &                          DAR_timeave,nR,bi,bj,myThid)
172                CALL TIMEAVE_NORMALIZE(
173         &                       EkoverE_nlave(1-Olx,1-Oly,1,1,1,np,ii),
174         &                          DAR_timeave,nR,bi,bj,myThid)
175              enddo
176             enddo
177    #endif
178    
179  c ANNA_TAVE  c ANNA_TAVE
180  #ifdef WAVES_DIAG_PCHL  #ifdef WAVES_DIAG_PCHL
181           do np=1,npmax           do np=1,npmax
182             CALL TIMEAVE_NORMALIZ(Pchlave(1-Olx,1-Oly,1,1,1,np),             CALL TIMEAVE_NORMALIZE(Pchlave(1-Olx,1-Oly,1,1,1,np),
183       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
184       &                 bi,bj,myThid)       &                 bi,bj,myThid)
185  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),
186  c     &                          DAR_timeave, nR ,  c     &                          DAR_timeave, nR ,
187  c     &                 bi,bj,myThid)  c     &                 bi,bj,myThid)
188           enddo           enddo
# Line 140  c ANNA end TAVE Line 191  c ANNA end TAVE
191    
192  #ifdef DAR_DIAG_RSTAR  #ifdef DAR_DIAG_RSTAR
193           do np=1,npmax           do np=1,npmax
194             CALL TIMEAVE_NORMALIZ(Rstarave(1-Olx,1-Oly,1,1,1,np),             CALL TIMEAVE_NORMALIZE(Rstarave(1-Olx,1-Oly,1,1,1,np),
195       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
196       &                 bi,bj,myThid)       &                 bi,bj,myThid)
197             CALL TIMEAVE_NORMALIZ(RNstarave(1-Olx,1-Oly,1,1,1,np),             CALL TIMEAVE_NORMALIZE(RNstarave(1-Olx,1-Oly,1,1,1,np),
198       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
199       &                 bi,bj,myThid)       &                 bi,bj,myThid)
200           enddo           enddo
201  #endif  #endif
202  #ifdef DAR_DIAG_DIVER  #ifdef DAR_DIAG_DIVER
203            CALL TIMEAVE_NORMALIZ(Diver1ave(1-Olx,1-Oly,1,1,1),            CALL TIMEAVE_NORMALIZE(Diver1ave(1-Olx,1-Oly,1,1,1),
204       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
205       &                 bi,bj,myThid)       &                 bi,bj,myThid)
206            CALL TIMEAVE_NORMALIZ(Diver2ave(1-Olx,1-Oly,1,1,1),            CALL TIMEAVE_NORMALIZE(Diver2ave(1-Olx,1-Oly,1,1,1),
207       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
208       &                 bi,bj,myThid)       &                 bi,bj,myThid)
209            CALL TIMEAVE_NORMALIZ(Diver3ave(1-Olx,1-Oly,1,1,1),            CALL TIMEAVE_NORMALIZE(Diver3ave(1-Olx,1-Oly,1,1,1),
210       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
211       &                 bi,bj,myThid)       &                 bi,bj,myThid)
212            CALL TIMEAVE_NORMALIZ(Diver4ave(1-Olx,1-Oly,1,1,1),            CALL TIMEAVE_NORMALIZE(Diver4ave(1-Olx,1-Oly,1,1,1),
213       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
214       &                 bi,bj,myThid)       &                 bi,bj,myThid)
215  #endif  #endif
216  c  c
217  #ifdef DAR_DIAG_GROW  #ifdef DAR_DIAG_GROW
218           do np=1,npmax           do np=1,npmax
219             CALL TIMEAVE_NORMALIZ(Growave(1-Olx,1-Oly,1,1,1,np),             CALL TIMEAVE_NORMALIZE(Growave(1-Olx,1-Oly,1,1,1,np),
220       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
221       &                 bi,bj,myThid)       &                 bi,bj,myThid)
222             CALL TIMEAVE_NORMALIZ(Growsqave(1-Olx,1-Oly,1,1,1,np),             CALL TIMEAVE_NORMALIZE(Growsqave(1-Olx,1-Oly,1,1,1,np),
223       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
224       &                 bi,bj,myThid)       &                 bi,bj,myThid)
225           enddo           enddo
# Line 177  c Line 228  c
228  c  c
229  #ifdef DAR_DIAG_NFIXP  #ifdef DAR_DIAG_NFIXP
230           do np=1,npmax           do np=1,npmax
231             CALL TIMEAVE_NORMALIZ(NfixPave(1-Olx,1-Oly,1,1,1,np),             CALL TIMEAVE_NORMALIZE(NfixPave(1-Olx,1-Oly,1,1,1,np),
232       &                          DAR_timeave, nR ,       &                          DAR_timeave, nR ,
233       &                 bi,bj,myThid)       &                 bi,bj,myThid)
234           enddo           enddo
# Line 190  c Line 241  c
241         call tracer_sumvol(         call tracer_sumvol(
242       &             PPave, nR, PPtot, myThid )       &             PPave, nR, PPtot, myThid )
243         PPtot=PPtot*120. _d 0 * 12. _d 0* 1. _d -3* 31536000. _d 0         PPtot=PPtot*120. _d 0 * 12. _d 0* 1. _d -3* 31536000. _d 0
244       &                      *1. _d-15       &                      *1. _d -15
245         call tracer_sumvol(         call tracer_sumvol(
246       &             Nfixave, nR, Nfixtot, myThid )       &             Nfixave, nR, Nfixtot, myThid )
247         Nfixtot=Nfixtot*14. _d 0 *1. _d -3* 31536000. _d 0         Nfixtot=Nfixtot*14. _d 0 *1. _d -3* 31536000. _d 0
# Line 199  c Line 250  c
250       &             Denitave, nR, Denittot, myThid )       &             Denitave, nR, Denittot, myThid )
251         Denittot=Denittot*14. _d 0 *1. _d -3* 31536000. _d 0         Denittot=Denittot*14. _d 0 *1. _d -3* 31536000. _d 0
252       &                      *1 _d -12       &                      *1 _d -12
253           _BEGIN_MASTER( myThid )
254         write(dar_ratetot_unit,*) myIter, PPtot, Nfixtot, Denittot         write(dar_ratetot_unit,*) myIter, PPtot, Nfixtot, Denittot
255           _END_MASTER( myThid )
256  #endif  #endif
257    
258          
259         IF ( PTRACERS_timeave_mdsio )  THEN         IF ( PTRACERS_timeave_mdsio )  THEN
260           WRITE(suff,'(I10.10)') myIter           WRITE(suff,'(I10.10)') myIter
261           CALL WRITE_FLD_XYZ_RL('DAR_PARtave.',suff,PARave,                   CALL WRITE_FLD_XYZ_RL('DAR_PARtave.',suff,PARave,
262       &        myIter,myThid)       &        myIter,myThid)
263           CALL WRITE_FLD_XYZ_RL('DAR_PPtave.',suff,PPave,           CALL WRITE_FLD_XYZ_RL('DAR_PPtave.',suff,PPave,
264       &        myIter,myThid)       &        myIter,myThid)
# Line 217  c Line 270  c
270       &        myIter,myThid)       &        myIter,myThid)
271           CALL WRITE_FLD_XYZ_RL('DAR_Denittave.',suff,Denitave,           CALL WRITE_FLD_XYZ_RL('DAR_Denittave.',suff,Denitave,
272       &        myIter,myThid)       &        myIter,myThid)
273    #ifdef DAR_DIAG_PARW
274             do ii = 1,tlam
275               WRITE(suff,'(A8,I2.2,A1,I10.10)') 'PARwave',
276         &           ii,'.',myIter
277               CALL WRITE_FLD_XYZ_RL(suff,' ',
278         &           PARwave(1-Olx,1-Oly,1,1,1,ii),
279         &           myIter,myThid)
280             enddo
281             do np = 1,npmax
282               WRITE(suff,'(A8,I2.2,A1,I10.10)') 'chl2cave',
283         &           np,'.',myIter
284               CALL WRITE_FLD_XYZ_RL(suff,' ',
285         &           chl2cave(1-Olx,1-Oly,1,1,1,np),
286         &           myIter,myThid)
287             enddo
288    #endif
289  #ifdef DAR_DIAG_ACDOM  #ifdef DAR_DIAG_ACDOM
290           CALL WRITE_FLD_XYZ_RL('DAR_aCDOMtave',suff,aCDOMave,           CALL WRITE_FLD_XYZ_RL('DAR_aCDOMtave',suff,aCDOMave,
291       &        myIter,myThid)       &        myIter,myThid)
# Line 238  c Line 307  c
307             CALL WRITE_FLD_XYZ_RL(suff,' ',             CALL WRITE_FLD_XYZ_RL(suff,' ',
308       &           Euave(1-Olx,1-Oly,1,1,1,ii),       &           Euave(1-Olx,1-Oly,1,1,1,ii),
309       &           myIter,myThid)       &           myIter,myThid)
310               WRITE(suff,'(A8,I2.2,A1,I10.10)') 'Estopave',
311         &           ii,'.',myIter
312               CALL WRITE_FLD_XYZ_RL(suff,' ',
313         &           Estave(1-Olx,1-Oly,1,1,1,ii),
314         &           myIter,myThid)
315             WRITE(suff,'(A8,I2.2,A1,I10.10)') 'Eutopave',             WRITE(suff,'(A8,I2.2,A1,I10.10)') 'Eutopave',
316       &           ii,'.',myIter       &           ii,'.',myIter
317             CALL WRITE_FLD_XYZ_RL(suff,' ',             CALL WRITE_FLD_XYZ_RL(suff,' ',
# Line 245  c Line 319  c
319       &           myIter,myThid)       &           myIter,myThid)
320           enddo           enddo
321  #endif  #endif
322    #ifdef DAR_DIAG_IRR_AMPS
323             do ii = 1,tlam
324               WRITE(suff,'(A5,I2.2,A1,I10.10)') 'amp1ave',
325         &           ii,'.',myIter
326               CALL WRITE_FLD_XYZ_RL(suff,' ',
327         &           amp1ave(1-Olx,1-Oly,1,1,1,ii),
328         &           myIter,myThid)
329               WRITE(suff,'(A5,I2.2,A1,I10.10)') 'amp2ave',
330         &           ii,'.',myIter
331               CALL WRITE_FLD_XYZ_RL(suff,' ',
332         &           amp2ave(1-Olx,1-Oly,1,1,1,ii),
333         &           myIter,myThid)
334             enddo
335    #endif
336  #ifdef DAR_DIAG_ABSORP  #ifdef DAR_DIAG_ABSORP
337           do ii = 1,tlam           do ii = 1,tlam
338             WRITE(suff,'(A4,I2.2,A1,I10.10)') 'aave',             WRITE(suff,'(A4,I2.2,A1,I10.10)') 'aave',
# Line 287  c Line 375  c
375       &           myIter,myThid)       &           myIter,myThid)
376           enddo           enddo
377  #endif  #endif
378    #ifdef DAR_RADTRANS
379             WRITE(suff,'(A7,A1,I10.10)') 'rmudave',
380         &         '.',myIter
381             CALL WRITE_FLD_XY_RL(suff,' ',
382         &         rmudave(1-Olx,1-Oly,1,1),
383         &         myIter,myThid)
384    #endif
385    #ifdef DAR_DIAG_EK
386             do np = 1,npmax
387               WRITE(suff,'(A8,I2.2,A1,I10.10)') 'Ekave',
388         &           np,'.',myIter
389               CALL WRITE_FLD_XYZ_RL(suff,' ',
390         &           Ekave(1-Olx,1-Oly,1,1,1,np),
391         &           myIter,myThid)
392               WRITE(suff,'(A8,I2.2,A1,I10.10)') 'EkoverEave',
393         &           np,'.',myIter
394               CALL WRITE_FLD_XYZ_RL(suff,' ',
395         &           EkoverEave(1-Olx,1-Oly,1,1,1,np),
396         &           myIter,myThid)
397               WRITE(suff,'(A8,I2.2,A1,I10.10)') 'acclimave',
398         &           np,'.',myIter
399               CALL WRITE_FLD_XYZ_RL(suff,' ',
400         &           acclimave(1-Olx,1-Oly,1,1,1,np),
401         &           myIter,myThid)
402               do ii = 1,tlam
403                 WRITE(suff,'(A8,I2.2,A1,I2.2,A1,I10.10)') 'Ek_nl_ave',
404         &             np,'_',ii,'.',myIter
405                 CALL WRITE_FLD_XYZ_RL(suff,' ',
406         &             Ek_nlave(1-Olx,1-Oly,1,1,1,np,ii),
407         &             myIter,myThid)
408                 WRITE(suff,'(A8,I2.2,A1,I2.2,A1,I10.10)') 'EkoverE_nl_ave',
409         &             np,'_',ii,'.',myIter
410                 CALL WRITE_FLD_XYZ_RL(suff,' ',
411         &             EkoverE_nlave(1-Olx,1-Oly,1,1,1,np,ii),
412         &             myIter,myThid)
413               enddo
414             enddo
415    #endif
416    
417  c ANNA_TAVE  c ANNA_TAVE
418  #ifdef WAVES_DIAG_PCHL  #ifdef WAVES_DIAG_PCHL
# Line 367  c Line 493  c
493       &        pf,'dar_tave',0,0,'dar_PAR_ave',PARave,myThid)       &        pf,'dar_tave',0,0,'dar_PAR_ave',PARave,myThid)
494           CALL MNC_CW_RL_W(           CALL MNC_CW_RL_W(
495       &        pf,'dar_tave',0,0,'dar_PP_ave',PPave,myThid)       &        pf,'dar_tave',0,0,'dar_PP_ave',PPave,myThid)
496    #ifdef QUOTA_DIAG_LIMIT
497             do np=1,npmax
498                WRITE(name,'(a,i3.3)') 'dar_Nlim_',np
499                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
500         &           Nlimave(1-OLx,1-OLy,1,1,1,np),myThid)
501                WRITE(name,'(a,i3.3)') 'dar_Flim_',np
502                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
503         &           Flimave(1-OLx,1-OLy,1,1,1,np),myThid)
504                WRITE(name,'(a,i3.3)') 'dar_Ilim_',np
505                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
506         &           Ilimave(1-OLx,1-OLy,1,1,1,np),myThid)
507             enddo
508             CALL MNC_CW_RL_W(
509         &        pf,'dar_tave',0,0,'dar_Tlim_ave',Tlimave,myThid)
510    #endif
511  #ifndef ALLOW_QUOTA  #ifndef ALLOW_QUOTA
512           CALL MNC_CW_RL_W(           CALL MNC_CW_RL_W(
513       &        pf,'dar_tave',0,0,'dar_Chl_ave',Chlave,myThid)       &        pf,'dar_tave',0,0,'dar_Chl_ave',Chlave,myThid)
# Line 375  c Line 516  c
516       &        pf,'dar_tave',0,0,'dar_Nfix_ave',Nfixave,myThid)       &        pf,'dar_tave',0,0,'dar_Nfix_ave',Nfixave,myThid)
517           CALL MNC_CW_RL_W(           CALL MNC_CW_RL_W(
518       &        pf,'dar_tave',0,0,'dar_Denit_ave',Denitave,myThid)       &        pf,'dar_tave',0,0,'dar_Denit_ave',Denitave,myThid)
519    #ifdef DAR_DIAG_PARW
520             do ii = 1,tlam
521                DO jj = 1,MAX_LEN_FNAM
522                  name(jj:jj) = ' '
523                ENDDO
524                WRITE(name,'(a,i2.2)') 'dar_PARw',ii
525                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
526         &           PARwave(1-OLx,1-OLy,1,1,1,ii),myThid)
527            enddo
528             do np = 1,npmax
529                DO jj = 1,MAX_LEN_FNAM
530                  name(jj:jj) = ' '
531                ENDDO
532                WRITE(name,'(a,i2.2)') 'dar_chl2c',np
533                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
534         &           chl2cave(1-OLx,1-OLy,1,1,1,np),myThid)
535            enddo
536    #endif
537  #ifdef DAR_DIAG_ACDOM  #ifdef DAR_DIAG_ACDOM
538           CALL MNC_CW_RL_W(           CALL MNC_CW_RL_W(
539       &        pf,'dar_tave',0,0,'dar_aCDOM_ave',aCDOMave,myThid)       &        pf,'dar_tave',0,0,'dar_aCDOM_ave',aCDOMave,myThid)
# Line 393  c Line 552  c
552              WRITE(name,'(a,i2.2)') 'dar_Eu',ii              WRITE(name,'(a,i2.2)') 'dar_Eu',ii
553              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
554       &           Euave(1-OLx,1-OLy,1,1,1,ii),myThid)       &           Euave(1-OLx,1-OLy,1,1,1,ii),myThid)
555                WRITE(name,'(a,i2.2)') 'dar_Estop',ii
556                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
557         &           Estave(1-OLx,1-OLy,1,1,1,ii),myThid)
558              WRITE(name,'(a,i2.2)') 'dar_Eutop',ii              WRITE(name,'(a,i2.2)') 'dar_Eutop',ii
559              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
560       &           Eutave(1-OLx,1-OLy,1,1,1,ii),myThid)       &           Eutave(1-OLx,1-OLy,1,1,1,ii),myThid)
561           enddo           enddo
562  #endif  #endif
563    #ifdef DAR_DIAG_IRR_AMPS
564             do ii = 1,tlam
565                DO jj = 1,MAX_LEN_FNAM
566                  name(jj:jj) = ' '
567                ENDDO
568                WRITE(name,'(a,i2.2)') 'dar_amp1_',ii
569                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
570         &           amp1ave(1-OLx,1-OLy,1,1,1,ii),myThid)
571                WRITE(name,'(a,i2.2)') 'dar_amp2_',ii
572                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
573         &           amp2ave(1-OLx,1-OLy,1,1,1,ii),myThid)
574             enddo
575    #endif
576  #ifdef DAR_DIAG_ABSORP  #ifdef DAR_DIAG_ABSORP
577           do ii = 1,tlam           do ii = 1,tlam
578              DO jj = 1,MAX_LEN_FNAM              DO jj = 1,MAX_LEN_FNAM
# Line 437  c Line 612  c
612       &           bbpartave(1-OLx,1-OLy,1,1,1,ii),myThid)       &           bbpartave(1-OLx,1-OLy,1,1,1,ii),myThid)
613           enddo           enddo
614  #endif  #endif
615    #ifdef DAR_RADTRANS
616             DO jj = 1,MAX_LEN_FNAM
617               name(jj:jj) = ' '
618             ENDDO
619             WRITE(name,'(a)') 'dar_rmud'
620             CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
621         &        rmudave(1-OLx,1-OLy,1,1),myThid)
622    #endif
623    #ifdef DAR_DIAG_EK
624             do np = 1,npmax
625                DO jj = 1,MAX_LEN_FNAM
626                  name(jj:jj) = ' '
627                ENDDO
628                WRITE(name,'(a,i2.2)') 'dar_Ek',np
629                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
630         &           Ekave(1-OLx,1-OLy,1,1,1,np),myThid)
631                WRITE(name,'(a,i2.2)') 'dar_EkoverE',np
632                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
633         &           EkoverEave(1-OLx,1-OLy,1,1,1,np),myThid)
634                WRITE(name,'(a,i2.2)') 'dar_acclim',np
635                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
636         &           acclimave(1-OLx,1-OLy,1,1,1,np),myThid)
637                do ii = 1,tlam
638                  WRITE(name,'(a,i2.2,a,i2.2)') 'dar_Ek_nl',np,'_',ii
639                  CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
640         &             Ek_nlave(1-OLx,1-OLy,1,1,1,np,ii),myThid)
641                  WRITE(name,'(a,i2.2,a,i2.2)') 'dar_EkoverE_nl',np,'_',ii
642                  CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
643         &             EkoverE_nlave(1-OLx,1-OLy,1,1,1,np,ii),myThid)
644                enddo
645            enddo
646    #endif
647  c ANNA_TAVE  c ANNA_TAVE
648  #ifdef WAVES_DIAG_PCHL  #ifdef WAVES_DIAG_PCHL
649           do np=1,npmax           do np=1,npmax
# Line 529  C      Reset averages to zero Line 736  C      Reset averages to zero
736          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
737           CALL TIMEAVE_RESET(PARave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(PARave,Nr,bi,bj,myThid)
738           CALL TIMEAVE_RESET(PPave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(PPave,Nr,bi,bj,myThid)
739    #ifdef QUOTA_DIAG_LIMIT
740             CALL TIMEAVE_RESET(Nlimave,Nr,bi,bj,myThid)
741             CALL TIMEAVE_RESET(Flimave,Nr,bi,bj,myThid)
742             CALL TIMEAVE_RESET(Ilimave,Nr,bi,bj,myThid)
743             CALL TIMEAVE_RESET(Tlimave,Nr,bi,bj,myThid)
744    #endif
745  #ifndef ALLOW_QUOTA  #ifndef ALLOW_QUOTA
746           CALL TIMEAVE_RESET(Chlave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(Chlave,Nr,bi,bj,myThid)
747  #endif  #endif
748           CALL TIMEAVE_RESET(Nfixave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(Nfixave,Nr,bi,bj,myThid)
749           CALL TIMEAVE_RESET(Denitave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(Denitave,Nr,bi,bj,myThid)
750    #ifdef DAR_DIAG_PARW
751             do ii = 1,tlam
752               CALL TIMEAVE_RESET(PARwave(1-OLx,1-OLy,1,1,1,ii),
753         &                                    Nr,bi,bj,myThid)
754             enddo
755             do np = 1,npmax
756               CALL TIMEAVE_RESET(chl2cave(1-OLx,1-OLy,1,1,1,np),
757         &                                    Nr,bi,bj,myThid)
758             enddo
759    #endif
760  #ifdef DAR_DIAG_ACDOM  #ifdef DAR_DIAG_ACDOM
761           CALL TIMEAVE_RESET(aCDOMave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(aCDOMave,Nr,bi,bj,myThid)
762  #endif  #endif
# Line 545  C      Reset averages to zero Line 768  C      Reset averages to zero
768       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
769             CALL TIMEAVE_RESET(Euave(1-OLx,1-OLy,1,1,1,ii),             CALL TIMEAVE_RESET(Euave(1-OLx,1-OLy,1,1,1,ii),
770       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
771               CALL TIMEAVE_RESET(Estave(1-OLx,1-OLy,1,1,1,ii),
772         &                                    Nr,bi,bj,myThid)
773             CALL TIMEAVE_RESET(Eutave(1-OLx,1-OLy,1,1,1,ii),             CALL TIMEAVE_RESET(Eutave(1-OLx,1-OLy,1,1,1,ii),
774       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
775           enddo           enddo
776  #endif  #endif
777    #ifdef DAR_DIAG_IRR_AMPS
778             do ii = 1,tlam
779               CALL TIMEAVE_RESET(amp1ave(1-OLx,1-OLy,1,1,1,ii),
780         &                                    Nr,bi,bj,myThid)
781               CALL TIMEAVE_RESET(amp2ave(1-OLx,1-OLy,1,1,1,ii),
782         &                                    Nr,bi,bj,myThid)
783             enddo
784    #endif
785  #ifdef DAR_DIAG_ABSORP  #ifdef DAR_DIAG_ABSORP
786           do ii = 1,tlam           do ii = 1,tlam
787             CALL TIMEAVE_RESET(aave(1-OLx,1-OLy,1,1,1,ii),             CALL TIMEAVE_RESET(aave(1-OLx,1-OLy,1,1,1,ii),
# Line 573  C      Reset averages to zero Line 806  C      Reset averages to zero
806       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
807           enddo           enddo
808  #endif  #endif
809  c ANNA_TAVE  #ifdef DAR_RADTRANS
810             CALL TIMEAVE_RESET(rmudave(1-OLx,1-OLy,1,1),
811         &                                  1,bi,bj,myThid)
812    #endif
813    #ifdef DAR_DIAG_EK
814             do np = 1,npmax
815               CALL TIMEAVE_RESET(Ekave(1-OLx,1-OLy,1,1,1,np),
816         &                                    Nr,bi,bj,myThid)
817               CALL TIMEAVE_RESET(EkoverEave(1-OLx,1-OLy,1,1,1,np),
818         &                                    Nr,bi,bj,myThid)
819               CALL TIMEAVE_RESET(acclimave(1-OLx,1-OLy,1,1,1,np),
820         &                                    Nr,bi,bj,myThid)
821               do ii = 1,tlam
822                 CALL TIMEAVE_RESET(Ek_nlave(1-OLx,1-OLy,1,1,1,np,ii),
823         &                                    Nr,bi,bj,myThid)
824                 CALL TIMEAVE_RESET(EkoverE_nlave(1-OLx,1-OLy,1,1,1,np,ii),
825         &                                    Nr,bi,bj,myThid)
826               enddo
827             enddo
828    #endif
829    c ANNA_TAVE
830  #ifdef WAVES_DIAG_PCHL  #ifdef WAVES_DIAG_PCHL
831           do np=1,npmax           do np=1,npmax
832             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 617  c Line 870  c
870           enddo           enddo
871  #endif  #endif
872  #endif  #endif
873           do k=1,Nr           DAR_TimeAve(bi,bj)=0.d0
            DAR_TimeAve(bi,bj,k)=0.d0  
          enddo  
874          ENDDO          ENDDO
875         ENDDO         ENDDO
876            
877        ENDIF        ENDIF
878    
879  #endif /* ALLOW_TIMEAVE */  #endif /* ALLOW_TIMEAVE */
         
880  #endif /* ALLOW_DARWIN */  #endif /* ALLOW_DARWIN */
881          
882        RETURN        RETURN
883        END        END
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22