Changeset 7abd262 in flexpart.git
- Timestamp:
- Sep 12, 2019, 7:14:16 AM (5 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug
- Children:
- c56da45
- Parents:
- 74d96de
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/advance.f90
r5184a7c r7abd262 123 123 real :: rhoa,rhograd,delz,dtf,rhoaux,dtftlw,uxscale,wpscale 124 124 integer(kind=2) :: icbt 125 real,parameter :: eps=nxmax/3.e5,eps2=1.e-9 125 real,parameter :: eps=nxmax/3.e5,eps2=1.e-9,eps3=tiny(1.0) 126 126 real :: ptot_lhh,Q_lhh,phi_lhh,ath,bth !modified by mc 127 127 real :: old_wp_buf,dcas,dcas1,del_test !added by mc … … 537 537 if (lsettling) then 538 538 do nsp=1,nspec 539 if (xmass(nrelpoint,nsp).gt.eps 2) exit539 if (xmass(nrelpoint,nsp).gt.eps3) exit 540 540 end do 541 541 if (nsp.gt.nspec) then 542 ! This should never happen 543 write(*,*) 'advance.f90: ERROR: could not find releasepoint' 544 stop 542 nsp=nspec 545 543 end if 546 544 if (density(nsp).gt.0.) then … … 707 705 if (lsettling) then 708 706 do nsp=1,nspec 709 if (xmass(nrelpoint,nsp).gt.eps 2) exit707 if (xmass(nrelpoint,nsp).gt.eps3) exit 710 708 end do 711 709 if (nsp.gt.nspec) then 712 ! This should never happen 713 write(*,*) 'advance.f90: ERROR: could not find releasepoint' 714 stop 710 nsp=nspec 715 711 end if 716 712 if (density(nsp).gt.0.) then … … 917 913 if (lsettling) then 918 914 do nsp=1,nspec 919 if (xmass(nrelpoint,nsp).gt.eps 2) exit915 if (xmass(nrelpoint,nsp).gt.eps3) exit 920 916 end do 921 917 if (nsp.gt.nspec) then 922 ! This should never happen 923 write(*,*) 'advance.f90: ERROR: could not find releasepoint' 924 stop 918 nsp=nspec 925 919 end if 926 920 if (density(nsp).gt.0.) then
Note: See TracChangeset
for help on using the changeset viewer.