Changeset 92a74b2 in flexpart.git
- Timestamp:
- Nov 21, 2016, 10:05:39 AM (7 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
- Children:
- 3d6060e
- Parents:
- 7062fab
- Location:
- src
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
src/conccalc.f90
rdeaac29 r92a74b2 207 207 endif 208 208 209 else ! attribution via uniform kernel 209 else ! attribution via uniform kernel 210 210 211 211 ddx=xl-real(ix) ! distance to left cell border … … 326 326 if ((itage.lt.10800).or.(xl.lt.0.5).or.(yl.lt.0.5).or. & 327 327 (xl.gt.real(numxgridn-1)-0.5).or. & 328 (yl.gt.real(numygridn-1)-0.5) ) then ! no kernel, direct attribution to grid cell328 (yl.gt.real(numygridn-1)-0.5).or.(usekernel.eq.0)) then ! no kernel, direct attribution to grid cell 329 329 if ((ix.ge.0).and.(jy.ge.0).and.(ix.le.numxgridn-1).and. & 330 330 (jy.le.numygridn-1)) then -
src/get_vdep_prob.f90
r7062fab r92a74b2 20 20 !********************************************************************** 21 21 22 subroutine advance_rec(itime,xt,yt,zt,prob)22 subroutine get_vdep_prob(itime,xt,yt,zt,prob) 23 23 ! i i i i o 24 24 !***************************************************************************** 25 25 ! * 26 ! Calculation of turbulent particle trajectories utilizing a * 27 ! zero-acceleration scheme, which is corrected by a numerically more * 28 ! accurate Petterssen scheme whenever possible. * 26 ! Calculation of the probability for dyr deposition * 29 27 ! * 30 ! Particle positions are read in, incremented, and returned to the calling * 31 ! program. * 32 ! * 33 ! In different regions of the atmosphere (PBL vs. free troposphere), * 34 ! different parameters are needed for advection, parameterizing turbulent * 35 ! velocities, etc. For efficiency, different interpolation routines have * 36 ! been written for these different cases, with the disadvantage that there * 37 ! exist several routines doing almost the same. They all share the * 38 ! included file 'interpol_mod'. The following * 39 ! interpolation routines are used: * 40 ! * 41 ! interpol_all(_nests) interpolates everything (called inside the PBL) * 42 ! interpol_misslev(_nests) if a particle moves vertically in the PBL, * 43 ! additional parameters are interpolated if it * 44 ! crosses a model level * 45 ! interpol_wind(_nests) interpolates the wind and determines the * 46 ! standard deviation of the wind (called outside * 47 ! PBL) also interpolates potential vorticity * 48 ! interpol_wind_short(_nests) only interpolates the wind (needed for the * 49 ! Petterssen scheme) * 50 ! interpol_vdep(_nests) interpolates deposition velocities * 51 ! * 52 ! * 53 ! Author: A. Stohl * 54 ! * 55 ! 16 December 1997 * 56 ! * 57 ! Changes: * 58 ! * 59 ! 8 April 2000: Deep convection parameterization * 60 ! * 61 ! May 2002: Petterssen scheme introduced * 28 ! Particle positions are read in - prob returned * 62 29 ! * 63 30 !***************************************************************************** … … 175 142 176 143 177 end subroutine advance_rec144 end subroutine get_vdep_prob 178 145 -
src/get_wetscav.f90
rc9cf570 r92a74b2 20 20 !********************************************************************** 21 21 22 subroutine get_wetscav(itime,ltsample,loutnext,jpart, grfraction,inc_count,blc_count,wetscav)22 subroutine get_wetscav(itime,ltsample,loutnext,jpart,ks,grfraction,inc_count,blc_count,wetscav) 23 23 ! i i i i o o 24 24 !***************************************************************************** … … 218 218 !********************************************************** 219 219 220 do ks=1,nspec ! loop over species221 220 wetdeposit(ks)=0. 222 221 wetscav=0. … … 338 337 endif !incloud 339 338 340 end do ! loop over species341 339 20 continue 342 340 -
src/makefile
r9669e1e r92a74b2 152 152 erf.o readavailable.o \ 153 153 ew.o readreleases.o \ 154 readdepo.o advance_rec.o \ 154 readdepo.o get_vdep_prob.o \ 155 get_wetscav.o \ 155 156 psim.o outgrid_init.o \ 156 157 outgrid_init_nest.o \ … … 271 272 272 273 ## DEPENDENCIES 273 advance_rec.o: cmapf_mod.o com_mod.o hanna_mod.o interpol_mod.o par_mod.o \274 get_vdep_prob.o: cmapf_mod.o com_mod.o hanna_mod.o interpol_mod.o par_mod.o \ 274 275 point_mod.o random_mod.o 275 276 advance.o: cmapf_mod.o com_mod.o hanna_mod.o interpol_mod.o par_mod.o \ … … 426 427 verttransform_gfs.o: cmapf_mod.o com_mod.o par_mod.o 427 428 verttransform_nests.o: com_mod.o par_mod.o 429 get_wetscav.o: com_mod.o par_mod.o point_mod.o 428 430 wetdepo.o: com_mod.o par_mod.o point_mod.o 429 431 wetdepokernel.o: com_mod.o par_mod.o unc_mod.o -
src/timemanager.f90
rc9cf570 r92a74b2 567 567 do ks=1,nspec 568 568 if ((xscav_frac1(j,ks).lt.0)) then 569 call get_wetscav(itime,lsynctime,loutnext,j, grfraction,idummy,idummy,wetscav)569 call get_wetscav(itime,lsynctime,loutnext,j,ks,grfraction,idummy,idummy,wetscav) 570 570 if (wetscav(ks).gt.0) then 571 xscav_frac1(j,ks)=wetscav(ks)*(zpoint2( j)-zpoint1(j))571 xscav_frac1(j,ks)=wetscav(ks)*(zpoint2(npoint(j))-zpoint1(npoint(j))) 572 572 else 573 573 xmass1(j,ks)=0. -
src/wetdepo.f90
rc9cf570 r92a74b2 103 103 33 continue 104 104 105 do ks=1,nspec ! loop over species 105 106 106 107 !************************************************** 107 108 ! CALCULATE DEPOSITION 108 109 !************************************************** 109 call get_wetscav(itime,ltsample,loutnext,jpart, grfraction,inc_count,blc_count,wetscav)110 call get_wetscav(itime,ltsample,loutnext,jpart,ks,grfraction,inc_count,blc_count,wetscav) 110 111 111 112 if (wetscav.gt.0.) then … … 136 137 endif 137 138 139 end do ! loop over species 140 138 141 ! Sabine Eckhardt, June 2008 create deposition runs only for forward runs 139 142 ! Add the wet deposition to accumulated amount on output grid and nested output grid
Note: See TracChangeset
for help on using the changeset viewer.