pkg/thsice/thsice_cost_final.F

6c747cb1b2 Patr*0001 #include "THSICE_OPTIONS.h"
                0002
e4eb5b7d6f Jean*0003 CBOP
                0004 C     !ROUTINE: THSICE_COST_FINAL
                0005 C     !INTERFACE:
                0006       SUBROUTINE THSICE_COST_FINAL( myThid )
6c747cb1b2 Patr*0007
e4eb5b7d6f Jean*0008 C     !DESCRIPTION:
                0009 C     *==========================================================*
                0010 C     | SUBROUTINE THSICE_COST_FINAL
                0011 C     *==========================================================*
6c747cb1b2 Patr*0012
e4eb5b7d6f Jean*0013 C     !USES:
                0014       IMPLICIT NONE
6c747cb1b2 Patr*0015
e4eb5b7d6f Jean*0016 C     == global variables ==
6c747cb1b2 Patr*0017 #include "EEPARAMS.h"
                0018 #include "SIZE.h"
                0019 #include "PARAMS.h"
                0020 #include "DYNVARS.h"
                0021 #include "THSICE_PARAMS.h"
                0022 #ifdef ALLOW_COST
65754df434 Mart*0023 # include "THSICE_COST.h"
                0024 # include "cost.h"
6c747cb1b2 Patr*0025 #endif
                0026
e4eb5b7d6f Jean*0027 C     !INPUT/OUTPUT PARAMETERS:
                0028       INTEGER myThid
6c747cb1b2 Patr*0029
                0030 #ifdef ALLOW_COST
e4eb5b7d6f Jean*0031 C     !LOCAL VARIABLES:
                0032       INTEGER bi,bj
6c747cb1b2 Patr*0033       _RL f_thsice
                0034       _RL no_thsice
e4eb5b7d6f Jean*0035 CEOP
6c747cb1b2 Patr*0036
e4eb5b7d6f Jean*0037 c     print *, 'ph-2 in thsice_cost_final'
6c747cb1b2 Patr*0038
                0039       f_thsice      = 0. _d 0
                0040       no_thsice     = 0. _d 0
                0041
e4eb5b7d6f Jean*0042 C--   Sum up all contributions.
                0043       DO bj = myByLo(myThid), myByHi(myThid)
                0044        DO bi = myBxLo(myThid), myBxHi(myThid)
6c747cb1b2 Patr*0045
e4eb5b7d6f Jean*0046           tile_fc(bi,bj) = tile_fc(bi,bj)
6c747cb1b2 Patr*0047      &          + mult_thsice        * objf_thsice(bi,bj)
                0048
                0049           f_thsice = f_thsice + objf_thsice(bi,bj)
                0050
e4eb5b7d6f Jean*0051        ENDDO
                0052       ENDDO
6c747cb1b2 Patr*0053
e4eb5b7d6f Jean*0054 C--   Do global summation for each part of the cost function
                0055 c     CALL GLOBAL_SUM_TILE_RL( objf_thsice, f_thsice, myThid )
                0056 c     CALL GLOBAL_SUM_TILE_RL( num_thsice, no_thsice, myThid )
6c747cb1b2 Patr*0057
e4eb5b7d6f Jean*0058       WRITE(standardMessageUnit,'(A,D22.15)')
6c747cb1b2 Patr*0059      &     ' --> f_thsice     = ', f_thsice
                0060
                0061 #endif /* ALLOW_COST */
                0062
e4eb5b7d6f Jean*0063       RETURN
                0064       END