modified: src/post.f90

src/post.f90

@@ -33,8 +33,8 @@
       use var
       real(kind=pr) :: ucl,avu,avv,avw,tke,gradu,graduy,graduz, &
                        Rec,Reb,avvp,avwp,avyp,avzp,avvf,avwf,avxp,avqp, &
-                       syncAv,syncSum
-      integer :: i,j,l,m
+                       syncAv,syncSum,C01,A01
+      integer :: i,j,l,m,N01,syncSumI
       logical :: file_ex
       character*70 :: filename2
 
@@ -81,11 +81,24 @@
        avyp=syncSum(sum(yp(1:np)))/npTot 
        avzp=syncSum(sum(zp(1:np)))/npTot 
        avqp=syncSum(sum(q_el(1:np)))/npTot
+       if (bcy.eq.'w'.and.bcz.eq.'w') then
+        N01=count((min((yp(1:np)-ymin),(ymax-yp(1:np))).le.(dimy*0.01_pr)).or. &
+                  (min((zp(1:np)-zmin),(zmax-zp(1:np))).le.(dimz*0.01_pr)))
+        A01=(dimy*dimz)-(dimy*0.98_pr)*(dimz*0.98_pr)
+       elseif (bcy.eq.'w'.and.bcz.eq.'p') then
+        N01=count(min((yp(1:np)-ymin),(ymax-yp(1:np))).le.(dimy*0.01_pr))
+        A01=(dimy*dimz)-(dimy*0.98_pr)*dimz
+       elseif (bcy.eq.'p'.and.bcz.eq.'w') then
+        N01=count(min((zp(1:np)-zmin),(zmax-zp(1:np))).le.(dimz*0.01_pr))
+        A01=(dimy*dimz)-dimy*(dimz*0.98_pr)
+       endif
+       C01=syncSumI(N01)/(A01*dimxTot)/pnd
       else
        npTot=0
        avyp=0._pr
        avzp=0._pr
        avqp=0._pr
+       C01=0._pr
       endif
 
 
@@ -119,12 +132,12 @@
         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','av(qp)'
+        'C01','av(qp)'
       endif
 
       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,avqp
+      tau_w,u_tau,Re_tau,C01,avqp
 
 
       close(10)