Changeset 8ee24a5 in flexpart.git for src

Timestamp:

Nov 24, 2016, 9:23:13 AM (7 years ago)

Author:

Sabine <sabine.eckhardt@…>

Branches:

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

Children:

1c0d5e6

Parents:

ceb7e4d

Message:

force releaseheights for bkwd depo runs

Location:

src

Files:

• get_wetscav.f90 (modified) (2 diffs)
• readcommand.f90 (modified) (2 diffs)
• readreleases.f90 (modified) (1 diff)

src/get_wetscav.f90

@@ -21 +21 @@
 
 subroutine get_wetscav(itime,ltsample,loutnext,jpart,ks,grfraction,inc_count,blc_count,wetscav)
-!                          i      i        i     i     o           o
+!                          i      i        i     i   i    o           o          o       o
 !*****************************************************************************
 !                                                                            *
@@ -139 +139 @@
     endif
 
+!  If total precipitation is less than 0.01 mm/h - no scavenging occurs
+    if ((lsp.lt.0.01).and.(convp.lt.0.01)) goto 20
 
 ! get the level were the actual particle is in

src/readcommand.f90

@@ -330 +330 @@
         ind_rel = 3
          write(*,*) ' #### FLEXPART WET DEPOSITION BACKWARD MODE    #### '
+         write(*,*) ' #### Releaseheight is forced to 0 - 20km      #### '
+         write(*,*) ' #### Release is performed above ground lev    #### '
          WETBKDEP=.true.
          allocate(xscav_frac1(maxpart,maxspec))
@@ -335 +337 @@
          ind_rel = 4
          write(*,*) ' #### FLEXPART DRY DEPOSITION BACKWARD MODE    #### '
+         write(*,*) ' #### Releaseheight is forced to 0 - 2*href    #### '
+         write(*,*) ' #### Release is performed above ground lev    #### '
          DRYBKDEP=.true.
          allocate(xscav_frac1(maxpart,maxspec))

src/readreleases.f90

@@ -516 +516 @@
   endif
 
+  ! If FLEXPART is run for backward deposition force zpoint
+  !*********************************************************************
+  if (WETBKDEP) then
+      zpoint1(numpoint)=0.
+      zpoint2(numpoint)=20000.
+      kindz(numpoint)=1
+  endif
+  if (DRYBKDEP) then
+      zpoint1(numpoint)=0.
+      zpoint2(numpoint)=2.*href
+      kindz(numpoint)=1
+  endif
+
+
   ! Check whether x coordinates of release point are within model domain
   !*********************************************************************