Changeset 7abd262 in flexpart.git


Ignore:
Timestamp:
Sep 12, 2019, 7:14:16 AM (5 years ago)
Author:
Espen Sollum ATMOS <eso@…>
Branches:
master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug
Children:
c56da45
Parents:
74d96de
Message:

Changed an error condition in advance.f90

File:
1 edited

Legend:

Unmodified
Added
Removed
  • src/advance.f90

    r5184a7c r7abd262  
    123123  real :: rhoa,rhograd,delz,dtf,rhoaux,dtftlw,uxscale,wpscale
    124124  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)
    126126  real :: ptot_lhh,Q_lhh,phi_lhh,ath,bth !modified by mc
    127127  real :: old_wp_buf,dcas,dcas1,del_test !added by mc
     
    537537        if (lsettling) then
    538538          do nsp=1,nspec
    539             if (xmass(nrelpoint,nsp).gt.eps2) exit
     539            if (xmass(nrelpoint,nsp).gt.eps3) exit
    540540          end do
    541541          if (nsp.gt.nspec) then
    542   ! This should never happen         
    543             write(*,*) 'advance.f90: ERROR: could not find releasepoint'
    544             stop
     542            nsp=nspec
    545543          end if
    546544          if (density(nsp).gt.0.) then
     
    707705    if (lsettling) then
    708706      do nsp=1,nspec
    709         if (xmass(nrelpoint,nsp).gt.eps2) exit
     707        if (xmass(nrelpoint,nsp).gt.eps3) exit
    710708      end do
    711709      if (nsp.gt.nspec) then
    712   ! This should never happen         
    713         write(*,*) 'advance.f90: ERROR: could not find releasepoint'
    714         stop
     710        nsp=nspec
    715711      end if
    716712      if (density(nsp).gt.0.) then
     
    917913    if (lsettling) then
    918914      do nsp=1,nspec
    919         if (xmass(nrelpoint,nsp).gt.eps2) exit
     915        if (xmass(nrelpoint,nsp).gt.eps3) exit
    920916      end do
    921917      if (nsp.gt.nspec) then
    922   ! This should never happen         
    923         write(*,*) 'advance.f90: ERROR: could not find releasepoint'
    924         stop
     918        nsp=nspec
    925919      end if
    926920      if (density(nsp).gt.0.) then
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG