Changeset 3c1da52 in flexpart.git for src/wetdepokernel_nest.f90
- Timestamp:
- Jul 15, 2015, 9:45:31 AM (9 years ago)
- Branches:
- svn-petra
- Children:
- c07b09e
- Parents:
- 31674b5
- git-author:
- Ignacio Pisso <Ignacio.Pisso@…> (07/15/15 09:23:56)
- git-committer:
- Ignacio Pisso <Ignacio.Pisso@…> (07/15/15 09:45:31)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/wetdepokernel_nest.f90
re200b7a r3c1da52 1 1 !********************************************************************** 2 ! Copyright 1998 ,1999,2000,2001,2002,2005,2007,2008,2009,2010*2 ! Copyright 1998-2015 * 3 3 ! Andreas Stohl, Petra Seibert, A. Frank, Gerhard Wotawa, * 4 4 ! Caroline Forster, Sabine Eckhardt, John Burkhart, Harald Sodemann * … … 20 20 !********************************************************************** 21 21 22 subroutine wetdepokernel_nest & 23 (nunc,deposit,x,y,nage,kp) 24 ! i i i i i i 22 subroutine wetdepokernel_nest (nunc,deposit,x,y,itage,nage,kp) 23 ! i i i i i i i 25 24 !***************************************************************************** 26 25 ! * … … 35 34 ! 2 September 2004: Adaptation from wetdepokernel. * 36 35 ! * 36 ! 37 ! PS, 2/2015: do not use kernel for itage < 3 h 38 ! same as for concentration in conccalc.f90 37 39 ! * 38 40 !***************************************************************************** … … 53 55 54 56 real :: x,y,deposit(maxspec),ddx,ddy,xl,yl,wx,wy,w 55 integer :: ix,jy,ixp,jyp,ks,kp,nunc,nage 57 integer :: ix,jy,ixp,jyp,ks,kp,nunc,nage,itage 56 58 57 59 … … 80 82 endif 81 83 84 if (itage .lt. itagekernmin) then 85 ! no kernel, direct attribution to grid cell 86 87 do ks=1,nspec 88 if (ix.ge.0 .and. jy.ge.0 .and. & 89 ix.le.numxgridn-1 .and. jy.le.numygridn-1) & 90 wetgriduncn(ix,jy,ks,kp,nunc,nage)= & 91 wetgriduncn(ix,jy,ks,kp,nunc,nage)+deposit(ks) 92 enddo 93 94 else 95 ! Determine mass fractions for four grid points & distribute 82 96 83 ! Determine mass fractions for four grid points 84 !********************************************** 97 do ks=1,nspec 85 98 86 do ks=1,nspec 99 if ((ix.ge.0).and.(jy.ge.0).and.(ix.le.numxgridn-1).and. & 100 (jy.le.numygridn-1)) then 101 w=wx*wy 102 wetgriduncn(ix,jy,ks,kp,nunc,nage)= & 103 wetgriduncn(ix,jy,ks,kp,nunc,nage)+deposit(ks)*w 104 endif 87 105 88 if ((ix.ge.0).and.(jy.ge.0).and.(ix.le.numxgridn-1).and. & 89 (jy.le.numygridn-1)) then 90 w=wx*wy 91 wetgriduncn(ix,jy,ks,kp,nunc,nage)= & 92 wetgriduncn(ix,jy,ks,kp,nunc,nage)+deposit(ks)*w 106 if ((ixp.ge.0).and.(jyp.ge.0).and.(ixp.le.numxgridn-1).and. & 107 (jyp.le.numygridn-1)) then 108 w=(1.-wx)*(1.-wy) 109 wetgriduncn(ixp,jyp,ks,kp,nunc,nage)= & 110 wetgriduncn(ixp,jyp,ks,kp,nunc,nage)+deposit(ks)*w 111 endif 112 113 if ((ixp.ge.0).and.(jy.ge.0).and.(ixp.le.numxgridn-1).and. & 114 (jy.le.numygridn-1)) then 115 w=(1.-wx)*wy 116 wetgriduncn(ixp,jy,ks,kp,nunc,nage)= & 117 wetgriduncn(ixp,jy,ks,kp,nunc,nage)+deposit(ks)*w 118 endif 119 120 if ((ix.ge.0).and.(jyp.ge.0).and.(ix.le.numxgridn-1).and. & 121 (jyp.le.numygridn-1)) then 122 w=wx*(1.-wy) 123 wetgriduncn(ix,jyp,ks,kp,nunc,nage)= & 124 wetgriduncn(ix,jyp,ks,kp,nunc,nage)+deposit(ks)*w 125 endif 126 127 end do 128 93 129 endif 94 95 if ((ixp.ge.0).and.(jyp.ge.0).and.(ixp.le.numxgridn-1).and. &96 (jyp.le.numygridn-1)) then97 w=(1.-wx)*(1.-wy)98 wetgriduncn(ixp,jyp,ks,kp,nunc,nage)= &99 wetgriduncn(ixp,jyp,ks,kp,nunc,nage)+deposit(ks)*w100 endif101 102 if ((ixp.ge.0).and.(jy.ge.0).and.(ixp.le.numxgridn-1).and. &103 (jy.le.numygridn-1)) then104 w=(1.-wx)*wy105 wetgriduncn(ixp,jy,ks,kp,nunc,nage)= &106 wetgriduncn(ixp,jy,ks,kp,nunc,nage)+deposit(ks)*w107 endif108 109 if ((ix.ge.0).and.(jyp.ge.0).and.(ix.le.numxgridn-1).and. &110 (jyp.le.numygridn-1)) then111 w=wx*(1.-wy)112 wetgriduncn(ix,jyp,ks,kp,nunc,nage)= &113 wetgriduncn(ix,jyp,ks,kp,nunc,nage)+deposit(ks)*w114 endif115 116 end do117 130 end subroutine wetdepokernel_nest
Note: See TracChangeset
for help on using the changeset viewer.