pkg/seaice/seaice_cost_init_varia.F

869864d4b6 Patr*0001 #include "SEAICE_OPTIONS.h"
                0002
2b959ba38e Mart*0003       SUBROUTINE SEAICE_COST_INIT_VARIA( myThid )
869864d4b6 Patr*0004
2b959ba38e Mart*0005 C     ==================================================================
                0006 C     SUBROUTINE SEAICE_COST_INIT_VARIA
                0007 C     ==================================================================
                0008 C     o Initialise the variable cost function part.
                0009 C     added sea-ice term: menemenlis@jpl.nasa.gov 26-Feb-2003
                0010 C
                0011 C     ==================================================================
869864d4b6 Patr*0012
2b959ba38e Mart*0013       IMPLICIT NONE
869864d4b6 Patr*0014
2b959ba38e Mart*0015 C     == global variables ==
869864d4b6 Patr*0016 #include "EEPARAMS.h"
                0017 #include "SIZE.h"
dcf0f3bd0f Patr*0018 #ifdef ALLOW_COST
03c669d1ab Jean*0019 # include "SEAICE_SIZE.h"
dcf0f3bd0f Patr*0020 # include "SEAICE_COST.h"
                0021 #endif
869864d4b6 Patr*0022
2b959ba38e Mart*0023 C     == routine arguments ==
                0024       INTEGER myThid
869864d4b6 Patr*0025
dcf0f3bd0f Patr*0026 #ifdef ALLOW_COST
2b959ba38e Mart*0027 C     == local variables ==
                0028       INTEGER bi,bj
                0029 #ifdef ALLOW_SEAICE_COST_EXPORT
                0030       INTEGER i,j
                0031 #endif
                0032 C     == end of interface ==
869864d4b6 Patr*0033
2b959ba38e Mart*0034 C--   Initialize the tiled cost function contributions.
                0035       DO bj = myByLo(myThid), myByHi(myThid)
                0036         DO bi = myBxLo(myThid), myBxHi(myThid)
                0037           objf_ice(bi,bj)        = 0. _d 0
1e22f5fc71 Patr*0038           objf_ice_export(bi,bj) = 0. _d 0
2b959ba38e Mart*0039           num_ice(bi,bj)         = 0. _d 0
1e22f5fc71 Patr*0040 #ifdef ALLOW_SEAICE_COST_EXPORT
2b959ba38e Mart*0041           DO j=1-OLy,sNy+OLy
                0042             DO i=1-OLx,sNx+OLx
1e22f5fc71 Patr*0043                uHeffExportCell(i,j,bi,bj) = 0. _d 0
                0044                vHeffExportCell(i,j,bi,bj) = 0. _d 0
2b959ba38e Mart*0045             ENDDO
                0046           ENDDO
869864d4b6 Patr*0047 #endif
2b959ba38e Mart*0048         ENDDO
                0049       ENDDO
869864d4b6 Patr*0050
9f5240b52a Jean*0051 #endif /* ALLOW_COST */
dcf0f3bd0f Patr*0052
2b959ba38e Mart*0053       RETURN
                0054       END