Changeset 657fa36 in flexpart.git
- Timestamp:
- Oct 18, 2016, 11:13:34 AM (7 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
- Children:
- 15fdb5c
- Parents:
- 5deb48c
- Location:
- src
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
src/conccalc.f90
r54cbd6c r657fa36 63 63 integer :: usekernel 64 64 65 usekernel= 065 usekernel=1 66 66 if (usekernel.ne.1) then 67 67 write (*,*) 'NOT USING THE KERNEL!' -
src/timemanager.f90
r5deb48c r657fa36 115 115 real :: xold,yold,zold,xmassfract 116 116 real, parameter :: e_inv = 1.0/exp(1.0) 117 logical :: firstdepocalc 117 118 118 !double precision xm(maxspec,maxpointspec_act), 119 119 ! + xm_depw(maxspec,maxpointspec_act), … … 377 377 378 378 if ((itime.eq.loutend).and.(outnum.gt.0.)) then 379 if ((iout.le.3.).or.(iout.eq.5) .or.(iout.eq.6)) then379 if ((iout.le.3.).or.(iout.eq.5)) then 380 380 if (surf_only.ne.1) then 381 381 if (lnetcdfout.eq.1) then … … 552 552 do ks=1,nspec 553 553 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) 555 555 if (decay(ks).gt.0.) then ! radioactive decay 556 556 decfact=exp(-real(abs(lsynctime))*decay(ks)) 557 557 else 558 558 decfact=1. 559 559 endif 560 560 if (DRYDEPSPEC(ks)) then ! dry deposition … … 562 562 xscav_frac1(j,ks)=xscav_frac1(j,ks)*(-1.)* & 563 563 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 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 569 569 else 570 570 xmass1(j,ks)=0
Note: See TracChangeset
for help on using the changeset viewer.