49       integer nu,nv,nvv,ngr  local
192       call db_data_get_int(sname,dname,1,ngr,result)
276        call rism(rtdb,rdffile,nu,nv,nvv,ngr,icl,
299        subroutine potcreat(nvv,nu,ngr,icr,kgrid,rgrid,tau,solperm,  argument
302        integer  nu,nv, ngr,nvv 
306        real*8 r1, r2, dr, dk, rgrid(1:ngr),kgrid(1:ngr),pi 
308        real*8 ul(1:nu,1:nvv,1:ngr),plj(1:nu,1:nvv,1:ngr)
324        do i=1,ngr 
370        subroutine ssoz(nvv,nu,ngr,nv,ims,mv,wu,c3,chi,hu)  argument
372        integer  nu, nvv, nv, ngr, i,j1,j,k1,ims(1:nv),mv(1:nv)          local
373        real*8 rgrid(1:ngr),kgrid(1:ngr)  
375        real*8 c3(1:nu,1:nvv,1:ngr),hu(1:nu,1:nvv,1:ngr)
376        real*8 ctt(1:nu,1:nvv,1:ngr),h1(1:nu,1:nvv,1:ngr)
377        real*8 wu(1:nu,1:nu,1:ngr), chi(1:nvv,1:nvv,1:ngr)  
379        do i=1,ngr
394        do i=1,ngr
406        do i=1,ngr
421        subroutine clos(nvv,nu,ngr,t,plj,cr,gam,icl)  argument
423        integer  nu, nvv, ngr, i,j1,j,icl  local
424        real*8 plj(1:nu,1:nvv,1:ngr)          
426        real*8 cr(1:nu,1:nvv,1:ngr),gam(1:nu,1:nvv,1:ngr)           
428        do i=1,ngr
457        subroutine rism(rtdb,rdffile,nu,nv,nvv,ngr,  argument
469        integer  nu, nv,nd,nvv, ngr, i,j1,j,k1,icl,k        
470        real*8 rgrid(1:ngr),kgrid(1:ngr)  
473        real*8 wu(1:nu,1:nu,1:ngr), chi(1:nvv,1:nvv,1:ngr)   
478        real*8 wv(1:nvv,1:nvv,1:ngr)
483        real*8 plj(1:nu,1:nvv,1:ngr), ul(1:nu,1:nvv,1:ngr)
485        real*8 cr(1:nu,1:nvv,1:ngr),c2(1:ngr),tt(1:nu,1:nvv)  
486        real*8 ht(1:ngr),c3(1:nu,1:nvv,1:ngr),cf3(1:nu,1:nvv,1:ngr)
487        real*8 cold(1:nu,1:nvv,1:ngr),cnew(1:nu,1:nvv,1:ngr)
488        real*8 gfold(1:nu,1:nvv,1:ngr),hu(1:nu,1:nvv,1:ngr) 
489        real*8 gold(1:nu,1:nvv,1:ngr), gnew(1:nu,1:nvv,1:ngr)   
490        real*8 hold(1:nu,1:nvv,1:ngr), hnew(1:nu,1:nvv,1:ngr)             
493        real*8 ggo(1:dd,1:nu,1:nvv,1:ngr),dgg(1:dd,1:nu,1:nvv,1:ngr)
495        write(76,*) nu,nv,nvv,ngr,icl,icr,tau,solperm
505        call chicreat(rdffile,nd,nv,ngr,ims,nvv,
515        call potcreat(nvv,nu,ngr,icr,kgrid,rgrid,tau,solperm,
520        do i=1,ngr
534        call clos(nvv,nu,ngr,t,plj,cnew,gold,icl)
537          do i=1,ngr
547          do i=1,ngr-1  
550          call sinft(c2,ngr)
552          do i=1,ngr-1
555          cf3(j,j1,ngr)=cf3(j,j1,ngr-1) 
561        do i=1,ngr
570        call ssoz(nvv,nu,ngr,nv,ims,mv,wu,cf3,chi,hu)
575          do i=1,ngr-1 
578          call sinft(ht,ngr)
579          do i=1,ngr-1
583          gnew(j,j1,ngr)=gnew(j,j1,ngr-1)
584          hnew(j,j1,ngr)=hnew(j,j1,ngr-1) 
598          do i=1,ngr
606        del0=(mmaxi/ngr/nu/nvv)**(1.0/2)
608        call  diis(nu,nvv,ngr,lambd,kd,dd,gold,gnew,ggo,dgg)
612          do i=1,ngr
623        call thermo(nu,nv,nvv,ngr,icl,ims,den,t,rgrid,qqu,qvv,
639        subroutine diis(nu,nvv,ngr,lam,k1,dd,hold,hnew,ggo,dgg)  argument
642        real*8 ggo(1:dd,1:nu,1:nvv,1:ngr),dgg(1:dd,1:nu,1:nvv,1:ngr)
643        real*8 ss(1:dd+1,1:dd+1),hold(1:nu,1:nvv,1:ngr)
644        real*8 s0(1:dd+1),s1,lam,hnew(1:nu,1:nvv,1:ngr)
645        integer j,j1,i,nu,nvv,ngr  local
661          do i=1,ngr     
668              ss(jd,jd1)=(dgg(jd,j,j1,i)*dgg(jd1,j,j1,i))/ngr+ss(jd,jd1)
682         do i=1,ngr
717        subroutine thermo(nu,nv,nvv,ngr,icl,ims,den,t,rgrid,qqu,qvv,  argument
722        integer  nu, nv,nvv, ngr, i,j1,j,icl,nd        
723        real*8 rgrid(1:ngr),dentv(1:nvv),cu(1:nu*nvv,1:ngr) 
725        real*8 den(1:nv),gu(1:nu*nvv,1:ngr),nc(1:nu*nvv,1:ngr)
727        real*8 cnew(1:nu,1:nvv,1:ngr), hnew(1:nu,1:nvv,1:ngr)
728        real*8 ul(1:nu,1:nvv,1:ngr),hq(1:nu,1:nvv,1:ngr)
729        real*8 ulr(1:nu,1:nvv,1:ngr),qvv(1:nvv),qqu(1:nu)
745          do i=2,ngr 
754         do i=1,ngr        
759         do i=1,ngr        
764         do i=1,ngr        
783          do i=2,ngr 
792          do i=2,ngr