/[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.4 by jahn, Mon Jul 30 15:21:51 2012 UTC revision 1.7 by stephd, Tue Oct 23 17:55:30 2012 UTC
# Line 96  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 108  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  #ifdef DAR_DIAG_IRR_AMPS
128           do ii = 1,tlam           do ii = 1,tlam
129            CALL TIMEAVE_NORMALIZ(c1ave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZ(amp1ave(1-Olx,1-Oly,1,1,1,ii),
130       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
131            CALL TIMEAVE_NORMALIZ(c2ave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZ(amp2ave(1-Olx,1-Oly,1,1,1,ii),
132       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
133           enddo           enddo
134  #endif  #endif
# Line 148  C      Normalize by integrated time Line 160  C      Normalize by integrated time
160           CALL TIMEAVE_NORMALIZ(rmudave(1-Olx,1-Oly,1,1),           CALL TIMEAVE_NORMALIZ(rmudave(1-Olx,1-Oly,1,1),
161       &                         DAR_timeave,1,bi,bj,myThid)       &                         DAR_timeave,1,bi,bj,myThid)
162  #endif  #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 214  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 241  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 262  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 271  c Line 321  c
321  #endif  #endif
322  #ifdef DAR_DIAG_IRR_AMPS  #ifdef DAR_DIAG_IRR_AMPS
323           do ii = 1,tlam           do ii = 1,tlam
324             WRITE(suff,'(A5,I2.2,A1,I10.10)') 'c1ave',             WRITE(suff,'(A5,I2.2,A1,I10.10)') 'amp1ave',
325       &           ii,'.',myIter       &           ii,'.',myIter
326             CALL WRITE_FLD_XYZ_RL(suff,' ',             CALL WRITE_FLD_XYZ_RL(suff,' ',
327       &           c1ave(1-Olx,1-Oly,1,1,1,ii),       &           amp1ave(1-Olx,1-Oly,1,1,1,ii),
328       &           myIter,myThid)       &           myIter,myThid)
329             WRITE(suff,'(A5,I2.2,A1,I10.10)') 'c2ave',             WRITE(suff,'(A5,I2.2,A1,I10.10)') 'amp2ave',
330       &           ii,'.',myIter       &           ii,'.',myIter
331             CALL WRITE_FLD_XYZ_RL(suff,' ',             CALL WRITE_FLD_XYZ_RL(suff,' ',
332       &           c2ave(1-Olx,1-Oly,1,1,1,ii),       &           amp2ave(1-Olx,1-Oly,1,1,1,ii),
333       &           myIter,myThid)       &           myIter,myThid)
334           enddo           enddo
335  #endif  #endif
# Line 332  c Line 382  c
382       &         rmudave(1-Olx,1-Oly,1,1),       &         rmudave(1-Olx,1-Oly,1,1),
383       &         myIter,myThid)       &         myIter,myThid)
384  #endif  #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 435  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 453  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)
# Line 463  c Line 565  c
565              DO jj = 1,MAX_LEN_FNAM              DO jj = 1,MAX_LEN_FNAM
566                name(jj:jj) = ' '                name(jj:jj) = ' '
567              ENDDO              ENDDO
568              WRITE(name,'(a,i2.2)') 'dar_c1_',ii              WRITE(name,'(a,i2.2)') 'dar_amp1_',ii
569              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
570       &           c1ave(1-OLx,1-OLy,1,1,1,ii),myThid)       &           amp1ave(1-OLx,1-OLy,1,1,1,ii),myThid)
571              WRITE(name,'(a,i2.2)') 'dar_c2_',ii              WRITE(name,'(a,i2.2)') 'dar_amp2_',ii
572              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
573       &           c2ave(1-OLx,1-OLy,1,1,1,ii),myThid)       &           amp2ave(1-OLx,1-OLy,1,1,1,ii),myThid)
574           enddo           enddo
575  #endif  #endif
576  #ifdef DAR_DIAG_ABSORP  #ifdef DAR_DIAG_ABSORP
# Line 518  c Line 620  c
620           CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,           CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
621       &        rmudave(1-OLx,1-OLy,1,1),myThid)       &        rmudave(1-OLx,1-OLy,1,1),myThid)
622  #endif  #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 621  C      Reset averages to zero Line 747  C      Reset averages to zero
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 632  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  #ifdef DAR_DIAG_IRR_AMPS
778           do ii = 1,tlam           do ii = 1,tlam
779             CALL TIMEAVE_RESET(c1ave(1-OLx,1-OLy,1,1,1,ii),             CALL TIMEAVE_RESET(amp1ave(1-OLx,1-OLy,1,1,1,ii),
780       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
781             CALL TIMEAVE_RESET(c2ave(1-OLx,1-OLy,1,1,1,ii),             CALL TIMEAVE_RESET(amp2ave(1-OLx,1-OLy,1,1,1,ii),
782       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
783           enddo           enddo
784  #endif  #endif
# Line 672  C      Reset averages to zero Line 810  C      Reset averages to zero
810           CALL TIMEAVE_RESET(rmudave(1-OLx,1-OLy,1,1),           CALL TIMEAVE_RESET(rmudave(1-OLx,1-OLy,1,1),
811       &                                  1,bi,bj,myThid)       &                                  1,bi,bj,myThid)
812  #endif  #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  c ANNA_TAVE
830  #ifdef WAVES_DIAG_PCHL  #ifdef WAVES_DIAG_PCHL
831           do np=1,npmax           do np=1,npmax
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22