osse/EnKF/distobs.F

        real function distobs(a,b,iobsloc,mobs,n)

        implicit none
        integer, intent(in) ::  mobs, n
        integer,   intent(in)  :: iobsloc(mobs) 
        real*8, intent(in) ::   a(mobs), b(n)
        integer i,j

        distobs=(a(1)-b(iobsloc(1)))**2
        do i=2,mobs
           j = iobsloc(i)
         distobs=distobs+(a(i)-b(j))**2
c         write (*,*) j, a(i), b(j)
        enddo
        distobs=sqrt(distobs)
        
        return
        end


1	afe	1.1	real function distobs(a,b,iobsloc,mobs,n)
2
3			implicit none
4			integer, intent(in) :: mobs, n
5			integer, intent(in) :: iobsloc(mobs)
6			real*8, intent(in) :: a(mobs), b(n)
7			integer i,j
8
9			distobs=(a(1)-b(iobsloc(1)))**2
10			do i=2,mobs
11			j = iobsloc(i)
12			distobs=distobs+(a(i)-b(j))**2
13			c write (,) j, a(i), b(j)
14			enddo
15			distobs=sqrt(distobs)
16
17			return
18			end
19
20