Changeset 1c0d5e6 in flexpart.git for src/wetdepokernel.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/wetdepokernel.f90

    re200b7a r1c0d5e6  
    7979  do ks=1,nspec
    8080
    81   if ((ix.ge.0).and.(jy.ge.0).and.(ix.le.numxgrid-1).and. &
     81  if (.not.usekernel) then
     82      wetgridunc(ix,jy,ks,kp,nunc,nage)= &
     83           wetgridunc(ix,jy,ks,kp,nunc,nage)+deposit(ks)
     84  else
     85    if ((ix.ge.0).and.(jy.ge.0).and.(ix.le.numxgrid-1).and. &
    8286       (jy.le.numygrid-1)) then
    83     w=wx*wy
     87      w=wx*wy
    8488      wetgridunc(ix,jy,ks,kp,nunc,nage)= &
    8589           wetgridunc(ix,jy,ks,kp,nunc,nage)+deposit(ks)*w
    86   endif
     90    endif
    8791
    88   if ((ixp.ge.0).and.(jyp.ge.0).and.(ixp.le.numxgrid-1).and. &
     92    if ((ixp.ge.0).and.(jyp.ge.0).and.(ixp.le.numxgrid-1).and. &
    8993       (jyp.le.numygrid-1)) then
    90     w=(1.-wx)*(1.-wy)
     94      w=(1.-wx)*(1.-wy)
    9195      wetgridunc(ixp,jyp,ks,kp,nunc,nage)= &
    9296           wetgridunc(ixp,jyp,ks,kp,nunc,nage)+deposit(ks)*w
    93   endif
     97    endif
    9498
    95   if ((ixp.ge.0).and.(jy.ge.0).and.(ixp.le.numxgrid-1).and. &
     99    if ((ixp.ge.0).and.(jy.ge.0).and.(ixp.le.numxgrid-1).and. &
    96100       (jy.le.numygrid-1)) then
    97     w=(1.-wx)*wy
     101      w=(1.-wx)*wy
    98102      wetgridunc(ixp,jy,ks,kp,nunc,nage)= &
    99103           wetgridunc(ixp,jy,ks,kp,nunc,nage)+deposit(ks)*w
    100   endif
     104    endif
    101105
    102   if ((ix.ge.0).and.(jyp.ge.0).and.(ix.le.numxgrid-1).and. &
     106    if ((ix.ge.0).and.(jyp.ge.0).and.(ix.le.numxgrid-1).and. &
    103107       (jyp.le.numygrid-1)) then
    104     w=wx*(1.-wy)
     108      w=wx*(1.-wy)
    105109      wetgridunc(ix,jyp,ks,kp,nunc,nage)= &
    106110           wetgridunc(ix,jyp,ks,kp,nunc,nage)+deposit(ks)*w
    107   endif
     111    endif
     112  endif
    108113  end do
    109114
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG