dyn3d/Read_reanalyse/correctbid.f



!===========================================================================
      subroutine correctbid(iim,nl,x)
!===========================================================================
      integer iim,nl
      real x(iim+1,nl)
      integer i,l
      real zz

      do l=1,nl
         do i=2,iim-1
            if(abs(x(i,l)).gt.1.e10) then
               zz=0.5*(x(i-1,l)+x(i+1,l))
!              print*,'correction ',i,l,x(i,l),zz
               x(i,l)=zz
            endif
         enddo
      enddo
      return
      end
1	guez	37
2
3
4			!===========================================================================
5			subroutine correctbid(iim,nl,x)
6			!===========================================================================
7			integer iim,nl
8			real x(iim+1,nl)
9			integer i,l
10			real zz
11
12			do l=1,nl
13			do i=2,iim-1
14			if(abs(x(i,l)).gt.1.e10) then
15			zz=0.5*(x(i-1,l)+x(i+1,l))
16			! print*,'correction ',i,l,x(i,l),zz
17			x(i,l)=zz
18			endif
19			enddo
20			enddo
21			return
22			end