Changeset 5f9d14a in flexpart.git for src/concoutput_nest.f90
- Timestamp:
- Apr 8, 2015, 2:23:27 PM (9 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
- Children:
- 1585284
- Parents:
- cd85138
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/concoutput_nest.f90
re200b7a r5f9d14a 95 95 character :: adate*8,atime*6 96 96 character(len=3) :: anspec 97 integer :: mind 98 ! mind eso:added to ensure identical results between 2&3-fields versions 97 99 98 100 … … 134 136 !******************************************************************* 135 137 138 mind=memind(2) 136 139 do kz=1,numzgrid 137 140 if (kz.eq.1) then … … 156 159 iix=max(min(nint(xl),nxmin1),0) 157 160 jjy=max(min(nint(yl),nymin1),0) 158 densityoutgrid(ix,jy,kz)=(rho(iix,jjy,kzz,2)*dz1+ & 159 rho(iix,jjy,kzz-1,2)*dz2)/dz 161 ! densityoutgrid(ix,jy,kz)=(rho(iix,jjy,kzz,2)*dz1+ & 162 ! rho(iix,jjy,kzz-1,2)*dz2)/dz 163 densityoutgrid(ix,jy,kz)=(rho(iix,jjy,kzz,mind)*dz1+ & 164 rho(iix,jjy,kzz-1,mind)*dz2)/dz 160 165 end do 161 166 end do 162 167 end do 163 168 164 do i=1,numreceptor 165 xl=xreceptor(i) 166 yl=yreceptor(i) 167 iix=max(min(nint(xl),nxmin1),0) 168 jjy=max(min(nint(yl),nymin1),0) 169 densityoutrecept(i)=rho(iix,jjy,1,2) 170 end do 169 do i=1,numreceptor 170 xl=xreceptor(i) 171 yl=yreceptor(i) 172 iix=max(min(nint(xl),nxmin1),0) 173 jjy=max(min(nint(yl),nymin1),0) 174 !densityoutrecept(i)=rho(iix,jjy,1,2) 175 densityoutrecept(i)=rho(iix,jjy,1,mind) 176 end do 171 177 172 178
Note: See TracChangeset
for help on using the changeset viewer.