/[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.9 by jahn, Fri Dec 27 17:29:00 2013 UTC revision 1.10 by benw, Tue May 19 14:30:15 2015 UTC
# Line 75  C      Normalize by integrated time Line 75  C      Normalize by integrated time
75       &                 bi,bj,myThid)       &                 bi,bj,myThid)
76  #ifdef QUOTA_DIAG_LIMIT  #ifdef QUOTA_DIAG_LIMIT
77           do np=1,npmax           do np=1,npmax
78    ! Carbon
79               CALL TIMEAVE_NORMALIZE(AP_C_ave(1-Olx,1-Oly,1,1,1,np),
80         &                           DAR_timeave, nR ,bi,bj,myThid)
81               CALL TIMEAVE_NORMALIZE(HP_C_ave(1-Olx,1-Oly,1,1,1,np),
82         &                           DAR_timeave, nR ,bi,bj,myThid)
83    ! Nitrogen
84               CALL TIMEAVE_NORMALIZE(AP_N_ave(1-Olx,1-Oly,1,1,1,np),
85         &                           DAR_timeave, nR ,bi,bj,myThid)
86               CALL TIMEAVE_NORMALIZE(HP_N_ave(1-Olx,1-Oly,1,1,1,np),
87         &                           DAR_timeave, nR ,bi,bj,myThid)
88             CALL TIMEAVE_NORMALIZE(Nlimave(1-Olx,1-Oly,1,1,1,np),             CALL TIMEAVE_NORMALIZE(Nlimave(1-Olx,1-Oly,1,1,1,np),
89       &                          DAR_timeave, nR ,bi,bj,myThid)       &                           DAR_timeave, nR ,bi,bj,myThid)
90    ! Phosphorus
91    #ifdef PQUOTA
92               CALL TIMEAVE_NORMALIZE(AP_P_ave(1-Olx,1-Oly,1,1,1,np),
93         &                           DAR_timeave, nR ,bi,bj,myThid)
94               CALL TIMEAVE_NORMALIZE(HP_P_ave(1-Olx,1-Oly,1,1,1,np),
95         &                           DAR_timeave, nR ,bi,bj,myThid)
96               CALL TIMEAVE_NORMALIZE(Plimave(1-Olx,1-Oly,1,1,1,np),
97         &                           DAR_timeave, nR ,bi,bj,myThid)
98    #endif
99    ! Iron
100    #ifdef FQUOTA
101               CALL TIMEAVE_NORMALIZE(AP_F_ave(1-Olx,1-Oly,1,1,1,np),
102         &                           DAR_timeave, nR ,bi,bj,myThid)
103               CALL TIMEAVE_NORMALIZE(HP_F_ave(1-Olx,1-Oly,1,1,1,np),
104         &                           DAR_timeave, nR ,bi,bj,myThid)
105             CALL TIMEAVE_NORMALIZE(Flimave(1-Olx,1-Oly,1,1,1,np),             CALL TIMEAVE_NORMALIZE(Flimave(1-Olx,1-Oly,1,1,1,np),
106       &                          DAR_timeave, nR ,bi,bj,myThid)       &                           DAR_timeave, nR ,bi,bj,myThid)
107    #endif
108    ! Light
109             CALL TIMEAVE_NORMALIZE(Ilimave(1-Olx,1-Oly,1,1,1,np),             CALL TIMEAVE_NORMALIZE(Ilimave(1-Olx,1-Oly,1,1,1,np),
110       &                          DAR_timeave, nR ,bi,bj,myThid)       &                           DAR_timeave, nR ,bi,bj,myThid)
111           enddo           enddo
112    ! Temperature
113           CALL TIMEAVE_NORMALIZE(Tlimave,DAR_timeave, nR ,           CALL TIMEAVE_NORMALIZE(Tlimave,DAR_timeave, nR ,
114       &                                           bi,bj,myThid)       &                                           bi,bj,myThid)
115  #endif  #endif
# Line 495  c Line 523  c
523       &        pf,'dar_tave',0,0,'dar_PP_ave',PPave,myThid)       &        pf,'dar_tave',0,0,'dar_PP_ave',PPave,myThid)
524  #ifdef QUOTA_DIAG_LIMIT  #ifdef QUOTA_DIAG_LIMIT
525           do np=1,npmax           do np=1,npmax
526              WRITE(name,'(a,i3.3)') 'dar_Nlim_',np  ! Carbon
527              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,             WRITE(name,'(a,i3.3)') 'dar_AP_C_',np
528       &           Nlimave(1-OLx,1-OLy,1,1,1,np),myThid)             CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
529              WRITE(name,'(a,i3.3)') 'dar_Flim_',np       &          AP_C_ave(1-OLx,1-OLy,1,1,1,np),myThid)
530              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,             WRITE(name,'(a,i3.3)') 'dar_HP_C_',np
531       &           Flimave(1-OLx,1-OLy,1,1,1,np),myThid)             CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
532              WRITE(name,'(a,i3.3)') 'dar_Ilim_',np       &          HP_C_ave(1-OLx,1-OLy,1,1,1,np),myThid)
533              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,  ! Nitrogen
534       &           Ilimave(1-OLx,1-OLy,1,1,1,np),myThid)             WRITE(name,'(a,i3.3)') 'dar_Nlim_',np
535               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
536         &          Nlimave(1-OLx,1-OLy,1,1,1,np),myThid)
537               WRITE(name,'(a,i3.3)') 'dar_AP_N_',np
538               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
539         &          AP_N_ave(1-OLx,1-OLy,1,1,1,np),myThid)
540               WRITE(name,'(a,i3.3)') 'dar_HP_N_',np
541               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
542         &          HP_N_ave(1-OLx,1-OLy,1,1,1,np),myThid)
543    ! Phosphorus
544    #ifdef PQUOTA
545               WRITE(name,'(a,i3.3)') 'dar_Plim_',np
546               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
547         &          Plimave(1-OLx,1-OLy,1,1,1,np),myThid)
548               WRITE(name,'(a,i3.3)') 'dar_AP_P_',np
549               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
550         &          AP_P_ave(1-OLx,1-OLy,1,1,1,np),myThid)
551               WRITE(name,'(a,i3.3)') 'dar_HP_P_',np
552               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
553         &          HP_P_ave(1-OLx,1-OLy,1,1,1,np),myThid)
554    #endif
555    ! Iron
556    #ifdef FQUOTA
557               WRITE(name,'(a,i3.3)') 'dar_Flim_',np
558               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
559         &          Flimave(1-OLx,1-OLy,1,1,1,np),myThid)
560               WRITE(name,'(a,i3.3)') 'dar_AP_F_',np
561               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
562         &          AP_F_ave(1-OLx,1-OLy,1,1,1,np),myThid)
563               WRITE(name,'(a,i3.3)') 'dar_HP_F_',np
564               CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
565         &          HP_F_ave(1-OLx,1-OLy,1,1,1,np),myThid)
566    #endif
567    ! Light
568              WRITE(name,'(a,i3.3)') 'dar_Ilim_',np
569              CALL MNC_CW_RL_W(pf,'dar_tave',0,0,name,
570         &         Ilimave(1-OLx,1-OLy,1,1,1,np),myThid)
571           enddo           enddo
572    ! Temperature
573           CALL MNC_CW_RL_W(           CALL MNC_CW_RL_W(
574       &        pf,'dar_tave',0,0,'dar_Tlim_ave',Tlimave,myThid)       &        pf,'dar_tave',0,0,'dar_Tlim_ave',Tlimave,myThid)
575  #endif  #endif
# Line 737  C      Reset averages to zero Line 802  C      Reset averages to zero
802           CALL TIMEAVE_RESET(PARave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(PARave,Nr,bi,bj,myThid)
803           CALL TIMEAVE_RESET(PPave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(PPave,Nr,bi,bj,myThid)
804  #ifdef QUOTA_DIAG_LIMIT  #ifdef QUOTA_DIAG_LIMIT
805           CALL TIMEAVE_RESET(Nlimave,Nr,bi,bj,myThid)  ! Carbon
806           CALL TIMEAVE_RESET(Flimave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(AP_C_ave,Nr,bi,bj,myThid)
807             CALL TIMEAVE_RESET(HP_C_ave,Nr,bi,bj,myThid)
808    ! Nitrogen
809             CALL TIMEAVE_RESET(AP_N_ave,Nr,bi,bj,myThid)
810             CALL TIMEAVE_RESET(HP_N_ave,Nr,bi,bj,myThid)
811             CALL TIMEAVE_RESET(Nlimave ,Nr,bi,bj,myThid)
812    ! Phosphorus
813    #ifdef PQUOTA
814             CALL TIMEAVE_RESET(AP_P_ave,Nr,bi,bj,myThid)
815             CALL TIMEAVE_RESET(HP_P_ave,Nr,bi,bj,myThid)
816             CALL TIMEAVE_RESET(Plimave ,Nr,bi,bj,myThid)
817    #endif
818    ! Iron
819    #ifdef FQUOTA
820             CALL TIMEAVE_RESET(AP_F_ave,Nr,bi,bj,myThid)
821             CALL TIMEAVE_RESET(HP_F_ave,Nr,bi,bj,myThid)
822             CALL TIMEAVE_RESET(Flimave ,Nr,bi,bj,myThid)
823    #endif
824           CALL TIMEAVE_RESET(Ilimave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(Ilimave,Nr,bi,bj,myThid)
825           CALL TIMEAVE_RESET(Tlimave,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(Tlimave,Nr,bi,bj,myThid)
826  #endif  #endif
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22