- Timestamp:
- Oct 24, 2016, 1:03:58 PM (8 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
- Children:
- deaac29
- Parents:
- d4a4250
- Location:
- src
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
src/advance.f90
r5844866 r0ee2193 569 569 prob(ks)=1.+(prob(ks)-1.)* & 570 570 exp(-vdepo(ks)*abs(dt)/(2.*href)) 571 ! write(*,*) 'Prob calc: ',zt,href,prob(ks)572 571 endif 573 572 end do -
src/advance_rec.f90
r5844866 r0ee2193 168 168 ! where f(n) is the exponential term 169 169 prob(ks)=1.+(prob(ks)-1.)* & 170 exp(-vdepo(ks)*abs( ideltas)/(2.*href))170 exp(-vdepo(ks)*abs(lsynctime)/(2.*href)) 171 171 endif 172 172 end do -
src/timemanager.f90
r15fdb5c r0ee2193 553 553 if ((xscav_frac1(j,ks).lt.0)) then 554 554 call advance_rec(itime,xtra1(j),ytra1(j),ztra1(j),prob_rec) 555 555 556 if (decay(ks).gt.0.) then ! radioactive decay 556 557 decfact=exp(-real(abs(lsynctime))*decay(ks)) … … 562 563 xscav_frac1(j,ks)=xscav_frac1(j,ks)*(-1.)* & 563 564 drydeposit(ks)/xmass1(j,ks) 564 !if (xscav_frac1(j,ks).eq.1) then565 ! write (*,*) 'xscav 1: ',j,prob_rec(ks),xmass1(j,ks), &566 ! ztra1(j),drydeposit(ks),decfact567 ! stop568 ! endif569 565 else 570 566 xmass1(j,ks)=0
Note: See TracChangeset
for help on using the changeset viewer.