Changeset 1c0d5e6 in flexpart.git for src/drydepokernel.f90


Ignore:
Timestamp:
Dec 1, 2016, 12:47:22 PM (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:
56a2e54
Parents:
8ee24a5
Message:

adding logical usekernel var-also for dry and wet

File:
1 edited

Legend:

Unmodified
Added
Removed
  • src/drydepokernel.f90

    re200b7a r1c0d5e6  
    8181    if ((abs(deposit(ks)).gt.0).and.DRYDEPSPEC(ks)) then
    8282
    83    if ((ix.ge.0).and.(jy.ge.0).and.(ix.le.numxgrid-1).and. &
     83    if (.not.usekernel) then
     84       drygridunc(ix,jy,ks,kp,nunc,nage)= &
     85           drygridunc(ix,jy,ks,kp,nunc,nage)+deposit(ks)
     86    else
     87      if ((ix.ge.0).and.(jy.ge.0).and.(ix.le.numxgrid-1).and. &
    8488        (jy.le.numygrid-1)) then
    85     w=wx*wy
     89       w=wx*wy
    8690      drygridunc(ix,jy,ks,kp,nunc,nage)= &
    8791           drygridunc(ix,jy,ks,kp,nunc,nage)+deposit(ks)*w
    88       continue
    89   endif
     92     endif
    9093
    91   if ((ixp.ge.0).and.(jyp.ge.0).and.(ixp.le.numxgrid-1).and. &
     94    if ((ixp.ge.0).and.(jyp.ge.0).and.(ixp.le.numxgrid-1).and. &
    9295       (jyp.le.numygrid-1)) then
    9396    w=(1.-wx)*(1.-wy)
    9497      drygridunc(ixp,jyp,ks,kp,nunc,nage)= &
    9598           drygridunc(ixp,jyp,ks,kp,nunc,nage)+deposit(ks)*w
    96   endif
     99    endif
    97100
    98   if ((ixp.ge.0).and.(jy.ge.0).and.(ixp.le.numxgrid-1).and. &
     101    if ((ixp.ge.0).and.(jy.ge.0).and.(ixp.le.numxgrid-1).and. &
    99102       (jy.le.numygrid-1)) then
    100     w=(1.-wx)*wy
     103      w=(1.-wx)*wy
    101104      drygridunc(ixp,jy,ks,kp,nunc,nage)= &
    102105           drygridunc(ixp,jy,ks,kp,nunc,nage)+deposit(ks)*w
    103   endif
     106    endif
    104107
    105   if ((ix.ge.0).and.(jyp.ge.0).and.(ix.le.numxgrid-1).and. &
     108    if ((ix.ge.0).and.(jyp.ge.0).and.(ix.le.numxgrid-1).and. &
    106109       (jyp.le.numygrid-1)) then
    107     w=wx*(1.-wy)
     110      w=wx*(1.-wy)
    108111      drygridunc(ix,jyp,ks,kp,nunc,nage)= &
    109112           drygridunc(ix,jyp,ks,kp,nunc,nage)+deposit(ks)*w
    110   endif
     113    endif
    111114
    112   endif
     115    endif ! kernel
     116    endif ! deposit>0
    113117
    114118    end do
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG