modified: src/post.f90

	modified:   src/pressure.f90

src/post.f90

@@ -116,15 +116,15 @@
         open(10,file=filename2)
         write(10,'(a,a)') '# ',version
         write(10,'(a,i10,i14,15(i10))') '#',(m,m=1,16)
-        write(10,'(a,a10,a14,15(a10))') &
-        '#','nt','t','dt','ucl','av(u)','av(v)','av(w)', &
+        write(10,'(a,a10,a14,14(a10))') &
+        '#','nt','t','dt','ucl','av(u)','av(v)','av(w)','rms(v)','rms(w)', &
         'Rec','Reb','tau_w','u_tau','Re_tau', &
-        'delta_v','gradu','av(qp)'
+        'delta_v','gradu'
       endif
 
-      write(10,'(x,i10,es14.6e2,21(es10.2e2))') &
-      nt,t,dt,ucl,avu,avv,avw,Rec,Reb, &
-      tau_w,u_tau,Re_tau,delta_v,gradu,avqp
+      write(10,'(x,i10,es14.6e2,20(es10.2e2))') &
+      nt,t,dt,ucl,avu,avv,avw,rmsv,rmsw,Rec,Reb, &
+      tau_w,u_tau,Re_tau,delta_v,gradu
 
 
       close(10)

src/pressure.f90

@@ -28,12 +28,12 @@
       if (celltype(i,j,l).ne.active) cycle
 
       h2sum=0._pr
-      if (celltype(i-1,j,l).ne.wall) h2sum=h2sum+1._pr/(xf(i)-xf(i-1))**2 
-      if (celltype(i+1,j,l).ne.wall) h2sum=h2sum+1._pr/(xf(i)-xf(i-1))**2
-      if (celltype(i,j-1,l).ne.wall) h2sum=h2sum+1._pr/(yf(j)-yf(j-1))**2
+      if (celltype(i+1,j,l).ne.wall) h2sum=h2sum+1._pr/(xf(i)-xf(i-1))**2 
+      if (celltype(i-1,j,l).ne.wall) h2sum=h2sum+1._pr/(xf(i)-xf(i-1))**2
       if (celltype(i,j+1,l).ne.wall) h2sum=h2sum+1._pr/(yf(j)-yf(j-1))**2
-      if (celltype(i,j,l-1).ne.wall) h2sum=h2sum+1._pr/(zf(l)-zf(l-1))**2
+      if (celltype(i,j-1,l).ne.wall) h2sum=h2sum+1._pr/(yf(j)-yf(j-1))**2
       if (celltype(i,j,l+1).ne.wall) h2sum=h2sum+1._pr/(zf(l)-zf(l-1))**2
+      if (celltype(i,j,l-1).ne.wall) h2sum=h2sum+1._pr/(zf(l)-zf(l-1))**2
                
 
       call mass2(ra,i,j,l)