/[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.3 by benw, Mon Jul 2 09:37:48 2012 UTC revision 1.6 by jahn, Thu Sep 20 14:13:29 2012 UTC
# Line 108  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 136  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 202  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 250  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 257  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 299  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 420  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 464  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 578  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 606  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.3  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22