Changeset 6985a98 in flexpart.git for src/wetdepokernel.f90
- Timestamp:
- May 8, 2017, 5:28:32 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:
- 2bec33e
- Parents:
- d9f0585 (diff), d1a8707 (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent. - File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/wetdepokernel.f90
r4c64400 r6985a98 95 95 do ks=1,nspec 96 96 97 if ((ix.ge.0).and.(jy.ge.0).and.(ix.le.numxgrid-1).and. & 97 if (.not.usekernel) then 98 wetgridunc(ix,jy,ks,kp,nunc,nage)= & 99 wetgridunc(ix,jy,ks,kp,nunc,nage)+deposit(ks) 100 else 101 if ((ix.ge.0).and.(jy.ge.0).and.(ix.le.numxgrid-1).and. & 98 102 (jy.le.numygrid-1)) then 99 w=wx*wy103 w=wx*wy 100 104 wetgridunc(ix,jy,ks,kp,nunc,nage)= & 101 105 wetgridunc(ix,jy,ks,kp,nunc,nage)+deposit(ks)*w 102 endif106 endif 103 107 104 if ((ixp.ge.0).and.(jyp.ge.0).and.(ixp.le.numxgrid-1).and. &108 if ((ixp.ge.0).and.(jyp.ge.0).and.(ixp.le.numxgrid-1).and. & 105 109 (jyp.le.numygrid-1)) then 106 w=(1.-wx)*(1.-wy)110 w=(1.-wx)*(1.-wy) 107 111 wetgridunc(ixp,jyp,ks,kp,nunc,nage)= & 108 112 wetgridunc(ixp,jyp,ks,kp,nunc,nage)+deposit(ks)*w 109 endif113 endif 110 114 111 if ((ixp.ge.0).and.(jy.ge.0).and.(ixp.le.numxgrid-1).and. &115 if ((ixp.ge.0).and.(jy.ge.0).and.(ixp.le.numxgrid-1).and. & 112 116 (jy.le.numygrid-1)) then 113 w=(1.-wx)*wy117 w=(1.-wx)*wy 114 118 wetgridunc(ixp,jy,ks,kp,nunc,nage)= & 115 119 wetgridunc(ixp,jy,ks,kp,nunc,nage)+deposit(ks)*w 116 endif120 endif 117 121 118 if ((ix.ge.0).and.(jyp.ge.0).and.(ix.le.numxgrid-1).and. &122 if ((ix.ge.0).and.(jyp.ge.0).and.(ix.le.numxgrid-1).and. & 119 123 (jyp.le.numygrid-1)) then 120 w=wx*(1.-wy)124 w=wx*(1.-wy) 121 125 wetgridunc(ix,jyp,ks,kp,nunc,nage)= & 122 126 wetgridunc(ix,jyp,ks,kp,nunc,nage)+deposit(ks)*w 123 endif 127 endif 128 endif 124 129 end do 125 130 end if
Note: See TracChangeset
for help on using the changeset viewer.