pkg/dic/dic_set_control.F

e7fb6342b1 Patr*0001 #include "DIC_OPTIONS.h"
e2c9aaa8c1 Jean*0002 #ifdef ALLOW_CTRL
                0003 # include "CTRL_OPTIONS.h"
                0004 #endif
e7fb6342b1 Patr*0005
989b1a2fcf Jean*0006 cphc$taf COMMON DIC_XX        adname = addic_xx
                0007 cphc$taf COMMON DIC_COST_CTRL adname = ADDIC_COST_CTRL
e7fb6342b1 Patr*0008
                0009 C !INTERFACE: ==========================================================
989b1a2fcf Jean*0010       SUBROUTINE DIC_SET_CONTROL( myThid )
e7fb6342b1 Patr*0011
                0012 C !DESCRIPTION:
                0013
                0014 C !USES: ===============================================================
989b1a2fcf Jean*0015       IMPLICIT NONE
e7fb6342b1 Patr*0016
                0017 C     == GLobal variables ==
                0018 #include "SIZE.h"
                0019 #include "EEPARAMS.h"
                0020 #include "PARAMS.h"
                0021 #ifdef DIC_BIOTIC
                0022 # include "DIC_VARS.h"
                0023 # include "DIC_CTRL.h"
                0024 #endif
211f542c8c Patr*0025 #ifdef ALLOW_CTRL
37d5136090 Jean*0026 #include "CTRL_SIZE.h"
4d72283393 Mart*0027 #include "CTRL.h"
edcd27be69 Mart*0028 #include "CTRL_DUMMY.h"
211f542c8c Patr*0029 #endif
e7fb6342b1 Patr*0030 C     == Routine arguments ==
                0031       INTEGER myThid
                0032
c845fbfeae Jean*0033 #ifdef ALLOW_CTRL
e7fb6342b1 Patr*0034 cph#ifdef DIC_BIOTIC
                0035 C     == Local arguments ==
cf705a6c8e Mart*0036 #ifdef ALLOW_DIC_CONTROL
e7fb6342b1 Patr*0037       INTEGER bi, bj
                0038       INTEGER i, j
11c3150c71 Mart*0039 #endif
e7fb6342b1 Patr*0040
                0041 c     == end of interface ==
                0042 CEOP
                0043
a70f44bc89 Patr*0044 #ifdef ALLOW_DIC_CONTROL
                0045
989b1a2fcf Jean*0046       DO bj=myByLo(myThid),myByHi(myThid)
                0047        DO bi=myBxLo(myThid),myBxHi(myThid)
                0048          DO j = 1, sNy
                0049           DO i = 1, sNx
6acab690ae Jona*0050 c           feload(i,j,bi,bj) = feload(i,j,bi,bj)*(1. _d 0 +xx_dic(1))
989b1a2fcf Jean*0051             rain_ratio(i,j,bi,bj) =
                0052      &                      rain_ratio(i,j,bi,bj)*(1. _d 0 +xx_dic(2))
                0053           ENDDO
                0054          ENDDO
                0055        ENDDO
                0056       ENDDO
                0057
                0058       _EXCH_XY_RL( alpha, myThid )
                0059       _EXCH_XY_RL( rain_ratio, myThid )
6acab690ae Jona*0060 c     _EXCH_XY_RL( feload, myThid )
989b1a2fcf Jean*0061
                0062       KScav       = KScav       * ( 1. _d 0 + xx_dic(3)*1. _d 6 )
                0063       ligand_stab = ligand_stab * ( 1. _d 0 + xx_dic(4)*1. _d 6 )
                0064       ligand_tot  = ligand_tot  * ( 1. _d 0 + xx_dic(5)*1. _d 6 )
e7fb6342b1 Patr*0065
                0066       print *,'COST KScav = ', KScav
                0067       print *,'COST ligand_stab = ', ligand_stab
                0068       print *,'COST ligand_tot = ', ligand_tot
                0069
989b1a2fcf Jean*0070 #endif /* ALLOW_DIC_CONTROL */
211f542c8c Patr*0071
e7fb6342b1 Patr*0072 cph#endif /* DIC_BIOTIC */
6acab690ae Jona*0073 #endif /* ALLOW_CTRL */
e7fb6342b1 Patr*0074
989b1a2fcf Jean*0075       RETURN
                0076       END