Changeset 5f9d14a in flexpart.git for src/concoutput_surf_nest_mpi.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_surf_nest_mpi.f90
- Property mode changed from 100755 to 100644
r8a65cb0 r5f9d14a 98 98 character :: adate*8,atime*6 99 99 character(len=3) :: anspec 100 integer :: mind 101 ! mind eso:added to get consistent results between 2&3-fields versions 100 102 101 103 ! Measure execution time … … 148 150 !******************************************************************* 149 151 152 mind=memind(2) 150 153 do kz=1,numzgrid 151 154 if (kz.eq.1) then … … 170 173 iix=max(min(nint(xl),nxmin1),0) 171 174 jjy=max(min(nint(yl),nymin1),0) 172 densityoutgrid(ix,jy,kz)=(rho(iix,jjy,kzz,2)*dz1+ & 173 rho(iix,jjy,kzz-1,2)*dz2)/dz 175 ! densityoutgrid(ix,jy,kz)=(rho(iix,jjy,kzz,2)*dz1+ & 176 ! rho(iix,jjy,kzz-1,2)*dz2)/dz 177 densityoutgrid(ix,jy,kz)=(rho(iix,jjy,kzz,mind)*dz1+ & 178 rho(iix,jjy,kzz-1,mind)*dz2)/dz 174 179 end do 175 180 end do 176 181 end do 177 182 178 do i=1,numreceptor 179 xl=xreceptor(i) 180 yl=yreceptor(i) 181 iix=max(min(nint(xl),nxmin1),0) 182 jjy=max(min(nint(yl),nymin1),0) 183 densityoutrecept(i)=rho(iix,jjy,1,2) 184 end do 183 do i=1,numreceptor 184 xl=xreceptor(i) 185 yl=yreceptor(i) 186 iix=max(min(nint(xl),nxmin1),0) 187 jjy=max(min(nint(yl),nymin1),0) 188 !densityoutrecept(i)=rho(iix,jjy,1,2) 189 densityoutrecept(i)=rho(iix,jjy,1,mind) 190 end do 185 191 186 192
Note: See TracChangeset
for help on using the changeset viewer.