Changeset 1d072c0 in flexpart.git
- Timestamp:
- Aug 7, 2017, 2:43:38 PM (7 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
- Children:
- aa8c34a
- Parents:
- b5127f9
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/coordtrafo.f90
re200b7a r1d072c0 46 46 47 47 integer :: i,j,k 48 real :: yrspc ! small real number relative to x 48 49 49 50 if (numpoint.eq.0) goto 30 … … 64 65 ! CHECK IF RELEASE POINTS ARE WITHIN DOMAIN 65 66 !****************************************** 66 67 68 yrspc = spacing(real(nymin1,kind=sp)) 69 67 70 do i=1,numpoint 68 71 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-571 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) &72 if (nglobal.and.(ypoint2(i).gt.real(nymin1,kind=dp)-1.e-5)) & 73 ypoint2(i)=real(nymin1,kind=dp)-10*yrspc 74 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) & 73 76 .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)))) then77 (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)))) then 76 79 write(*,*) ' NOTICE: RELEASE POINT OUT OF DOMAIN DETECTED.' 77 80 write(*,*) ' IT IS REMOVED NOW ... ' 78 if (i. ge.1000) then81 if (i.le.1000) then 79 82 write(*,*) ' COMMENT: ',compoint(i) 80 83 else
Note: See TracChangeset
for help on using the changeset viewer.