openacc code cannot be parallelized by 19.10

I have to reinstall the PGI community 19.10 instead of 19.4 which licence is out of date. However, when I compiled one of my openacc code which worked well by 19.4, it cannot be parallelized any more.

main:
80, Generating copyin(fsm(:)) [if not already present]
Generating copy(obs(:)) [if not already present]
81, Generating implicit copy(v,x,z,y) [if not already present]
82, Generating Tesla code
82, !$acc loop seq
90, !$acc loop seq
90, Inner sequential loop scheduled on accelerator
116, Reference argument passing prevents parallelization: z
Reference argument passing prevents parallelization: y
Reference argument passing prevents parallelization: x
Reference argument passing prevents parallelization: v
Reference argument passing prevents parallelization: u2
Reference argument passing prevents parallelization: gy
Reference argument passing prevents parallelization: gx
Reference argument passing prevents parallelization: dep
Reference argument passing prevents parallelization: gz
Reference argument passing prevents parallelization: u1
Reference argument passing prevents parallelization: t7

The part of paralleled code is

!$acc routine(indisp) seq

    !$acc data copy(obs) copyin(fsm)
    !$acc kernels loop gang vector private(latf,lonf,lats,lons,t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,u1,u2,dep,gx,gy,gz) reduction(+:sum1,sum2,sum3)
    do i=1,nob
        
        sum1 = 0.d0
        sum2 = 0.d0
        sum3 = 0.d0
        lats = obs(i)%lat*pi/180.d0
        lons = obs(i)%lon*pi/180.d0

        do j=1,nf
            latf = fsm(j)%lat*pi/180.d0
            lonf = fsm(j)%lon*pi/180.d0
            t6 = dsin(lats)*dsin(latf)+dcos(lats)*dcos(latf)*dcos(lons-lonf)
		    if(lons.lt.lonf) then 
			    t6=-dacos(t6)
		    else
			    t6=dacos(t6)
		    endif
		    t3 = t6*re                      
		    t12 = dcos(lats)*dsin(lons-lonf)/dsin(t6)       
		    t13 = (dcos(latf)*dsin(lats)-dsin(latf)*dcos(lats)*dcos(lons-lonf))/dsin(t6)               
		    t1 = t3*t12
		    t2 = t3*t13
            t8 = fsm(j)%strike*pi/180.d0
            t7 = fsm(j)%dip*pi/180.d0
		    t4 =  t1*dsin(t8)+t2*dcos(t8)
		    t5 = -t1*dcos(t8)+t2*dsin(t8)
		    x = t4+0.5d0*fsm(j)%length                   
		    y = t5+0.5d0*fsm(j)%width*dcos(t7)
		    dep = fsm(j)%depth+0.5d0*fsm(j)%width*dsin(t7)  
            
				t9 = fsm(j)%rake*pi/180.d0
				u1 = fsm(j)%slip*dcos(t9)
				u2 = fsm(j)%slip*dsin(t9)
				
				call indisp(u1,u2,0.d0,v,fsm(j)%length,fsm(j)%width,t7,dep,x,y,z,gx,gy,gz)

				t10 = gx*dsin(t8)-gy*dcos(t8)
				t11 = gx*dcos(t8)+gy*dsin(t8)
				sum1 = sum1+t10
				sum2 = sum2+t11
				sum3 = sum3+gz
        enddo
        obs(i)%ux = sum1
        obs(i)%uy = sum2
        obs(i)%uz = sum3
    enddo
    !$acc end data

Did someone encounter the similar problem?

Hi XinZhou,

Unfortunately since the code snip-it isn’t complete, I can only guess as to the issue and can’t determine the difference in behavior between 19.4 and 19.10.

116, Reference argument passing prevents parallelization: z

By default in OpenACC, scalars are private. Though when scalars are passed by reference (the default in Fortran), the address of the scalar escapes analysis and the compiler must presume another reference to the scalar can be made, hence preventing parallelization.

The preferred solution is to use the F2003 “VALUE” attribute when declaring the passed in scalars in the called routine so the scalars are passed by value. However, this does require an interface for the routine, which it doesn’t appear that you have given you’re using “!$acc routine(indisp) seq”.

The second solution is to manually add the scalars to a “private” clause. While you have most of the scalars in a private clause, you are missing v, x, z, and y. Hopefully adding these scalars to the private clause will fix the issue.

If not, please provide a full reproducing example and I’ll be able to look more in depth as to what’s going on.

-Mat

Thank you, Mat. I tried to add the scales x, y, z and v to the ‘private’ clause, the code indeed can be parallelized, but the result is wrong and I also obtained

116, Reference argument passing prevents parallelization: gy
Reference argument passing prevents parallelization: gx
Reference argument passing prevents parallelization: dep
Reference argument passing prevents parallelization: z
Reference argument passing prevents parallelization: y
Reference argument passing prevents parallelization: x
Reference argument passing prevents parallelization: v
Reference argument passing prevents parallelization: gz
Reference argument passing prevents parallelization: u2
Reference argument passing prevents parallelization: u1
Reference argument passing prevents parallelization: t7

The full code is attached here.
Main.f90

module para
		real*8,parameter ::  re=6371.d0, mu=4.41d10, lambda=7.53d10
		real*8 :: v, z
    real*8 :: pi
    integer :: nf, nob
    type station
        real*8 lon, lat, ux, uy, uz
    end type
    type fault
        real*8 lat, lon, depth, slip, rake, strike, dip, length, width
    end type
end module
    
program main
    use para
    use openacc
    implicit none
    integer :: i,j,k,n,m
    type(fault),target,allocatable :: fsm(:)
    type(station),target,allocatable :: obs(:)
    real*8 x,y,dep,u1,u2,u3,gx,gy,gz,sum1,sum2,sum3
    real*8 lons,lats,lonf,latf
    real*8 t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,t14,t15
    real :: start,finish
!$acc routine(indisp) seq

		call CPU_TIME(start)

		pi=4.d0*datan(1.d0)

		v=(lambda+mu)/(lambda+2.d0*mu)
    
    open(unit=1,file="fault.dat",status='old',err=101)
    read(1,*)
    read(1,*) nf
    
    allocate(fsm(nf))
    
    do i=1,nf
        read(1,*) fsm(i)
    enddo
    close(1)
    
    open(unit=2,file="observation point.dat",status='old',err=102)
    read(2,*)
    k=0
30  read(2,*,end=40)
    k=k+1
    goto 30
40  continue
    
    nob=k
    allocate(obs(nob))
    
    rewind(2)
    read(2,*)
    do i=1,nob
        read(2,*) obs(i)%lat, obs(i)%lon
    enddo
    close(2)
    
    
    !$acc data copy(obs) copyin(re,fsm)
    !$acc kernels loop gang vector private(latf,lonf,lats,lons,t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,u1,u2,dep,gx,gy,gz,v,x,y,z) reduction(+:sum1,sum2,sum3)
    do i=1,nob
        
        sum1 = 0.d0
        sum2 = 0.d0
        sum3 = 0.d0
        lats = obs(i)%lat*pi/180.d0
        lons = obs(i)%lon*pi/180.d0

        do j=1,nf
            latf = fsm(j)%lat*pi/180.d0
            lonf = fsm(j)%lon*pi/180.d0
            t6 = dsin(lats)*dsin(latf)+dcos(lats)*dcos(latf)*dcos(lons-lonf)
		    if(lons.lt.lonf) then 
			    t6=-dacos(t6)
		    else
			    t6=dacos(t6)
		    endif
		    t3 = t6*re                            
		    t12 = dcos(lats)*dsin(lons-lonf)/dsin(t6)       
		    t13 = (dcos(latf)*dsin(lats)-dsin(latf)*dcos(lats)*dcos(lons-lonf))/dsin(t6)                    
		    t1 = t3*t12
		    t2 = t3*t13
            t8 = fsm(j)%strike*pi/180.d0
            t7 = fsm(j)%dip*pi/180.d0
		    t4 =  t1*dsin(t8)+t2*dcos(t8)
		    t5 = -t1*dcos(t8)+t2*dsin(t8)
		    x = t4+0.5d0*fsm(j)%length          
		    y = t5+0.5d0*fsm(j)%width*dcos(t7)
		    dep = fsm(j)%depth+0.5d0*fsm(j)%width*dsin(t7)  
            
				t9 = fsm(j)%rake*pi/180.d0
				u1 = fsm(j)%slip*dcos(t9)
				u2 = fsm(j)%slip*dsin(t9)
				
				call indisp(u1,u2,0.d0,v,fsm(j)%length,fsm(j)%width,t7,dep,x,y,z,gx,gy,gz)

				t10 = gx*dsin(t8)-gy*dcos(t8)
				t11 = gx*dcos(t8)+gy*dsin(t8)
				sum1 = sum1+t10
				sum2 = sum2+t11
				sum3 = sum3+gz

        enddo
        obs(i)%ux = sum1
        obs(i)%uy = sum2
        obs(i)%uz = sum3
    enddo
    !$acc end data
        
        open(22,file="displacement_Okada.dat")
        write(22,'(A)') " lon.   lat.   Ue(m),   Un(m),   Uz(m)"
        do k=1,nob
            write(22,55) obs(k)
        enddo
        close(22)
55	format(2f12.7,1x,3E15.5)    
    
    deallocate(fsm)
    deallocate(obs)
    
    call CPU_TIME(finish)
    write(*,*) 'it takes',finish-start,'s'

    goto 103
101	write(*,*) "fault model file is not found !"
	goto 103
102   write(*,*) "observation point file is not found !"
103	continue

end program

The subroutines in fortran77 are

	subroutine indisp(u1,u2,u3,alpha,L,W,dip,c,x,y,z,ux,uy,uz)
!$acc routine seq
	implicit none
	integer*4 i,j
	real*8 u1,u2,u3,alpha,L,W,dip,c,x,y,z,ux,uy,uz,
     .       uas(3),uas_v(3),ubs(3),ucs(3),       
     .       uad(3),uad_v(3),ubd(3),ucd(3),     
     .       uat(3),uat_v(3),ubt(3),uct(3),    
     .       fas(3,4),fas_v(3,4),fbs(3,4),fcs(3,4),   
     .       fad(3,4),fad_v(3,4),fbd(3,4),fcd(3,4), 
     .       fat(3,4),fat_v(3,4),fbt(3,4),fct(3,4),  
     .       d,p,kesai,ita,pi,
     .       temp1(3),temp2(3),temp3(3),
     .       temp4(3),temp5(3),temp6(3),
     .       temp7(3),temp8(3),temp9(3),
     .       uxs,uys,uzs,                           
     .       uxd,uyd,uzd,                
     .       uxt,uyt,uzt          
