osse_MkII/filter/distobs.F

        real function distobs(a,b,iobsloc,mobs,currobs,n)
        implicit none
c#include osse_parms.h
        integer, intent(in) ::  currobs, n,mobs
        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,currobs
           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,currobs,n)
2	implicit none
3	c#include osse_parms.h
4	integer, intent(in) :: currobs, n,mobs
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,currobs
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