pkg/atm2d/yr_end_diags.F

8e101fde6e Jeff*0001 #include "ctrparam.h"
                0002 #include "ATM2D_OPTIONS.h"
                0003
                0004 C     !INTERFACE:
                0005       SUBROUTINE YR_END_DIAGS(  iyr, myTime, myIter, myThid )
                0006 C     *==========================================================*
                0007 C     | Normalize and dump year-end annual mean diagnostics      |
                0008 c     |                                                          |
                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 "ATM2D_VARS.h"
                0018
                0019
                0020 C     !INPUT/OUTPUT PARAMETERS:
                0021 C     === Routine arguments ===
                0022 C     iyr    - current year of simulation
                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       INTEGER iyr
                0027       _RL     myTime
                0028       INTEGER myIter
                0029       INTEGER myThid
                0030
                0031 C     LOCAL VARIABLES:
                0032       _RL secYr
9274434acc Jean*0033       DATA secYr /31536000.D0/
8e101fde6e Jeff*0034       _RL norm_factor
                0035
                0036       norm_factor = globalArea * secYr / dtcouplo
                0037       sst_tave = sst_tave / norm_factor
                0038       sss_tave = sss_tave / norm_factor
                0039       HF2ocn_tave = HF2ocn_tave / norm_factor
                0040       FW2ocn_tave = FW2ocn_tave / norm_factor
678e639be8 Jeff*0041       CO2flx_tave = CO2flx_tave / norm_factor
8e101fde6e Jeff*0042
                0043       OPEN(25,ACCESS='APPEND',FILE='resocean.dat',STATUS='old')
                0044
                0045       WRITE(25,*) '      '
                0046       WRITE(25,*) '~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~'
                0047       WRITE(25,*) 'YEAR MEAN', iyr
                0048       WRITE(25,*) 'SST SSS', sst_tave, sss_tave
                0049       WRITE(25,*) 'HFLUX FWFLUX',  HF2ocn_tave,  FW2ocn_tave
                0050       WRITE(25,*) 'SH ICE', SHice_max*1. _d -12, SHice_min*1. _d -12
                0051       WRITE(25,*) 'NH ICE', NHice_max*1. _d -12, NHice_min*1. _d -12
9274434acc Jean*0052       IF  (CO2flx_tave .NE. 0. _d 0)
678e639be8 Jeff*0053      &   WRITE(25,*) 'OCEAN CO2FLUX', CO2flx_tave
8e101fde6e Jeff*0054       WRITE(25,*) '~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~'
                0055       WRITE(25,*) '      '
                0056
                0057       CLOSE(25)
                0058
                0059       SHice_min = 1. _d 18
                0060       NHice_min = 1. _d 18
                0061       SHice_max = 0. _d 0
                0062       NHice_max = 0. _d 0
                0063       sst_tave =  0. _d 0
                0064       sss_tave =  0. _d 0
                0065       HF2ocn_tave =  0. _d 0
                0066       FW2ocn_tave =  0. _d 0
678e639be8 Jeff*0067       CO2flx_tave =  0. _d 0
8e101fde6e Jeff*0068
                0069       RETURN
                0070       END
                0071