pkg/autodiff/autodiff_init_varia.F

6d54cf9ca1 Ed H*0001 #include "AUTODIFF_OPTIONS.h"
72a14826a0 Gael*0002 #ifdef ALLOW_CTRL
                0003 # include "CTRL_OPTIONS.h"
                0004 #endif
09d2417e9e Patr*0005
                0006 CBOP
6cb1079f9d Jean*0007 C     !ROUTINE: AUTODIFF_INIT_VARIA
09d2417e9e Patr*0008 C     !INTERFACE:
6cb1079f9d Jean*0009       SUBROUTINE AUTODIFF_INIT_VARIA( myThid )
                0010
09d2417e9e Patr*0011 C     !DESCRIPTION: \bv
                0012 C     *==========================================================*
6cb1079f9d Jean*0013 C     | SUBROUTINE AUTODIFF_INIT_VARIA
                0014 C     | o Initialise to zero some active arrays
09d2417e9e Patr*0015 C     *==========================================================*
                0016 C     \ev
                0017
                0018 C     !USES:
                0019       IMPLICIT NONE
                0020 C     === Global variables ===
                0021 #include "SIZE.h"
                0022 #include "EEPARAMS.h"
                0023 #include "PARAMS.h"
c88a41bb8a Patr*0024 #include "FFIELDS.h"
bef5cbe0ec Patr*0025 #include "SURFACE.h"
09d2417e9e Patr*0026
                0027 C     !INPUT/OUTPUT PARAMETERS:
                0028 C     == Routine arguments ==
6cb1079f9d Jean*0029 C     myThid :: My Thread Id number
09d2417e9e Patr*0030       INTEGER myThid
                0031
6cb1079f9d Jean*0032 #ifdef ALLOW_AUTODIFF
09d2417e9e Patr*0033 C     !LOCAL VARIABLES:
                0034 C     == Local variables ==
6cb1079f9d Jean*0035 C     bi,bj  :: tile indices
                0036 C     i,j,k  :: Loop counters
09d2417e9e Patr*0037       INTEGER bi, bj
cf705a6c8e Mart*0038 #ifdef ALLOW_EP_FLUX
6cb1079f9d Jean*0039       INTEGER i, j, k
9f5240b52a Jean*0040 #endif
09d2417e9e Patr*0041 CEOP
                0042
bef5cbe0ec Patr*0043 C--   Scalar fields
6cb1079f9d Jean*0044       TsurfCor = 0. _d 0
                0045       SsurfCor = 0. _d 0
bef5cbe0ec Patr*0046
09d2417e9e Patr*0047 C--   Over all tiles
                0048       DO bj = myByLo(myThid), myByHi(myThid)
                0049        DO bi = myBxLo(myThid), myBxHi(myThid)
                0050
                0051 C-      3D arrays
c88a41bb8a Patr*0052 #ifdef ALLOW_EP_FLUX
6cb1079f9d Jean*0053         DO k=1,Nr
a19fa93b10 Jean*0054          DO j=1-OLy,sNy+OLy
                0055           DO i=1-OLx,sNx+OLx
6cb1079f9d Jean*0056            EfluxY(i,j,k,bi,bj) = 0.
                0057            EfluxP(i,j,k,bi,bj) = 0.
c88a41bb8a Patr*0058           ENDDO
                0059          ENDDO
                0060         ENDDO
6cb1079f9d Jean*0061 #endif
09d2417e9e Patr*0062
                0063        ENDDO
                0064       ENDDO
                0065
6cb1079f9d Jean*0066 #endif /* ALLOW_AUTODIFF */
                0067       RETURN
09d2417e9e Patr*0068       END