!$acc routine(defor) seq
	
	d=c-z
	p=y*dcos(dip)+d*dsin(dip)

	kesai=x
	ita=p
	call defor(c,kesai,ita,z,y,dip,alpha,temp1,temp2,temp3,temp4,
     .           temp5,temp6,temp7,temp8,temp9)
	do i=1,3
		fas(i,1)=temp1(i)
		fbs(i,1)=temp2(i)
		fcs(i,1)=temp3(i)
		fad(i,1)=temp4(i)
		fbd(i,1)=temp5(i)
		fcd(i,1)=temp6(i)
		fat(i,1)=temp7(i)
		fbt(i,1)=temp8(i)
		fct(i,1)=temp9(i)
	enddo
	
	kesai=x
	ita=p-W
	call defor(c,kesai,ita,z,y,dip,alpha,temp1,temp2,temp3,temp4,
     .           temp5,temp6,temp7,temp8,temp9)
	do i=1,3
		fas(i,2)=temp1(i)
		fbs(i,2)=temp2(i)
		fcs(i,2)=temp3(i)
		fad(i,2)=temp4(i)
		fbd(i,2)=temp5(i)
		fcd(i,2)=temp6(i)
		fat(i,2)=temp7(i)
		fbt(i,2)=temp8(i)
		fct(i,2)=temp9(i)
	enddo

	kesai=x-L
	ita=p
	call defor(c,kesai,ita,z,y,dip,alpha,temp1,temp2,temp3,temp4,
     .           temp5,temp6,temp7,temp8,temp9)
	do i=1,3
		fas(i,3)=temp1(i)
		fbs(i,3)=temp2(i)
		fcs(i,3)=temp3(i)
		fad(i,3)=temp4(i)
		fbd(i,3)=temp5(i)
		fcd(i,3)=temp6(i)
		fat(i,3)=temp7(i)
		fbt(i,3)=temp8(i)
		fct(i,3)=temp9(i)
	enddo

	kesai=x-L
	ita=p-W
	call defor(c,kesai,ita,z,y,dip,alpha,temp1,temp2,temp3,temp4,
     .           temp5,temp6,temp7,temp8,temp9)
	do i=1,3
		fas(i,4)=temp1(i)
		fbs(i,4)=temp2(i)
		fcs(i,4)=temp3(i)
		fad(i,4)=temp4(i)
		fbd(i,4)=temp5(i)
		fcd(i,4)=temp6(i)
		fat(i,4)=temp7(i)
		fbt(i,4)=temp8(i)
		fct(i,4)=temp9(i)
	enddo

	d=c+z
	p=y*dcos(dip)+d*dsin(dip)

	kesai=x
	ita=p
	call defor(c,kesai,ita,-z,y,dip,alpha,temp1,temp2,temp3,temp4,
     .           temp5,temp6,temp7,temp8,temp9)
	do i=1,3
		fas_v(i,1)=temp1(i)
		fad_v(i,1)=temp4(i)
		fat_v(i,1)=temp7(i)
	enddo
	
	kesai=x
	ita=p-W
	call defor(c,kesai,ita,-z,y,dip,alpha,temp1,temp2,temp3,temp4,
     .           temp5,temp6,temp7,temp8,temp9)
	do i=1,3
		fas_v(i,2)=temp1(i)
		fad_v(i,2)=temp4(i)
		fat_v(i,2)=temp7(i)
	enddo

	kesai=x-L
	ita=p
	call defor(c,kesai,ita,-z,y,dip,alpha,temp1,temp2,temp3,temp4,
     .           temp5,temp6,temp7,temp8,temp9)
	do i=1,3
		fas_v(i,3)=temp1(i)
		fad_v(i,3)=temp4(i)
		fat_v(i,3)=temp7(i)
	enddo

	kesai=x-L
	ita=p-W
	call defor(c,kesai,ita,-z,y,dip,alpha,temp1,temp2,temp3,temp4,
     .           temp5,temp6,temp7,temp8,temp9)
	do i=1,3
		fas_v(i,4)=temp1(i)
		fad_v(i,4)=temp4(i)
		fat_v(i,4)=temp7(i)
	enddo

	do i=1,3
		uas(i)=0.d0
		uas_v(i)=0.d0
		ubs(i)=0.d0
		ucs(i)=0.d0
		uad(i)=0.d0
		uad_v(i)=0.d0
		ubd(i)=0.d0
		ucd(i)=0.d0
		uat(i)=0.d0
		uat_v(i)=0.d0
		ubt(i)=0.d0
		uct(i)=0.d0
	enddo
	do i=1,3
		uas(i)=fas(i,1)-fas(i,2)-fas(i,3)+fas(i,4)
		uas_v(i)=fas_v(i,1)-fas_v(i,2)-fas_v(i,3)+fas_v(i,4)
		ubs(i)=fbs(i,1)-fbs(i,2)-fbs(i,3)+fbs(i,4)
		ucs(i)=fcs(i,1)-fcs(i,2)-fcs(i,3)+fcs(i,4)
		uad(i)=fad(i,1)-fad(i,2)-fad(i,3)+fad(i,4)
		uad_v(i)=fad_v(i,1)-fad_v(i,2)-fad_v(i,3)+fad_v(i,4)
		ubd(i)=fbd(i,1)-fbd(i,2)-fbd(i,3)+fbd(i,4)
		ucd(i)=fcd(i,1)-fcd(i,2)-fcd(i,3)+fcd(i,4)
		uat(i)=fat(i,1)-fat(i,2)-fat(i,3)+fat(i,4)
		uat_v(i)=fat_v(i,1)-fat_v(i,2)-fat_v(i,3)+fat_v(i,4)
		ubt(i)=fbt(i,1)-fbt(i,2)-fbt(i,3)+fbt(i,4)
		uct(i)=fct(i,1)-fct(i,2)-fct(i,3)+fct(i,4)
	enddo
	
	pi=4.d0*datan(1.d0)
	
	uxs=u1*(uas(1)-uas_v(1)+ubs(1)+z*ucs(1))/2.d0/pi
	uys=u1*((uas(2)-uas_v(2)+ubs(2)+z*ucs(2))*dcos(dip)-
     .	(uas(3)-uas_v(3)+ubs(3)+z*ucs(3))*dsin(dip))/2.d0/pi
	uzs=u1*((uas(2)-uas_v(2)+ubs(2)-z*ucs(2))*dsin(dip)+
     .	(uas(3)-uas_v(3)+ubs(3)-z*ucs(3))*dcos(dip))/2.d0/pi
	
	uxd=u2*(uad(1)-uad_v(1)+ubd(1)+z*ucd(1))/2.d0/pi
	uyd=u2*((uad(2)-uad_v(2)+ubd(2)+z*ucd(2))*dcos(dip)-
     .	(uad(3)-uad_v(3)+ubd(3)+z*ucd(3))*dsin(dip))/2.d0/pi
	uzd=u2*((uad(2)-uad_v(2)+ubd(2)-z*ucd(2))*dsin(dip)+
     .	(uad(3)-uad_v(3)+ubd(3)-z*ucd(3))*dcos(dip))/2.d0/pi
	
	uxt=u3*(uat(1)-uat_v(1)+ubt(1)+z*uct(1))/2.d0/pi
	uyt=u3*((uat(2)-uat_v(2)+ubt(2)+z*uct(2))*dcos(dip)-
     .	(uat(3)-uat_v(3)+ubt(3)+z*uct(3))*dsin(dip))/2.d0/pi
	uzt=u3*((uat(2)-uat_v(2)+ubt(2)-z*uct(2))*dsin(dip)+
     .	(uat(3)-uat_v(3)+ubt(3)-z*uct(3))*dcos(dip))/2.d0/pi


	
	ux=uxs+uxd+uxt
	uy=uys+uyd+uyt
	uz=uzs+uzd+uzt

	end subroutine

	subroutine defor(c,kesai,ita,z,y,dip,alpha,fas,fbs,fcs,
     .       fad,fbd,fcd,fat,fbt,fct)
!$acc routine seq
	implicit none
	real*8 c,kesai,ita,z,y,dip,alpha,fas(3),fas_v(3),fbs(3),fcs(3),
     .       fad(3),fad_v(3),fbd(3),fcd(3),
     .       fat(3),fat_v(3),fbt(3),fct(3),
     .       d,p,q,r,yy,dd,cc,theta,
     .       x,xx,I1,I2,I3,I4,y11,y32,z32,x11,x32,h,temp1,temp2,temp3
	real*8 sheld
	parameter(sheld=1.d-25)

	d=c-z
	p=y*dcos(dip)+d*dsin(dip)
	q=y*dsin(dip)-d*dcos(dip)
	r=dsqrt(kesai**2+ita**2+q**2)
	yy=ita*dcos(dip)+q*dsin(dip)
	dd=ita*dsin(dip)-q*dcos(dip)
	cc=dd+z

	if(dabs(q).lt.sheld) then
		theta=0.d0
	else
		theta=datan(kesai*ita/q/r)
	endif

	xx=dsqrt(kesai**2+q**2)
	
	if(dabs(dcos(dip)).lt.sheld) then
		if(dabs(r+ita).lt.sheld) then
	I3=0.5d0*(ita/(r+dd)+yy*q/((r+dd)**2)+dlog(r-ita))
		else
	I3=0.5d0*(ita/(r+dd)+yy*q/((r+dd)**2)-dlog(r+ita))
		endif
		if(abs(kesai).lt.(10.d0**(-16))) then
	I4=0.d0
		else
	I4=0.5d0*kesai*yy/((r+dd)**2)
		endif
	else
		if(dabs(r+ita).lt.sheld) then
	I3=yy/dcos(dip)/(r+dd)-(-dlog(r-ita)-dsin(dip)*dlog(r+dd))
     .	/dcos(dip)/dcos(dip)
		else
	I3=yy/dcos(dip)/(r+dd)-(dlog(r+ita)-dsin(dip)*dlog(r+dd))
     .	/dcos(dip)/dcos(dip)
		endif
		if(abs(kesai).lt.(10.d0**(-16))) then
	I4=0.d0
		else
	I4=dsin(dip)*kesai/dcos(dip)/(r+dd)+2.d0*
     .	datan((ita*(xx+q*dcos(dip))+xx*(r+xx)*dsin(dip))/kesai/(r+xx)
     .    /dcos(dip))/dcos(dip)/dcos(dip)
		endif
	endif
	
	I1=-kesai/(r+dd)*dcos(dip)-I4*dsin(dip)
	I2=dlog(r+dd)+I3*dsin(dip)

	if(dabs(r+kesai).lt.sheld) then
		x11=0.d0
		x32=0.d0
	else
		x11=1.d0/r/(r+kesai)
		x32=(2.d0*r+kesai)/(r**3)/((r+kesai)**2)
	endif

	if(dabs(r+ita).lt.sheld) then
		y11=0.d0
		y32=0.d0
	else
		y11=1.d0/r/(r+ita)
		y32=(2.d0*r+ita)/(r**3)/((r+ita)**2)
	endif
	h=q*dcos(dip)-z
	z32=dsin(dip)/(r**3)-h*y32

	
	fas(1)=(theta+alpha*kesai*q*y11)/2.d0
	fas(2)=alpha*q/2.d0/r
	if(dabs(r+ita).lt.sheld) then
		fas(3)=-(1.d0-alpha)*dlog(r-ita)/2.d0-alpha*q*q*y11/2.d0
	else
		fas(3)=(1.d0-alpha)*dlog(r+ita)/2.d0-alpha*q*q*y11/2.d0
	endif

	fbs(1)=-kesai*q*y11-theta-(1.d0-alpha)/alpha*I1*dsin(dip)
	fbs(2)=-q/r+(1.d0-alpha)/alpha*yy*dsin(dip)/(r+dd)
	fbs(3)=q*q*y11-(1.d0-alpha)/alpha*I2*dsin(dip)

	fcs(1)=(1.d0-alpha)*kesai*y11*dcos(dip)-alpha*kesai*q*z32
	fcs(2)=(1.d0-alpha)*(dcos(dip)/r+2.d0*q*y11*dsin(dip))
     .	-alpha*cc*q/(r**3)
	fcs(3)=(1.d0-alpha)*q*y11*dcos(dip)
     .	-alpha*(cc*ita/(r**3)-z*y11+kesai**2*z32)

	
	fad(1)=alpha*q/2.d0/r
	fad(2)=(theta+alpha*ita*q*x11)/2.d0
	if(dabs(r+kesai).lt.sheld) then
		fad(3)=-(1.d0-alpha)*dlog(r-kesai)/2.d0-alpha*q*q*x11/2.d0
	else
		fad(3)=(1.d0-alpha)*dlog(r+kesai)/2.d0-alpha*q*q*x11/2.d0
	endif

	fbd(1)=-q/r+(1.d0-alpha)/alpha*I3*dsin(dip)*dcos(dip)
	fbd(2)=-ita*q*x11-theta-(1.d0-alpha)/alpha*kesai*dsin(dip)*
     .	dcos(dip)/(r+dd)
	fbd(3)=q*q*x11+(1.d0-alpha)/alpha*I4*dsin(dip)*dcos(dip)

	fcd(1)=(1.d0-alpha)*dcos(dip)/r-q*y11*dsin(dip)-alpha*cc*q/(r**3)
	fcd(2)=(1.d0-alpha)*yy*x11-alpha*cc*ita*q*x32
	fcd(3)=-dd*x11-kesai*y11*dsin(dip)-alpha*cc*(x11-q*q*x32)

	
	if(dabs(r+ita).lt.sheld) then
		fat(1)=(1.d0-alpha)*dlog(r-ita)/2.d0-alpha*q*q*y11/2.d0
	else
		fat(1)=-(1.d0-alpha)*dlog(r+ita)/2.d0-alpha*q*q*y11/2.d0
	endif
	if(dabs(r+kesai).lt.sheld) then
		fat(2)=(1.d0-alpha)*dlog(r-kesai)/2.d0-alpha*q*q*x11/2.d0
	else
		fat(2)=-(1.d0-alpha)*dlog(r+kesai)/2.d0-alpha*q*q*x11/2.d0
	endif
	fat(3)=(theta-alpha*q*(ita*x11+kesai*y11))/2.d0

	fbt(1)=q*q*y11-(1.d0-alpha)/alpha*I3*dsin(dip)*dsin(dip)
	fbt(2)=q*q*x11+(1.d0-alpha)/alpha*kesai*dsin(dip)*
     .	dsin(dip)/(r+dd)
	fbt(3)=q*(ita*x11+kesai*y11)-theta-(1.d0-alpha)/alpha*I4*
     .	dsin(dip)*dsin(dip)

	fct(1)=-(1.d0-alpha)*(dsin(dip)/r+q*Y11*dcos(dip))
     .       -alpha*(z*y11-q*q*z32)
	fct(2)=(1.d0-alpha)*2.d0*kesai*y11*dsin(dip)+dd*x11
     .       -alpha*cc*(x11-q*q*x32)
	fct(3)=(1.d0-alpha)*(yy*x11+kesai*y11*dcos(dip))
     .       +alpha*q*(cc*ita*x32+kesai*z32)

	end subroutine

Good, glad I was correct about the parallelization issue.

Though to look at a verification issue, I’ll need to be able to run the program. Can you please provide the data files as well?

Thanks,
Mat

Hi Mat, there are two input data files.

  1. ‘fault.dat’
