File indexing completed on 2018-03-02 18:41:47 UTC
view on githubraw file Latest commit add29e06 on 2018-01-31 20:35:05 UTC
fda3710353 Oliv*0001 #include "LONGSTEP_OPTIONS.h"
0002
0003
0004
0005
0006 SUBROUTINE LONGSTEP_FILL_3D_FAC(
0007 U lsCnt, lsFld,
0008 I fld, fac, myNr, myThid)
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018 IMPLICIT NONE
0019
0020 #include "SIZE.h"
0021 #include "EEPARAMS.h"
0022 #include "PARAMS.h"
0023 #include "GRID.h"
0024
0025
0026
0027
0028
0029
0030
0031
0032 INTEGER myNr
0033 _RL lsFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
0034 _RL Fld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
9637a20806 Jean*0035 _RS Fac(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
fda3710353 Oliv*0036 INTEGER lsCnt(nSx,nSy)
0037 INTEGER myThid
0038
0039 #ifdef ALLOW_LONGSTEP
0040
0041
0042 INTEGER i,j,k,bi,bj
0043
0044
0045 DO bj=myByLo(myThid), myByHi(myThid)
0046 DO bi=myBxLo(myThid), myBxHi(myThid)
0047 lsCnt(bi,bj) = lsCnt(bi,bj) + 1
0048 DO k = 1,myNr
0049 DO j = 1-OLy,sNy+OLy
0050 DO i = 1-OLx,sNx+OLx
0051 lsFld(i,j,k,bi,bj) = lsFld(i,j,k,bi,bj)
0052 & + fld(i,j,k,bi,bj)*fac(i,j,k,bi,bj)
0053 ENDDO
0054 ENDDO
0055 ENDDO
0056 ENDDO
0057 ENDDO
0058 #endif /* ALLOW_LONGSTEP */
0059
0060 RETURN
0061 END
0062