pkg/pp81/pp81_init_varia.F

08be60903a Mart*0001 #include "PP81_OPTIONS.h"
                0002
8861629416 Jean*0003       SUBROUTINE PP81_INIT_VARIA( myThid )
                0004 C     *==========================================================*
                0005 C     | SUBROUTINE PP81_INIT_VARIA
                0006 C     | o Routine to initialize PP81 parameters and variables.
                0007 C     *==========================================================*
                0008 C     | Initialize PP81 parameters and variables.
                0009 C     *==========================================================*
08be60903a Mart*0010       IMPLICIT NONE
                0011
                0012 C     === Global variables ===
                0013 #include "SIZE.h"
                0014 #include "EEPARAMS.h"
                0015 #include "PARAMS.h"
5e48dccc42 Jean*0016 c#include "GRID.h"
08be60903a Mart*0017 #include "PP81.h"
                0018
                0019 C     === Routine arguments ===
                0020 C     myThid -  Number of this instance of PP81_INIT
                0021       INTEGER myThid
                0022
                0023 #ifdef ALLOW_PP81
                0024
                0025 C     === Local variables ===
                0026 C     i,j,k,bi,bj - Loop counters
                0027
                0028       INTEGER i, j, k, bi, bj
                0029
                0030 C-----------------------------------------------------------------------
                0031 C Initialize PP variables PPviscAr, PPdiffKrT
                0032 C-----------------------------------------------------------------------
                0033
5e48dccc42 Jean*0034       DO bj = myByLo(myThid), myByHi(myThid)
                0035        DO bi = myBxLo(myThid), myBxHi(myThid)
                0036         DO k = 1, Nr
                0037          DO j = 1-OLy, sNy+OLy
                0038           DO i = 1-OLx, sNx+OLx
                0039            PPviscAr (i,j,k,bi,bj) = viscArNr(k)
78524d1402 Jean*0040            PPdiffKr (i,j,k,bi,bj) = diffKrNrS(k)
5e48dccc42 Jean*0041           ENDDO
                0042          ENDDO
                0043         ENDDO
                0044        ENDDO
                0045       ENDDO
08be60903a Mart*0046
                0047 #endif /* ALLOW_PP81 */
                0048
5e48dccc42 Jean*0049       RETURN
                0050       END