/[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.5 by jahn, Thu Aug 23 21:48:24 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  #ifdef DAR_DIAG_IRR_AMPS
118           do ii = 1,tlam           do ii = 1,tlam
119            CALL TIMEAVE_NORMALIZ(c1ave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZ(amp1ave(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(c2ave(1-Olx,1-Oly,1,1,1,ii),            CALL TIMEAVE_NORMALIZ(amp2ave(1-Olx,1-Oly,1,1,1,ii),
122       &                          DAR_timeave,nR,bi,bj,myThid)       &                          DAR_timeave,nR,bi,bj,myThid)
123           enddo           enddo
124  #endif  #endif
# Line 262  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 271  c Line 278  c
278  #endif  #endif
279  #ifdef DAR_DIAG_IRR_AMPS  #ifdef DAR_DIAG_IRR_AMPS
280           do ii = 1,tlam           do ii = 1,tlam
281             WRITE(suff,'(A5,I2.2,A1,I10.10)') 'c1ave',             WRITE(suff,'(A5,I2.2,A1,I10.10)') 'amp1ave',
282       &           ii,'.',myIter       &           ii,'.',myIter
283             CALL WRITE_FLD_XYZ_RL(suff,' ',             CALL WRITE_FLD_XYZ_RL(suff,' ',
284       &           c1ave(1-Olx,1-Oly,1,1,1,ii),       &           amp1ave(1-Olx,1-Oly,1,1,1,ii),
285       &           myIter,myThid)       &           myIter,myThid)
286             WRITE(suff,'(A5,I2.2,A1,I10.10)') 'c2ave',             WRITE(suff,'(A5,I2.2,A1,I10.10)') 'amp2ave',
287       &           ii,'.',myIter       &           ii,'.',myIter
288             CALL WRITE_FLD_XYZ_RL(suff,' ',             CALL WRITE_FLD_XYZ_RL(suff,' ',
289       &           c2ave(1-Olx,1-Oly,1,1,1,ii),       &           amp2ave(1-Olx,1-Oly,1,1,1,ii),
290       &           myIter,myThid)       &           myIter,myThid)
291           enddo           enddo
292  #endif  #endif
# Line 453  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)
# Line 463  c Line 473  c
473              DO jj = 1,MAX_LEN_FNAM              DO jj = 1,MAX_LEN_FNAM
474                name(jj:jj) = ' '                name(jj:jj) = ' '
475              ENDDO              ENDDO
476              WRITE(name,'(a,i2.2)') 'dar_c1_',ii              WRITE(name,'(a,i2.2)') 'dar_amp1_',ii
477              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
478       &           c1ave(1-OLx,1-OLy,1,1,1,ii),myThid)       &           amp1ave(1-OLx,1-OLy,1,1,1,ii),myThid)
479              WRITE(name,'(a,i2.2)') 'dar_c2_',ii              WRITE(name,'(a,i2.2)') 'dar_amp2_',ii
480              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
481       &           c2ave(1-OLx,1-OLy,1,1,1,ii),myThid)       &           amp2ave(1-OLx,1-OLy,1,1,1,ii),myThid)
482           enddo           enddo
483  #endif  #endif
484  #ifdef DAR_DIAG_ABSORP  #ifdef DAR_DIAG_ABSORP
# Line 632  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  #ifdef DAR_DIAG_IRR_AMPS
652           do ii = 1,tlam           do ii = 1,tlam
653             CALL TIMEAVE_RESET(c1ave(1-OLx,1-OLy,1,1,1,ii),             CALL TIMEAVE_RESET(amp1ave(1-OLx,1-OLy,1,1,1,ii),
654       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
655             CALL TIMEAVE_RESET(c2ave(1-OLx,1-OLy,1,1,1,ii),             CALL TIMEAVE_RESET(amp2ave(1-OLx,1-OLy,1,1,1,ii),
656       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
657           enddo           enddo
658  #endif  #endif

Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5

  ViewVC Help
Powered by ViewVC 1.1.22