pkg/atm2d/sum_yr_end_diags.F

8e101fde6e Jeff*0001 #include "ctrparam.h"
                0002 #include "ATM2D_OPTIONS.h"
                0003
                0004 C     !INTERFACE:
                0005       SUBROUTINE SUM_YR_END_DIAGS(  myTime, myIter, myThid )
                0006 C     *==========================================================*
                0007 C     | Tabulate year-end annual mean diagnostics                |
                0008 c     | Called at top of coupled period, after seaice growth/melt|
                0009 C     *==========================================================*
                0010         IMPLICIT NONE
                0011
                0012 C     === Global Atmosphere Variables ===
                0013 #include "ATMSIZE.h"
                0014 #include "SIZE.h"
                0015 #include "GRID.h"
                0016 #include "EEPARAMS.h"
                0017 #include "THSICE_VARS.h"
                0018 #include "ATM2D_VARS.h"
                0019
                0020
                0021 C     !INPUT/OUTPUT PARAMETERS:
                0022 C     === Routine arguments ===
                0023 C     myTime - current simulation time (ocean model time)
                0024 C     myIter - iteration number (ocean model)
                0025 C     myThid - Thread no. that called this routine.
                0026       _RL     myTime
                0027       INTEGER myIter
                0028       INTEGER myThid
                0029
                0030 C     LOCAL VARIABLES:
                0031       INTEGER i,j
                0032       _RL iceArea
                0033
                0034       iceArea = 0. _d 0
                0035       DO j=1,int(sNy/2)
                0036         DO i=1,sNx
9274434acc Jean*0037           IF (iceMask(i,j,1,1).GT. 0. _d 0)
8e101fde6e Jeff*0038      &        iceArea = iceArea + rA(i,j,1,1)*iceMask(i,j,1,1)
                0039         ENDDO
                0040       ENDDO
9274434acc Jean*0041       SHice_min=min(SHice_min,iceArea)
                0042       SHice_max=max(SHice_max,iceArea)
8e101fde6e Jeff*0043
                0044       iceArea = 0. _d 0
                0045       DO j=int(sNy/2),sNy
                0046         DO i=1,sNx
9274434acc Jean*0047           IF (iceMask(i,j,1,1).GT. 0. _d 0)
8e101fde6e Jeff*0048      &        iceArea = iceArea + rA(i,j,1,1)*iceMask(i,j,1,1)
                0049         ENDDO
                0050       ENDDO
9274434acc Jean*0051       NHice_min=min(NHice_min,iceArea)
                0052       NHice_max=max(NHice_max,iceArea)
8e101fde6e Jeff*0053
                0054       DO j=1,sNy
                0055         DO i=1,sNx
                0056
                0057           IF (maskC(i,j,1,1,1) .EQ. 1.) THEN
                0058             sst_tave= sst_tave + rA(i,j,1,1)*sstFromOcn(i,j)
                0059             sss_tave= sss_tave + rA(i,j,1,1)*sssFromOcn(i,j)
                0060             HF2ocn_tave = HF2ocn_tave + rA(i,j,1,1)*pass_qnet(i,j)
                0061             FW2ocn_tave = FW2ocn_tave + rA(i,j,1,1)*
                0062      &       (pass_precip(i,j) + pass_evap(i,j) + pass_runoff(i,j))
678e639be8 Jeff*0063
                0064             CO2flx_tave = CO2flx_tave + rA(i,j,1,1)*oFluxCO2(i,j)
8e101fde6e Jeff*0065           ENDIF
                0066
                0067         ENDDO
                0068       ENDDO
                0069
                0070       RETURN
                0071       END
                0072