ETS  \$Id: Doxyfile 2162 2020-02-26 14:16:09Z g2dpc $
 All Classes Files Functions Variables Pages
_eqa_m_r.f
Go to the documentation of this file.
1 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
2  subroutine eqa(
3  * keyctr,igdf,nstep,platok, psax,i_betp,betplx,
4  * rax,zax, rxpnt,zxpnt, psbo, psdel,
5  * rk,zk,nk, pcequi,ncequi, psitok,
6  * rloop,zloop,nloop,
7  * rprob,zprob,nprob,
8  * zli3,betpol,betful,
9  * necon,wecon,ntipe , nflag, errarr)
10 
11  include 'double.inc'
12  include 'parevo.inc'
13  parameter(nkp=njlim)
14  include 'dim.inc'
15  ! include 'urs.inc'
16  !parameter(nursp4=nursp+4,nursp6=nursp4*6)
17  include 'compol.inc'
18  include 'compol_add.inc'
19 c -----------------------------
20 
21  real*8 rk(*),zk(*),psitok(*),wecon(*)
22  real*8 rloop(*),zloop(*),rprob(*),zprob(*)
23  real*8 pcequi(*)
24 
25  real*8 errarr(*)
26  integer ntipe(*),necon(*)
27  integer nk,nloop,nprob,nstep,keyctr
28 
29  real*8 alf0,alf1,alf2,bet0,bet1,bet2
30  dimension alm(4,4),blm(4),xlm(4),iwrk(4)
31  real time_beg,time_end,time_b,time_e,dtim1,dtim2,dtim3
32  !dimension pstab(nursp),qtab(nursp)
33  !real*8 RRK(nursp4),CCK(nursp4),WRK(nursp6)
34  !real*8 CWK(4)
35  common
36  * /c_kpr/kpr
37 
38  abs(xx)=dabs(xx)
39  sqrt(xx)=dsqrt(xx)
40 
41 
42 
43  cnor=1.d0
44  iplasm=iplas
45  nroi=nr
46  ntetj=nt
47  ngav=keyctr
48 
49  epscrz=5.d-6
50  epspsm=5.d-7
51  epsfpv=5.d-4
52  epsro =1.d-7
53 
54  rolim=ro(iplas,jrolim)
55 
56  if(ngav.ne.2) then
57  fpv=0.d0
58  endif
59 
60 
61  !pause ' '
62 
63  !tok=tokf
64  !qcen=qc
65 
66 
67  itout=0
68  iter=0
69  itin=0
70  ich=0
71 
72 
73 
74  call f_bndmat(rk,zk,nk,rloop,zloop,nloop,rprob,zprob,nprob)
75  !write(*,*)'f_bndmat'
76 
77  1000 continue
78 
79  iter=iter+1
80  itin=itin+1
81 
82  !write(*,*)'iter=',iter,itin
83  !write(*,*)'ich=',ich
84 
85  call f_metric
86  !write(*,*)'metric'
87  call f_matcof
88  call matpla
89  !write(*,*)'matpla'
90  call extpol
91  !write(*,*)'ext '
92 
93  call rigext
94  !write(*,*)'rigext '
95 
96  call solext
97 
98  !write(*,*)'solext '
99  call f_matrix
100  !write(*,*)'matrix'
101 
102  call f_rightg
103  !write(*,*)'rightg'
104 !%%%%%%%%%%%%%%%%%%% test %%%%%%%%%%%%%%%%%%%%%
105 ! if(nstep.eq.1) then
106 ! call rightg_test
107 ! call rightp_test
108 ! call f_solve(0,g)
109 ! call g_test
110 ! stop
111 ! endif
112 !%%%%%%%%%%%%%%%%%%% test %%%%%%%%%%%%%%%%%%%%%
113 
114  if(i_betp.eq.1) then
115  if(iter.gt.4) call skbetp(betplx,betpol)
116  endif
117 
118  call f_solve(0,g)
119  !write(*,*)'solve(g)'
120 
121  call f_rightp
122  !write(*,*)'rightp'
123 
124  call f_solve(1,psii)
125 
126  !write(*,*)'solve(psii)'
127 
128  call f_psiful
129  !call f_wrd
130  !call out
131  !write(*,*)'psiful '
132  call artfil
133  if(kpr.eq.1) then
134  write(*,*)'artfil ',clr,clz
135  write(*,*)'psiax,psim',psiax,psim
136  endif
137  !if(ngav.gt.0) then
138  ! call f_wrd
139  ! pause 'wrd'
140  !endif
141 
142 
143  if(ngav.eq.0 .AND. igdf.eq.2) then
144  call qst_b
145  call grdef(igdf)
146  endif
147 
148 
149  call ada(erro)
150  ! call f_wrd
151  !pause 'wrd'
152  if(kpr.eq.1) then
153  write(*,*)'ada '
154  write(*,*)'erro=',erro
155  endif
156 !!!!!!!!!!!! accurasy test parameters !!!!!!!!!!!!!!!
157 
158  cab=abs(clr)+abs(clz)
159 
160  if(ngav/10*10.eq.ngav) then
161 
162  errpsm=0.d0
163  errfpv=0.d0
164 
165  elseif(ngav.eq.1) then
166 
167  errpsm=abs((psiax-psim)/psipla)
168  errpsb=abs((psip-psibon)/psipla)
169  fvv=sqrt(f(iplas)**2+fpv)
170  !errfpv=abs((fvac-fvv)/(f(1)-f(iplas)))
171  errfpv=0.d0
172 
173  elseif(ngav.eq.2) then
174 
175  errpsm=abs((psiax-psim)/psipla)
176  errpsb=abs((psip-psibon)/psipla)
177 
178  fvv=sqrt(f(iplas)**2+fpv)
179  errfpv=abs((fvac-fvv)/(f(1)-f(iplas)))
180 
181  endif
182 
183 
184 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
185 
186  if(erro.lt.epsro .OR. itin.gt.itrmax) then
187 
188  if(ich.ne.0) go to 2000
189 
190  if(kpr.eq.1) then
191  write(*,*) 'itin',itin
192  write(*,*) 'crz',cab
193  write(*,*) 'errpsm',errpsm
194  write(*,*) 'errpsb',errpsb
195  write(*,*) 'errfpv',errfpv
196  endif
197  if( (cab.lt.epscrz) .AND.
198  * (errpsm.lt.epspsm) .AND.
199  * (errfpv.lt.epsfpv) ) go to 3000
200 
201  if( itout.gt.nitmax ) go to 3000
202 
203  go to 2000
204 
205  endif
206 
207  go to 1000
208 
209  2000 continue
210  itin=0
211 
212  if(ich.eq.0) then
213 
214  cr0=clr
215  cz0=clz
216  pm0=psim
217  fv0=f(iplas)**2+fpv
218 
219  rm0=rm
220  zm0=zm
221  rolim0=rolim
222  fpv0=fpv
223 
224  ich=1
225 
226  drm=0.101*(ro(2,1)-ro(1,1))
227  zm=zm+drm
228 
229  !call f_wrd
230  call reform
231  !write(*,*)'reform1 '
232  !call f_wrd
233 
234  !pause ' '
235  go to 1000
236 
237  elseif(ich.eq.1) then
238 
239  cr1=clr
240  cz1=clz
241  pm1=psim
242  fv1=f(iplas)**2+fpv
243 
244 
245  ich=2
246  rm=rm0+drm
247  zm=zm0
248  rolim=rolim0
249  fpv=fpv0
250 
251  call reform
252  ! write(6,*)'reform2 '
253  !pause ' '
254  go to 1000
255 
256  elseif(ich.eq.2 .AND. ngav/10*10.eq.ngav) then
257 
258  cr2=clr
259  cz2=clz
260 
261  dcrdr=(cr2-cr0)/drm
262  dczdr=(cz2-cz0)/drm
263 
264  dcrdz=(cr1-cr0)/drm
265  dczdz=(cz1-cz0)/drm
266 
267  det=dcrdr*dczdz-dczdr*dcrdz
268 
269  delr= (cz0*dcrdz-cr0*dczdz)/det
270  delz= (cr0*dczdr-cz0*dcrdr)/det
271 
272 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
273  dell=sqrt(delr**2+delz**2)
274  if(dell.gt.drm*10.d0) then
275  nshift=dell/(drm*10.d0)
276 
277  delr=drm*10.d0*delr/dell
278  delz=drm*10.d0*delz/dell
279 
280  if(nshift.gt.10) nshift=10
281 
282  do 6793 ish=1,nshift
283 
284  !write(6,*)'slow shift',ish,nshift
285  rm=rm0+ delr*ish
286  zm=zm0+ delz*ish
287 
288  !write(6,*) 'delr',delr
289  !write(6,*) 'delz',delz
290  !pause ' '
291 
292  call reform
293  !write(6,*)'reform3 '
294  !!!pause ' '
295  call f_metric
296  call extpol
297  call f_matcof
298  call matpla
299  call rigext
300  call solext
301  call f_matrix
302  call f_rightg
303  call f_solve(0,g)
304  call f_rightp
305  call f_solve(1,psii)
306  call f_psiful
307  call artfil
308  !write(6,*)'artfil ',clr,clz
309  call ada(erro)
310  !write(6,*)'ada '
311  !write(6,*)'erro=',erro
312  !write(6,*)' '
313 
314  6793 continue
315 
316 
317 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
318  else
319 
320  rm=rm0+ delr
321  zm=zm0+ delz
322 
323  !write(6,*) 'delr',delr
324  !write(6,*) 'delz',delz
325  !pause ' '
326 
327  call reform
328  !write(6,*)'reform3 '
329  !!!pause ' '
330  endif
331 
332  ich=0
333  itout=itout+1
334  go to 1000
335 
336  elseif(ich.eq.2 .AND. ngav.gt.0) then
337 
338  cr2=clr
339  cz2=clz
340  pm2=psim
341  fv2=f(iplas)**2+fpv
342 
343  ich=3
344  rm=rm0
345  zm=zm0
346  drolim=-drm
347  rolim=rolim0+drolim
348  fpv=fpv0
349 
350  call reform
351  ! write(6,*)'reform3 '
352  !pause ' '
353  go to 1000
354 
355  elseif(ich.eq.3 .AnD. ngav.eq.1) then
356 
357  cr3=clr
358  cz3=clz
359  pm3=psim
360 
361  dcrdr=(cr2-cr0)/drm
362  dczdr=(cz2-cz0)/drm
363  dpmdr=(pm2-pm0)/drm
364 
365  dcrdz=(cr1-cr0)/drm
366  dczdz=(cz1-cz0)/drm
367  dpmdz=(pm1-pm0)/drm
368 
369  dcrdro=(cr3-cr0)/drolim
370  dczdro=(cz3-cz0)/drolim
371  dpmdro=(pm3-pm0)/drolim
372 
373  !write(6,*) 'deriv'
374 
375 
376  alm(1,1)=dcrdr
377  alm(1,2)=dcrdz
378  alm(1,3)=dcrdro
379 
380  alm(2,1)=dczdr
381  alm(2,2)=dczdz
382  alm(2,3)=dczdro
383 
384  alm(3,1)=dpmdr
385  alm(3,2)=dpmdz
386  alm(3,3)=dpmdro
387  do ic=1,3
388  do jc=1,3
389  !write(6,*)'i,j',ic,jc
390  !write(6,*)'a(i,j)',Alm(ic,jc)
391  enddo
392  enddo
393 
394  blm(1)=-cr0
395  blm(2)=-cz0
396  blm(3)= psiax-pm0
397 
398  !write(6,*) 'mat a'
399 
400  call ge(3,4,alm,blm,xlm,iwrk)
401 
402  !write(6,*) 'ge '
403 
404  delr=xlm(1)
405  delz=xlm(2)
406  delro=xlm(3)
407 
408  !write(6,*) 'delr',delr
409  !write(6,*) 'delz',delz
410  !write(6,*) 'delro',delro
411  !!!pause ' '
412 
413  rm=rm0+ delr
414  zm=zm0+ delz
415  rolim=rolim0+delro
416 
417  call reform
418  !write(6,*)'reform*** '
419  !!!pause ' '
420  ich=0
421  itout=itout+1
422  go to 1000
423 
424  elseif(ich.eq.3 .AnD. ngav.gt.1) then
425 
426  cr3=clr
427  cz3=clz
428  pm3=psim
429  fv3=f(iplas)**2+fpv
430 
431  ich=4
432  rm=rm0
433  zm=zm0
434  rolim=rolim0
435  !dfpv=0.005d0
436  dfpv=(f(1)**2-fvac**2)*1.d-3
437  fpv=fpv0+dfpv
438 
439  call reform
440  ! write(6,*)'reform4 '
441  !!! pause ' '
442  go to 1000
443 
444  elseif(ich.eq.4) then
445 
446  cr4=clr
447  cz4=clz
448  pm4=psim
449  fv4=f(iplas)**2+fpv
450 
451  dcrdr=(cr2-cr0)/drm
452  dczdr=(cz2-cz0)/drm
453  dpmdr=(pm2-pm0)/drm
454  dfvdr=(fv2-fv0)/drm
455 
456  dcrdz=(cr1-cr0)/drm
457  dczdz=(cz1-cz0)/drm
458  dpmdz=(pm1-pm0)/drm
459  dfvdz=(fv1-fv0)/drm
460 
461  dcrdro=(cr3-cr0)/drolim
462  dczdro=(cz3-cz0)/drolim
463  dpmdro=(pm3-pm0)/drolim
464  dfvdro=(fv3-fv0)/drolim
465 
466  dcrdf=(cr4-cr0)/dfpv
467  dczdf=(cz4-cz0)/dfpv
468  dpmdf=(pm4-pm0)/dfpv
469  dfvdf=(fv4-fv0)/dfpv
470 
471  ! write(6,*) 'deriv'
472 
473 
474  alm(1,1)=dcrdr
475  alm(1,2)=dcrdz
476  alm(1,3)=dcrdro
477  alm(1,4)=dcrdf
478 
479  alm(2,1)=dczdr
480  alm(2,2)=dczdz
481  alm(2,3)=dczdro
482  alm(2,4)=dczdf
483 
484  alm(3,1)=dpmdr
485  alm(3,2)=dpmdz
486  alm(3,3)=dpmdro
487  alm(3,4)=dpmdf
488 
489  alm(4,1)=dfvdr
490  alm(4,2)=dfvdz
491  alm(4,3)=dfvdro
492  alm(4,4)=dfvdf
493  do ic=1,4
494  do jc=1,4
495  !write(6,*)'i,j',ic,jc
496  !write(6,*)'a(i,j)',Alm(ic,jc)
497  enddo
498  enddo
499 
500  blm(1)=-cr0
501  blm(2)=-cz0
502  blm(3)= psiax-pm0
503  blm(4)= fvac**2-fv0
504 
505  !write(6,*) 'mat a'
506 
507  call ge(4,4,alm,blm,xlm,iwrk)
508 
509  !write(6,*) 'ge '
510 
511  delr=xlm(1)
512  delz=xlm(2)
513  delro=xlm(3)
514  delfpv=xlm(4)
515 
516  !write(6,*) 'delr',delr
517  !write(6,*) 'delz',delz
518  !write(6,*) 'delro',delro
519  !write(6,*) 'delfv',delfpv
520  !!!pause ' '
521 
522  rm=rm0+ delr
523  zm=zm0+ delz
524  rolim=rolim0+delro
525  fpv=fpv0+delfpv
526 
527  call reform
528  ! write(6,*)'reform*** '
529  !!!pause ' '
530  ich=0
531  itout=itout+1
532  go to 1000
533 
534  endif
535 
536  3000 continue
537 
538 
539  nflag=0
540 
541  if(itout.gt.nitmax) nflag=1
542 
543  errarr(1)=erro ! out of accuracy tests
544  errarr(2)=cab !
545  errarr(3)=errpsm !
546  errarr(4)=errfpv !
547 
548 
549 
550  ! if(ngav.gt.0 .AND. itin.lt.15) go to 1000
551 
552  ! if(ngav.eq.1 .AND. itout.lt.2) then
553  ! call f_wrd
554  ! go to 2000
555  ! endif
556 
557  !if(ngav.eq.2 .AND. itout.lt.2) then
558  ! f_call wrd
559  ! go to 2000
560  !endif
561 
562  if( itout.lt.1) then
563  ! f_call wrd
564  go to 2000
565  endif
566 
567  !if(ngav/10*10.eq.ngav) then
568  if(ngav.eq.0) then
569  call qst_b
570  psiax=psim
571  psipla=psim-psip
572  fvac = f(iplas)
573  endif
574 
575 ! write(fname,'(a,a)') path(1:kname),'ddps0.pr'
576 ! open(1,file='ddps0.pr')
577 ! do 567 i=1,iplas
578 ! 567 write(1,*) dfdpsi(i),f(i),i
579 ! close(1)
580  if(kpr.eq.1) then
581  write(*,*) 'iter itout',iter,itout
582  endif
583 c----------------------------------------------------------
584  !call f_flux(psitok,rk,zk,nk)
585  !call flux_p(psitok,rk,zk,nk)
586  !call flux_g(psitok,rk,zk,nk)
587  call flux_r(psitok,ncequi)
588 
589  platok = tokp
590  rax = rm
591  zax = zm
592  rxpnt = rx0
593  zxpnt = zx0
594  psax = psim
595  psbo = psip
596  psdel = psim - psip
597 
598  call bt_pol(betpol)
599  call f_wrd
600  ! call out(rbnd,zbnd,zli3,betpol,betful,parpla)
601 c write(6,*)'wrd '
602 
603  iter=0
604  itin=0
605 
606 
607  return
608  end
609 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
610 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
611 
612  subroutine eqa_ax( dt,time,
613  * keyctr,igdf,nstep,platok, psax,i_betp,betplx,
614  * rax,zax, rxpnt,zxpnt, psbo, psdel,
615  * rk,zk,nk, pcequi,ncequi, psitok,
616  * rloop,zloop,nloop, rprob,zprob,nprob,
617  * zli3,betpol,betful,
618  * necon,wecon,ntipe , nflag, errarr)
619 
620 
621  include 'double.inc'
622  include 'parevo.inc'
623  parameter(nkp=njlim)
624  include 'dim.inc'
625  include 'compol.inc'
626  include 'compol_add.inc'
627 c -----------------------------
628 
629  real*8 rk(*),zk(*),psitok(*),wecon(*)
630  real*8 pcequi(*)
631  real*8 rloop(*),zloop(*),rprob(*),zprob(*)
632 
633  real*8 errarr(*)
634  integer ntipe(*),necon(*)
635  integer nk,nstep,keyctr
636 
637  real*8 alf0,alf1,alf2,bet0,bet1,bet2
638  real time_beg,time_end,time_b,time_e,dtim1,dtim2,dtim3
639  common
640  * /c_kpr/kpr
641 
642 ! save nstepO
643  integer nstepo
644 
645  abs(xx)=dabs(xx)
646  sqrt(xx)=dsqrt(xx)
647 
648  ngav=keyctr
649  kstep=nstep
650  dtim=dt
651  ctim=time
652 
653  itrmax=50
654  nitmax=5
655  nitdel=7
656  nitbeg=5
657 
658  cnor=1.d0
659 
660  !pause ' '
661 
662  call f_ext_fil(pcequi,ncequi)
663 
664 c write(*,*)'extrec'
665 
666  if(nstep.ne.nstepo) then
667  itin=0
668  erru=0.d0
669  call renet
670  !call f_wrd
671  call f_bndmat(rk,zk,nk,rloop,zloop,nloop,rprob,zprob,nprob)
672  !call f_wrd
673  !pause 'pause'
674  endif
675 
676  1000 continue
677 
678  iter=iter+1
679  itin=itin+1
680 
681  if(kpr.eq.1) then
682  write(*,*)' '
683  write(*,*)'iter=',iter,itin
684  endif
685  call f_metric
686 c write(*,*)'metric'
687  call f_matcof
688 c write(*,*)'matcof'
689  call matpla
690 c ! write(*,*)'matpla'
691  call extpol
692 c write(*,*)'ext '
693  call rigext
694 c ! write(*,*)'rigext '
695  call solext
696 c ! write(*,*)'solext '
697  call f_matrix
698 c write(*,*)'matrix'
699 
700  call f_rightg
701 
702 c write(*,*)'rightg'
703 
704  if(i_betp.eq.1) then
705  if(iter.gt.4) call skbetp(betplx,betpol)
706  endif
707 
708 
709  call f_solve(0,g)
710  !f_call wrd
711  !call out
712  !stop
713 c write(*,*)'solve(g)'
714  call f_rightp
715 c write(*,*)'rightp'
716 
717  call f_solve(1,psii)
718 
719 c write(*,*)'solve(psii)'
720 
721  !call flux_g(psitok,rk,zk,nk)
722  !call flux_r(psitok,ncequi)
723  ! write(6,*)'flux'
724 
725  call f_psiful
726  !f_call wrd
727  !call out
728 c write(*,*)'psiful '
729 
730  if(ngav.le.0 .AND. igdf.eq.2) then
731  call qst_b
732  call grdef(igdf)
733  endif
734 
735  call f_remesh(erro)
736  erru=erro
737  !call flux_p(psitok,rk,zk,nk)
738  call flux_r(psitok,ncequi)
739 
740  if(ngav.eq.0) then
741 
742  errpsm=0.d0
743 
744  elseif(ngav.eq.1) then
745 
746  errpsm=abs((psiax-psim)/psipla)
747 
748  endif
749 
750 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
751 
752 
753 c write(*,*) 'itin',itin
754 c write(*,*) 'errpsm',errpsm
755 
756  3000 continue
757 
758 
759 c write(*,*) 'iter, itout ==', iter, itout
760 c----------------------------------------------------------
761 
762  platok = tokp
763  rax = rm
764  zax = zm
765  rxpnt = rx0
766  zxpnt = zx0
767  psax = psim
768  psbo = psip
769  psdel = psim-psip
770 
771  if(ngav.le.0) then
772  psiax = psim
773  psipla = psim-psip
774  !Fvac = f(iplas)
775  endif
776  call bt_pol(betpol)
777  !call f_wrd
778  !call out(rbnd,zbnd,zli3,betpol,betful,parpla)
779  !write(*,*)'f_wrd '
780 
781  nflag=0
782 
783  errarr(1) = erro ! out of accuracy tests
784  !errarr(2) = cab !
785  !errarr(3) = errpsm !
786  !errarr(4) = errfpv !
787 
788  nstepo=nstep
789 
790  if(ngav.lt.0) call retab_l
791 
792  return
793  end
794 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
795 
eqa_ax
subroutine eqa_ax(dt, time,
Definition: _eqa_m_r.f:612
renet
subroutine renet
Definition: _mesh.f:6
f_rightp
subroutine f_rightp
Definition: _rig.f:260
f_wrd
subroutine f_wrd
Definition: _wrd.f:107
solext
subroutine solext
Definition: _sol.f:410
extpol
subroutine extpol
Definition: _extp.f:480
reform
subroutine reform
Definition: _mesh.f:2919
flux_r
subroutine flux_r(psitok, ncequi)
Definition: _fluxt.f:169
f_rightg
subroutine f_rightg
Definition: _rig.f:12
artfil
subroutine artfil
Definition: _extp.f:617
grdef
subroutine grdef(igdf)
Definition: com_sub.f:7
f_remesh
subroutine f_remesh(erro)
Definition: _mesh.f:617
skbetp
subroutine skbetp(betplx, betpol)
Definition: com_sub.f:781
matpla
subroutine matpla
Definition: com_sub.f:1061
ada
subroutine ada(erro)
Definition: _mesh.f:2616
bt_pol
subroutine bt_pol(betpol)
Definition: com_sub.f:747
f_ext_fil
subroutine f_ext_fil(pcequi, ncequi)
Definition: _ext_m.f:4
qst_b
subroutine qst_b
Definition: com_sub.f:576
f_metric
subroutine f_metric
Definition: _metric.f:16
f_matcof
subroutine f_matcof
Definition: _metric.f:279
eqa
subroutine eqa(
Definition: _eqa_m_r.f:2
ge
subroutine ge(N, NZ, A, X, Y, IP)
Definition: scu.f:420
retab_l
subroutine retab_l
Definition: B_eqb.f:541
f_bndmat
subroutine f_bndmat(rk, zk, nk, rlop, zlop, nlop, rprob, zprob, nprob)
Definition: _bnd.f:1
f_matrix
subroutine f_matrix
Definition: _metric.f:344
rigext
subroutine rigext
Definition: _rig.f:344
f_psiful
subroutine f_psiful
Definition: _extp.f:572
f_solve
subroutine f_solve(isol, wdm)
Definition: _sol.f:1