program f90_MANNING_CIRC implicit none real(8),parameter::pi=3.14159265358979323846D0 real(8)::nn,DD,Q,ii real(8)::v,R,A,S,theta,h,Qmax,hm character(len=50)::dummy integer j call getarg(1,dummy);read(dummy,*) DD call getarg(2,dummy);read(dummy,*) ii call getarg(3,dummy);read(dummy,*) nn write(6,'(i4,3f10.5)') 0,DD,ii,nn Qmax=0.0D0 hm=0.0D0 do j=1,180 theta=dble(j)/180.0D0*pi A=0.25*DD*DD*(theta-sin(theta)*cos(theta)) R=0.25*DD*(1.0D0-sin(theta)*cos(theta)/theta) v=1.0D0/nn*R**(2.0D0/3.0D0)*sqrt(ii) Q=A*v h=0.5D0*DD*(1.0D0-cos(theta)) write(6,'(i4,4f10.5)') j,h/DD,h,Q,v if(Qmax