Changeset 6473ad3 in flexpart.git for src/conccalc.f90
- Timestamp:
- Feb 9, 2017, 5:02:33 PM (7 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
- Children:
- 33279f7
- Parents:
- 75a4ded
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/conccalc.f90
r278f4ed r6473ad3 385 385 griduncn(ix,jyp,kz,ks,nrelpointer,nclass(i),nage)= & 386 386 griduncn(ix,jyp,kz,ks,nrelpointer,nclass(i),nage)+ & 387 xmass1(i,ks)/rhoi*weight*w* xscav_frac1(i,ks)387 xmass1(i,ks)/rhoi*weight*w*max(xscav_frac1(i,ks),0.0) 388 388 end do 389 389 else … … 405 405 griduncn(ixp,jyp,kz,ks,nrelpointer,nclass(i),nage)= & 406 406 griduncn(ixp,jyp,kz,ks,nrelpointer,nclass(i),nage)+ & 407 xmass1(i,ks)/rhoi*weight*w* xscav_frac1(i,ks)407 xmass1(i,ks)/rhoi*weight*w*max(xscav_frac1(i,ks),0.0) 408 408 end do 409 409 else … … 422 422 griduncn(ixp,jy,kz,ks,nrelpointer,nclass(i),nage)= & 423 423 griduncn(ixp,jy,kz,ks,nrelpointer,nclass(i),nage)+ & 424 xmass1(i,ks)/rhoi*weight*w* xscav_frac1(i,ks)424 xmass1(i,ks)/rhoi*weight*w*max(xscav_frac1(i,ks),0.0) 425 425 end do 426 426 else
Note: See TracChangeset
for help on using the changeset viewer.