lat.     lon.    depth       slip(m)     rake   strike   dip   length   width
400
  38.353  44.401   1.894    0.679E-05   -70.000 348.000  40.000   1.000   1.000
  38.354  44.410   2.536    0.497E-04  -154.074 348.000  40.000   1.000   1.000
  38.356  44.418   3.179    0.329E-03  -160.000 348.000  40.000   1.000   1.000
  38.357  44.427   3.822    0.788E-03  -160.000 348.000  40.000   1.000   1.000
  38.359  44.435   4.465    0.132E-02  -158.726 348.000  40.000   1.000   1.000
  38.360  44.444   5.107    0.182E-02  -157.161 348.000  40.000   1.000   1.000
  38.361  44.453   5.750    0.217E-02  -155.945 348.000  40.000   1.000   1.000
  38.363  44.461   6.393    0.234E-02  -155.354 348.000  40.000   1.000   1.000
  38.364  44.470   7.036    0.231E-02  -155.487 348.000  40.000   1.000   1.000
  38.366  44.478   7.679    0.212E-02  -156.329 348.000  40.000   1.000   1.000
  38.367  44.487   8.321    0.181E-02  -157.718 348.000  40.000   1.000   1.000
  38.369  44.496   8.964    0.144E-02  -159.218 348.000  40.000   1.000   1.000
  38.370  44.504   9.607    0.106E-02  -160.000 348.000  40.000   1.000   1.000
  38.371  44.513  10.250    0.708E-03  -160.000 348.000  40.000   1.000   1.000
  38.373  44.521  10.893    0.419E-03  -160.000 348.000  40.000   1.000   1.000
  38.374  44.530  11.535    0.205E-03  -150.504 348.000  40.000   1.000   1.000
  38.376  44.539  12.178    0.123E-03  -101.228 348.000  40.000   1.000   1.000
  38.377  44.547  12.821    0.195E-03   -70.000 348.000  40.000   1.000   1.000
  38.379  44.556  13.464    0.281E-03   -70.000 348.000  40.000   1.000   1.000
  38.380  44.564  14.106    0.349E-03   -77.240 348.000  40.000   1.000   1.000
  38.362  44.399   1.894    0.203E-05   -70.000 348.000  40.000   1.000   1.000
  38.363  44.407   2.536    0.121E-03  -159.921 348.000  40.000   1.000   1.000
  38.364  44.416   3.179    0.777E-03  -160.000 348.000  40.000   1.000   1.000
  38.366  44.424   3.822    0.184E-02  -159.507 348.000  40.000   1.000   1.000
  38.367  44.433   4.465    0.307E-02  -157.396 348.000  40.000   1.000   1.000
  38.369  44.442   5.107    0.421E-02  -155.178 348.000  40.000   1.000   1.000
  38.370  44.450   5.750    0.506E-02  -153.376 348.000  40.000   1.000   1.000
  38.372  44.459   6.393    0.550E-02  -152.161 348.000  40.000   1.000   1.000
  38.373  44.467   7.036    0.550E-02  -151.564 348.000  40.000   1.000   1.000
  38.375  44.476   7.679    0.511E-02  -151.536 348.000  40.000   1.000   1.000
  38.376  44.485   8.321    0.445E-02  -151.952 348.000  40.000   1.000   1.000
  38.377  44.493   8.964    0.364E-02  -152.572 348.000  40.000   1.000   1.000
  38.379  44.502   9.607    0.280E-02  -153.005 348.000  40.000   1.000   1.000
  38.380  44.510  10.250    0.201E-02  -152.697 348.000  40.000   1.000   1.000
  38.382  44.519  10.893    0.134E-02  -150.363 348.000  40.000   1.000   1.000
  38.383  44.528  11.535    0.837E-03  -141.811 348.000  40.000   1.000   1.000
  38.385  44.536  12.178    0.561E-03  -120.036 348.000  40.000   1.000   1.000
  38.386  44.545  12.821    0.538E-03   -93.840 348.000  40.000   1.000   1.000
  38.387  44.553  13.464    0.632E-03   -82.753 348.000  40.000   1.000   1.000
  38.389  44.562  14.106    0.707E-03   -84.239 348.000  40.000   1.000   1.000
  38.370  44.396   1.894    0.000E+00     0.000 348.000  40.000   1.000   1.000
  38.372  44.405   2.536    0.239E-03  -160.000 348.000  40.000   1.000   1.000
  38.373  44.413   3.179    0.144E-02  -160.000 348.000  40.000   1.000   1.000
  38.375  44.422   3.822    0.335E-02  -157.972 348.000  40.000   1.000   1.000
  38.376  44.431   4.465    0.557E-02  -154.945 348.000  40.000   1.000   1.000
  38.378  44.439   5.107    0.765E-02  -152.097 348.000  40.000   1.000   1.000
  38.379  44.448   5.750    0.924E-02  -149.750 348.000  40.000   1.000   1.000
  38.380  44.456   6.393    0.101E-01  -147.963 348.000  40.000   1.000   1.000
  38.382  44.465   7.036    0.103E-01  -146.704 348.000  40.000   1.000   1.000
  38.383  44.474   7.679    0.969E-02  -145.892 348.000  40.000   1.000   1.000
  38.385  44.482   8.321    0.862E-02  -145.403 348.000  40.000   1.000   1.000
  38.386  44.491   8.964    0.726E-02  -145.043 348.000  40.000   1.000   1.000
  38.388  44.499   9.607    0.580E-02  -144.508 348.000  40.000   1.000   1.000
  38.389  44.508  10.250    0.441E-02  -143.286 348.000  40.000   1.000   1.000
  38.390  44.517  10.893    0.321E-02  -140.443 348.000  40.000   1.000   1.000
  38.392  44.525  11.535    0.229E-02  -134.336 348.000  40.000   1.000   1.000
  38.393  44.534  12.178    0.169E-02  -123.383 348.000  40.000   1.000   1.000
  38.395  44.542  12.821    0.143E-02  -109.715 348.000  40.000   1.000   1.000
  38.396  44.551  13.464    0.140E-02  -100.339 348.000  40.000   1.000   1.000
  38.398  44.560  14.106    0.148E-02   -99.383 348.000  40.000   1.000   1.000
  38.379  44.394   1.894    0.000E+00     0.000 348.000  40.000   1.000   1.000
  38.381  44.402   2.536    0.417E-03  -160.000 348.000  40.000   1.000   1.000
  38.382  44.411   3.179    0.239E-02  -159.524 348.000  40.000   1.000   1.000
  38.384  44.420   3.822    0.548E-02  -155.625 348.000  40.000   1.000   1.000
  38.385  44.428   4.465    0.907E-02  -151.780 348.000  40.000   1.000   1.000
  38.386  44.437   5.107    0.125E-01  -148.445 348.000  40.000   1.000   1.000
  38.388  44.445   5.750    0.152E-01  -145.695 348.000  40.000   1.000   1.000
  38.389  44.454   6.393    0.167E-01  -143.489 348.000  40.000   1.000   1.000
  38.391  44.463   7.036    0.171E-01  -141.753 348.000  40.000   1.000   1.000
  38.392  44.471   7.679    0.164E-01  -140.390 348.000  40.000   1.000   1.000
  38.394  44.480   8.321    0.148E-01  -139.285 348.000  40.000   1.000   1.000
  38.395  44.488   8.964    0.128E-01  -138.281 348.000  40.000   1.000   1.000
  38.396  44.497   9.607    0.105E-01  -137.150 348.000  40.000   1.000   1.000
  38.398  44.506  10.250    0.830E-02  -135.544 348.000  40.000   1.000   1.000
  38.399  44.514  10.893    0.637E-02  -132.934 348.000  40.000   1.000   1.000
  38.401  44.523  11.535    0.483E-02  -128.626 348.000  40.000   1.000   1.000
  38.402  44.531  12.178    0.374E-02  -122.189 348.000  40.000   1.000   1.000
  38.404  44.540  12.821    0.312E-02  -114.580 348.000  40.000   1.000   1.000
  38.405  44.549  13.464    0.288E-02  -108.724 348.000  40.000   1.000   1.000
  38.406  44.557  14.106    0.293E-02  -107.443 348.000  40.000   1.000   1.000
  38.388  44.391   1.894    0.106E-04  -160.000 348.000  40.000   1.000   1.000
  38.389  44.400   2.536    0.657E-03  -160.000 348.000  40.000   1.000   1.000
  38.391  44.409   3.179    0.361E-02  -157.078 348.000  40.000   1.000   1.000
  38.392  44.417   3.822    0.824E-02  -152.391 348.000  40.000   1.000   1.000
  38.394  44.426   4.465    0.137E-01  -148.136 348.000  40.000   1.000   1.000
  38.395  44.434   5.107    0.189E-01  -144.513 348.000  40.000   1.000   1.000
  38.397  44.443   5.750    0.230E-01  -141.495 348.000  40.000   1.000   1.000
  38.398  44.452   6.393    0.256E-01  -139.001 348.000  40.000   1.000   1.000
  38.399  44.460   7.036    0.263E-01  -136.939 348.000  40.000   1.000   1.000
  38.401  44.469   7.679    0.254E-01  -135.216 348.000  40.000   1.000   1.000
  38.402  44.477   8.321    0.233E-01  -133.731 348.000  40.000   1.000   1.000
  38.404  44.486   8.964    0.203E-01  -132.361 348.000  40.000   1.000   1.000
  38.405  44.495   9.607    0.170E-01  -130.947 348.000  40.000   1.000   1.000
  38.407  44.503  10.250    0.138E-01  -129.261 348.000  40.000   1.000   1.000
  38.408  44.512  10.893    0.109E-01  -127.003 348.000  40.000   1.000   1.000
  38.409  44.520  11.535    0.852E-02  -123.838 348.000  40.000   1.000   1.000
  38.411  44.529  12.178    0.677E-02  -119.630 348.000  40.000   1.000   1.000
  38.412  44.538  12.821    0.566E-02  -114.915 348.000  40.000   1.000   1.000
  38.414  44.546  13.464    0.513E-02  -111.194 348.000  40.000   1.000   1.000
  38.415  44.555  14.106    0.510E-02  -110.143 348.000  40.000   1.000   1.000
  38.397  44.389   1.894    0.140E-04  -160.000 348.000  40.000   1.000   1.000
  38.398  44.398   2.536    0.915E-03  -160.000 348.000  40.000   1.000   1.000
  38.400  44.406   3.179    0.502E-02  -154.054 348.000  40.000   1.000   1.000
  38.401  44.415   3.822    0.115E-01  -148.775 348.000  40.000   1.000   1.000
  38.403  44.423   4.465    0.193E-01  -144.244 348.000  40.000   1.000   1.000
  38.404  44.432   5.107    0.268E-01  -140.413 348.000  40.000   1.000   1.000
  38.405  44.441   5.750    0.329E-01  -137.193 348.000  40.000   1.000   1.000
  38.407  44.449   6.393    0.366E-01  -134.483 348.000  40.000   1.000   1.000
  38.408  44.458   7.036    0.379E-01  -132.190 348.000  40.000   1.000   1.000
  38.410  44.466   7.679    0.368E-01  -130.230 348.000  40.000   1.000   1.000
  38.411  44.475   8.321    0.339E-01  -128.517 348.000  40.000   1.000   1.000
  38.413  44.484   8.964    0.298E-01  -126.957 348.000  40.000   1.000   1.000
  38.414  44.492   9.607    0.253E-01  -125.434 348.000  40.000   1.000   1.000
  38.415  44.501  10.250    0.207E-01  -123.800 348.000  40.000   1.000   1.000
  38.417  44.510  10.893    0.166E-01  -121.874 348.000  40.000   1.000   1.000
  38.418  44.518  11.535    0.132E-01  -119.474 348.000  40.000   1.000   1.000
  38.420  44.527  12.178    0.106E-01  -116.550 348.000  40.000   1.000   1.000
  38.421  44.535  12.821    0.891E-02  -113.425 348.000  40.000   1.000   1.000
  38.423  44.544  13.464    0.802E-02  -110.962 348.000  40.000   1.000   1.000
  38.424  44.553  14.106    0.789E-02  -110.184 348.000  40.000   1.000   1.000
  38.406  44.387   1.894    0.000E+00     0.000 348.000  40.000   1.000   1.000
  38.407  44.395   2.536    0.113E-02  -160.000 348.000  40.000   1.000   1.000
  38.408  44.404   3.179    0.643E-02  -150.762 348.000  40.000   1.000   1.000
  38.410  44.413   3.822    0.151E-01  -144.900 348.000  40.000   1.000   1.000
  38.411  44.421   4.465    0.255E-01  -140.107 348.000  40.000   1.000   1.000
  38.413  44.430   5.107    0.358E-01  -136.079 348.000  40.000   1.000   1.000
  38.414  44.438   5.750    0.442E-01  -132.674 348.000  40.000   1.000   1.000
  38.416  44.447   6.393    0.495E-01  -129.787 348.000  40.000   1.000   1.000
  38.417  44.456   7.036    0.514E-01  -127.327 348.000  40.000   1.000   1.000
  38.418  44.464   7.679    0.500E-01  -125.218 348.000  40.000   1.000   1.000
  38.420  44.473   8.321    0.462E-01  -123.384 348.000  40.000   1.000   1.000
  38.421  44.481   8.964    0.408E-01  -121.750 348.000  40.000   1.000   1.000
  38.423  44.490   9.607    0.348E-01  -120.230 348.000  40.000   1.000   1.000
  38.424  44.499  10.250    0.287E-01  -118.722 348.000  40.000   1.000   1.000
  38.426  44.507  10.893    0.232E-01  -117.112 348.000  40.000   1.000   1.000
  38.427  44.516  11.535    0.186E-01  -115.292 348.000  40.000   1.000   1.000
  38.429  44.524  12.178    0.150E-01  -113.235 348.000  40.000   1.000   1.000
  38.430  44.533  12.821    0.126E-01  -111.137 348.000  40.000   1.000   1.000
  38.431  44.542  13.464    0.113E-01  -109.507 348.000  40.000   1.000   1.000
  38.433  44.550  14.106    0.110E-01  -108.969 348.000  40.000   1.000   1.000
  38.414  44.384   1.894    0.000E+00     0.000 348.000  40.000   1.000   1.000
  38.416  44.393   2.536    0.120E-02  -159.329 348.000  40.000   1.000   1.000
  38.417  44.402   3.179    0.747E-02  -146.984 348.000  40.000   1.000   1.000
  38.419  44.410   3.822    0.182E-01  -140.606 348.000  40.000   1.000   1.000
  38.420  44.419   4.465    0.317E-01  -135.549 348.000  40.000   1.000   1.000
  38.422  44.427   5.107    0.452E-01  -131.310 348.000  40.000   1.000   1.000
  38.423  44.436   5.750    0.564E-01  -127.727 348.000  40.000   1.000   1.000
  38.424  44.445   6.393    0.634E-01  -124.701 348.000  40.000   1.000   1.000
  38.426  44.453   7.036    0.659E-01  -122.142 348.000  40.000   1.000   1.000
  38.427  44.462   7.679    0.643E-01  -119.972 348.000  40.000   1.000   1.000
  38.429  44.470   8.321    0.594E-01  -118.118 348.000  40.000   1.000   1.000
  38.430  44.479   8.964    0.525E-01  -116.514 348.000  40.000   1.000   1.000
  38.432  44.488   9.607    0.447E-01  -115.091 348.000  40.000   1.000   1.000
  38.433  44.496  10.250    0.370E-01  -113.778 348.000  40.000   1.000   1.000
  38.434  44.505  10.893    0.300E-01  -112.502 348.000  40.000   1.000   1.000
  38.436  44.513  11.535    0.240E-01  -111.195 348.000  40.000   1.000   1.000
  38.437  44.522  12.178    0.194E-01  -109.837 348.000  40.000   1.000   1.000
  38.439  44.531  12.821    0.163E-01  -108.528 348.000  40.000   1.000   1.000
  38.440  44.539  13.464    0.145E-01  -107.543 348.000  40.000   1.000   1.000
  38.442  44.548  14.106    0.141E-01  -107.227 348.000  40.000   1.000   1.000
  38.423  44.382   1.894    0.000E+00     0.000 348.000  40.000   1.000   1.000
  38.425  44.391   2.536    0.963E-03  -157.211 348.000  40.000   1.000   1.000
  38.426  44.399   3.179    0.761E-02  -142.419 348.000  40.000   1.000   1.000
  38.427  44.408   3.822    0.202E-01  -135.665 348.000  40.000   1.000   1.000
  38.429  44.416   4.465    0.368E-01  -130.299 348.000  40.000   1.000   1.000
  38.430  44.425   5.107    0.539E-01  -125.817 348.000  40.000   1.000   1.000
  38.432  44.434   5.750    0.681E-01  -122.087 348.000  40.000   1.000   1.000
  38.433  44.442   6.393    0.771E-01  -119.000 348.000  40.000   1.000   1.000
  38.435  44.451   7.036    0.803E-01  -116.447 348.000  40.000   1.000   1.000
  38.436  44.459   7.679    0.782E-01  -114.334 348.000  40.000   1.000   1.000
  38.438  44.468   8.321    0.722E-01  -112.579 348.000  40.000   1.000   1.000
  38.439  44.477   8.964    0.637E-01  -111.117 348.000  40.000   1.000   1.000
  38.440  44.485   9.607    0.542E-01  -109.894 348.000  40.000   1.000   1.000
  38.442  44.494  10.250    0.447E-01  -108.860 348.000  40.000   1.000   1.000
  38.443  44.502  10.893    0.361E-01  -107.972 348.000  40.000   1.000   1.000
  38.445  44.511  11.535    0.289E-01  -107.193 348.000  40.000   1.000   1.000
  38.446  44.520  12.178    0.233E-01  -106.502 348.000  40.000   1.000   1.000
  38.448  44.528  12.821    0.194E-01  -105.920 348.000  40.000   1.000   1.000
  38.449  44.537  13.464    0.173E-01  -105.527 348.000  40.000   1.000   1.000
  38.450  44.545  14.106    0.167E-01  -105.436 348.000  40.000   1.000   1.000
  38.432  44.380   1.894    0.000E+00     0.000 348.000  40.000   1.000   1.000
  38.433  44.388   2.536    0.356E-03  -157.986 348.000  40.000   1.000   1.000
  38.435  44.397   3.179    0.648E-02  -137.694 348.000  40.000   1.000   1.000
  38.436  44.405   3.822    0.203E-01  -130.025 348.000  40.000   1.000   1.000
  38.438  44.414   4.465    0.399E-01  -124.035 348.000  40.000   1.000   1.000
  38.439  44.423   5.107    0.607E-01  -119.265 348.000  40.000   1.000   1.000
  38.441  44.431   5.750    0.780E-01  -115.490 348.000  40.000   1.000   1.000
  38.442  44.440   6.393    0.890E-01  -112.501 348.000  40.000   1.000   1.000
  38.443  44.448   7.036    0.929E-01  -110.125 348.000  40.000   1.000   1.000
  38.445  44.457   7.679    0.903E-01  -108.229 348.000  40.000   1.000   1.000
  38.446  44.466   8.321    0.831E-01  -106.720 348.000  40.000   1.000   1.000
  38.448  44.474   8.964    0.731E-01  -105.533 348.000  40.000   1.000   1.000
  38.449  44.483   9.607    0.619E-01  -104.627 348.000  40.000   1.000   1.000
  38.451  44.491  10.250    0.509E-01  -103.976 348.000  40.000   1.000   1.000
  38.452  44.500  10.893    0.409E-01  -103.569 348.000  40.000   1.000   1.000
  38.453  44.509  11.535    0.326E-01  -103.397 348.000  40.000   1.000   1.000
  38.455  44.517  12.178    0.261E-01  -103.440 348.000  40.000   1.000   1.000
  38.456  44.526  12.821    0.216E-01  -103.641 348.000  40.000   1.000   1.000
  38.458  44.534  13.464    0.191E-01  -103.888 348.000  40.000   1.000   1.000
  38.459  44.543  14.106    0.185E-01  -104.044 348.000  40.000   1.000   1.000
  38.441  44.377   1.894    0.000E+00     0.000 348.000  40.000   1.000   1.000
  38.442  44.386   2.536    0.000E+00     0.000 348.000  40.000   1.000   1.000
  38.444  44.394   3.179    0.456E-02  -136.022 348.000  40.000   1.000   1.000
  38.445  44.403   3.822    0.187E-01  -123.548 348.000  40.000   1.000   1.000
  38.447  44.412   4.465    0.408E-01  -116.242 348.000  40.000   1.000   1.000
  38.448  44.420   5.107    0.650E-01  -111.300 348.000  40.000   1.000   1.000
  38.449  44.429   5.750    0.853E-01  -107.765 348.000  40.000   1.000   1.000
  38.451  44.437   6.393    0.979E-01  -105.157 348.000  40.000   1.000   1.000
  38.452  44.446   7.036    0.102E+00  -103.197 348.000  40.000   1.000   1.000
  38.454  44.455   7.679    0.991E-01  -101.718 348.000  40.000   1.000   1.000
  38.455  44.463   8.321    0.908E-01  -100.621 348.000  40.000   1.000   1.000
  38.457  44.472   8.964    0.794E-01   -99.854 348.000  40.000   1.000   1.000
  38.458  44.480   9.607    0.668E-01   -99.396 348.000  40.000   1.000   1.000
  38.459  44.489  10.250    0.545E-01   -99.257 348.000  40.000   1.000   1.000
  38.461  44.498  10.893    0.435E-01   -99.464 348.000  40.000   1.000   1.000
  38.462  44.506  11.535    0.344E-01  -100.046 348.000  40.000   1.000   1.000
  38.464  44.515  12.178    0.273E-01  -100.991 348.000  40.000   1.000   1.000
  38.465  44.523  12.821    0.224E-01  -102.155 348.000  40.000   1.000   1.000
  38.467  44.532  13.464    0.197E-01  -103.185 348.000  40.000   1.000   1.000
  38.468  44.541  14.106    0.190E-01  -103.635 348.000  40.000   1.000   1.000
  38.450  44.375   1.894    0.159E-03   -70.000 348.000  40.000   1.000   1.000
  38.451  44.383   2.536    0.000E+00     0.000 348.000  40.000   1.000   1.000
  38.452  44.392   3.179    0.256E-02  -142.985 348.000  40.000   1.000   1.000
  38.454  44.401   3.822    0.161E-01  -114.212 348.000  40.000   1.000   1.000
  38.455  44.409   4.465    0.401E-01  -106.017 348.000  40.000   1.000   1.000
  38.457  44.418   5.107    0.669E-01  -101.721 348.000  40.000   1.000   1.000
  38.458  44.426   5.750    0.892E-01   -99.013 348.000  40.000   1.000   1.000
  38.460  44.435   6.393    0.103E+00   -97.179 348.000  40.000   1.000   1.000
  38.461  44.444   7.036    0.107E+00   -95.908 348.000  40.000   1.000   1.000
  38.462  44.452   7.679    0.103E+00   -95.046 348.000  40.000   1.000   1.000
  38.464  44.461   8.321    0.941E-01   -94.523 348.000  40.000   1.000   1.000
  38.465  44.469   8.964    0.816E-01   -94.318 348.000  40.000   1.000   1.000
  38.467  44.478   9.607    0.681E-01   -94.452 348.000  40.000   1.000   1.000
  38.468  44.487  10.250    0.551E-01   -94.982 348.000  40.000   1.000   1.000
  38.470  44.495  10.893    0.435E-01   -95.995 348.000  40.000   1.000   1.000
  38.471  44.504  11.535    0.340E-01   -97.581 348.000  40.000   1.000   1.000
  38.473  44.512  12.178    0.267E-01   -99.741 348.000  40.000   1.000   1.000
  38.474  44.521  12.821    0.218E-01  -102.210 348.000  40.000   1.000   1.000
  38.475  44.530  13.464    0.191E-01  -104.293 348.000  40.000   1.000   1.000
  38.477  44.538  14.106    0.184E-01  -105.105 348.000  40.000   1.000   1.000
  38.458  44.372   1.894    0.224E-03   -70.000 348.000  40.000   1.000   1.000
  38.460  44.381   2.536    0.000E+00     0.000 348.000  40.000   1.000   1.000
  38.461  44.390   3.179    0.820E-03  -160.000 348.000  40.000   1.000   1.000
  38.463  44.398   3.822    0.137E-01   -97.958 348.000  40.000   1.000   1.000
  38.464  44.407   4.465    0.388E-01   -93.207 348.000  40.000   1.000   1.000
  38.466  44.415   5.107    0.668E-01   -91.109 348.000  40.000   1.000   1.000
  38.467  44.424   5.750    0.896E-01   -89.894 348.000  40.000   1.000   1.000
  38.468  44.433   6.393    0.103E+00   -89.167 348.000  40.000   1.000   1.000
  38.470  44.441   7.036    0.107E+00   -88.779 348.000  40.000   1.000   1.000
  38.471  44.450   7.679    0.102E+00   -88.673 348.000  40.000   1.000   1.000
  38.473  44.458   8.321    0.922E-01   -88.845 348.000  40.000   1.000   1.000
  38.474  44.467   8.964    0.793E-01   -89.333 348.000  40.000   1.000   1.000
  38.476  44.476   9.607    0.655E-01   -90.220 348.000  40.000   1.000   1.000
  38.477  44.484  10.250    0.523E-01   -91.637 348.000  40.000   1.000   1.000
  38.478  44.493  10.893    0.408E-01   -93.763 348.000  40.000   1.000   1.000
  38.480  44.501  11.535    0.315E-01   -96.780 348.000  40.000   1.000   1.000
  38.481  44.510  12.178    0.245E-01  -100.709 348.000  40.000   1.000   1.000
  38.483  44.519  12.821    0.199E-01  -105.071 348.000  40.000   1.000   1.000
  38.484  44.527  13.464    0.174E-01  -108.620 348.000  40.000   1.000   1.000
  38.486  44.536  14.106    0.169E-01  -109.846 348.000  40.000   1.000   1.000
  38.467  44.370   1.894    0.243E-03  -139.736 348.000  40.000   1.000   1.000
  38.469  44.379   2.536    0.000E+00     0.000 348.000  40.000   1.000   1.000
  38.470  44.387   3.179    0.000E+00     0.000 348.000  40.000   1.000   1.000
  38.471  44.396   3.822    0.129E-01   -79.787 348.000  40.000   1.000   1.000
  38.473  44.404   4.465    0.376E-01   -80.921 348.000  40.000   1.000   1.000
  38.474  44.413   5.107    0.647E-01   -81.388 348.000  40.000   1.000   1.000
  38.476  44.422   5.750    0.861E-01   -81.715 348.000  40.000   1.000   1.000
  38.477  44.430   6.393    0.981E-01   -82.084 348.000  40.000   1.000   1.000
  38.479  44.439   7.036    0.101E+00   -82.577 348.000  40.000   1.000   1.000
  38.480  44.447   7.679    0.955E-01   -83.256 348.000  40.000   1.000   1.000
  38.482  44.456   8.321    0.853E-01   -84.197 348.000  40.000   1.000   1.000
  38.483  44.465   8.964    0.724E-01   -85.515 348.000  40.000   1.000   1.000
  38.484  44.473   9.607    0.590E-01   -87.380 348.000  40.000   1.000   1.000
  38.486  44.482  10.250    0.465E-01   -90.040 348.000  40.000   1.000   1.000
  38.487  44.490  10.893    0.358E-01   -93.819 348.000  40.000   1.000   1.000
  38.489  44.499  11.535    0.273E-01   -99.025 348.000  40.000   1.000   1.000
  38.490  44.508  12.178    0.213E-01  -105.624 348.000  40.000   1.000   1.000
  38.492  44.516  12.821    0.175E-01  -112.643 348.000  40.000   1.000   1.000
  38.493  44.525  13.464    0.156E-01  -117.971 348.000  40.000   1.000   1.000
  38.494  44.534  14.106    0.153E-01  -119.505 348.000  40.000   1.000   1.000
  38.476  44.368   1.894    0.270E-02  -160.000 348.000  40.000   1.000   1.000
  38.477  44.376   2.536    0.207E-02  -160.000 348.000  40.000   1.000   1.000
  38.479  44.385   3.179    0.836E-03  -160.000 348.000  40.000   1.000   1.000
  38.480  44.393   3.822    0.127E-01   -73.026 348.000  40.000   1.000   1.000
  38.482  44.402   4.465    0.358E-01   -73.705 348.000  40.000   1.000   1.000
  38.483  44.411   5.107    0.604E-01   -74.854 348.000  40.000   1.000   1.000
  38.485  44.419   5.750    0.791E-01   -75.910 348.000  40.000   1.000   1.000
  38.486  44.428   6.393    0.889E-01   -76.977 348.000  40.000   1.000   1.000
  38.487  44.436   7.036    0.900E-01   -78.166 348.000  40.000   1.000   1.000
  38.489  44.445   7.679    0.844E-01   -79.588 348.000  40.000   1.000   1.000
  38.490  44.454   8.321    0.743E-01   -81.383 348.000  40.000   1.000   1.000
  38.492  44.462   8.964    0.622E-01   -83.754 348.000  40.000   1.000   1.000
  38.493  44.471   9.607    0.499E-01   -87.006 348.000  40.000   1.000   1.000
  38.495  44.480  10.250    0.388E-01   -91.570 348.000  40.000   1.000   1.000
  38.496  44.488  10.893    0.295E-01   -97.970 348.000  40.000   1.000   1.000
  38.497  44.497  11.535    0.227E-01  -106.531 348.000  40.000   1.000   1.000
  38.499  44.505  12.178    0.182E-01  -116.661 348.000  40.000   1.000   1.000
  38.500  44.514  12.821    0.158E-01  -126.234 348.000  40.000   1.000   1.000
  38.502  44.523  13.464    0.147E-01  -132.510 348.000  40.000   1.000   1.000
  38.503  44.531  14.106    0.146E-01  -133.843 348.000  40.000   1.000   1.000
  38.485  44.365   1.894    0.756E-02  -160.000 348.000  40.000   1.000   1.000
  38.486  44.374   2.536    0.718E-02  -160.000 348.000  40.000   1.000   1.000
  38.488  44.382   3.179    0.505E-02  -151.181 348.000  40.000   1.000   1.000
  38.489  44.391   3.822    0.134E-01   -81.738 348.000  40.000   1.000   1.000
  38.491  44.400   4.465    0.338E-01   -72.926 348.000  40.000   1.000   1.000
  38.492  44.408   5.107    0.546E-01   -72.351 348.000  40.000   1.000   1.000
  38.493  44.417   5.750    0.697E-01   -73.196 348.000  40.000   1.000   1.000
  38.495  44.425   6.393    0.769E-01   -74.552 348.000  40.000   1.000   1.000
  38.496  44.434   7.036    0.767E-01   -76.298 348.000  40.000   1.000   1.000
  38.498  44.443   7.679    0.707E-01   -78.516 348.000  40.000   1.000   1.000
  38.499  44.451   8.321    0.613E-01   -81.407 348.000  40.000   1.000   1.000
  38.501  44.460   8.964    0.505E-01   -85.313 348.000  40.000   1.000   1.000
  38.502  44.469   9.607    0.399E-01   -90.749 348.000  40.000   1.000   1.000
  38.503  44.477  10.250    0.308E-01   -98.378 348.000  40.000   1.000   1.000
  38.505  44.486  10.893    0.239E-01  -108.687 348.000  40.000   1.000   1.000
  38.506  44.494  11.535    0.194E-01  -121.096 348.000  40.000   1.000   1.000
  38.508  44.503  12.178    0.170E-01  -133.262 348.000  40.000   1.000   1.000
  38.509  44.512  12.821    0.160E-01  -142.428 348.000  40.000   1.000   1.000
  38.511  44.520  13.464    0.157E-01  -147.325 348.000  40.000   1.000   1.000
  38.512  44.529  14.106    0.157E-01  -147.927 348.000  40.000   1.000   1.000
  38.494  44.363   1.894    0.105E-01  -160.000 348.000  40.000   1.000   1.000
  38.495  44.371   2.536    0.102E-01  -160.000 348.000  40.000   1.000   1.000
  38.496  44.380   3.179    0.781E-02  -139.545 348.000  40.000   1.000   1.000
  38.498  44.389   3.822    0.152E-01   -85.548 348.000  40.000   1.000   1.000
  38.499  44.397   4.465    0.322E-01   -73.729 348.000  40.000   1.000   1.000
  38.501  44.406   5.107    0.487E-01   -72.214 348.000  40.000   1.000   1.000
  38.502  44.414   5.750    0.600E-01   -73.055 348.000  40.000   1.000   1.000
  38.504  44.423   6.393    0.646E-01   -74.856 348.000  40.000   1.000   1.000
  38.505  44.432   7.036    0.631E-01   -77.384 348.000  40.000   1.000   1.000
  38.506  44.440   7.679    0.572E-01   -80.761 348.000  40.000   1.000   1.000
  38.508  44.449   8.321    0.488E-01   -85.328 348.000  40.000   1.000   1.000
  38.509  44.458   8.964    0.397E-01   -91.648 348.000  40.000   1.000   1.000
  38.511  44.466   9.607    0.314E-01  -100.444 348.000  40.000   1.000   1.000
  38.512  44.475  10.250    0.250E-01  -112.166 348.000  40.000   1.000   1.000
  38.514  44.483  10.893    0.209E-01  -125.911 348.000  40.000   1.000   1.000
  38.515  44.492  11.535    0.189E-01  -138.910 348.000  40.000   1.000   1.000
  38.516  44.501  12.178    0.182E-01  -148.527 348.000  40.000   1.000   1.000
  38.518  44.509  12.821    0.180E-01  -154.179 348.000  40.000   1.000   1.000
  38.519  44.518  13.464    0.180E-01  -156.649 348.000  40.000   1.000   1.000
  38.521  44.526  14.106    0.180E-01  -156.649 348.000  40.000   1.000   1.000
  38.502  44.360   1.894    0.955E-02  -160.000 348.000  40.000   1.000   1.000
  38.504  44.369   2.536    0.882E-02  -160.000 348.000  40.000   1.000   1.000
  38.505  44.378   3.179    0.722E-02  -122.078 348.000  40.000   1.000   1.000
  38.507  44.386   3.822    0.160E-01   -79.532 348.000  40.000   1.000   1.000
  38.508  44.395   4.465    0.303E-01   -72.343 348.000  40.000   1.000   1.000
  38.510  44.403   5.107    0.432E-01   -72.362 348.000  40.000   1.000   1.000
  38.511  44.412   5.750    0.515E-01   -74.304 348.000  40.000   1.000   1.000
  38.512  44.421   6.393    0.542E-01   -77.230 348.000  40.000   1.000   1.000
  38.514  44.429   7.036    0.520E-01   -81.104 348.000  40.000   1.000   1.000
  38.515  44.438   7.679    0.465E-01   -86.234 348.000  40.000   1.000   1.000
  38.517  44.447   8.321    0.394E-01   -93.175 348.000  40.000   1.000   1.000
  38.518  44.455   8.964    0.324E-01  -102.623 348.000  40.000   1.000   1.000
  38.520  44.464   9.607    0.266E-01  -114.941 348.000  40.000   1.000   1.000
  38.521  44.472  10.250    0.230E-01  -129.076 348.000  40.000   1.000   1.000
  38.522  44.481  10.893    0.212E-01  -142.191 348.000  40.000   1.000   1.000
  38.524  44.490  11.535    0.207E-01  -151.697 348.000  40.000   1.000   1.000
  38.525  44.498  12.178    0.206E-01  -157.062 348.000  40.000   1.000   1.000
  38.527  44.507  12.821    0.206E-01  -159.312 348.000  40.000   1.000   1.000
  38.528  44.515  13.464    0.205E-01  -160.000 348.000  40.000   1.000   1.000
  38.530  44.524  14.106    0.204E-01  -159.835 348.000  40.000   1.000   1.000
  38.511  44.358   1.894    0.680E-02  -160.000 348.000  40.000   1.000   1.000
  38.513  44.367   2.536    0.573E-02  -159.257 348.000  40.000   1.000   1.000
  38.514  44.375   3.179    0.594E-02   -98.021 348.000  40.000   1.000   1.000
  38.515  44.384   3.822    0.158E-01   -71.454 348.000  40.000   1.000   1.000
  38.517  44.392   4.465    0.283E-01   -70.000 348.000  40.000   1.000   1.000
  38.518  44.401   5.107    0.389E-01   -72.363 348.000  40.000   1.000   1.000
  38.520  44.410   5.750    0.454E-01   -75.891 348.000  40.000   1.000   1.000
  38.521  44.418   6.393    0.472E-01   -80.230 348.000  40.000   1.000   1.000
  38.523  44.427   7.036    0.450E-01   -85.596 348.000  40.000   1.000   1.000
  38.524  44.436   7.679    0.403E-01   -92.472 348.000  40.000   1.000   1.000
  38.525  44.444   8.321    0.345E-01  -101.492 348.000  40.000   1.000   1.000
  38.527  44.453   8.964    0.292E-01  -113.103 348.000  40.000   1.000   1.000
  38.528  44.461   9.607    0.254E-01  -126.746 348.000  40.000   1.000   1.000
  38.530  44.470  10.250    0.234E-01  -140.200 348.000  40.000   1.000   1.000
  38.531  44.479  10.893    0.227E-01  -150.770 348.000  40.000   1.000   1.000
  38.533  44.487  11.535    0.226E-01  -157.203 348.000  40.000   1.000   1.000
  38.534  44.496  12.178    0.225E-01  -159.848 348.000  40.000   1.000   1.000
  38.536  44.504  12.821    0.224E-01  -160.000 348.000  40.000   1.000   1.000
  38.537  44.513  13.464    0.223E-01  -160.000 348.000  40.000   1.000   1.000
  38.538  44.522  14.106    0.221E-01  -160.000 348.000  40.000   1.000   1.000
  38.520  44.356   1.894    0.545E-02  -160.000 348.000  40.000   1.000   1.000
  38.521  44.364   2.536    0.447E-02  -150.095 348.000  40.000   1.000   1.000
  38.523  44.373   3.179    0.651E-02   -87.259 348.000  40.000   1.000   1.000
  38.524  44.381   3.822    0.160E-01   -70.208 348.000  40.000   1.000   1.000
  38.526  44.390   4.465    0.272E-01   -70.056 348.000  40.000   1.000   1.000
  38.527  44.399   5.107    0.368E-01   -72.957 348.000  40.000   1.000   1.000
  38.529  44.407   5.750    0.427E-01   -76.868 348.000  40.000   1.000   1.000
  38.530  44.416   6.393    0.444E-01   -81.555 348.000  40.000   1.000   1.000
  38.531  44.425   7.036    0.426E-01   -87.267 348.000  40.000   1.000   1.000
  38.533  44.433   7.679    0.385E-01   -94.484 348.000  40.000   1.000   1.000
  38.534  44.442   8.321    0.334E-01  -103.775 348.000  40.000   1.000   1.000
  38.536  44.450   8.964    0.288E-01  -115.425 348.000  40.000   1.000   1.000
  38.537  44.459   9.607    0.255E-01  -128.680 348.000  40.000   1.000   1.000
  38.539  44.468  10.250    0.238E-01  -141.380 348.000  40.000   1.000   1.000
  38.540  44.476  10.893    0.232E-01  -151.194 348.000  40.000   1.000   1.000
  38.541  44.485  11.535    0.231E-01  -157.151 348.000  40.000   1.000   1.000
  38.543  44.493  12.178    0.230E-01  -159.667 348.000  40.000   1.000   1.000
  38.544  44.502  12.821    0.229E-01  -160.000 348.000  40.000   1.000   1.000
  38.546  44.511  13.464    0.227E-01  -160.000 348.000  40.000   1.000   1.000
  38.547  44.519  14.106    0.226E-01  -159.998 348.000  40.000   1.000   1.000
  1. ‘observation point.dat’
 LAT        LON      
 38.561937 44.256844 
 38.580687 44.340892 
 38.580689 44.397142 
 38.543187 44.340892 
 38.552865 44.382917 
 38.552865 44.411370 
 38.534112 44.382916 
 38.534112 44.411366 
 38.486938 44.256844 
 38.505688 44.340892 
 38.515365 44.382917 
 38.519905 44.404257 
 38.519905 44.418483 
 38.510826 44.404257 
 38.510826 44.418484 
 38.496612 44.382916 
 38.496612 44.411366 
 38.468188 44.340892 
 38.477865 44.382917 
 38.477865 44.411371 
 38.459112 44.382916 
 38.459112 44.411366 
 38.580687 44.453391 
 38.590365 44.495417 
 38.590365 44.523870 
 38.571612 44.495416 
 38.571612 44.523865 
 38.552863 44.439166 
 38.552865 44.466964 
 38.539251 44.432054 
 38.539247 44.445624 
 38.529572 44.432053 
 38.529574 44.445634 
 38.534112 44.466966 
 38.543189 44.509642 
 38.590363 44.551665 
 38.590365 44.579463 
 38.571609 44.551664 
 38.576747 44.572350 
 38.576747 44.586576 
 38.567073 44.572356 
 38.569643 44.583347 
 38.569646 44.590453 
 38.564508 44.583347 
 38.564497 44.590453 
 38.594904 44.600803 
 38.594904 44.615029 
 38.585826 44.600803 
 38.585826 44.615030 
 38.594904 44.629255 
 38.594905 44.643481 
 38.585826 44.629256 
 38.585826 44.643482 
 38.576747 44.600803 
 38.576747 44.615029 
 38.569646 44.596931 
 38.569643 44.604033 
 38.564497 44.596932 
 38.564508 44.604034 
 38.569646 44.611151 
 38.569643 44.618259 
 38.564497 44.611152 
 38.564508 44.618259 
 38.576747 44.629256 
 38.576747 44.643482 
 38.567073 44.629248 
 38.567073 44.643471 
 38.557398 44.544553 
 38.557405 44.558124 
 38.548324 44.544553 
 38.548326 44.558123 
 38.552865 44.579463 
 38.534110 44.551664 
 38.534112 44.579465 
 38.552865 44.607916 
 38.552865 44.636369 
 38.534112 44.607915 
 38.534112 44.636364 
 38.505688 44.453391 
 38.515365 44.495417 
 38.515365 44.523870 
 38.496612 44.495416 
 38.496612 44.523865 
 38.477863 44.439167 
 38.482405 44.459851 
 38.482405 44.474077 
 38.475888 44.456621 
 38.475891 44.463736 
 38.471355 44.456622 
 38.471363 44.463733 
 38.473326 44.474077 
 38.464251 44.432054 
 38.466220 44.442397 
 38.467502 44.447894 
 38.467502 44.451130 
 38.464932 44.447893 
 38.465859 44.450271 
 38.465859 44.451848 
 38.465026 44.449962 
 38.465026 44.450796 
 38.464193 44.449962 
 38.464193 44.450796 
 38.464586 44.451973 
 38.461683 44.442388 
 38.462958 44.447898 
 38.463359 44.450283 
 38.463359 44.451873 
 38.462526 44.449962 
 38.462526 44.450796 
 38.462526 44.451915 
 38.460997 44.447893 
 38.460996 44.451136 
 38.454572 44.432053 
 38.457143 44.442397 
 38.458431 44.447893 
 38.458431 44.451132 
 38.455860 44.447894 
 38.456693 44.449962 
 38.456693 44.450796 
 38.455860 44.449962 
 38.455860 44.450796 
 38.456693 44.451629 
 38.456693 44.452462 
 38.455860 44.451629 
 38.455860 44.452462 
 38.455026 44.449962 
 38.455026 44.450796 
 38.455026 44.451925 
 38.453301 44.440780 
 38.453289 44.444654 
 38.450747 44.440781 
 38.451693 44.443953 
 38.451693 44.445536 
 38.450360 44.443803 
 38.450860 44.444962 
 38.450860 44.445796 
 38.450026 44.444962 
 38.450026 44.445796 
 38.453288 44.447893 
 38.454193 44.450278 
 38.454193 44.451863 
 38.453360 44.449962 
 38.453360 44.450796 
 38.452526 44.449962 
 38.452526 44.450796 
 38.452951 44.451973 
 38.451693 44.447115 
 38.451693 44.448695 
 38.450860 44.446629 
 38.450860 44.447462 
 38.450026 44.446629 
 38.450026 44.447462 
 38.450860 44.448296 
 38.450860 44.449129 
 38.450026 44.448296 
 38.450026 44.449129 
 38.451693 44.450274 
 38.451693 44.451855 
 38.450860 44.449962 
 38.450860 44.450796 
 38.450026 44.449962 
 38.450026 44.450796 
 38.450860 44.451629 
 38.450860 44.452462 
 38.450026 44.451629 
 38.450026 44.452462 
 38.468359 44.453522 
 38.468359 44.455823 
 38.467526 44.453296 
 38.467526 44.454129 
 38.466693 44.453296 
 38.466693 44.454129 
 38.467526 44.454962 
 38.467526 44.456219 
 38.466693 44.454962 
 38.466693 44.456223 
 38.467502 44.458876 
 38.465859 44.453429 
 38.465859 44.455828 
 38.464587 44.453614 
 38.464592 44.455814 
 38.464934 44.458869 
 38.466213 44.463733 
 38.463359 44.453296 
 38.463359 44.454129 
 38.463359 44.455825 
 38.462526 44.453296 
 38.462526 44.454129 
 38.462526 44.455822 
 38.462958 44.458920 
 38.461693 44.453296 
 38.461693 44.454129 
 38.460859 44.453296 
 38.460859 44.454129 
 38.461338 44.455814 
 38.460026 44.453441 
 38.460026 44.455827 
 38.460994 44.458863 
 38.461682 44.463736 
 38.466220 44.470847 
 38.466213 44.477953 
 38.461683 44.470854 
 38.461682 44.477964 
 38.459193 44.453296 
 38.459193 44.454129 
 38.458359 44.453296 
 38.458359 44.454129 
 38.459193 44.454962 
 38.459193 44.456222 
 38.458359 44.454962 
 38.458359 44.456221 
 38.457526 44.453296 
 38.457526 44.454129 
 38.457526 44.454962 
 38.457526 44.456220 
 38.458429 44.458871 
 38.456180 44.453611 
 38.456693 44.454962 
 38.456693 44.456223 
 38.455860 44.454962 
 38.455860 44.456218 
 38.455026 44.453550 
 38.455026 44.455820 
 38.455860 44.458877 
 38.457147 44.463733 
 38.454193 44.453456 
 38.454193 44.455826 
 38.453360 44.453296 
 38.453360 44.454129 
 38.452526 44.453296 
 38.452526 44.454129 
 38.452949 44.455814 
 38.453290 44.458871 
 38.451693 44.453441 
 38.451693 44.454962 
 38.451693 44.455795 
 38.451693 44.456629 
 38.450361 44.453612 
 38.450381 44.455814 
 38.450726 44.458863 
 38.451998 44.463733 
 38.457143 44.470847 
 38.457147 44.477953 
 38.452009 44.470847 
 38.453289 44.476333 
 38.453288 44.479573 
 38.450724 44.476327 
 38.451693 44.478807 
 38.451693 44.480386 
 38.450357 44.478750 
 38.450357 44.480381 
 38.477865 44.495417 
 38.477865 44.523870 
 38.466214 44.484432 
 38.467506 44.489915 
 38.467502 44.493796 
 38.465859 44.489098 
 38.465859 44.491498 
 38.464592 44.489111 
 38.465026 44.490795 
 38.465026 44.491629 
 38.464193 44.490795 
 38.464193 44.491629 
 38.465859 44.493078 
 38.465859 44.494656 
 38.464586 44.492951 
 38.465026 44.494129 
 38.465026 44.494962 
 38.464193 44.494129 
 38.464193 44.494962 
 38.461682 44.484415 
 38.463359 44.489100 
 38.463359 44.490795 
 38.463359 44.491629 
 38.462526 44.489104 
 38.462526 44.491390 
 38.462958 44.493820 
 38.460972 44.489917 
 38.460996 44.493789 
 38.464248 44.502530 
 38.457147 44.484432 
 38.457143 44.491534 
 38.453288 44.482812 
 38.453290 44.486054 
 38.451693 44.481965 
 38.451693 44.483545 
 38.450357 44.482014 
 38.450860 44.483295 
 38.450860 44.484129 
 38.450026 44.483295 
 38.450026 44.484129 
 38.450726 44.486062 
 38.452009 44.491534 
 38.458431 44.497032 
 38.458430 44.500272 
 38.455860 44.497032 
 38.455860 44.500269 
 38.457143 44.505759 
 38.451998 44.498652 
 38.452009 44.505759 
 38.464248 44.516757 
 38.464248 44.530983 
 38.454574 44.516749 
 38.454574 44.530972 
 38.505688 44.565890 
 38.505690 44.622141 
 38.477863 44.551665 
 38.477865 44.579463 
 38.459110 44.551664 
 38.459112 44.579465 
 38.468190 44.622141 
 38.411938 44.256844 
 38.430688 44.340892 
 38.440366 44.382917 
 38.440366 44.411371 
 38.426748 44.375804 
 38.426748 44.390031 
 38.419647 44.371932 
 38.419644 44.379034 
 38.414498 44.371932 
 38.414509 44.379035 
 38.417074 44.390027 
 38.426748 44.404257 
 38.426748 44.418484 
 38.417074 44.404250 
 38.419647 44.414594 
 38.420931 44.419452 
 38.421693 44.421629 
 38.421693 44.422462 
 38.420860 44.421629 
 38.420860 44.422462 
 38.421693 44.423296 
 38.421693 44.424532 
 38.420860 44.423296 
 38.420860 44.424533 
 38.420026 44.421629 
 38.420026 44.422462 
 38.420026 44.423296 
 38.420026 44.424535 
 38.418360 44.419450 
 38.418681 44.421937 
 38.419193 44.423296 
 38.419193 44.424531 
 38.418360 44.423296 
 38.418360 44.424129 
 38.418360 44.424962 
 38.417526 44.421873 
 38.417526 44.424135 
 38.414498 44.414593 
 38.415789 44.419452 
 38.415801 44.423328 
 38.413223 44.419457 
 38.414193 44.421753 
 38.414193 44.424136 
 38.412862 44.421938 
 38.413360 44.423296 
 38.413360 44.424531 
 38.412526 44.423296 
 38.412526 44.424129 
 38.412526 44.424962 
 38.402864 44.326667 
 38.409980 44.344122 
 38.411261 44.349637 
 38.411693 44.352043 
 38.411693 44.353647 
 38.410860 44.352042 
 38.410860 44.353628 
 38.409287 44.349614 
 38.409287 44.352853 
 38.405432 44.344123 
 38.405431 44.351234 
 38.409980 44.358355 
 38.409974 44.365465 
 38.407040 44.355311 
 38.407042 44.356949 
 38.405860 44.355247 
 38.405860 44.356861 
 38.407045 44.358594 
 38.407059 44.360772 
 38.405860 44.358492 
 38.405860 44.360761 
 38.404146 44.356091 
 38.404154 44.359964 
 38.405431 44.365455 
 38.400888 44.344122 
 38.400891 44.351237 
 38.397621 44.342510 
 38.397640 44.346367 
 38.395066 44.342511 
 38.395074 44.346375 
 38.396364 44.351234 
 38.400888 44.358355 
 38.400891 44.365465 
 38.397641 44.356097 
 38.397622 44.359963 
 38.395074 44.356091 
 38.395860 44.358296 
 38.395860 44.359129 
 38.395026 44.358296 
 38.395026 44.359129 
 38.395415 44.360772 
 38.394193 44.358454 
 38.394193 44.360758 
 38.397641 44.363830 
 38.397643 44.367073 
 38.395074 44.363836 
 38.395075 44.367075 
 38.384111 44.326666 
 38.384113 44.354467 
 38.411261 44.370291 
 38.411693 44.372714 
 38.411693 44.374345 
 38.410860 44.372657 
 38.410860 44.374271 
 38.409287 44.370313 
 38.409285 44.373557 
 38.409980 44.379022 
 38.405431 44.371932 
 38.405432 44.379034 
 38.407406 44.390031 
 38.400891 44.371915 
 38.402526 44.376208 
 38.402526 44.377462 
 38.402526 44.379006 
 38.401237 44.376612 
 38.401224 44.378813 
 38.401577 44.381291 
 38.399610 44.377404 
 38.399610 44.381320 
 38.396364 44.371932 
 38.396355 44.379034 
 38.398327 44.390031 
 38.409974 44.400372 
 38.409980 44.407488 
 38.406717 44.398754 
 38.407038 44.401191 
 38.407039 44.402823 
 38.405860 44.401171 
 38.405860 44.402772 
 38.404145 44.398754 
 38.405026 44.401165 
 38.405026 44.402462 
 38.405026 44.403296 
 38.403817 44.401191 
 38.403816 44.402824 
 38.405432 44.407485 
 38.411261 44.413004 
 38.411693 44.415411 
 38.411693 44.417015 
 38.410860 44.415421 
 38.410860 44.416629 
 38.410860 44.417462 
 38.409287 44.412973 
 38.409634 44.415397 
 38.410026 44.416629 
 38.410026 44.417462 
 38.409193 44.416629 
 38.409193 44.417462 
 38.408360 44.415425 
 38.408360 44.417002 
 38.411693 44.418296 
 38.411693 44.419129 
 38.411693 44.419962 
 38.411693 44.420796 
 38.410860 44.418594 
 38.410860 44.420186 
 38.411693 44.421629 
 38.411693 44.422462 
 38.411693 44.423296 
 38.411693 44.424531 
 38.410860 44.421629 
 38.410860 44.422462 
 38.410860 44.424136 
 38.410026 44.418296 
 38.410026 44.419129 
 38.409193 44.418296 
 38.409193 44.419129 
 38.409633 44.420298 
 38.408360 44.418579 
 38.408360 44.420157 
 38.409270 44.423328 
 38.407526 44.411629 
 38.407526 44.412462 
 38.406693 44.411629 
 38.406693 44.412462 
 38.407526 44.413296 
 38.407526 44.414129 
 38.406693 44.413296 
 38.406693 44.414129 
 38.405860 44.412204 
 38.405860 44.413812 
 38.406717 44.416213 
 38.404145 44.412976 
 38.404145 44.416213 
 38.407039 44.418662 
 38.407040 44.420297 
 38.405860 44.418615 
 38.405860 44.419962 
 38.405860 44.420796 
 38.407043 44.421939 
 38.407059 44.424133 
 38.405860 44.421843 
 38.405860 44.423296 
 38.405860 44.424129 
 38.405860 44.424962 
 38.405026 44.418601 
 38.405026 44.420198 
 38.403816 44.418662 
 38.403816 44.420298 
 38.405026 44.421629 
 38.405026 44.422462 
 38.405026 44.424136 
 38.404193 44.421629 
 38.404193 44.422462 
 38.403360 44.421629 
 38.403360 44.422462 
 38.404193 44.423296 
 38.404193 44.424534 
 38.403360 44.423296 
 38.403360 44.424129 
 38.403360 44.424962 
 38.400891 44.400372 
 38.401579 44.405238 
 38.402526 44.408290 
 38.402526 44.410698 
 38.401693 44.407896 
 38.401693 44.409129 
 38.400860 44.407891 
 38.400860 44.409129 
 38.401693 44.409962 
 38.401693 44.410796 
 38.400860 44.409962 
 38.400860 44.410796 
 38.399610 44.405188 
 38.399610 44.409090 
 38.396364 44.400373 
 38.396355 44.407485 
 38.402526 44.411629 
 38.402526 44.412462 
 38.402526 44.413851 
 38.401223 44.412127 
 38.401221 44.413764 
 38.401577 44.416213 
 38.399610 44.413004 
 38.399610 44.416213 
 38.402526 44.418578 
 38.402526 44.419962 
 38.402526 44.420796 
 38.401221 44.418662 
 38.401693 44.419962 
 38.401693 44.420796 
 38.400860 44.419962 
 38.400860 44.420796 
 38.402526 44.421629 
 38.402526 44.422462 
 38.402526 44.423296 
 38.402526 44.424129 
 38.402526 44.424962 
 38.401225 44.421939 
 38.401237 44.424133 
 38.400026 44.418608 
 38.400026 44.419962 
 38.400026 44.420796 
 38.399193 44.418608 
 38.399193 44.419962 
 38.399193 44.420796 
 38.400026 44.421629 
 38.400026 44.422462 
 38.400026 44.424136 
 38.399193 44.421826 
 38.399193 44.424136 
 38.396364 44.414593 
 38.397642 44.419455 
 38.397621 44.423327 
 38.395074 44.419451 
 38.395860 44.421629 
 38.395860 44.422462 
 38.395026 44.421629 
 38.395026 44.422462 
 38.395415 44.424133 
 38.394193 44.421629 
 38.394193 44.422462 
 38.394193 44.424136 
 38.384113 44.382916 
 38.391214 44.400373 
 38.392503 44.405229 
 38.393360 44.408291 
 38.393360 44.410582 
 38.392161 44.408293 
 38.392176 44.410487 
 38.389934 44.405236 
 38.389950 44.409099 
 38.386683 44.400372 
 38.386684 44.407488 
 38.393360 44.412203 
 38.393360 44.413812 
 38.392179 44.412128 
 38.392526 44.413296 
 38.392526 44.414129 
 38.391693 44.413296 
 38.391693 44.414129 
 38.392502 44.416213 
 38.389933 44.412973 
 38.389932 44.416213 
 38.391220 44.421711 
 38.386683 44.414601 
 38.386684 44.421721 
 38.379574 44.404250 
 38.379574 44.418473 
 38.336939 44.256844 
 38.336939 44.369343 
 38.444899 44.432054 
 38.448764 44.440777 
 38.449193 44.443872 
 38.449193 44.444962 
 38.449193 44.445796 
 38.448360 44.443927 
 38.448360 44.444962 
 38.448360 44.445796 
 38.447526 44.439129 
 38.447526 44.440378 
 38.446693 44.439129 
 38.446693 44.440377 
 38.447132 44.442157 
 38.445860 44.439983 
 38.445860 44.441629 
 38.445860 44.442462 
 38.447526 44.443296 
 38.447526 44.444129 
 38.446693 44.443296 
 38.446693 44.444129 
 38.447133 44.445439 
 38.445860 44.443962 
 38.445860 44.445541 
 38.448761 44.447898 
 38.448761 44.451106 
 38.446787 44.447893 
 38.446787 44.451133 
 38.444213 44.440780 
 38.444217 44.444657 
 38.441654 44.440780 
 38.442526 44.443913 
 38.442526 44.445514 
 38.441315 44.443802 
 38.441316 44.445440 
 38.445026 44.446629 
 38.445026 44.447462 
 38.444193 44.446629 
 38.444193 44.447462 
 38.445026 44.448296 
 38.445026 44.449129 
 38.444193 44.448296 
 38.444193 44.449129 
 38.443360 44.447099 
 38.443360 44.448296 
 38.443360 44.449129 
 38.444217 44.451130 
 38.442526 44.447106 
 38.442526 44.448697 
 38.441693 44.446629 
 38.441693 44.447462 
 38.440860 44.446629 
 38.440860 44.447462 
 38.441317 44.448706 
 38.441645 44.451130 
 38.438391 44.428829 
 38.438391 44.435279 
 38.433864 44.428814 
 38.435143 44.433673 
 38.435141 44.436918 
 38.432903 44.432867 
 38.433360 44.434129 
 38.433360 44.434962 
 38.432526 44.434129 
 38.432526 44.434962 
 38.431693 44.432853 
 38.431693 44.434444 
 38.432574 44.436911 
 38.440026 44.439983 
 38.440026 44.442390 
 38.438737 44.439974 
 38.439193 44.441629 
 38.439193 44.442462 
 38.438360 44.441629 
 38.438360 44.442462 
 38.440026 44.443967 
 38.440026 44.445543 
 38.438723 44.443802 
 38.439193 44.444962 
 38.439193 44.445796 
 38.438360 44.444962 
 38.438360 44.445796 
 38.437526 44.439980 
 38.437526 44.441629 
 38.437526 44.442462 
 38.436693 44.439980 
 38.436693 44.442276 
 38.437526 44.443899 
 38.437526 44.445505 
 38.436693 44.443899 
 38.436693 44.445505 
 38.440026 44.447118 
 38.440026 44.448695 
 38.439193 44.446629 
 38.439193 44.447462 
 38.438360 44.446629 
 38.438360 44.447462 
 38.438721 44.448706 
 38.439077 44.451134 
 38.437526 44.447103 
 38.437526 44.448296 
 38.437526 44.449129 
 38.436693 44.446629 
 38.436693 44.447462 
 38.436693 44.448698 
 38.437110 44.451106 
 38.435483 44.439974 
 38.435495 44.442158 
 38.434193 44.439129 
 38.434193 44.439962 
 38.434193 44.440796 
 38.434193 44.442390 
 38.435142 44.444651 
 38.433360 44.439129 
 38.433360 44.440377 
 38.432526 44.439129 
 38.432526 44.439962 
 38.432526 44.440796 
 38.432906 44.442158 
 38.431693 44.439129 
 38.431693 44.440377 
 38.431693 44.441629 
 38.431693 44.442462 
 38.432905 44.443802 
 38.433360 44.444962 
 38.433360 44.445796 
 38.432526 44.444962 
 38.432526 44.445796 
 38.431693 44.443296 
 38.431693 44.444129 
 38.431693 44.444962 
 38.431693 44.445796 
 38.435499 44.447073 
 38.435860 44.448296 
 38.435860 44.449129 
 38.435026 44.448296 
 38.435026 44.449129 
 38.434193 44.446629 
 38.434193 44.447462 
 38.434193 44.448695 
 38.435499 44.450338 
 38.435860 44.451629 
 38.435860 44.452462 
 38.435026 44.451629 
 38.435026 44.452462 
 38.434193 44.450270 
 38.434193 44.451629 
 38.434193 44.452462 
 38.432575 44.447894 
 38.432574 44.451131 
 38.447480 44.456621 
 38.447474 44.463736 
 38.442932 44.456622 
 38.442931 44.463733 
 38.447480 44.470854 
 38.447474 44.477964 
 38.442932 44.470847 
 38.442931 44.477953 
 38.439098 44.455009 
 38.439080 44.458866 
 38.437110 44.455023 
 38.437526 44.458114 
 38.437526 44.459739 
 38.436693 44.458114 
 38.436693 44.459739 
 38.439077 44.462111 
 38.439077 44.465353 
 38.437110 44.462135 
 38.437110 44.465344 
 38.435496 44.453613 
 38.435483 44.455814 
 38.434193 44.453422 
 38.434193 44.455829 
 38.435140 44.458866 
 38.432905 44.453614 
 38.433360 44.454962 
 38.433360 44.456222 
 38.432526 44.454962 
 38.432526 44.456220 
 38.431693 44.453490 
 38.431693 44.455825 
 38.432574 44.458874 
 38.435143 44.462111 
 38.435143 44.465353 
 38.432574 44.462114 
 38.432903 44.464543 
 38.432903 44.466176 
 38.431693 44.464541 
 38.431693 44.466132 
 38.438388 44.470854 
 38.438391 44.477964 
 38.433855 44.470847 
 38.433864 44.477953 
 38.428714 44.428814 
 38.430002 44.433671 
 38.430860 44.435796 
 38.430860 44.436629 
 38.430860 44.437682 
 38.430026 44.435796 
 38.430026 44.436629 
 38.429193 44.435796 
 38.429193 44.436629 
 38.429677 44.437775 
 38.428360 44.432848 
 38.428360 44.434129 
 38.428360 44.434962 
 38.427526 44.432462 
 38.427526 44.433296 
 38.426693 44.432462 
 38.426693 44.433296 
 38.427526 44.434129 
 38.427526 44.434962 
 38.426693 44.434129 
 38.426693 44.434962 
 38.427434 44.436915 
 38.424183 44.428829 
 38.425459 44.433659 
 38.425459 44.436872 
 38.423497 44.433674 
 38.424193 44.435796 
 38.424193 44.436629 
 38.423360 44.435796 
 38.423360 44.436629 
 38.423859 44.437775 
 38.422526 44.436010 
 38.422526 44.437462 
 38.422526 44.438296 
 38.430860 44.439129 
 38.430860 44.439962 
 38.430860 44.440796 
 38.430860 44.442258 
 38.430026 44.439129 
 38.430026 44.439962 
 38.430026 44.440796 
 38.429193 44.439129 
 38.429193 44.440377 
 38.429675 44.442159 
 38.430003 44.444657 
 38.428360 44.439983 
 38.428360 44.442377 
 38.427526 44.439129 
 38.427526 44.439962 
 38.427526 44.440796 
 38.426693 44.439129 
 38.426693 44.440378 
 38.427526 44.441629 
 38.427526 44.442462 
 38.426693 44.441629 
 38.426693 44.442462 
 38.428360 44.443962 
 38.428360 44.445540 
 38.427086 44.443802 
 38.427526 44.444962 
 38.427526 44.445796 
 38.426693 44.444962 
 38.426693 44.445796 
 38.430002 44.447894 
 38.430003 44.451130 
 38.428360 44.447117 
 38.428360 44.448695 
 38.427526 44.446629 
 38.427526 44.447462 
 38.426693 44.446629 
 38.426693 44.447462 
 38.427085 44.448706 
 38.427433 44.451133 
 38.425455 44.440777 
 38.425459 44.444688 
 38.424193 44.439129 
 38.424193 44.440378 
 38.423360 44.439129 
 38.423360 44.440377 
 38.423857 44.442160 
 38.422526 44.439982 
 38.422526 44.442359 
 38.423860 44.443803 
 38.424193 44.444962 
 38.424193 44.445796 
 38.423360 44.444962 
 38.423360 44.445796 
 38.422526 44.443953 
 38.422526 44.445536 
 38.425459 44.447898 
 38.425459 44.451106 
 38.424193 44.446629 
 38.424193 44.447462 
 38.423360 44.446629 
 38.423360 44.447462 
 38.424193 44.448296 
 38.424193 44.449129 
 38.423360 44.448296 
 38.423360 44.449129 
 38.422526 44.446629 
 38.422526 44.447462 
 38.422526 44.448296 
 38.422526 44.449129 
 38.423496 44.451136 
 38.420930 44.427193 
 38.420931 44.430434 
 38.419193 44.425796 
 38.419193 44.426629 
 38.418360 44.425796 
 38.418360 44.426629 
 38.418679 44.427966 
 38.417526 44.426390 
 38.417526 44.428020 
 38.418360 44.430435 
 38.420931 44.433672 
 38.421269 44.436136 
 38.421693 44.437462 
 38.421693 44.438296 
 38.420860 44.437462 
 38.420860 44.438296 
 38.420026 44.436019 
 38.420026 44.437613 
 38.418360 44.433671 
 38.418360 44.436909 
 38.416693 44.426497 
 38.416693 44.428097 
 38.415451 44.426321 
 38.415860 44.427462 
 38.415860 44.428296 
 38.415026 44.427462 
 38.415026 44.428296 
 38.416693 44.429129 
 38.416693 44.429962 
 38.416693 44.431268 
 38.415860 44.429129 
 38.415860 44.429962 
 38.415026 44.429129 
 38.415026 44.429962 
 38.415451 44.431235 
 38.413225 44.427185 
 38.414193 44.429129 
 38.414193 44.429962 
 38.414193 44.430796 
 38.414193 44.431629 
 38.413360 44.429129 
 38.413360 44.429962 
 38.412526 44.429129 
 38.412526 44.429962 
 38.413360 44.430796 
 38.413360 44.431629 
 38.412526 44.430796 
 38.412526 44.431629 
 38.415788 44.433672 
 38.416693 44.436019 
 38.416693 44.437613 
 38.415860 44.435796 
 38.415860 44.436629 
 38.415026 44.435796 
 38.415026 44.436629 
 38.415451 44.437777 
 38.414193 44.432462 
 38.414193 44.433296 
 38.414193 44.434428 
 38.412857 44.432867 
 38.412857 44.434499 
 38.413224 44.436920 
 38.419644 44.442397 
 38.421693 44.446629 
 38.421693 44.447462 
 38.420860 44.446629 
 38.420860 44.447462 
 38.421269 44.448706 
 38.420026 44.447113 
 38.420026 44.448296 
 38.420026 44.449129 
 38.421269 44.450338 
 38.421693 44.451629 
 38.421693 44.452462 
 38.420860 44.451629 
 38.420860 44.452462 
 38.420026 44.450278 
 38.420026 44.451863 
 38.418360 44.447894 
 38.419193 44.449962 
 38.419193 44.450796 
 38.418360 44.449962 
 38.418360 44.450796 
 38.418678 44.451971 
 38.417526 44.449962 
 38.417526 44.450796 
 38.417526 44.451925 
 38.414509 44.442397 
 38.415788 44.447893 
 38.416693 44.450278 
 38.416693 44.451863 
 38.415860 44.449962 
 38.415860 44.450796 
 38.415026 44.449962 
 38.415026 44.450796 
 38.415860 44.451629 
 38.415860 44.452462 
 38.415026 44.451629 
 38.415026 44.452462 
 38.413222 44.447893 
 38.414193 44.450274 
 38.414193 44.451855 
 38.412857 44.450337 
 38.413360 44.451629 
 38.413360 44.452462 
 38.412526 44.451629 
 38.412526 44.452462 
 38.430860 44.453522 
 38.430860 44.454962 
 38.430860 44.456222 
 38.429677 44.453612 
 38.429660 44.455814 
 38.430860 44.458095 
 38.430860 44.459724 
 38.429675 44.457995 
 38.429678 44.459639 
 38.428360 44.453429 
 38.428360 44.455828 
 38.427526 44.453296 
 38.427526 44.454129 
 38.426693 44.453296 
 38.426693 44.454129 
 38.427092 44.455814 
 38.427435 44.458869 
 38.430860 44.461337 
 38.430860 44.462941 
 38.429680 44.461277 
 38.429681 44.462911 
 38.430002 44.465352 
 38.428360 44.460795 
 38.428360 44.461629 
 38.428360 44.462963 
 38.427086 44.461277 
 38.427086 44.462911 
 38.427432 44.465353 
 38.425455 44.455023 
 38.425860 44.458155 
 38.425860 44.459770 
 38.425026 44.458079 
 38.425026 44.459711 
 38.423472 44.455009 
 38.423494 44.458863 
 38.425459 44.462135 
 38.425860 44.464541 
 38.425860 44.466127 
 38.425026 44.464542 
 38.425026 44.466147 
 38.424193 44.460795 
 38.424193 44.461629 
 38.423360 44.460795 
 38.423360 44.461629 
 38.423863 44.462911 
 38.422526 44.461381 
 38.422526 44.462961 
 38.424193 44.464129 
 38.424193 44.464962 
 38.423360 44.464129 
 38.423360 44.464962 
 38.423863 44.466175 
 38.422526 44.464540 
 38.422526 44.466119 
 38.428720 44.470847 
 38.428714 44.477953 
 38.425459 44.468554 
 38.425455 44.472463 
 38.423862 44.467809 
 38.423860 44.469448 
 38.422526 44.467699 
 38.422526 44.469129 
 38.422526 44.469962 
 38.423472 44.472463 
 38.424183 44.477964 
 38.421693 44.453296 
 38.421693 44.454129 
 38.420860 44.453296 
 38.420860 44.454129 
 38.421693 44.454962 
 38.421693 44.455795 
 38.421693 44.456629 
 38.420860 44.454962 
 38.420860 44.456221 
 38.420026 44.453456 
 38.420026 44.455826 
 38.420929 44.458871 
 38.419193 44.453296 
 38.419193 44.454129 
 38.418360 44.453296 
 38.418360 44.454129 
 38.418704 44.455814 
 38.417526 44.453550 
 38.417526 44.455820 
 38.419193 44.457462 
 38.419193 44.458295 
 38.418360 44.457462 
 38.418360 44.458295 
 38.419193 44.459129 
 38.419193 44.459962 
 38.418360 44.459129 
 38.418360 44.459962 
 38.417526 44.458064 
 38.417526 44.459698 
 38.419647 44.463733 
 38.415801 44.455011 
 38.415790 44.458871 
 38.414193 44.453441 
 38.414193 44.455827 
 38.413360 44.453296 
 38.413360 44.454129 
 38.412526 44.453296 
 38.412526 44.454129 
 38.413360 44.454962 
 38.413360 44.455795 
 38.413360 44.456629 
 38.412526 44.454962 
 38.412526 44.456218 
 38.414193 44.458200 
 38.414193 44.459797 
 38.413360 44.457462 
 38.413360 44.458295 
 38.412526 44.457462 
 38.412526 44.458295 
 38.412860 44.459640 
 38.414498 44.463733 
 38.420930 44.468592 
 38.420918 44.472463 
 38.418677 44.467809 
 38.419193 44.469129 
 38.419193 44.469962 
 38.418360 44.469129 
 38.418360 44.469962 
 38.417526 44.467765 
 38.417526 44.469388 
 38.418360 44.472463 
 38.419647 44.477953 
 38.414509 44.470847 
 38.414498 44.477953 
 38.448761 44.482790 
 38.448761 44.486005 
 38.446787 44.482813 
 38.447526 44.484962 
 38.447526 44.485795 
 38.446693 44.484962 
 38.446693 44.485795 
 38.447132 44.486931 
 38.445860 44.485118 
 38.445860 44.486705 
 38.447480 44.491522 
 38.442931 44.484432 
 38.442932 44.491534 
 38.444906 44.502530 
 38.435827 44.488304 
 38.435827 44.502530 
 38.440366 44.523870 
 38.426748 44.488304 
 38.426748 44.502530 
 38.419647 44.484432 
 38.419644 44.491534 
 38.416693 44.481968 
 38.416693 44.483552 
 38.415451 44.482014 
 38.415451 44.483648 
 38.415790 44.486054 
 38.413223 44.482815 
 38.413226 44.486062 
 38.414509 44.491534 
 38.417074 44.502526 
 38.428714 44.512872 
 38.428720 44.519985 
 38.424183 44.512871 
 38.424184 44.519988 
 38.426748 44.530983 
 38.417074 44.516749 
 38.419647 44.527093 
 38.419644 44.534210 
 38.414498 44.527093 
 38.414509 44.534210 
 38.411261 44.427238 
 38.411693 44.429646 
 38.411693 44.431254 
 38.410860 44.429677 
 38.410860 44.430796 
 38.410860 44.431629 
 38.409285 44.427190 
 38.409633 44.429601 
 38.410026 44.430796 
 38.410026 44.431629 
 38.409193 44.430796 
 38.409193 44.431629 
 38.408360 44.429694 
 38.408360 44.431272 
 38.411693 44.432858 
 38.411693 44.434463 
 38.410860 44.432462 
 38.410860 44.433296 
 38.410860 44.434438 
 38.411261 44.436872 
 38.410026 44.432462 
 38.410026 44.433296 
 38.409193 44.432462 
 38.409193 44.433296 
 38.409634 44.434500 
 38.408360 44.432848 
 38.408360 44.434425 
 38.409286 44.436915 
 38.407043 44.426323 
 38.407041 44.427965 
 38.405860 44.425796 
 38.405860 44.426629 
 38.405860 44.428045 
 38.406717 44.430435 
 38.405026 44.426457 
 38.405026 44.428071 
 38.404193 44.425796 
 38.404193 44.426629 
 38.403360 44.425796 
 38.403360 44.426629 
 38.404193 44.427462 
 38.404193 44.428296 
 38.403360 44.427462 
 38.403360 44.428296 
 38.405026 44.429670 
 38.405026 44.431263 
 38.404193 44.429129 
 38.404193 44.429962

