Changes in src/coordtrafo.f90 [1d072c0:e200b7a] in flexpart.git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/coordtrafo.f90
r1d072c0 re200b7a 46 46 47 47 integer :: i,j,k 48 real :: yrspc ! small real number relative to x49 48 50 49 if (numpoint.eq.0) goto 30 … … 65 64 ! CHECK IF RELEASE POINTS ARE WITHIN DOMAIN 66 65 !****************************************** 67 68 yrspc = spacing(real(nymin1,kind=sp)) 69 66 70 67 do i=1,numpoint 71 68 if (sglobal.and.(ypoint1(i).lt.1.e-6)) ypoint1(i)=1.e-6 72 if (nglobal.and.(ypoint2(i).gt.real(nymin1 ,kind=dp)-1.e-5)) &73 ypoint2(i)=real(nymin1 ,kind=dp)-10*yrspc74 if ((ypoint1(i).lt.1.e-6).or.(ypoint1(i).ge.real(nymin1,kind=dp)-1.e-6) &75 .or.(ypoint2(i).lt.1.e-6).or.(ypoint2(i).ge.real(nymin1 ,kind=dp)-yrspc) &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) & 76 73 .or.((.not.xglobal).and.((xpoint1(i).lt.1.e-6).or. & 77 (xpoint1(i).ge.real(nxmin1 ,kind=dp)-1.e-6).or.(xpoint2(i).lt.1.e-6).or. &78 (xpoint2(i).ge.real(nxmin1 ,kind=dp)-1.e-6)))) then74 (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 79 76 write(*,*) ' NOTICE: RELEASE POINT OUT OF DOMAIN DETECTED.' 80 77 write(*,*) ' IT IS REMOVED NOW ... ' 81 if (i. le.1000) then78 if (i.ge.1000) then 82 79 write(*,*) ' COMMENT: ',compoint(i) 83 80 else
Note: See TracChangeset
for help on using the changeset viewer.