/[MITgcm]/MITgcm_contrib/jscott/pkg_atm2d/month_end_diags.F
ViewVC logotype

Diff of /MITgcm_contrib/jscott/pkg_atm2d/month_end_diags.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by jscott, Fri Aug 11 18:55:50 2006 UTC revision 1.2 by jscott, Tue Aug 22 20:21:39 2006 UTC
# Line 33  C     LOCAL VARIABLES: Line 33  C     LOCAL VARIABLES:
33        INTEGER i,j        INTEGER i,j
34        _RS norm_factor        _RS norm_factor
35    
36          norm_factor = ndmonth(monid)*86400.0
37    
38        DO j=1,sNy        DO j=1,sNy
39          DO i=1,sNx          DO i=1,sNx
# Line 49  C     LOCAL VARIABLES: Line 50  C     LOCAL VARIABLES:
50       &                            sum_qrel(i,j)       &                            sum_qrel(i,j)
51            sum_frel_ta(i,j,monid)= sum_frel_ta(i,j,monid) +            sum_frel_ta(i,j,monid)= sum_frel_ta(i,j,monid) +
52       &                            sum_frel(i,j)       &                            sum_frel(i,j)
53              sum_iceMask_ta(i,j,monid)= sum_iceMask_ta(i,j,monid) +
54         &                               sum_iceMask(i,j)
55              sum_iceHeight_ta(i,j,monid)= sum_iceHeight_ta(i,j,monid) +
56         &                                 sum_iceHeight(i,j)
57              sum_iceTime_ta(i,j,monid)= sum_iceTime_ta(i,j,monid) +
58         &                               sum_iceTime(i,j)
59              sum_oceMxLT_ta(i,j,monid)= sum_oceMxLT_ta(i,j,monid) +
60         &                               sum_oceMxLT(i,j)
61              sum_oceMxLS_ta(i,j,monid)= sum_oceMxLS_ta(i,j,monid) +
62         &                               sum_oceMxLS(i,j)
63    
           norm_factor = ndmonth(monid)*86400.0  
64            qnet_atm(i,j)=qnet_atm(i,j)/ norm_factor            qnet_atm(i,j)=qnet_atm(i,j)/ norm_factor
65            evap_atm(i,j)=evap_atm(i,j)/ norm_factor            evap_atm(i,j)=evap_atm(i,j)/ norm_factor
66            precip_atm(i,j)= precip_atm(i,j)/ norm_factor            precip_atm(i,j)= precip_atm(i,j)/ norm_factor
67            runoff_atm(i,j)=  runoff_atm(i,j)/ norm_factor            runoff_atm(i,j)=  runoff_atm(i,j)/ norm_factor
68            sum_qrel(i,j)= sum_qrel(i,j)/ norm_factor            sum_qrel(i,j)= sum_qrel(i,j)/ norm_factor
69            sum_frel(i,j)= sum_frel(i,j)/ norm_factor            sum_frel(i,j)= sum_frel(i,j)/ norm_factor
70              sum_iceMask(i,j)= sum_iceMask(i,j)/ norm_factor
71              sum_iceHeight(i,j)= sum_iceHeight(i,j)/ norm_factor
72              sum_iceTime(i,j)= sum_iceTime(i,j)/ norm_factor
73              sum_oceMxLT(i,j)= sum_oceMxLT(i,j)/ norm_factor
74              sum_oceMxLS(i,j)= sum_oceMxLS(i,j)/ norm_factor
75    
76          ENDDO          ENDDO
77        ENDDO        ENDDO
# Line 74  C     LOCAL VARIABLES: Line 89  C     LOCAL VARIABLES:
89       &                     myIter, myThid)       &                     myIter, myThid)
90        CALL WRITE_FLD_XY_RL('mFrelfluxtave.', suff, sum_frel,        CALL WRITE_FLD_XY_RL('mFrelfluxtave.', suff, sum_frel,
91       &                     myIter, myThid)       &                     myIter, myThid)
92          CALL WRITE_FLD_XY_RL('mIceMasktave.', suff, sum_iceMask,
93         &                     myIter, myThid)
94          CALL WRITE_FLD_XY_RL('mIceHeighttave.', suff, sum_iceHeight,
95         &                     myIter, myThid)
96          CALL WRITE_FLD_XY_RL('mIceTimetave.', suff, sum_iceTime,
97         &                     myIter, myThid)
98          CALL WRITE_FLD_XY_RL('mOceMxLTtave.', suff, sum_oceMxLT,
99         &                     myIter, myThid)
100          CALL WRITE_FLD_XY_RL('mOceMxLStave.', suff, sum_oceMxLS,
101         &                     myIter, myThid)
102    
103        DO j=1,sNy        DO j=1,sNy
104          DO i=1,sNx          DO i=1,sNx
# Line 84  C     LOCAL VARIABLES: Line 109  C     LOCAL VARIABLES:
109            runoff_atm(i,j)= 0. _d 0            runoff_atm(i,j)= 0. _d 0
110            sum_qrel(i,j)= 0. _d 0            sum_qrel(i,j)= 0. _d 0
111            sum_frel(i,j)= 0. _d 0            sum_frel(i,j)= 0. _d 0
112              sum_iceMask(i,j)= 0. _d 0
113              sum_iceHeight(i,j)= 0. _d 0
114              sum_iceTime(i,j)= 0. _d 0
115              sum_oceMxLT(i,j)= 0. _d 0
116              sum_oceMxLS(i,j)= 0. _d 0
117    
118          ENDDO          ENDDO
119        ENDDO        ENDDO
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22