/[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.6 by jahn, Thu Sep 20 14:13:29 2012 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 96  C      Normalize by integrated time Line 108  C      Normalize by integrated time
108       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
109            CALL TIMEAVE_NORMALIZ(Euave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZ(Euave(1-Olx,1-Oly,1,1,1,ii),
110       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
111              CALL TIMEAVE_NORMALIZ(Estave(1-Olx,1-Oly,1,1,1,ii),
112         &                          DAR_timeave,nR,bi,bj,myThid)
113            CALL TIMEAVE_NORMALIZ(Eutave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZ(Eutave(1-Olx,1-Oly,1,1,1,ii),
114       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
115           enddo           enddo
116  #endif  #endif
117    #ifdef DAR_DIAG_IRR_AMPS
118             do ii = 1,tlam
119              CALL TIMEAVE_NORMALIZ(amp1ave(1-Olx,1-Oly,1,1,1,ii),
120         &                          DAR_timeave,nR,bi,bj,myThid)
121              CALL TIMEAVE_NORMALIZ(amp2ave(1-Olx,1-Oly,1,1,1,ii),
122         &                          DAR_timeave,nR,bi,bj,myThid)
123             enddo
124    #endif
125  #ifdef DAR_DIAG_ABSORP  #ifdef DAR_DIAG_ABSORP
126           do ii = 1,tlam           do ii = 1,tlam
127            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 146  C      Normalize by integrated time
146       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
147           enddo           enddo
148  #endif  #endif
149    #ifdef DAR_RADTRANS
150             CALL TIMEAVE_NORMALIZ(rmudave(1-Olx,1-Oly,1,1),
151         &                         DAR_timeave,1,bi,bj,myThid)
152    #endif
153    
154  c ANNA_TAVE  c ANNA_TAVE
155  #ifdef WAVES_DIAG_PCHL  #ifdef WAVES_DIAG_PCHL
# Line 190  c Line 216  c
216         call tracer_sumvol(         call tracer_sumvol(
217       &             PPave, nR, PPtot, myThid )       &             PPave, nR, PPtot, myThid )
218         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
219       &                      *1. _d-15       &                      *1. _d -15
220         call tracer_sumvol(         call tracer_sumvol(
221       &             Nfixave, nR, Nfixtot, myThid )       &             Nfixave, nR, Nfixtot, myThid )
222         Nfixtot=Nfixtot*14. _d 0 *1. _d -3* 31536000. _d 0         Nfixtot=Nfixtot*14. _d 0 *1. _d -3* 31536000. _d 0
# Line 238  c Line 264  c
264             CALL WRITE_FLD_XYZ_RL(suff,' ',             CALL WRITE_FLD_XYZ_RL(suff,' ',
265       &           Euave(1-Olx,1-Oly,1,1,1,ii),       &           Euave(1-Olx,1-Oly,1,1,1,ii),
266       &           myIter,myThid)       &           myIter,myThid)
267               WRITE(suff,'(A8,I2.2,A1,I10.10)') 'Estopave',
268         &           ii,'.',myIter
269               CALL WRITE_FLD_XYZ_RL(suff,' ',
270         &           Estave(1-Olx,1-Oly,1,1,1,ii),
271         &           myIter,myThid)
272             WRITE(suff,'(A8,I2.2,A1,I10.10)') 'Eutopave',             WRITE(suff,'(A8,I2.2,A1,I10.10)') 'Eutopave',
273       &           ii,'.',myIter       &           ii,'.',myIter
274             CALL WRITE_FLD_XYZ_RL(suff,' ',             CALL WRITE_FLD_XYZ_RL(suff,' ',
# Line 245  c Line 276  c
276       &           myIter,myThid)       &           myIter,myThid)
277           enddo           enddo
278  #endif  #endif
279    #ifdef DAR_DIAG_IRR_AMPS
280             do ii = 1,tlam
281               WRITE(suff,'(A5,I2.2,A1,I10.10)') 'amp1ave',
282         &           ii,'.',myIter
283               CALL WRITE_FLD_XYZ_RL(suff,' ',
284         &           amp1ave(1-Olx,1-Oly,1,1,1,ii),
285         &           myIter,myThid)
286               WRITE(suff,'(A5,I2.2,A1,I10.10)') 'amp2ave',
287         &           ii,'.',myIter
288               CALL WRITE_FLD_XYZ_RL(suff,' ',
289         &           amp2ave(1-Olx,1-Oly,1,1,1,ii),
290         &           myIter,myThid)
291             enddo
292    #endif
293  #ifdef DAR_DIAG_ABSORP  #ifdef DAR_DIAG_ABSORP
294           do ii = 1,tlam           do ii = 1,tlam
295             WRITE(suff,'(A4,I2.2,A1,I10.10)') 'aave',             WRITE(suff,'(A4,I2.2,A1,I10.10)') 'aave',
# Line 287  c Line 332  c
332       &           myIter,myThid)       &           myIter,myThid)
333           enddo           enddo
334  #endif  #endif
335    #ifdef DAR_RADTRANS
336             WRITE(suff,'(A7,A1,I10.10)') 'rmudave',
337         &         '.',myIter
338             CALL WRITE_FLD_XY_RL(suff,' ',
339         &         rmudave(1-Olx,1-Oly,1,1),
340         &         myIter,myThid)
341    #endif
342    
343  c ANNA_TAVE  c ANNA_TAVE
344  #ifdef WAVES_DIAG_PCHL  #ifdef WAVES_DIAG_PCHL
# Line 367  c Line 419  c
419       &        pf,'dar_tave',0,0,'dar_PAR_ave',PARave,myThid)       &        pf,'dar_tave',0,0,'dar_PAR_ave',PARave,myThid)
420           CALL MNC_CW_RL_W(           CALL MNC_CW_RL_W(
421       &        pf,'dar_tave',0,0,'dar_PP_ave',PPave,myThid)       &        pf,'dar_tave',0,0,'dar_PP_ave',PPave,myThid)
422    #ifdef QUOTA_DIAG_LIMIT
423             do np=1,npmax
424                WRITE(name,'(a,i3.3)') 'dar_Nlim_',np
425                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
426         &           Nlimave(1-OLx,1-OLy,1,1,1,np),myThid)
427                WRITE(name,'(a,i3.3)') 'dar_Flim_',np
428                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
429         &           Flimave(1-OLx,1-OLy,1,1,1,np),myThid)
430                WRITE(name,'(a,i3.3)') 'dar_Ilim_',np
431                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
432         &           Ilimave(1-OLx,1-OLy,1,1,1,np),myThid)
433             enddo
434             CALL MNC_CW_RL_W(
435         &        pf,'dar_tave',0,0,'dar_Tlim_ave',Tlimave,myThid)
436    #endif
437  #ifndef ALLOW_QUOTA  #ifndef ALLOW_QUOTA
438           CALL MNC_CW_RL_W(           CALL MNC_CW_RL_W(
439       &        pf,'dar_tave',0,0,'dar_Chl_ave',Chlave,myThid)       &        pf,'dar_tave',0,0,'dar_Chl_ave',Chlave,myThid)
# Line 393  c Line 460  c
460              WRITE(name,'(a,i2.2)') 'dar_Eu',ii              WRITE(name,'(a,i2.2)') 'dar_Eu',ii
461              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
462       &           Euave(1-OLx,1-OLy,1,1,1,ii),myThid)       &           Euave(1-OLx,1-OLy,1,1,1,ii),myThid)
463                WRITE(name,'(a,i2.2)') 'dar_Estop',ii
464                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
465         &           Estave(1-OLx,1-OLy,1,1,1,ii),myThid)
466              WRITE(name,'(a,i2.2)') 'dar_Eutop',ii              WRITE(name,'(a,i2.2)') 'dar_Eutop',ii
467              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
468       &           Eutave(1-OLx,1-OLy,1,1,1,ii),myThid)       &           Eutave(1-OLx,1-OLy,1,1,1,ii),myThid)
469           enddo           enddo
470  #endif  #endif
471    #ifdef DAR_DIAG_IRR_AMPS
472             do ii = 1,tlam
473                DO jj = 1,MAX_LEN_FNAM
474                  name(jj:jj) = ' '
475                ENDDO
476                WRITE(name,'(a,i2.2)') 'dar_amp1_',ii
477                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
478         &           amp1ave(1-OLx,1-OLy,1,1,1,ii),myThid)
479                WRITE(name,'(a,i2.2)') 'dar_amp2_',ii
480                CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
481         &           amp2ave(1-OLx,1-OLy,1,1,1,ii),myThid)
482             enddo
483    #endif
484  #ifdef DAR_DIAG_ABSORP  #ifdef DAR_DIAG_ABSORP
485           do ii = 1,tlam           do ii = 1,tlam
486              DO jj = 1,MAX_LEN_FNAM              DO jj = 1,MAX_LEN_FNAM
# Line 437  c Line 520  c
520       &           bbpartave(1-OLx,1-OLy,1,1,1,ii),myThid)       &           bbpartave(1-OLx,1-OLy,1,1,1,ii),myThid)
521           enddo           enddo
522  #endif  #endif
523    #ifdef DAR_RADTRANS
524             DO jj = 1,MAX_LEN_FNAM
525               name(jj:jj) = ' '
526             ENDDO
527             WRITE(name,'(a)') 'dar_rmud'
528             CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
529         &        rmudave(1-OLx,1-OLy,1,1),myThid)
530    #endif
531  c ANNA_TAVE  c ANNA_TAVE
532  #ifdef WAVES_DIAG_PCHL  #ifdef WAVES_DIAG_PCHL
533           do np=1,npmax           do np=1,npmax
# Line 529  C      Reset averages to zero Line 620  C      Reset averages to zero
620          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
621           CALL TIMEAVE_RESET(PARave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(PARave,Nr,bi,bj,myThid)
622           CALL TIMEAVE_RESET(PPave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(PPave,Nr,bi,bj,myThid)
623    #ifdef QUOTA_DIAG_LIMIT
624             CALL TIMEAVE_RESET(Nlimave,Nr,bi,bj,myThid)
625             CALL TIMEAVE_RESET(Flimave,Nr,bi,bj,myThid)
626             CALL TIMEAVE_RESET(Ilimave,Nr,bi,bj,myThid)
627             CALL TIMEAVE_RESET(Tlimave,Nr,bi,bj,myThid)
628    #endif
629  #ifndef ALLOW_QUOTA  #ifndef ALLOW_QUOTA
630           CALL TIMEAVE_RESET(Chlave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(Chlave,Nr,bi,bj,myThid)
631  #endif  #endif
# Line 545  C      Reset averages to zero Line 642  C      Reset averages to zero
642       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
643             CALL TIMEAVE_RESET(Euave(1-OLx,1-OLy,1,1,1,ii),             CALL TIMEAVE_RESET(Euave(1-OLx,1-OLy,1,1,1,ii),
644       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
645               CALL TIMEAVE_RESET(Estave(1-OLx,1-OLy,1,1,1,ii),
646         &                                    Nr,bi,bj,myThid)
647             CALL TIMEAVE_RESET(Eutave(1-OLx,1-OLy,1,1,1,ii),             CALL TIMEAVE_RESET(Eutave(1-OLx,1-OLy,1,1,1,ii),
648       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
649           enddo           enddo
650  #endif  #endif
651    #ifdef DAR_DIAG_IRR_AMPS
652             do ii = 1,tlam
653               CALL TIMEAVE_RESET(amp1ave(1-OLx,1-OLy,1,1,1,ii),
654         &                                    Nr,bi,bj,myThid)
655               CALL TIMEAVE_RESET(amp2ave(1-OLx,1-OLy,1,1,1,ii),
656         &                                    Nr,bi,bj,myThid)
657             enddo
658    #endif
659  #ifdef DAR_DIAG_ABSORP  #ifdef DAR_DIAG_ABSORP
660           do ii = 1,tlam           do ii = 1,tlam
661             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 680  C      Reset averages to zero
680       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
681           enddo           enddo
682  #endif  #endif
683    #ifdef DAR_RADTRANS
684             CALL TIMEAVE_RESET(rmudave(1-OLx,1-OLy,1,1),
685         &                                  1,bi,bj,myThid)
686    #endif
687  c ANNA_TAVE  c ANNA_TAVE
688  #ifdef WAVES_DIAG_PCHL  #ifdef WAVES_DIAG_PCHL
689           do np=1,npmax           do np=1,npmax
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22