Skip to content
Snippets Groups Projects
Commit 46880916 authored by holger's avatar holger
Browse files

modified: src/momentum.f90

parent 8e77bb23
No related branches found
No related tags found
No related merge requests found
...@@ -106,7 +106,7 @@ ...@@ -106,7 +106,7 @@
if (celltype(i,j,l).ne.active) cycle if (celltype(i,j,l).ne.active) cycle
if (((bcx.eq.'p').or.((bcx.eq.'w').and.(i-2.ge.imin-1).and.(i+1.le.imax))).and. & if (((bcx.eq.'p'.or.bcx.eq.'i').or.((bcx.eq.'w').and.(i-2.ge.imin-1).and.(i+1.le.imax))).and. &
((bcy.eq.'p').or.((bcy.eq.'w').and.(j-2.ge.jmin-1).and.(j+1.le.jmax))).and. & ((bcy.eq.'p').or.((bcy.eq.'w').and.(j-2.ge.jmin-1).and.(j+1.le.jmax))).and. &
((bcz.eq.'p').or.((bcz.eq.'w').and.(l-2.ge.lmin-1).and.(l+1.le.lmax)))) then ((bcz.eq.'p').or.((bcz.eq.'w').and.(l-2.ge.lmin-1).and.(l+1.le.lmax)))) then
call mass2(ra,i,j,l) call mass2(ra,i,j,l)
...@@ -116,7 +116,7 @@ ...@@ -116,7 +116,7 @@
defect_c(i,j,l)=0._pr defect_c(i,j,l)=0._pr
endif endif
if (((bcx.eq.'p').or.((bcx.eq.'w').and.(i-3.ge.imin).and.(i+3.le.imax-1))).and. & if (((bcx.eq.'p'.or.bcx.eq.'i').or.((bcx.eq.'w').and.(i-3.ge.imin).and.(i+3.le.imax-1))).and. &
((bcy.eq.'p').or.((bcy.eq.'w').and.(j-3.ge.jmin).and.(j+3.le.jmax))).and. & ((bcy.eq.'p').or.((bcy.eq.'w').and.(j-3.ge.jmin).and.(j+3.le.jmax))).and. &
((bcz.eq.'p').or.((bcz.eq.'w').and.(l-3.ge.lmin).and.(l+3.le.lmax)))) then ((bcz.eq.'p').or.((bcz.eq.'w').and.(l-3.ge.lmin).and.(l+3.le.lmax)))) then
call momx1(tu,qu,i,j,l) call momx1(tu,qu,i,j,l)
...@@ -127,7 +127,7 @@ ...@@ -127,7 +127,7 @@
defect_u(i,j,l)=0._pr defect_u(i,j,l)=0._pr
endif endif
if (((bcx.eq.'p').or.((bcx.eq.'w').and.(i-3.ge.imin).and.(i+3.le.imax))).and. & if (((bcx.eq.'p'.or.bcx.eq.'i').or.((bcx.eq.'w').and.(i-3.ge.imin).and.(i+3.le.imax))).and. &
((bcy.eq.'p').or.((bcy.eq.'w').and.(j-3.ge.jmin).and.(j+3.le.jmax-1))).and. & ((bcy.eq.'p').or.((bcy.eq.'w').and.(j-3.ge.jmin).and.(j+3.le.jmax-1))).and. &
((bcz.eq.'p').or.((bcz.eq.'w').and.(l-3.ge.lmin).and.(l+3.le.lmax)))) then ((bcz.eq.'p').or.((bcz.eq.'w').and.(l-3.ge.lmin).and.(l+3.le.lmax)))) then
call momy1(tv,qv,i,j,l) call momy1(tv,qv,i,j,l)
...@@ -138,7 +138,7 @@ ...@@ -138,7 +138,7 @@
defect_v(i,j,l)=0._pr defect_v(i,j,l)=0._pr
endif endif
if (((bcx.eq.'p').or.((bcx.eq.'w').and.(i-3.ge.imin).and.(i+3.le.imax))).and. & if (((bcx.eq.'p'.or.bcx.eq.'i').or.((bcx.eq.'w').and.(i-3.ge.imin).and.(i+3.le.imax))).and. &
((bcy.eq.'p').or.((bcy.eq.'w').and.(j-3.ge.jmin).and.(j+3.le.jmax))).and. & ((bcy.eq.'p').or.((bcy.eq.'w').and.(j-3.ge.jmin).and.(j+3.le.jmax))).and. &
((bcz.eq.'p').or.((bcz.eq.'w').and.(l-3.ge.lmin).and.(l+3.le.lmax-1)))) then ((bcz.eq.'p').or.((bcz.eq.'w').and.(l-3.ge.lmin).and.(l+3.le.lmax-1)))) then
call momz1(tw,qw,i,j,l) call momz1(tw,qw,i,j,l)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment