# awk script 
#
# the begin command is executed BEFORE anything else is done
BEGIN {FS=" ";
"date" | getline date;
print date
}
#
{N=N+1}
{if(N==1) {printf "Input file [ %12s ]\n",FILENAME ;
print "#----------------------Input---------------------#"}}
{if($1=="Length") {lgth=$2;print$0}}
{if($1=="Diameter") {diam=$2;print$0}}
{if($1=="Flow") {flow=$2;print$0}}
{if($1=="Rgas") {rgas=$2;print$0}}
{if($1=="kgas") {kgas=$2;print$0}}
{if($1=="Cpgas") {cpgas=$2*1000.;print$0}}
{if($1=="mugas") {mugas=$2;print$0}}
{if($1=="Tinlet") {tinlet=$2;print$0}}
{if($1=="Pinlet") {pinlet=$2;print$0}}
{if($1=="Tamb") {Tamb=$2;print$0}}
#
END {
# input dimension are mm's 
{gasdens=pinlet/rgas/(tinlet+273.15)}
{prndlt=cpgas*mugas/kgas}
# (note: diam in input in mm's)
{rend=4.*1000.*flow/3.14159/diam/mugas}
{n=.4}
{if(Tamb<tinlet) {n=.3}}
{hinf=.023*(rend^.8)*(prndlt^n)*kgas/(diam/1000.)}
# if flow is laminar.....
{if(rend<2300) {hinf=4.36*kgas/(diam/1000.)}}
{vmean=1000.*1000.*flow/gasdens/.7854/diam/diam}
# entrance effect using 1.0+1.2*(D/X)^m
{m=1.5}
{x1=lgth/20.;hx1=hinf*(1.0+1.2*(diam/x1)^m)}
{x2=x1+lgth/20.;hx2=hinf*(1.0+1.2*(diam/x2)^m)}
{x3=x2+lgth/20.;hx3=hinf*(1.0+1.2*(diam/x3)^m)}
{x4=x3+lgth/20.;hx4=hinf*(1.0+1.2*(diam/x4)^m)}
{x5=x4+lgth/20.;hx5=hinf*(1.0+1.2*(diam/x5)^m)}
{x6=x5+lgth/20.;hx6=hinf*(1.0+1.2*(diam/x6)^m)}
{x7=x6+lgth/20.;hx7=hinf*(1.0+1.2*(diam/x7)^m)}
{x8=x7+lgth/20.;hx8=hinf*(1.0+1.2*(diam/x8)^m)}
{x9=x8+lgth/20.;hx9=hinf*(1.0+1.2*(diam/x9)^m)}
{x10=x9+lgth/20.;hx10=hinf*(1.0+1.2*(diam/x10)^m)}
{x11=x10+lgth/20.;hx11=hinf*(1.0+1.2*(diam/x11)^m)}
{x12=x11+lgth/20.;hx12=hinf*(1.0+1.2*(diam/x12)^m)}
{x13=x12+lgth/20.;hx13=hinf*(1.0+1.2*(diam/x13)^m)}
{x14=x13+lgth/20.;hx14=hinf*(1.0+1.2*(diam/x14)^m)}
{x15=x14+lgth/20.;hx15=hinf*(1.0+1.2*(diam/x15)^m)}
{x16=x15+lgth/20.;hx16=hinf*(1.0+1.2*(diam/x16)^m)}
{x17=x16+lgth/20.;hx17=hinf*(1.0+1.2*(diam/x17)^m)}
{x18=x17+lgth/20.;hx18=hinf*(1.0+1.2*(diam/x18)^m)}
{x19=x18+lgth/20.;hx19=hinf*(1.0+1.2*(diam/x19)^m)}
{x20=x19+lgth/20.;hx20=hinf*(1.0+1.2*(diam/x20)^m)}
{havg=(hx1+hx2+hx3+hx4+hx5+hx6+hx7+hx8+hx9+hx10+hx11+hx12+hx13+hx14+hx15)}
{havg=(havg+hx16+hx17+hx18+hx19+hx20)/20.}
{print" "}
#
{print "#----------------------Output--------------------#"}
{printf("%-12s %12.3f\n","Gas Density", gasdens)}
{printf("%-12s %12.3f\n","Prandlt No.",prndlt)}
{printf("%-12s %12.0f\n","Reyn(D)",rend)}
{printf("%-12s %12.1f\n","Velocity",vmean)}
{printf("%-12s %12.1f\n","Hinfinity",hinf)}
{printf("%-12s %12.1f\n","Haverage",havg)}
{print" "}
{print "Entrance effect included: "}
{printf("%-10s %-10s\n","X(mm)","h(local)")}
{printf("%-10.2f %-10.1f\n",x1,hx1)}
{printf("%-10.2f %-10.1f\n",x2,hx2)}
{printf("%-10.2f %-10.1f\n",x3,hx3)}
{printf("%-10.2f %-10.1f\n",x4,hx4)}
{printf("%-10.2f %-10.1f\n",x5,hx5)}
{printf("%-10.2f %-10.1f\n",x6,hx6)}
{printf("%-10.2f %-10.1f\n",x7,hx7)}
{printf("%-10.2f %-10.1f\n",x8,hx8)}
{printf("%-10.2f %-10.1f\n",x9,hx9)}
{printf("%-10.2f %-10.1f\n",x10,hx10)}
{printf("%-10.2f %-10.1f\n",x11,hx11)}
{printf("%-10.2f %-10.1f\n",x12,hx12)}
{printf("%-10.2f %-10.1f\n",x13,hx13)}
{printf("%-10.2f %-10.1f\n",x14,hx14)}
{printf("%-10.2f %-10.1f\n",x15,hx15)}
{printf("%-10.2f %-10.1f\n",x16,hx16)}
{printf("%-10.2f %-10.1f\n",x17,hx17)}
{printf("%-10.2f %-10.1f\n",x18,hx18)}
{printf("%-10.2f %-10.1f\n",x19,hx19)}
{printf("%-10.2f %-10.1f\n",x20,hx20)}
#
{print " "}
{print "#.......................End of File..............#"}
}