49 if (numpoint.eq.0) goto 30
55 xpoint1(i)=(xpoint1(i)-xlon0)/dx
56 xpoint2(i)=(xpoint2(i)-xlon0)/dx
57 ypoint1(i)=(ypoint1(i)-ylat0)/dy
58 ypoint2(i)=(ypoint2(i)-ylat0)/dy
68 if (sglobal.and.(ypoint1(i).lt.1.e-6)) ypoint1(i)=1.e-6
69 if (nglobal.and.(ypoint2(i).gt.
real(nymin1)-1.e-5)) &
70 ypoint2(i)=
real(nymin1)-1.e-5
71 if ((ypoint1(i).lt.1.e-6).or.(ypoint1(i).ge.
real(nymin1)-1.e-6) &
72 .or.(ypoint2(i).lt.1.e-6).or.(ypoint2(i).ge.
real(nymin1)-1.e-6) &
73 .or.((.not.xglobal).and.((xpoint1(i).lt.1.e-6).or. &
74 (xpoint1(i).ge.
real(nxmin1)-1.e-6).or.(xpoint2(i).lt.1.e-6).or. &
75 (xpoint2(i).ge.
real(nxmin1)-1.e-6)))) then
76 write(*,*)
' NOTICE: RELEASE POINT OUT OF DOMAIN DETECTED.'
77 write(*,*)
' IT IS REMOVED NOW ... '
79 write(*,*)
' COMMENT: ',compoint(i)
81 write(*,*)
' COMMENT: ',compoint(1001)
83 if (i.lt.numpoint)
then
85 xpoint1(j-1)=xpoint1(j)
86 ypoint1(j-1)=ypoint1(j)
87 xpoint2(j-1)=xpoint2(j)
88 ypoint2(j-1)=ypoint2(j)
89 zpoint1(j-1)=zpoint1(j)
90 zpoint2(j-1)=zpoint2(j)
93 ireleasestart(j-1)=ireleasestart(j)
94 ireleaseend(j-1)=ireleaseend(j)
95 if (j.le.1000) compoint(j-1)=compoint(j)
97 xmass(j-1,k)=xmass(j,k)
103 if (numpoint.gt.0) goto 15
107 30
if(numpoint.eq.0)
then
108 write(*,*)
' FLEXPART MODEL SUBROUTINE COORDTRAFO: ERROR ! '
109 write(*,*)
' NO PARTICLE RELEASES ARE DEFINED!'
110 write(*,*)
' CHECK FILE RELEASES...'