Changeset 657fa36 in flexpart.git


Ignore:
Timestamp:
Oct 18, 2016, 11:13:34 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:
15fdb5c
Parents:
5deb48c
Message:

usekernel off, bkwd drydeposit calculation deleted

Location:
src
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • src/conccalc.f90

    r54cbd6c r657fa36  
    6363  integer :: usekernel
    6464
    65   usekernel=0
     65  usekernel=1
    6666  if (usekernel.ne.1) then
    6767     write (*,*) 'NOT USING THE KERNEL!'
  • src/timemanager.f90

    r5deb48c r657fa36  
    115115  real :: xold,yold,zold,xmassfract
    116116  real, parameter :: e_inv = 1.0/exp(1.0)
    117   logical :: firstdepocalc
     117
    118118  !double precision xm(maxspec,maxpointspec_act),
    119119  !    +                 xm_depw(maxspec,maxpointspec_act),
     
    377377
    378378      if ((itime.eq.loutend).and.(outnum.gt.0.)) then
    379         if ((iout.le.3.).or.(iout.eq.5).or.(iout.eq.6)) then
     379        if ((iout.le.3.).or.(iout.eq.5)) then
    380380          if (surf_only.ne.1) then
    381381            if (lnetcdfout.eq.1) then
     
    552552      do ks=1,nspec
    553553         if  ((xscav_frac1(j,ks).lt.0)) then
    554          call advance_rec(itime,xtra1(j),ytra1(j),ztra1(j),prob_rec)
     554            call advance_rec(itime,xtra1(j),ytra1(j),ztra1(j),prob_rec)
    555555            if (decay(ks).gt.0.) then             ! radioactive decay
    556556                decfact=exp(-real(abs(lsynctime))*decay(ks))
    557557            else
    558                  decfact=1.
     558                decfact=1.
    559559            endif
    560560            if (DRYDEPSPEC(ks)) then        ! dry deposition
     
    562562               xscav_frac1(j,ks)=xscav_frac1(j,ks)*(-1.)* &
    563563               drydeposit(ks)/xmass1(j,ks)
    564 !               write (*,*) 'notance: ',prob(ks),xmass1(j,ks),ztra1(j)
    565                if (decay(ks).gt.0.) then   ! correct for decay (see wetdepo)
    566                   drydeposit(ks)=drydeposit(ks)* &
    567                   exp(real(abs(ldeltat))*decay(ks))
    568                 endif
     564               !if (xscav_frac1(j,ks).eq.1) then
     565                  write (*,*) 'xscav 1: ',j,prob_rec(ks),xmass1(j,ks), &
     566                                ztra1(j),drydeposit(ks),decfact
     567              !     stop
     568              ! endif
    569569             else
    570570                xmass1(j,ks)=0
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG