boundcond_domainfill Subroutine

public subroutine boundcond_domainfill(itime, loutend)

                                                                        *

Particles are created by this subroutine continuously throughout the * simulation at the boundaries of the domain-filling box. * All particles carry the same amount of mass which alltogether comprises the* mass of air within the box, which remains (more or less) constant. * * Author: A. Stohl * * 16 October 2002 * * Changes * 2022, L. Bakels: OpenMP parallelisation * *

                                                                        *

Variables: * * nx_we(2) grid indices for western and eastern boundary of domain- * filling trajectory calculations * ny_sn(2) grid indices for southern and northern boundary of domain- * filling trajectory calculations * *

Western and eastern boundary condition

OMP doesn't work here yet because particles are being spawned inside of the loop !$OMP PARALLEL PRIVATE(i,jy,k,j,ii,deltaz,boundarea,indz,indzp,indexh,windl,rhol, & !$OMP windhl,rhohl,windx,rhox,fluxofmass,mmass,ixm,jym,ixp,jyp,ddx,ddy,rddx, & !$OMP rddy,p1,p2,p3,p4,indzm,mm,indzh,pvpart,ylat,ix,cosfact,ipart) & !$OMP REDUCTION(+:numactiveparticles,numparticlecount_tmp,accmasst)

Southern and northern boundary condition

!$OMP DO

Arguments