/[MITgcm]/MITgcm_contrib/gael/pkg/smooth2/smooth_init3D.F
ViewVC logotype

Diff of /MITgcm_contrib/gael/pkg/smooth2/smooth_init3D.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by gforget, Sun Oct 25 21:24:04 2009 UTC revision 1.2 by jmc, Fri Nov 13 19:52:32 2009 UTC
# Line 13  C     *================================= Line 13  C     *=================================
13    
14  cgf the choices of smooth3Dtype and smooth3Dsize need comments...  cgf the choices of smooth3Dtype and smooth3Dsize need comments...
15  cgf  cgf
16  cgf smooth3DtypeH=      1) HORIZONTAL ALONG GRID AXIS  cgf smooth3DtypeH= 1) HORIZONTAL ALONG GRID AXIS
17  cgf                                     2-3) GMREDI TYPES  cgf              2-3) GMREDI TYPES
18  cgf                                     4) HORIZONTAL BUT WITH ROTATED AXIS  cgf                4) HORIZONTAL BUT WITH ROTATED AXIS
19  cgf for now I focus on the simpler smooth3DtypeH=1 case  cgf for now I focus on the simpler smooth3DtypeH=1 case
20    
21    
# Line 92  c avoid excessive vertical smoothing: Line 92  c avoid excessive vertical smoothing:
92        if (smooth3DsizeZ(smoothOpNb).NE.3) then        if (smooth3DsizeZ(smoothOpNb).NE.3) then
93        DO bj=jtlo,jthi        DO bj=jtlo,jthi
94         DO bi=itlo,ithi         DO bi=itlo,ithi
95          DO k=1,Nr          DO k=1,Nr
96           DO j=1-OLy,sNy+OLy           DO j=1-OLy,sNy+OLy
97            DO i=1-OLx,sNx+OLx            DO i=1-OLx,sNx+OLx
98    
99         smooth3D_KzMax=drC(k)         smooth3D_KzMax=drC(k)
100         smooth3D_KzMax=smooth3D_KzMax*smooth3D_KzMax/smooth3DtotTime/2         smooth3D_KzMax=smooth3D_KzMax*smooth3D_KzMax/smooth3DtotTime/2
101         if (smooth3D_kappaR(i,j,k,bi,bj).GT.smooth3D_KzMax) then         if (smooth3D_kappaR(i,j,k,bi,bj).GT.smooth3D_KzMax) then
102         smooth3D_kappaR(i,j,k,bi,bj)=smooth3D_KzMax         smooth3D_kappaR(i,j,k,bi,bj)=smooth3D_KzMax
103         endif         endif
104            ENDDO            ENDDO
105           ENDDO           ENDDO
# Line 200  c write diffusion operator to file Line 200  c write diffusion operator to file
200       & nR,smooth3D_kappaR,10,1,mythid)       & nR,smooth3D_kappaR,10,1,mythid)
201    
202    
203        END        END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22