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
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