File indexing completed on 2018-03-02 18:41:48 UTC
view on githubraw file Latest commit add29e06 on 2018-01-31 20:35:05 UTC
5d0c92356b Jean*0001 #include "MATRIX_OPTIONS.h"
8f6dc1243d Jean*0002
2cd7c28080 Jean*0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013 SUBROUTINE MATRIX_STORE_TENDENCY_EXP(
0014 I iTracer, bi, bj,
66e282e41e Jean*0015 U tracNew,
2cd7c28080 Jean*0016 I myTime, myIter, myThid )
0017
0018
8385b3bb5f Ed H*0019
0020
0021
0022
a9544d2f65 Sama*0023
2cd7c28080 Jean*0024
a9544d2f65 Sama*0025 IMPLICIT NONE
8f6dc1243d Jean*0026 #include "SIZE.h"
a9544d2f65 Sama*0027 #include "EEPARAMS.h"
0028 #include "PARAMS.h"
8f6dc1243d Jean*0029 #include "PTRACERS_SIZE.h"
0030 #include "PTRACERS_PARAMS.h"
0031 #include "MATRIX.h"
0032
2cd7c28080 Jean*0033
66e282e41e Jean*0034
2cd7c28080 Jean*0035 INTEGER iTracer
a9544d2f65 Sama*0036 INTEGER bi, bj
66e282e41e Jean*0037 _RL myTime
0038 _RL tracNew(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
a9544d2f65 Sama*0039 INTEGER myIter
0040 INTEGER myThid
0041
0042 #ifdef ALLOW_MATRIX
2cd7c28080 Jean*0043
8f6dc1243d Jean*0044 INTEGER i,j,k
2cd7c28080 Jean*0045
a9544d2f65 Sama*0046
66e282e41e Jean*0047 DO k=1,Nr
0048 DO j=1-OLy, sNy+OLy
0049 DO i=1-OLx, sNx+OLx
0050 MATRIX(i,j,k,bi,bj,iTracer,1) = MATRIX(i,j,k,bi,bj,iTracer,1)
0051 & + ( tracNew(i,j,k) - PTRACERS_initial(i,j,k,bi,bj,iTracer) )
0052 tracNew(i,j,k) = PTRACERS_initial(i,j,k,bi,bj,iTracer)
a9544d2f65 Sama*0053 ENDDO
66e282e41e Jean*0054 ENDDO
0055 ENDDO
2cd7c28080 Jean*0056
0057
0058 IF ( iTracer.EQ.PTRACERS_numInUse ) THEN
0059 expMatrixCounter(bi,bj) = expMatrixCounter(bi,bj) + 1
0060 ENDIF
8f6dc1243d Jean*0061
a9544d2f65 Sama*0062 #endif /* ALLOW_MATRIX */
8f6dc1243d Jean*0063 RETURN
a9544d2f65 Sama*0064 END
8f6dc1243d Jean*0065
5d0c92356b Jean*0066
2cd7c28080 Jean*0067
0068
0069
0070
0071 SUBROUTINE MATRIX_STORE_TENDENCY_IMP(
0072 I bi, bj, myTime, myIter, myThid )
5d0c92356b Jean*0073
2cd7c28080 Jean*0074
8385b3bb5f Ed H*0075
0076
0077
a9544d2f65 Sama*0078
2cd7c28080 Jean*0079
a9544d2f65 Sama*0080 IMPLICIT NONE
8f6dc1243d Jean*0081 #include "SIZE.h"
a9544d2f65 Sama*0082 #include "EEPARAMS.h"
0083 #include "PARAMS.h"
8f6dc1243d Jean*0084 #include "PTRACERS_SIZE.h"
0085 #include "PTRACERS_PARAMS.h"
0086 #include "PTRACERS_FIELDS.h"
0087 #include "MATRIX.h"
0088
2cd7c28080 Jean*0089
a9544d2f65 Sama*0090 INTEGER bi, bj
0091 _RL myTime
0092 INTEGER myIter
0093 INTEGER myThid
0094
0095 #ifdef ALLOW_MATRIX
2cd7c28080 Jean*0096
8f6dc1243d Jean*0097 INTEGER iTracer
2cd7c28080 Jean*0098 INTEGER i,j,k
0099
8f6dc1243d Jean*0100
a9544d2f65 Sama*0101 DO iTracer=1,PTRACERS_numInUse
0102 DO k=1,Nr
2cd7c28080 Jean*0103 DO j=1-OLy, sNy+OLy
0104 DO i=1-OLx, sNx+OLx
5d0c92356b Jean*0105 MATRIX(i,j,k,bi,bj,iTracer,2) =
8f6dc1243d Jean*0106 & MATRIX(i,j,k,bi,bj,iTracer,2)
5d0c92356b Jean*0107 & + pTracer(i,j,k,bi,bj,iTracer)
0108 pTracer(i,j,k,bi,bj,iTracer)=
8f6dc1243d Jean*0109 & PTRACERS_initial(i,j,k,bi,bj,iTracer)
a9544d2f65 Sama*0110 ENDDO
5d0c92356b Jean*0111 ENDDO
a9544d2f65 Sama*0112 ENDDO
0113 ENDDO
5d0c92356b Jean*0114 impMatrixCounter(bi,bj) = impMatrixCounter(bi,bj) + 1
8f6dc1243d Jean*0115
a9544d2f65 Sama*0116 #endif /* ALLOW_MATRIX */
8f6dc1243d Jean*0117 RETURN
a9544d2f65 Sama*0118 END