Validated the trimmed curvature computations, reset mindist=3.0

This commit is contained in:
James McClure
2014-04-21 20:45:22 -04:00
parent 07fe4f65e9
commit 42884c5ae7

View File

@@ -1065,7 +1065,7 @@ int main(int argc, char **argv)
// Averaging variables
//...........................................................................
// local averages (to each MPI process)
double trimdist=0.0; // pixel distance to trim surface for specified averages
double trimdist=3.0; // pixel distance to trim surface for specified averages
double awn,ans,aws,lwns,nwp_volume;
double As;
double vol_w, vol_n; // volumes the exclude the interfacial region
@@ -2316,12 +2316,12 @@ int main(int argc, char **argv)
van_global(1) = van_global(1) / vol_n_global;
van_global(2) = van_global(2) / vol_n_global;
if(rank==0){
printf("awn_global = %f \n",awn_global);
printf("trawn_global = %f \n",trawn_global);
printf("Jwn_global = %f \n",Jwn_global);
printf("trJwn_global = %f \n",trJwn_global);
}
// if(rank==0){
// printf("awn_global = %f \n",awn_global);
// printf("trawn_global = %f \n",trawn_global);
// printf("Jwn_global = %f \n",Jwn_global);
// printf("trJwn_global = %f \n",trJwn_global);
// }
// Normalize surface averages by the interfacial area