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_AVERAGE_3D_FAC(
0007 I lsCnt,
0008 U lsFld,
0009 I fac,
0010 I myNr, myThid)
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020 IMPLICIT NONE
0021
0022 #include "SIZE.h"
0023 #include "EEPARAMS.h"
0024 #include "PARAMS.h"
0025 #include "GRID.h"
0026
0027
0028
0029
0030
0031
0032
0033 INTEGER myNr
0034 _RL lsFld(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 _RL factor
0044
0045
0046 DO bj=myByLo(myThid), myByHi(myThid)
0047 DO bi=myBxLo(myThid), myBxHi(myThid)
0048 DO k = 1,myNr
0049 DO j = 1-OLy,sNy+OLy
0050 DO i = 1-OLx,sNx+OLx
b4566ab81c Jean*0051 factor = lsCnt(bi,bj)
0052 factor = factor*fac(i,j,k,bi,bj)
0053 IF (factor.NE.zeroRL) factor = 1. _d 0 / factor
fda3710353 Oliv*0054 lsFld(i,j,k,bi,bj) = lsFld(i,j,k,bi,bj)*factor
0055 ENDDO
0056 ENDDO
0057 ENDDO
0058 ENDDO
0059 ENDDO
0060 #endif /* ALLOW_LONGSTEP */
0061
0062 RETURN
0063 END
0064