pkg/thsice/thsice_ave.F

fc7306ba7d Jean*0001 #include "THSICE_OPTIONS.h"
                0002
87ea84cac6 Jean*0003 CBOP
fc7306ba7d Jean*0004 C     !ROUTINE: THSICE_AVE
                0005 C     !INTERFACE:
7269783f6f Jean*0006       SUBROUTINE THSICE_AVE(
41e6ab0e43 Jean*0007      I                       bi, bj, myTime, myIter, myThid )
87ea84cac6 Jean*0008 C     !DESCRIPTION: \bv
fc7306ba7d Jean*0009 C     *==========================================================*
7269783f6f Jean*0010 C     | S/R  THSICE_AVE
b7411f1a84 Jean*0011 C     | o fill some thsice diagnostics
fc7306ba7d Jean*0012 C     *==========================================================*
                0013 C     \ev
                0014
                0015 C     !USES:
                0016       IMPLICIT NONE
87ea84cac6 Jean*0017
fc7306ba7d Jean*0018 C     == Global variables ==
                0019 #include "SIZE.h"
87ea84cac6 Jean*0020 #include "EEPARAMS.h"
                0021 #include "PARAMS.h"
                0022 #include "FFIELDS.h"
fc7306ba7d Jean*0023 #include "THSICE_PARAMS.h"
87ea84cac6 Jean*0024 #include "THSICE_VARS.h"
fc7306ba7d Jean*0025
                0026 C     !INPUT/OUTPUT PARAMETERS:
87ea84cac6 Jean*0027 C     frwAtm  :: net fresh-water flux (E-P) to the atmosphere [kg/m2/s]
                0028 C     flxAtm  :: net heat flux from the atmosphere (+=down) [W/m2]
                0029 C     bi,bj   :: tile indices
41e6ab0e43 Jean*0030 C     myTime  :: time counter for this thread
87ea84cac6 Jean*0031 C     myIter  :: iteration counter for this thread
                0032 C     myThid  :: thread number for this instance of the routine.
7269783f6f Jean*0033 c     _RL frwAtm(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
                0034 c     _RL flxAtm(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
fc7306ba7d Jean*0035       INTEGER bi,bj
87ea84cac6 Jean*0036       _RL  myTime
                0037       INTEGER myIter
fc7306ba7d Jean*0038       INTEGER myThid
87ea84cac6 Jean*0039 CEOP
fc7306ba7d Jean*0040
                0041 #ifdef ALLOW_THSICE
b7411f1a84 Jean*0042 C     !FUNCTIONS:
41e6ab0e43 Jean*0043 c     LOGICAL  DIFFERENT_MULTIPLE
                0044 c     EXTERNAL DIFFERENT_MULTIPLE
fc7306ba7d Jean*0045
b7411f1a84 Jean*0046 C     !LOCAL VARIABLES:
                0047 c     INTEGER wrtIter
87ea84cac6 Jean*0048
7269783f6f Jean*0049 c     IF ( DIFFERENT_MULTIPLE(thSIce_diagFreq,
                0050 c    &                        myTime+deltaTClock,deltaTClock)
                0051 c    &     .AND. fluidIsWater ) THEN
                0052 c       wrtIter = myIter+1
                0053 c       CALL WRITE_LOCAL_RL('ice_flxAtm','I10',1,flxAtm,
                0054 c    &                      bi,bj,1,wrtIter,myThid)
                0055 c       CALL WRITE_LOCAL_RL('ice_frwAtm','I10',1,frwAtm,
                0056 c    &                      bi,bj,1,wrtIter,myThid)
                0057 c     ENDIF
87ea84cac6 Jean*0058
040c28a424 Jean*0059 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0060
                0061 #ifdef ALLOW_DIAGNOSTICS
ddc7b8ee22 Jean*0062       IF ( useDiagnostics ) THEN
040c28a424 Jean*0063
                0064 C-- Atmospheric fluxes :
                0065        IF ( fluidIsWater ) THEN
7269783f6f Jean*0066         CALL DIAGNOSTICS_FILL(icFlxAtm,'SIflxAtm',0,1,1,bi,bj,myThid)
                0067         CALL DIAGNOSTICS_FILL(icFrwAtm,'SIfrwAtm',0,1,1,bi,bj,myThid)
040c28a424 Jean*0068        ENDIF
                0069
                0070 C-- Oceanic fluxes :
                0071 c #ifndef CHECK_ENERGY_CONSERV
41e6ab0e43 Jean*0072        CALL DIAGNOSTICS_FILL_RS(Qnet,    'SIflx2oc',0,1,1,bi,bj,myThid)
                0073        CALL DIAGNOSTICS_FILL_RS(EmPmR,   'SIfrw2oc',0,1,1,bi,bj,myThid)
                0074        CALL DIAGNOSTICS_FILL_RS(saltFlux,'SIsaltFx',0,1,1,bi,bj,myThid)
040c28a424 Jean*0075 c #endif /* CHECK_ENERGY_CONSERV */
                0076
ddc7b8ee22 Jean*0077       ENDIF
040c28a424 Jean*0078 #endif /* ALLOW_DIAGNOSTICS */
ddc7b8ee22 Jean*0079
fc7306ba7d Jean*0080 #endif  /* ALLOW_THSICE */
                0081
                0082       RETURN
                0083       END