Xin

Thanks, that helps.

This one is subtle in that “v” gets set before the parallel loop so should be used in a “firstprivate” clause. “private” variables are uninitialized, while “firstprivate” get initialize to the value from the host. Though, “firstprivate” is only available on “parallel” loops, so let switch your compute region to use:

 !$acc parallel loop gang vector firstprivate(v) private(latf,lonf,lats,lons,t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,u1,u2,dep,gx,gy,gz,x,y,z)

With this change, the result I computed on the CPU match the ones computed on the GPU.

% pgfortran func.f main.f90 -Minfo=accel -ta=tesla -o gpu1910.out -V19.10; ./gpu1910.out
func.f:
indisp:
      4, Generating acc routine seq
         Generating Tesla code
defor:
    193, Generating acc routine seq
         Generating Tesla code
main.f90:
PGF90-W-0155-Constant or Parameter used in data clause - re (main.f90: 64)
  0 inform,   1 warnings,   0 severes, 0 fatal for main
main:
     64, Generating copy(obs(:)) [if not already present]
         Generating copyin(fsm(:)) [if not already present]
     65, Generating Tesla code
         66, !$acc loop gang, vector(128) ! blockidx%x threadidx%x
         74, !$acc loop seq
    100, Reference argument passing prevents parallelization: gy
         Reference argument passing prevents parallelization: gx
         Reference argument passing prevents parallelization: dep
         Reference argument passing prevents parallelization: z
         Reference argument passing prevents parallelization: y
         Reference argument passing prevents parallelization: x
         Reference argument passing prevents parallelization: v
         Reference argument passing prevents parallelization: gz
         Reference argument passing prevents parallelization: u2
         Reference argument passing prevents parallelization: u1
         Reference argument passing prevents parallelization: t7
 it takes   0.2446420     s
% diff displacement_Okada.dat displacement_Okada.cpu.dat
%

-Mat

Thank you very much for your help, Mat. I found the important thing is using “$acc parallel loop …” to replace “$acc kernels loop …”. I can obtain a correct result evenif I delete “firstprivate(v)” and “private()”. I have not known why “$acc kernels loop …” cannot work by 19.10 in my case but it can work well in the following example.

program main
    use openacc
    implicit none
    integer :: i,j,k
    integer,parameter :: n=5000
    real*8 :: a(n,n), b(n,n), c(n,n)
    real*8 :: sum, start, finish, s

    call CPU_TIME(start)

    a = 1.d0
    forall(i=1:n) a(i,i) = dble(i)

    b = 1.d0
    c = 0.d0

    !$acc data copyout(c) copyin(a,b)
    !$acc kernels loop gang vector reduction(+:sum)
    do i=1,n
       do k=1,n
          sum = 0.d0
          do j=1,n
             sum = sum+a(i,j)*b(j,k)
          enddo
          c(i,k) = sum
       enddo
     enddo
     !$acc end data

    call CPU_TIME(finish)

    write(*,*) ((c(i,j),j=1,10),i=1,10)
    write(*,*) 'it takes',finish-start,'s'

end program