16 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
17 |
#include "PARAMS.h" |
#include "PARAMS.h" |
18 |
#include "GRID.h" |
#include "GRID.h" |
19 |
|
#include "SHAP_FILT.h" |
20 |
|
|
21 |
C == Routine arguments |
C == Routine arguments |
22 |
_RL uVel(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
_RL uVel(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
40 |
ENDDO |
ENDDO |
41 |
|
|
42 |
C Extract small-scale noise from tmpFldX (delta_ii^n) |
C Extract small-scale noise from tmpFldX (delta_ii^n) |
43 |
DO N=1,nShap |
DO N=1,nShapUV |
44 |
N1=1+mod(N+1,2) |
N1=1+mod(N+1,2) |
45 |
N2=1+mod( N ,2) |
N2=1+mod( N ,2) |
46 |
DO J=1-OLy,sNy+OLy |
DO J=1-OLy,sNy+OLy |
78 |
#endif /* SEQUENTIAL_2D_SHAP */ |
#endif /* SEQUENTIAL_2D_SHAP */ |
79 |
|
|
80 |
C Extract small-scale noise from tmpFldY (delta_jj^n) |
C Extract small-scale noise from tmpFldY (delta_jj^n) |
81 |
DO N=1,nShap |
DO N=1,nShapUV |
82 |
N1=1+mod(N+1,2) |
N1=1+mod(N+1,2) |
83 |
N2=1+mod( N ,2) |
N2=1+mod( N ,2) |
84 |
DO J=1-OLy+1,sNy+OLy-1 |
DO J=1-OLy+1,sNy+OLy-1 |