[Fortran] 纯文本查看 复制代码
subroutine divergent_point_eprgrah(root,n_root,ww_min,ww_max,e1,e2,qq,lmd_k,lmd_kq,eef,kt,uu_d,ss,vvo_q,gdelt,layers) !e1:横坐标精度 e2:纵坐标精度
use number
implicit none
real(8)::qq,lmd_k,lmd_kq,eef,kt,uu_d,ss,vvo_q,gdelt
real(8)::eprgrah,epigrah
real(8)::e1,e2,ww_max,ww_min
real(8)::cint,dint,fcint,fdint,a,b,c,fa,fb,fc !自带参数
real(8)::time_dep,layers,ww0
real(8),dimension(100)::root !存放根
integer::i,j,k,n,sign,n_root
eprgrah=0.0d0
epigrah=0.0d0
n_root=0
sign=0
cint=ww_min
dint=cint+1.0d-5/c_ev
ww0=500.0d0*0.02d0*0.001d0
dint=ww0/c_ev
call CPU_TIME(time_dep)
write(*,*) 'time_dep1',time_dep
call epgrah(fdint,epigrah,dint,qq,gdelt,lmd_k,lmd_kq,eef,kt,uu_d,ss,vvo_q,layers)
call CPU_TIME(time_dep)
write(*,*) 'time_dep2',time_dep
a=cint
! fa=fcint
fb=fdint
sign=1
end