pkg/smooth/smooth_init_fixed.F

0c3d35c9cd Gael*0001 #include "SMOOTH_OPTIONS.h"
                0002
9f5240b52a Jean*0003       SUBROUTINE SMOOTH_INIT_FIXED( myThid )
0c3d35c9cd Gael*0004
                0005 C     *==========================================================*
                0006 C     | SUBROUTINE smooth_init_fixed
                0007 C     | o Routine that initializes smoothing/correlation operators
                0008 C     *==========================================================*
                0009
                0010       IMPLICIT NONE
                0011 #include "SIZE.h"
                0012 #include "EEPARAMS.h"
                0013 #include "PARAMS.h"
                0014 #include "GRID.h"
                0015 #include "SMOOTH.h"
                0016
9f5240b52a Jean*0017       INTEGER myThid
7eed09f170 Gael*0018
9f5240b52a Jean*0019       INTEGER smoothOpNb
                0020       INTEGER i, j, k, bi, bj
7eed09f170 Gael*0021
9f5240b52a Jean*0022       DO bj=myByLo(myThid),myByHi(myThid)
                0023        DO bi=myBxLo(myThid),myBxHi(myThid)
7eed09f170 Gael*0024         DO k=1,Nr
                0025          DO j=1-OLy,sNy+OLy
                0026           DO i=1-OLx,sNx+OLx
9f5240b52a Jean*0027            smooth_recip_hFacC(i,j,k,bi,bj) = _recip_hFacC(i,j,k,bi,bj)
                0028            smooth_hFacW(i,j,k,bi,bj) = _hFacW(i,j,k,bi,bj)
                0029            smooth_hFacS(i,j,k,bi,bj) = _hFacS(i,j,k,bi,bj)
7eed09f170 Gael*0030           ENDDO
                0031          ENDDO
                0032         ENDDO
                0033        ENDDO
                0034       ENDDO
0c3d35c9cd Gael*0035
                0036       DO smoothOpNb=1,smoothOpNbMax
9f5240b52a Jean*0037        IF (smooth2Dtype(smoothOpNb).NE.0) THEN
                0038         CALL smooth_init2D( smoothOpNb, myThid )
                0039        ENDIF
                0040       ENDDO
0c3d35c9cd Gael*0041
                0042       DO smoothOpNb=1,smoothOpNbMax
9f5240b52a Jean*0043        IF (smooth2Dtype(smoothOpNb).NE.0) THEN
                0044         CALL smooth_filtervar2D( smoothOpNb, myThid )
                0045        ENDIF
0c3d35c9cd Gael*0046       ENDDO
                0047
                0048       DO smoothOpNb=1,smoothOpNbMax
9f5240b52a Jean*0049        IF ( (smooth3DtypeZ(smoothOpNb).NE.0).OR.
                0050      &      (smooth3DtypeH(smoothOpNb).NE.0) ) THEN
                0051         CALL smooth_init3D( smoothOpNb, myThid )
                0052        ENDIF
0c3d35c9cd Gael*0053       ENDDO
                0054
                0055       DO smoothOpNb=1,smoothOpNbMax
9f5240b52a Jean*0056        IF ( (smooth3DtypeZ(smoothOpNb).NE.0).OR.
                0057      &      (smooth3DtypeH(smoothOpNb).NE.0) ) THEN
                0058         CALL smooth_filtervar3D( smoothOpNb, myThid )
                0059        ENDIF
0c3d35c9cd Gael*0060       ENDDO
                0061
9f5240b52a Jean*0062       RETURN
                0063       END