Changeset 07c3e71 in flexpart.git for src/readspecies.f90

Timestamp:

Jul 24, 2019, 12:16:51 PM (5 years ago)

Author:

Ignacio Pisso <ip@…>

Branches:

master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug

Children:

25b4532

Parents:

95a45d3

Message:

check dsigma for dry deposition velocity

File:

• src/readspecies.f90 (modified) (2 diffs)

src/readspecies.f90

@@ -296 +296 @@
       if (density(pos_spec) .gt. 0) then
         write(*,'(a)') '  Dry deposition is turned         :   ON'
-      else
+        if (reldiff(pos_spec).gt.0) then
+           stop 'density>0 (SPECIES is a particle) implies reldiff <=0  '
+        endif
+      else
+        if (reldiff(pos_spec).le.0) then
+           stop 'density<=0 (SPECIES is a gas) implies reldiff >0  '
+        endif      
         write(*,'(a)') '  Dry deposition is (density<0)    :   OFF'
       end if
@@ -344 +350 @@
   end if
 
-  if (dsigma(i).eq.0.) dsigma(i)=1.0001   ! avoid floating exception
+  !  if (dsigma(i).eq.0.) dsigma(i)=1.0001   ! avoid floating exception
+  if (dsigma(i).le.1.) then !dsigma(i)=1.0001   ! avoid floating exception
+    !write(*,*) '#### FLEXPART MODEL ERROR!                      ####'
+    write(*,*) '#### FLEXPART MODEL WARNING                     ####'
+    write(*,*) '#### in SPECIES_',aspecnumb, '                             ####'
+    write(*,*) '#### from v10.4 dsigma has to be larger than 1  ####'  
+    write(*,*) '#### to adapt older SPECIES files,              ####' 
+    write(*,*) '#### if dsigma was < 1                          ####' 
+    write(*,*) '#### use the reciprocal of the old dsigma       ####'  
+    if (.not.debug_mode) then 
+       stop
+    else
+       write(*,*) 'debug mode: continue'
+    endif
+  endif
 
   if ((reldiff(i).gt.0.).and.(density(i).gt.0.)) then