pkg/smooth/smooth_ensemble3D.F

#include "CPP_OPTIONS.h"

      subroutine smooth_ensemble3D (mythid)


      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "GRID.h"
#include "PARAMS.h"
c#include "tamc.h"
#include "smooth.h"


      _RL fld_in(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)
      _RL fld_mean(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)
      _RL fld_meansq(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)
      _RL fld_norm(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)
      character*( 80) fnamegeneric
      integer i,j,k,bi,bj
      integer itlo,ithi
      integer jtlo,jthi
      integer myThid
      _RL port_rand, port_rand_norm
      integer iii,ii,nbRand

      jtlo = mybylo(mythid)
      jthi = mybyhi(mythid)
      itlo = mybxlo(mythid)
      ithi = mybxhi(mythid)

      
C initialize the random number generator:
      fld_in(1,1,1,1,1)=port_rand(1)
      nbRand=1000

C main loop:
      DO ii=1,nbRand
            WRITE(standardMessageUnit,'(A,I4)')
     & 'smooth_ensemble3D ii: ',ii

C define a perturbation:
      DO bj=jtlo,jthi
       DO bi=itlo,ithi
        DO k=1,Nr
         DO j=1-OLy,sNy+OLy
          DO i=1-OLx,sNx+OLx
           fld_in(i,j,k,bi,bj)=0.
           if (maskC(i,j,k,bi,bj).NE.0) then
           fld_in(i,j,k,bi,bj)=port_rand_norm()
           endif
          ENDDO
         ENDDO
        ENDDO
       ENDDO
      ENDDO
      _EXCH_XYZ_RL ( fld_in, myThid )

C the smoothing itself:
      call smooth_correl3D(fld_in,1,mythid)

      write(fnamegeneric(1:80),'(1a,i3.3)')
     &    'wc01_3Densemble.',ii
      call mdswritefield(fnamegeneric,32,.false.,'RL',
     & nR,fld_in,1,1,mythid)

      ENDDO


      end
1	gforget	1.1	#include "CPP_OPTIONS.h"
2
3			subroutine smooth_ensemble3D (mythid)
4
5
6			IMPLICIT NONE
7			#include "SIZE.h"
8			#include "EEPARAMS.h"
9			#include "GRID.h"
10			#include "PARAMS.h"
11			c#include "tamc.h"
12			#include "smooth.h"
13
14
15			_RL fld_in(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)
16			_RL fld_mean(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)
17			_RL fld_meansq(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)
18			_RL fld_norm(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)
19			character*( 80) fnamegeneric
20			integer i,j,k,bi,bj
21			integer itlo,ithi
22			integer jtlo,jthi
23			integer myThid
24			_RL port_rand, port_rand_norm
25			integer iii,ii,nbRand
26
27			jtlo = mybylo(mythid)
28			jthi = mybyhi(mythid)
29			itlo = mybxlo(mythid)
30			ithi = mybxhi(mythid)
31
32
33
34			C initialize the random number generator:
35			fld_in(1,1,1,1,1)=port_rand(1)
36			nbRand=1000
37
38			C main loop:
39			DO ii=1,nbRand
40			WRITE(standardMessageUnit,'(A,I4)')
41			& 'smooth_ensemble3D ii: ',ii
42
43			C define a perturbation:
44			DO bj=jtlo,jthi
45			DO bi=itlo,ithi
46			DO k=1,Nr
47			DO j=1-OLy,sNy+OLy
48			DO i=1-OLx,sNx+OLx
49			fld_in(i,j,k,bi,bj)=0.
50			if (maskC(i,j,k,bi,bj).NE.0) then
51			fld_in(i,j,k,bi,bj)=port_rand_norm()
52			endif
53			ENDDO
54			ENDDO
55			ENDDO
56			ENDDO
57			ENDDO
58	gforget	1.2	_EXCH_XYZ_RL ( fld_in, myThid )
59	gforget	1.1
60			C the smoothing itself:
61			call smooth_correl3D(fld_in,1,mythid)
62
63			write(fnamegeneric(1:80),'(1a,i3.3)')
64			& 'wc01_3Densemble.',ii
65			call mdswritefield(fnamegeneric,32,.false.,'RL',
66			& nR,fld_in,1,1,mythid)
67
68			ENDDO
69
70
71			end