pkg/seaice/seaice_cost_accumulate_mean.F

1e22f5fc71 Patr*0001 #include "SEAICE_OPTIONS.h"
                0002
2b959ba38e Mart*0003 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0004 CBOP
                0005 C     !ROUTINE: SEAICE_COST_ACCUMULATE_MEAN
                0006
                0007 C     !INTERFACE:
989b1a2fcf Jean*0008       SUBROUTINE SEAICE_COST_ACCUMULATE_MEAN(
                0009      &           myTime, myIter, myThid )
1e22f5fc71 Patr*0010
2b959ba38e Mart*0011 C     !DESCRIPTION:
                0012 C     Accumulate time-mean state for cost evalualtion
                0013
                0014 C     !USES:
                0015       IMPLICIT NONE
1e22f5fc71 Patr*0016 C     == Global variables ===
                0017 #include "SIZE.h"
                0018 #include "EEPARAMS.h"
                0019 #include "PARAMS.h"
                0020 #include "GRID.h"
03c669d1ab Jean*0021 #include "SEAICE_SIZE.h"
1e22f5fc71 Patr*0022 #include "SEAICE.h"
a7cac4d824 Patr*0023 #ifdef ALLOW_COST
2b959ba38e Mart*0024 # include "SEAICE_COST.h"
a7cac4d824 Patr*0025 # include "cost.h"
                0026 #endif
1e22f5fc71 Patr*0027
2b959ba38e Mart*0028 C     !INPUT/OUTPUT PARAMETERS:
                0029 C     myThid ::  my Thread Id number
57867cd15e Patr*0030       _RL myTime
989b1a2fcf Jean*0031       INTEGER myIter
                0032       INTEGER myThid
2b959ba38e Mart*0033 CEOP
1e22f5fc71 Patr*0034
                0035 #ifdef ALLOW_COST
46ba63dfb8 Gael*0036 #ifdef ALLOW_SEAICE_COST_EXPORT
                0037
2b959ba38e Mart*0038 C     !LOCAL VARIABLES:
                0039       INTEGER i, j, bi, bj
1e22f5fc71 Patr*0040
57867cd15e Patr*0041       IF ( myTime .GT. (EndTime - lastinterval) ) THEN
1736eae30e Patr*0042
1e22f5fc71 Patr*0043 C--   Calculate cost function on tile of this instance
2b959ba38e Mart*0044        DO bj = myByLo(myThid), myByHi(myThid)
                0045         DO bi = myBxLo(myThid), myBxHi(myThid)
989b1a2fcf Jean*0046          DO j=1,sNy
                0047           DO i=1,sNx
                0048             uHeffExportCell(i,j,bi,bj) = uHeffExportCell(i,j,bi,bj)
                0049      &      + ( 0.92 _d 0*(HEFF(i,j,bi,bj)+HEFF(i-1,j,bi,bj))*0.5 _d 0
                0050      &        + 0.33 _d 0*(HSNOW(i,j,bi,bj)+HSNOW(i-1,j,bi,bj))*0.5 _d 0
                0051      &        )*UICE(i,j,bi,bj)
                0052      &         *_dyG(i,j,bi,bj)/lastinterval*deltaTClock
                0053             vHeffExportCell(i,j,bi,bj) = vHeffExportCell(i,j,bi,bj)
                0054      &      + ( 0.92 _d 0*(HEFF(i,j,bi,bj)+HEFF(i,j-1,bi,bj))*0.5 _d 0
                0055      &        + 0.33 _d 0*(HSNOW(i,j,bi,bj)+HSNOW(i,j-1,bi,bj))*0.5 _d 0
                0056      &        )*VICE(i,j,bi,bj)
                0057      &         *_dxG(i,j,bi,bj)/lastinterval*deltaTClock
                0058           ENDDO
                0059          ENDDO
                0060         ENDDO
                0061        ENDDO
1e22f5fc71 Patr*0062
57867cd15e Patr*0063       ENDIF
                0064
989b1a2fcf Jean*0065 #endif /* ALLOW_SEAICE_COST_EXPORT */
                0066 #endif /* ALLOW_COST */
1e22f5fc71 Patr*0067
989b1a2fcf Jean*0068       RETURN
1e22f5fc71 Patr*0069       END