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