/[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.8 by jahn, Wed Dec 4 21:19:12 2013 UTC
# Line 75  C      Normalize by integrated time Line 75  C      Normalize by integrated time
75       &                 bi,bj,myThid)       &                 bi,bj,myThid)
76           CALL TIMEAVE_NORMALIZ(PPave, DAR_timeave, nR ,           CALL TIMEAVE_NORMALIZ(PPave, DAR_timeave, nR ,
77       &                 bi,bj,myThid)       &                 bi,bj,myThid)
78    #ifdef QUOTA_DIAG_LIMIT
79             do np=1,npmax
80               CALL TIMEAVE_NORMALIZ(Nlimave(1-Olx,1-Oly,1,1,1,np),
81         &                          DAR_timeave, nR ,bi,bj,myThid)
82               CALL TIMEAVE_NORMALIZ(Flimave(1-Olx,1-Oly,1,1,1,np),
83         &                          DAR_timeave, nR ,bi,bj,myThid)
84               CALL TIMEAVE_NORMALIZ(Ilimave(1-Olx,1-Oly,1,1,1,np),
85         &                          DAR_timeave, nR ,bi,bj,myThid)
86             enddo
87             CALL TIMEAVE_NORMALIZ(Tlimave,DAR_timeave, nR ,
88         &                                           bi,bj,myThid)
89    #endif
90  #ifndef ALLOW_QUOTA  #ifndef ALLOW_QUOTA
91           CALL TIMEAVE_NORMALIZ(Chlave, DAR_timeave, nR ,           CALL TIMEAVE_NORMALIZ(Chlave, DAR_timeave, nR ,
92       &                 bi,bj,myThid)       &                 bi,bj,myThid)
# Line 84  C      Normalize by integrated time Line 96  C      Normalize by integrated time
96           CALL TIMEAVE_NORMALIZ(Denitave, DAR_timeave, nR ,           CALL TIMEAVE_NORMALIZ(Denitave, DAR_timeave, nR ,
97       &                 bi,bj,myThid)       &                 bi,bj,myThid)
98    
99    #ifdef DAR_DIAG_PARW
100             do ii = 1,tlam
101              CALL TIMEAVE_NORMALIZ(PARwave(1-Olx,1-Oly,1,1,1,ii),
102         &                          DAR_timeave,nR,bi,bj,myThid)
103             enddo
104             do np = 1,npmax
105              CALL TIMEAVE_NORMALIZ(chl2cave(1-Olx,1-Oly,1,1,1,np),
106         &                          DAR_timeave,nR,bi,bj,myThid)
107             enddo
108    #endif
109  #ifdef DAR_DIAG_ACDOM  #ifdef DAR_DIAG_ACDOM
110           CALL TIMEAVE_NORMALIZ(aCDOMave, DAR_timeave, nR ,           CALL TIMEAVE_NORMALIZ(aCDOMave, DAR_timeave, nR ,
111       &                 bi,bj,myThid)       &                 bi,bj,myThid)
# Line 96  C      Normalize by integrated time Line 118  C      Normalize by integrated time
118       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
119            CALL TIMEAVE_NORMALIZ(Euave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZ(Euave(1-Olx,1-Oly,1,1,1,ii),
120       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
121              CALL TIMEAVE_NORMALIZ(Estave(1-Olx,1-Oly,1,1,1,ii),
122         &                          DAR_timeave,nR,bi,bj,myThid)
123            CALL TIMEAVE_NORMALIZ(Eutave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZ(Eutave(1-Olx,1-Oly,1,1,1,ii),
124       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
125           enddo           enddo
126  #endif  #endif
127    #ifdef DAR_DIAG_IRR_AMPS
128             do ii = 1,tlam
129              CALL TIMEAVE_NORMALIZ(amp1ave(1-Olx,1-Oly,1,1,1,ii),
130         &                          DAR_timeave,nR,bi,bj,myThid)
131              CALL TIMEAVE_NORMALIZ(amp2ave(1-Olx,1-Oly,1,1,1,ii),
132         &                          DAR_timeave,nR,bi,bj,myThid)
133             enddo
134    #endif
135  #ifdef DAR_DIAG_ABSORP  #ifdef DAR_DIAG_ABSORP
136           do ii = 1,tlam           do ii = 1,tlam
137            CALL TIMEAVE_NORMALIZ(aave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZ(aave(1-Olx,1-Oly,1,1,1,ii),
# Line 124  C      Normalize by integrated time Line 156  C      Normalize by integrated time
156       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
157           enddo           enddo
158  #endif  #endif
159    #ifdef DAR_RADTRANS
160             CALL TIMEAVE_NORMALIZ(rmudave(1-Olx,1-Oly,1,1),
161         &                         DAR_timeave,1,bi,bj,myThid)
162    #endif
163    #ifdef DAR_DIAG_EK
164             do np = 1,npmax
165              CALL TIMEAVE_NORMALIZ(Ekave(1-Olx,1-Oly,1,1,1,np),
166         &                          DAR_timeave,nR,bi,bj,myThid)
167              CALL TIMEAVE_NORMALIZ(EkoverEave(1-Olx,1-Oly,1,1,1,np),
168         &                          DAR_timeave,nR,bi,bj,myThid)
169              CALL TIMEAVE_NORMALIZ(acclimave(1-Olx,1-Oly,1,1,1,np),
170         &                          DAR_timeave,nR,bi,bj,myThid)
171              do ii = 1,tlam
172                CALL TIMEAVE_NORMALIZ(Ek_nlave(1-Olx,1-Oly,1,1,1,np,ii),
173         &                          DAR_timeave,nR,bi,bj,myThid)
174                CALL TIMEAVE_NORMALIZ(
175         &                       EkoverE_nlave(1-Olx,1-Oly,1,1,1,np,ii),
176         &                          DAR_timeave,nR,bi,bj,myThid)
177              enddo
178             enddo
179    #endif
180    
181  c ANNA_TAVE  c ANNA_TAVE
182  #ifdef WAVES_DIAG_PCHL  #ifdef WAVES_DIAG_PCHL
# Line 190  c Line 243  c
243         call tracer_sumvol(         call tracer_sumvol(
244       &             PPave, nR, PPtot, myThid )       &             PPave, nR, PPtot, myThid )
245         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
246       &                      *1. _d-15       &                      *1. _d -15
247         call tracer_sumvol(         call tracer_sumvol(
248       &             Nfixave, nR, Nfixtot, myThid )       &             Nfixave, nR, Nfixtot, myThid )
249         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 252  c
252       &             Denitave, nR, Denittot, myThid )       &             Denitave, nR, Denittot, myThid )
253         Denittot=Denittot*14. _d 0 *1. _d -3* 31536000. _d 0         Denittot=Denittot*14. _d 0 *1. _d -3* 31536000. _d 0
254       &                      *1 _d -12       &                      *1 _d -12
255           _BEGIN_MASTER( myThid )
256         write(dar_ratetot_unit,*) myIter, PPtot, Nfixtot, Denittot         write(dar_ratetot_unit,*) myIter, PPtot, Nfixtot, Denittot
257           _END_MASTER( myThid )
258  #endif  #endif
259    
260                
# Line 217  c Line 272  c
272       &        myIter,myThid)       &        myIter,myThid)
273           CALL WRITE_FLD_XYZ_RL('DAR_Denittave.',suff,Denitave,           CALL WRITE_FLD_XYZ_RL('DAR_Denittave.',suff,Denitave,
274       &        myIter,myThid)       &        myIter,myThid)
275    #ifdef DAR_DIAG_PARW
276             do ii = 1,tlam
277               WRITE(suff,'(A8,I2.2,A1,I10.10)') 'PARwave',
278         &           ii,'.',myIter
279               CALL WRITE_FLD_XYZ_RL(suff,' ',
280         &           PARwave(1-Olx,1-Oly,1,1,1,ii),
281         &           myIter,myThid)
282             enddo
283             do np = 1,npmax
284               WRITE(suff,'(A8,I2.2,A1,I10.10)') 'chl2cave',
285         &           np,'.',myIter
286               CALL WRITE_FLD_XYZ_RL(suff,' ',
287         &           chl2cave(1-Olx,1-Oly,1,1,1,np),
288         &           myIter,myThid)
289             enddo
290    #endif
291  #ifdef DAR_DIAG_ACDOM  #ifdef DAR_DIAG_ACDOM
292           CALL WRITE_FLD_XYZ_RL('DAR_aCDOMtave',suff,aCDOMave,           CALL WRITE_FLD_XYZ_RL('DAR_aCDOMtave',suff,aCDOMave,
293       &        myIter,myThid)       &        myIter,myThid)
# Line 238  c Line 309  c
309             CALL WRITE_FLD_XYZ_RL(suff,' ',             CALL WRITE_FLD_XYZ_RL(suff,' ',
310       &           Euave(1-Olx,1-Oly,1,1,1,ii),       &           Euave(1-Olx,1-Oly,1,1,1,ii),
311       &           myIter,myThid)       &           myIter,myThid)
312               WRITE(suff,'(A8,I2.2,A1,I10.10)') 'Estopave',
313         &           ii,'.',myIter
314               CALL WRITE_FLD_XYZ_RL(suff,' ',
315         &           Estave(1-Olx,1-Oly,1,1,1,ii),
316         &           myIter,myThid)
317             WRITE(suff,'(A8,I2.2,A1,I10.10)') 'Eutopave',             WRITE(suff,'(A8,I2.2,A1,I10.10)') 'Eutopave',
318       &           ii,'.',myIter       &           ii,'.',myIter
319             CALL WRITE_FLD_XYZ_RL(suff,' ',             CALL WRITE_FLD_XYZ_RL(suff,' ',
# Line 245  c Line 321  c
321       &           myIter,myThid)       &           myIter,myThid)
322           enddo           enddo
323  #endif  #endif
324    #ifdef DAR_DIAG_IRR_AMPS
325             do ii = 1,tlam
326               WRITE(suff,'(A5,I2.2,A1,I10.10)') 'amp1ave',
327         &           ii,'.',myIter
328               CALL WRITE_FLD_XYZ_RL(suff,' ',
329         &           amp1ave(1-Olx,1-Oly,1,1,1,ii),
330         &           myIter,myThid)
331               WRITE(suff,'(A5,I2.2,A1,I10.10)') 'amp2ave',
332         &           ii,'.',myIter
333               CALL WRITE_FLD_XYZ_RL(suff,' ',
334         &           amp2ave(1-Olx,1-Oly,1,1,1,ii),
335         &           myIter,myThid)
336             enddo
337    #endif
338  #ifdef DAR_DIAG_ABSORP  #ifdef DAR_DIAG_ABSORP
339           do ii = 1,tlam           do ii = 1,tlam
340             WRITE(suff,'(A4,I2.2,A1,I10.10)') 'aave',             WRITE(suff,'(A4,I2.2,A1,I10.10)') 'aave',
# Line 287  c Line 377  c
377       &           myIter,myThid)       &           myIter,myThid)
378           enddo           enddo
379  #endif  #endif
380    #ifdef DAR_RADTRANS
381             WRITE(suff,'(A7,A1,I10.10)') 'rmudave',
382         &         '.',myIter
383             CALL WRITE_FLD_XY_RL(suff,' ',
384         &         rmudave(1-Olx,1-Oly,1,1),
385         &         myIter,myThid)
386    #endif
387    #ifdef DAR_DIAG_EK
388             do np = 1,npmax
389               WRITE(suff,'(A8,I2.2,A1,I10.10)') 'Ekave',
390         &           np,'.',myIter
391               CALL WRITE_FLD_XYZ_RL(suff,' ',
392         &           Ekave(1-Olx,1-Oly,1,1,1,np),
393         &           myIter,myThid)
394               WRITE(suff,'(A8,I2.2,A1,I10.10)') 'EkoverEave',
395         &           np,'.',myIter
396               CALL WRITE_FLD_XYZ_RL(suff,' ',
397         &           EkoverEave(1-Olx,1-Oly,1,1,1,np),
398         &           myIter,myThid)
399               WRITE(suff,'(A8,I2.2,A1,I10.10)') 'acclimave',
400         &           np,'.',myIter
401               CALL WRITE_FLD_XYZ_RL(suff,' ',
402         &           acclimave(1-Olx,1-Oly,1,1,1,np),
403         &           myIter,myThid)
404               do ii = 1,tlam
405                 WRITE(suff,'(A8,I2.2,A1,I2.2,A1,I10.10)') 'Ek_nl_ave',
406         &             np,'_',ii,'.',myIter
407                 CALL WRITE_FLD_XYZ_RL(suff,' ',
408         &             Ek_nlave(1-Olx,1-Oly,1,1,1,np,ii),
409         &             myIter,myThid)
410                 WRITE(suff,'(A8,I2.2,A1,I2.2,A1,I10.10)') 'EkoverE_nl_ave',
411         &             np,'_',ii,'.',myIter
412                 CALL WRITE_FLD_XYZ_RL(suff,' ',
413         &             EkoverE_nlave(1-Olx,1-Oly,1,1,1,np,ii),
414         &             myIter,myThid)
415               enddo
416             enddo
417    #endif
418    
419  c ANNA_TAVE  c ANNA_TAVE
420  #ifdef WAVES_DIAG_PCHL  #ifdef WAVES_DIAG_PCHL
# Line 367  c Line 495  c
495       &        pf,'dar_tave',0,0,'dar_PAR_ave',PARave,myThid)       &        pf,'dar_tave',0,0,'dar_PAR_ave',PARave,myThid)
496           CALL MNC_CW_RL_W(           CALL MNC_CW_RL_W(
497       &        pf,'dar_tave',0,0,'dar_PP_ave',PPave,myThid)       &        pf,'dar_tave',0,0,'dar_PP_ave',PPave,myThid)
498    #ifdef QUOTA_DIAG_LIMIT
499             do np=1,npmax
500                WRITE(name,'(a,i3.3)') 'dar_Nlim_',np
501                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
502         &           Nlimave(1-OLx,1-OLy,1,1,1,np),myThid)
503                WRITE(name,'(a,i3.3)') 'dar_Flim_',np
504                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
505         &           Flimave(1-OLx,1-OLy,1,1,1,np),myThid)
506                WRITE(name,'(a,i3.3)') 'dar_Ilim_',np
507                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
508         &           Ilimave(1-OLx,1-OLy,1,1,1,np),myThid)
509             enddo
510             CALL MNC_CW_RL_W(
511         &        pf,'dar_tave',0,0,'dar_Tlim_ave',Tlimave,myThid)
512    #endif
513  #ifndef ALLOW_QUOTA  #ifndef ALLOW_QUOTA
514           CALL MNC_CW_RL_W(           CALL MNC_CW_RL_W(
515       &        pf,'dar_tave',0,0,'dar_Chl_ave',Chlave,myThid)       &        pf,'dar_tave',0,0,'dar_Chl_ave',Chlave,myThid)
# Line 375  c Line 518  c
518       &        pf,'dar_tave',0,0,'dar_Nfix_ave',Nfixave,myThid)       &        pf,'dar_tave',0,0,'dar_Nfix_ave',Nfixave,myThid)
519           CALL MNC_CW_RL_W(           CALL MNC_CW_RL_W(
520       &        pf,'dar_tave',0,0,'dar_Denit_ave',Denitave,myThid)       &        pf,'dar_tave',0,0,'dar_Denit_ave',Denitave,myThid)
521    #ifdef DAR_DIAG_PARW
522             do ii = 1,tlam
523                DO jj = 1,MAX_LEN_FNAM
524                  name(jj:jj) = ' '
525                ENDDO
526                WRITE(name,'(a,i2.2)') 'dar_PARw',ii
527                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
528         &           PARwave(1-OLx,1-OLy,1,1,1,ii),myThid)
529            enddo
530             do np = 1,npmax
531                DO jj = 1,MAX_LEN_FNAM
532                  name(jj:jj) = ' '
533                ENDDO
534                WRITE(name,'(a,i2.2)') 'dar_chl2c',np
535                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
536         &           chl2cave(1-OLx,1-OLy,1,1,1,np),myThid)
537            enddo
538    #endif
539  #ifdef DAR_DIAG_ACDOM  #ifdef DAR_DIAG_ACDOM
540           CALL MNC_CW_RL_W(           CALL MNC_CW_RL_W(
541       &        pf,'dar_tave',0,0,'dar_aCDOM_ave',aCDOMave,myThid)       &        pf,'dar_tave',0,0,'dar_aCDOM_ave',aCDOMave,myThid)
# Line 393  c Line 554  c
554              WRITE(name,'(a,i2.2)') 'dar_Eu',ii              WRITE(name,'(a,i2.2)') 'dar_Eu',ii
555              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
556       &           Euave(1-OLx,1-OLy,1,1,1,ii),myThid)       &           Euave(1-OLx,1-OLy,1,1,1,ii),myThid)
557                WRITE(name,'(a,i2.2)') 'dar_Estop',ii
558                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
559         &           Estave(1-OLx,1-OLy,1,1,1,ii),myThid)
560              WRITE(name,'(a,i2.2)') 'dar_Eutop',ii              WRITE(name,'(a,i2.2)') 'dar_Eutop',ii
561              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
562       &           Eutave(1-OLx,1-OLy,1,1,1,ii),myThid)       &           Eutave(1-OLx,1-OLy,1,1,1,ii),myThid)
563           enddo           enddo
564  #endif  #endif
565    #ifdef DAR_DIAG_IRR_AMPS
566             do ii = 1,tlam
567                DO jj = 1,MAX_LEN_FNAM
568                  name(jj:jj) = ' '
569                ENDDO
570                WRITE(name,'(a,i2.2)') 'dar_amp1_',ii
571                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
572         &           amp1ave(1-OLx,1-OLy,1,1,1,ii),myThid)
573                WRITE(name,'(a,i2.2)') 'dar_amp2_',ii
574                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
575         &           amp2ave(1-OLx,1-OLy,1,1,1,ii),myThid)
576             enddo
577    #endif
578  #ifdef DAR_DIAG_ABSORP  #ifdef DAR_DIAG_ABSORP
579           do ii = 1,tlam           do ii = 1,tlam
580              DO jj = 1,MAX_LEN_FNAM              DO jj = 1,MAX_LEN_FNAM
# Line 437  c Line 614  c
614       &           bbpartave(1-OLx,1-OLy,1,1,1,ii),myThid)       &           bbpartave(1-OLx,1-OLy,1,1,1,ii),myThid)
615           enddo           enddo
616  #endif  #endif
617    #ifdef DAR_RADTRANS
618             DO jj = 1,MAX_LEN_FNAM
619               name(jj:jj) = ' '
620             ENDDO
621             WRITE(name,'(a)') 'dar_rmud'
622             CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
623         &        rmudave(1-OLx,1-OLy,1,1),myThid)
624    #endif
625    #ifdef DAR_DIAG_EK
626             do np = 1,npmax
627                DO jj = 1,MAX_LEN_FNAM
628                  name(jj:jj) = ' '
629                ENDDO
630                WRITE(name,'(a,i2.2)') 'dar_Ek',np
631                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
632         &           Ekave(1-OLx,1-OLy,1,1,1,np),myThid)
633                WRITE(name,'(a,i2.2)') 'dar_EkoverE',np
634                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
635         &           EkoverEave(1-OLx,1-OLy,1,1,1,np),myThid)
636                WRITE(name,'(a,i2.2)') 'dar_acclim',np
637                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
638         &           acclimave(1-OLx,1-OLy,1,1,1,np),myThid)
639                do ii = 1,tlam
640                  WRITE(name,'(a,i2.2,a,i2.2)') 'dar_Ek_nl',np,'_',ii
641                  CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
642         &             Ek_nlave(1-OLx,1-OLy,1,1,1,np,ii),myThid)
643                  WRITE(name,'(a,i2.2,a,i2.2)') 'dar_EkoverE_nl',np,'_',ii
644                  CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
645         &             EkoverE_nlave(1-OLx,1-OLy,1,1,1,np,ii),myThid)
646                enddo
647            enddo
648    #endif
649  c ANNA_TAVE  c ANNA_TAVE
650  #ifdef WAVES_DIAG_PCHL  #ifdef WAVES_DIAG_PCHL
651           do np=1,npmax           do np=1,npmax
# Line 529  C      Reset averages to zero Line 738  C      Reset averages to zero
738          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
739           CALL TIMEAVE_RESET(PARave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(PARave,Nr,bi,bj,myThid)
740           CALL TIMEAVE_RESET(PPave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(PPave,Nr,bi,bj,myThid)
741    #ifdef QUOTA_DIAG_LIMIT
742             CALL TIMEAVE_RESET(Nlimave,Nr,bi,bj,myThid)
743             CALL TIMEAVE_RESET(Flimave,Nr,bi,bj,myThid)
744             CALL TIMEAVE_RESET(Ilimave,Nr,bi,bj,myThid)
745             CALL TIMEAVE_RESET(Tlimave,Nr,bi,bj,myThid)
746    #endif
747  #ifndef ALLOW_QUOTA  #ifndef ALLOW_QUOTA
748           CALL TIMEAVE_RESET(Chlave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(Chlave,Nr,bi,bj,myThid)
749  #endif  #endif
750           CALL TIMEAVE_RESET(Nfixave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(Nfixave,Nr,bi,bj,myThid)
751           CALL TIMEAVE_RESET(Denitave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(Denitave,Nr,bi,bj,myThid)
752    #ifdef DAR_DIAG_PARW
753             do ii = 1,tlam
754               CALL TIMEAVE_RESET(PARwave(1-OLx,1-OLy,1,1,1,ii),
755         &                                    Nr,bi,bj,myThid)
756             enddo
757             do np = 1,npmax
758               CALL TIMEAVE_RESET(chl2cave(1-OLx,1-OLy,1,1,1,np),
759         &                                    Nr,bi,bj,myThid)
760             enddo
761    #endif
762  #ifdef DAR_DIAG_ACDOM  #ifdef DAR_DIAG_ACDOM
763           CALL TIMEAVE_RESET(aCDOMave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(aCDOMave,Nr,bi,bj,myThid)
764  #endif  #endif
# Line 545  C      Reset averages to zero Line 770  C      Reset averages to zero
770       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
771             CALL TIMEAVE_RESET(Euave(1-OLx,1-OLy,1,1,1,ii),             CALL TIMEAVE_RESET(Euave(1-OLx,1-OLy,1,1,1,ii),
772       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
773               CALL TIMEAVE_RESET(Estave(1-OLx,1-OLy,1,1,1,ii),
774         &                                    Nr,bi,bj,myThid)
775             CALL TIMEAVE_RESET(Eutave(1-OLx,1-OLy,1,1,1,ii),             CALL TIMEAVE_RESET(Eutave(1-OLx,1-OLy,1,1,1,ii),
776       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
777           enddo           enddo
778  #endif  #endif
779    #ifdef DAR_DIAG_IRR_AMPS
780             do ii = 1,tlam
781               CALL TIMEAVE_RESET(amp1ave(1-OLx,1-OLy,1,1,1,ii),
782         &                                    Nr,bi,bj,myThid)
783               CALL TIMEAVE_RESET(amp2ave(1-OLx,1-OLy,1,1,1,ii),
784         &                                    Nr,bi,bj,myThid)
785             enddo
786    #endif
787  #ifdef DAR_DIAG_ABSORP  #ifdef DAR_DIAG_ABSORP
788           do ii = 1,tlam           do ii = 1,tlam
789             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 808  C      Reset averages to zero
808       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
809           enddo           enddo
810  #endif  #endif
811    #ifdef DAR_RADTRANS
812             CALL TIMEAVE_RESET(rmudave(1-OLx,1-OLy,1,1),
813         &                                  1,bi,bj,myThid)
814    #endif
815    #ifdef DAR_DIAG_EK
816             do np = 1,npmax
817               CALL TIMEAVE_RESET(Ekave(1-OLx,1-OLy,1,1,1,np),
818         &                                    Nr,bi,bj,myThid)
819               CALL TIMEAVE_RESET(EkoverEave(1-OLx,1-OLy,1,1,1,np),
820         &                                    Nr,bi,bj,myThid)
821               CALL TIMEAVE_RESET(acclimave(1-OLx,1-OLy,1,1,1,np),
822         &                                    Nr,bi,bj,myThid)
823               do ii = 1,tlam
824                 CALL TIMEAVE_RESET(Ek_nlave(1-OLx,1-OLy,1,1,1,np,ii),
825         &                                    Nr,bi,bj,myThid)
826                 CALL TIMEAVE_RESET(EkoverE_nlave(1-OLx,1-OLy,1,1,1,np,ii),
827         &                                    Nr,bi,bj,myThid)
828               enddo
829             enddo
830    #endif
831  c ANNA_TAVE  c ANNA_TAVE
832  #ifdef WAVES_DIAG_PCHL  #ifdef WAVES_DIAG_PCHL
833           do np=1,npmax           do np=1,npmax
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22