File indexing completed on 2023-03-03 06:09:35 UTC
view on githubraw file Latest commit 4125e77f on 2022-12-27 05:51:21 UTC
8fe5594931 Patr*0001 #include "CPP_OPTIONS.h"
4f7b6674b6 Patr*0002
9366854e02 Chri*0003
0004
0005
8fe5594931 Patr*0006 SUBROUTINE INI_MIXING( myThid )
9366854e02 Chri*0007
0008
67be4914b2 Jean*0009
0010
9366854e02 Chri*0011
0012
8fe5594931 Patr*0013
9366854e02 Chri*0014
0015 IMPLICIT NONE
8fe5594931 Patr*0016
0017 #include "SIZE.h"
0018 #include "EEPARAMS.h"
0019 #include "PARAMS.h"
0020 #include "GRID.h"
0021 #include "DYNVARS.h"
0022
9366854e02 Chri*0023
4125e77f3d Jean*0024
8fe5594931 Patr*0025 INTEGER myThid
0026
9366854e02 Chri*0027
4125e77f3d Jean*0028 #if (defined ALLOW_3D_DIFFKR) || (defined ALLOW_BL79_LAT_VARY)
239f048bdf Jean*0029
0030
8fe5594931 Patr*0031 INTEGER bi, bj
4125e77f3d Jean*0032 INTEGER i, j
0033 #endif
0034 #ifdef ALLOW_3D_DIFFKR
0035 INTEGER k
0036 #endif
9366854e02 Chri*0037
8fe5594931 Patr*0038
94a8024bbe Jean*0039 #ifdef ALLOW_3D_DIFFKR
8fe5594931 Patr*0040 DO bj = myByLo(myThid), myByHi(myThid)
0041 DO bi = myBxLo(myThid), myBxHi(myThid)
8fc5af62b8 Jean*0042 DO k=1,Nr
517dbdc414 Jean*0043 DO j = 1-OLy, sNy+OLy
0044 DO i = 1-OLx, sNx+OLx
8fc5af62b8 Jean*0045 diffKr(i,j,k,bi,bj) = diffKrNrS(k)
8fe5594931 Patr*0046 ENDDO
0047 ENDDO
0048 ENDDO
0049 ENDDO
0050 ENDDO
94a8024bbe Jean*0051 IF ( diffKrFile .NE. ' ' ) THEN
0052 CALL READ_FLD_XYZ_RL( diffKrFile, ' ', diffKr, 0, myThid )
0053 _EXCH_XYZ_RL( diffKr, myThid )
0054 ENDIF
0055 #endif /* ALLOW_3D_DIFFKR */
8fe5594931 Patr*0056
e40c34e398 Dimi*0057 #ifdef ALLOW_BL79_LAT_VARY
0058 DO bj = myByLo(myThid), myByHi(myThid)
0059 DO bi = myBxLo(myThid), myBxHi(myThid)
517dbdc414 Jean*0060 DO j = 1-OLy, sNy+OLy
0061 DO i = 1-OLx, sNx+OLx
e40c34e398 Dimi*0062 BL79LatArray(i,j,bi,bj) = 0.0 _d 0
86e9b00011 Jean*0063 IF ( ABS(YC(i,j,bi,bj)) .LT. BL79LatVary )
4125e77f3d Jean*0064 & BL79LatArray(i,j,bi,bj) = halfRL
0065 & *( 1. _d 0 + COS(YC(i,j,bi,bj)*PI/BL79LatVary) )
e40c34e398 Dimi*0066 ENDDO
0067 ENDDO
0068 ENDDO
0069 ENDDO
0070 #endif
0071
8fe5594931 Patr*0072 RETURN
0073 END