pkg/shelfice/shelfice_cost_accumulate.F

359d9ec3d7 Patr*0001 #include "SHELFICE_OPTIONS.h"
                0002
                0003       subroutine shelfice_cost_accumulate( myThid )
                0004 C     *==========================================================*
                0005 C     | subroutine shelfice_cost_accumulate                      |
                0006 C     | o accumulate mean state for cost evalualtion             |
                0007 C     *==========================================================*
                0008 C     |                                                          |
                0009 C     *==========================================================*
                0010       IMPLICIT NONE
                0011
                0012 C     == Global variables ===
                0013 #include "SIZE.h"
                0014 #include "EEPARAMS.h"
                0015 #include "PARAMS.h"
                0016 #include "DYNVARS.h"
                0017 #include "GRID.h"
                0018 #include "SHELFICE.h"
                0019 #include "SHELFICE_COST.h"
                0020
                0021 C     == Routine arguments ==
                0022 C     myThid - Thread number for this instance of the routine.
                0023       integer myThid
                0024
                0025 #ifdef ALLOW_COST
                0026 C     == Local variables
d28329b3d8 Jean*0027       integer bi, bj
                0028       integer i, j
359d9ec3d7 Patr*0029
                0030 C--   Calculate cost function on tile of this instance
9f5240b52a Jean*0031       do bj = myByLo(myThid), myByHi(myThid)
                0032        do bi = myBxLo(myThid), myBxHi(myThid)
                0033         do j=1,sNy
                0034          do i=1,sNx
                0035           cMeanSHIforT(i,j,bi,bj) = cMeanSHIforT(i,j,bi,bj)
359d9ec3d7 Patr*0036      &         + shelficeHeatFlux(i,j,bi,bj)*deltaTClock
9f5240b52a Jean*0037           cMeanSHIforS(i,j,bi,bj) = cMeanSHIforS(i,j,bi,bj)
359d9ec3d7 Patr*0038      &         + shelficeFreshWaterFlux(i,j,bi,bj)*deltaTClock
d28329b3d8 Jean*0039          enddo
9f5240b52a Jean*0040         enddo
d28329b3d8 Jean*0041        enddo
                0042       enddo
359d9ec3d7 Patr*0043
                0044 #endif /* ALLOW_COST */
                0045
d28329b3d8 Jean*0046       RETURN
359d9ec3d7 Patr*0047       END