Changeset 6473ad3 in flexpart.git


Ignore:
Timestamp:
Feb 9, 2017, 5:02:33 PM (2 years ago)
Author:
Sabine <sabine.eckhardt@…>
Branches:
dev, release-10, univie
Children:
33279f7
Parents:
75a4ded
Message:

BUGFIX: add max(xscav1) to all write statements

File:
1 edited

Legend:

Unmodified
Added
Removed
  • src/conccalc.f90

    r278f4ed r6473ad3  
    385385                   griduncn(ix,jyp,kz,ks,nrelpointer,nclass(i),nage)= &
    386386                     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)
    388388                 end do
    389389              else
     
    405405                   griduncn(ixp,jyp,kz,ks,nrelpointer,nclass(i),nage)= &
    406406                     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)
    408408                 end do
    409409              else
     
    422422                   griduncn(ixp,jy,kz,ks,nrelpointer,nclass(i),nage)= &
    423423                     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)
    425425                 end do
    426426              else
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG