copy

! using coefficient of fiting
open(14,file="output/"//filename(1:20)//"_output_fit_coefficient.dat",status="OLD")
read(14,*)ncoeffi
print*,"ncoeffi=",ncoeffi
do icoeffi=1,ncoeffi
read(14,*)
read(14,*) norder
print*,"norder=",norder
read(14,*) dmin, dmax
print*,"dmin=",dmin
do iorder=0,norder
read(14,*)cvars,cvars,coeffi(iorder)
print*,"coeffi=",iorder,coeffi(iorder)
enddo
if(icoeffi==1)then
open(20,file="output/"//filename(1:20)//"_output_dcpdn.dat")
write(20,"(A)")"# dis(mm) dcpdn"
endif
do i=1,np
if(d(i)>=dmin .and. d(i) dcpdn=0
do iorder=1,norder
dcpdn=dcpdn+iorder*coeffi(iorder)*d(i)**(iorder-1)
enddo
write(20,"(2(E19.10e3))")d(i),dcpdn
endif
enddo
if(icoeffi==ncoeffi)then
close(20)
endif
enddo
close(14)