ETS  \$Id: Doxyfile 2162 2020-02-26 14:16:09Z g2dpc $
 All Classes Files Functions Variables Pages
_rig.f
Go to the documentation of this file.
1 !! file contains following routines:
2 !!
3 !! funp
4 !! funf
5 !! f_rightg
6 !! f_rightp
7 !! rigext
8 !! toksur
9 
10 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
11 
12  subroutine f_rightg
13 
14 !! riht-hand side for problem L(g)=J
15 
16  include 'double.inc'
17  include 'parevo.inc'
18  parameter(nkp=njlim)
19  include 'dim.inc'
20  include 'compol.inc'
21  include 'compol_add.inc'
22  common
23  * /c_kpr/kpr
24  common/com_flag/kastr
25 ! save cur_mu
26  real*8 cur_mu
27 
28  tokp=0.d0
29 
30  do 1 il=1,neqp
31 
32  right(il)=0.d0
33 
34  1 continue
35 
36  do 2 i=1,nr
37  do 2 j=1,nt
38 
39  cur(i,j)=0.d0
40 
41  2 continue
42 
43  !call aprod
44 
45  if(kstep.eq.0 .OR. (ngav.eq.0 .AnD. kastr.eq.0)) then
46  !if(ngav.eq.0) then
47 ch4astra if(kstep.eq.0) then
48  !!!!!!!!!!!!!!!!!!!!
49 !! central point
50 
51  curp=tabp(1.d0)
52  curf=tabf(1.d0)
53 
54  dpdpsi(1)=curp
55  dfdpsi(1)=curf
56 
57  curcen=rm*curp+curf/rm
58 
59  sqcen=0.d0
60 
61  do 10 j=2,nt1
62 
63  sqcen=sqcen+sq1(1,j)+sq4(1,j)
64 
65  10 continue
66 
67  tokp=curcen*sqcen
68  tokff=curf*sqcen/rm
69  tokpp=curp*sqcen*rm
70 
71  il=numlin(1,1,nr,nt)
72 
73  right(il)=curcen*sqcen
74 
75 !! regular points
76 
77  do 20 i=2,iplas
78  do 25 j=2,nt1
79 
80  if(i.ne.iplas) then
81  sqk=sq1(i,j)+sq2(i-1,j)+sq3(i-1,j-1)+sq4(i,j-1)
82  r0=r(i,j)
83  psn=psin(i,j)
84  else
85  sqk=sq2(i-1,j)+sq3(i-1,j-1)
86  r0=(3.d0*r(i,j)+r(i-1,j))*0.25d0
87  psn=(3.d0*psin(i,j)+psin(i-1,j))*0.25d0
88  !psn=0.d0
89  endif
90 
91  curp=tabp(psn)
92  curf=tabf(psn)
93 
94  cur(i,j)=r0*curp+curf/r0
95 
96  tokp=tokp+cur(i,j)*sqk
97  tokff=tokff+curf*sqk/r0
98  tokpp=tokpp+curp*sqk*r0
99 
100  il=numlin(i,j,nr,nt)
101  right(il)=cur(i,j)*sqk
102 
103  25 continue
104 
105  dpdpsi(i)=curp
106  dfdpsi(i)=curf
107 
108  20 continue
109 
110  cnor=amu0*tok/tokp
111  cur_mu=amu0*tok
112 
113  if(kpr.eq.1) then
114  write(*,*) 'cnor,tok,tokp'
115  write(*,*) cnor,tok,tokp/amu0
116  endif
117 
118 
119  do 30 j=1,nt
120 
121  cur(1,j)=curcen*cnor
122 
123  30 continue
124 
125  do 40 i=2,iplas
126  do 40 j=2,nt1
127 
128  cur(i,j)=cur(i,j)*cnor
129 
130  40 continue
131 
132  do 45 i=1,nr
133 
134  cur(i,1)=cur(i,nt1)
135  cur(i,nt)=cur(i,2)
136 
137  45 continue
138 
139  do 50 il=1,neq
140 
141  right(il)=right(il)*cnor !-zpro(il)
142 
143  50 continue
144 
145  do 55 i=1,iplas
146 
147  dpdpsi(i)=dpdpsi(i)*cnor
148  dfdpsi(i)=dfdpsi(i)*cnor
149 
150  55 continue
151 
152  tokp=tok
153 
154  else!!!!!!!!!!!!!!!!!!!
155  !!!!!!!!!!!!!!!!!!!!!!!!
156 
157  if(kastr.eq.1) then
158  call pres_d_psi
159  endif
160 
161 ch4astra
162  call f_procof(1,cur_mu)
163  !if(ngav.ne.0) call f_procof(1,cur_mu)
164 
165 !! central point
166 
167  curcen=rm*dpdpsi(1)+dfdpsi(1)/rm
168 
169  sqcen=0.d0
170 
171  do 100 j=2,nt1
172 
173  sqcen=sqcen+sq1(1,j)+sq4(1,j)
174 
175  100 continue
176 
177  tokp=curcen*sqcen
178 
179  il=numlin(1,1,nr,nt)
180 
181  right(il)=curcen*sqcen
182 
183  do 300 j=1,nt
184 
185  cur(1,j)=curcen
186 
187  300 continue
188 
189 !! regular points
190 
191  do 200 i=2,iplas
192  do 200 j=2,nt1
193 
194  if(i.ne.iplas) then
195  sqk=sq1(i,j)+sq2(i-1,j)+sq3(i-1,j-1)+sq4(i,j-1)
196  r0=r(i,j)
197  else
198  sqk=sq2(i-1,j)+sq3(i-1,j-1)
199  r0=(3.d0*r(i,j)+r(i-1,j))*0.25d0
200  endif
201 
202  curp=dpdpsi(i)
203  curf=dfdpsi(i)
204 
205  cur(i,j)=r0*curp+curf/r0
206 
207  tokp=tokp+cur(i,j)*sqk
208 
209  il=numlin(i,j,nr,nt)
210  right(il)=cur(i,j)*sqk
211 
212  200 continue
213 
214  !write(*,*) 'rightp:tokp=',tokp
215 
216  !!!!!!!!!!!!!!!!!!!!!!!!!!
217  if(ngav.eq.2) call toksur!
218  !!!!!!!!!!!!!!!!!!!!!!!!!!
219 
220  cnor=cur_mu/tokp
221 
222 
223  do 450 i=1,nr
224 
225  cur(i,1)=cur(i,nt1)
226  cur(i,nt)=cur(i,2)
227 
228  450 continue
229 
230  tokp=tokp/amu0
231  write(*,*) 'rightp:tok tokp ',tok,tokp
232 
233  if(ngav.eq.0 ) then
234 
235  tokp=cur_mu/amu0
236 
237  do i=1,iplas
238  do j=1,nt
239  cur(i,j)=cur(i,j)*cnor
240  enddo
241  enddo
242 
243  do il=1,neq
244  right(il)=right(il)*cnor ! !-zpro(il)
245  enddo
246 
247  endif
248 
249  if(kpr.eq.1) then
250  write(*,*) 'cnor cur_mu',cnor,cur_mu
251  endif
252 
253  endif
254 
255  return
256  end
257 
258 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
259 
260  subroutine f_rightp
261 
262  include 'double.inc'
263  include 'parevo.inc'
264  parameter(nkp=njlim)
265  include 'dim.inc'
266  include 'compol.inc'
267  include 'compol_add.inc'
268 
269 
270  common/comaaa/ a12(nrp,ntp),a23(nrp,ntp),a34(nrp,ntp),
271  + a14(nrp,ntp),a13(nrp,ntp),a24(nrp,ntp)
272 
273  real*8 psib(ntp)
274 
275  i=nr !!!!!!!!!!!!!!!!!!!!!!
276 
277  do 10 j=2,nt1
278 
279  a1=a13(i-1,j-1)
280  a2=a34(i-1,j-1)+a12(i-1,j)
281  a3=a24(i-1,j)
282 
283  g1=g(i-1,j-1)
284  g2=g(i-1,j)
285  g3=g(i-1,j+1)
286 
287  dltk=(dlt(i,j-1)+dlt(i,j))*0.5d0
288  dgdnl=a1*g1+a2*g2+a3*g3
289  dgdn(j)=dgdnl/dltk
290 
291  10 continue
292 
293  dgdn(1)=dgdn(nt1)
294  dgdn(nt)=dgdn(2)
295 
296  do 20 j=2,nt1
297 
298  psb=0.d0
299 
300  do 22 jb=2,nt1
301 
302  psb=psb+binadg(jb,j)*(dgdn(jb)+dgdn(jb+1))*0.5d0
303 
304  22 continue
305 
306  psib(j)=psb
307  psii(i,j)=-psb
308 
309  20 continue
310 
311  psib(1)=psib(nt1)
312  psib(nt)=psib(2)
313 
314  psii(i,1)=-psib(nt1)
315  psii(i,nt)=-psib(2) !
316  !!
317  ! !
318  i=nr1 !!!!!!!!!!!!!!!!!!!!! !
319  ! !
320  do 100 j=2,nt1 !!
321  !
322  a7=a24(i,j-1)
323  a8=a34(i,j-1)+a12(i,j)
324  a9=a13(i,j)
325 
326  il=numlin(i,j,nr,nt)
327 
328  right(il)= (a7*psib(j-1)+a8*psib(j)+a9*psib(j+1))
329 
330  100 continue
331 
332 
333 ! do 200 il=1,neq
334 
335 ! right(il)=right(il)!-zpro(il)
336 
337 ! 200 continue
338 
339 
340  return
341  end
342 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
343 
344  subroutine rigext
345 
346  include 'double.inc'
347  include 'parevo.inc'
348  parameter(nkp=njlim)
349  include 'dim.inc'
350  include 'compol.inc'
351  include 'compol_add.inc'
352 
353  common/comaaa/ a12(nrp,ntp),a23(nrp,ntp),a34(nrp,ntp),
354  + a14(nrp,ntp),a13(nrp,ntp),a24(nrp,ntp)
355 
356  real*8 psib(ntp)
357 
358  do 1 il=1,neqp
359 
360  right(il)=0.d0
361 
362  1 continue
363 
364  do 20 j=2,nt1
365 
366  psib(j)=psie(iplas,j)
367 
368  20 continue
369 
370  psib(1)=psib(nt1)
371  psib(nt)=psib(2)
372  !
373  !!
374  ! !
375  i=iplas-1 !!!!!!!!!!!!!!!!! !
376  ! !
377  do 100 j=2,nt1 !!
378  !
379  a7=a24(i,j-1)
380  a8=a34(i,j-1)+a12(i,j)
381  a9=a13(i,j)
382 
383  il=numlin(i,j,nr,nt)
384 
385  right(il)=-(a7*psib(j-1)+a8*psib(j)+a9*psib(j+1))
386 
387  100 continue
388 
389  return
390  end
391 
392 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
393 
394  subroutine f_rigbon
395 
396  include 'double.inc'
397  include 'parevo.inc'
398  parameter(nkp=njlim)
399  include 'dim.inc'
400  include 'compol.inc'
401 
402  common/comaaa/ a12(nrp,ntp),a23(nrp,ntp),a34(nrp,ntp),
403  + a14(nrp,ntp),a13(nrp,ntp),a24(nrp,ntp)
404 
405  real*8 psib(ntp)
406 
407 
408  do 20 j=2,nt1
409 
410  !psib(j)=psi(iplas,j)
411  psib(j)=psip
412 
413  20 continue
414 
415  !open(3,file='boncon.wr')
416  !read(3,*) (psib(j),j=2,nt1)
417  !close(3)
418 
419 
420 
421 
422  psib(1)=psib(nt1)
423  psib(nt)=psib(2)
424  !
425  !!
426  ! !
427  i=iplas-1 !!!!!!!!!!!!!!!!! !
428  ! !
429  do 100 j=2,nt1 !!
430  !
431  a7=a24(i,j-1)
432  a8=a34(i,j-1)+a12(i,j)
433  a9=a13(i,j)
434 
435  il=numlin(i,j,nr,nt)
436 
437  right(il)=right(il)-(a7*psib(j-1)+a8*psib(j)+a9*psib(j+1))
438 
439  100 continue
440 
441  return
442  end
443 
444 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
445 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
446 
447  subroutine toksur
448 
449  include 'double.inc'
450  include 'parevo.inc'
451  parameter(nkp=njlim)
452  include 'dim.inc'
453  include 'compol.inc'
454  include 'compol_add.inc'
455 
456  common/comaaa/ a12(nrp,ntp),a23(nrp,ntp),a34(nrp,ntp),
457  + a14(nrp,ntp),a13(nrp,ntp),a24(nrp,ntp)
458 
459  sqrt(xx)=dsqrt(xx)
460 
461  i=iplas
462 
463  toksfi=0.d0
464 
465  do 10 j=2,nt1
466 
467  il=numlin(i,j,nr,nt)
468 
469  a1=a13(i-1,j-1)
470  a2=a34(i-1,j-1)+a12(i-1,j)
471  a3=a24(i-1,j)
472  a4=a23(i-1,j-1)
473  a6=a23(i-1,j)
474 
475  a5=-(a1+a2+a3+a4+a6)
476 
477  ps1=psii(i-1,j-1)
478  ps2=psii(i-1,j)
479  ps3=psii(i-1,j+1)
480 
481  ps4=psii(i,j-1)
482  ps5=psii(i,j)
483  ps6=psii(i,j+1)
484 
485  dpidni=a1*ps1+a2*ps2+a3*ps3+a4*ps4+a5*ps5+a6*ps6-right(il)
486 
487  ps1=psie(i-1,j-1)
488  ps2=psie(i-1,j)
489  ps3=psie(i-1,j+1)
490 
491  ps4=psie(i,j-1)
492  ps5=psie(i,j)
493  ps6=psie(i,j+1)
494 
495  dpedni=a1*ps1+a2*ps2+a3*ps3+a4*ps4+a5*ps5+a6*ps6
496 
497  a4=a14(i,j-1)
498  a6=a14(i,j)
499  a7=a24(i,j-1)
500  a8=a34(i,j-1)+a12(i,j)
501  a9=a13(i,j)
502 
503  a5=-(a4+a6+a7+a8+a9)
504 
505  ps4=psii(i,j-1)
506  ps5=psii(i,j)
507  ps6=psii(i,j+1)
508 
509  ps7=psii(i+1,j-1)
510  ps8=psii(i+1,j)
511  ps9=psii(i+1,j+1)
512 
513  dpidne=-(a4*ps4+a5*ps5+a6*ps6+a7*ps7+a8*ps8+a9*ps9)
514 
515  dpsdni=dpidni+dpedni
516  dpsdne=dpidne+dpedni
517 
518  !write(6,*) 'De,Di,j',dpsdne,dpsdni,j
519 
520  sqk=sq2(i-1,j)+sq3(i-1,j-1)
521 
522  dlt0=(dlt(i,j-1)+dlt(i,j))*0.5d0
523  r0=r(i,j)
524 
525  !fpv= Fvac**2-f(i)**2
526  !fpv=-7.10d0
527 
528  curs(j)=(dlt0/r0**2)*fpv/(dpsdni+dpsdne)
529  !write(6,*) 'cursurf',curs(j)
530 
531  right(il)=right(il)+curs(j)*dlt0
532 
533  cur(i,j)=curs(j)*dlt0/sqk+cur(i,j)
534  tokp=tokp!+curs(j)*dlt0
535  toksfi=toksfi+curs(j)*dlt0
536 
537  10 continue
538 
539  curs(1)=curs(nt1)
540  curs(nt)=curs(2)
541 
542  fvv=sqrt(f(iplas)**2+fpv)
543  write(6,*) 'Fp,Fvac,fv',f(iplas),fvac,fvv
544  return
545  end
546 
547 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
548 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
549 
550  subroutine rightp_test
551 
552  include 'double.inc'
553  include 'parevo.inc'
554  parameter(nkp=njlim)
555  include 'dim.inc'
556  include 'compol.inc'
557  include 'compol_add.inc'
558 
559  common/comaaa/ a12(nrp,ntp),a23(nrp,ntp),a34(nrp,ntp),
560  + a14(nrp,ntp),a13(nrp,ntp),a24(nrp,ntp)
561 
562  real*8 psib(ntp)
563 
564  i=nr !!<<<<<<<<<<<<<
565 
566  do 10 j=2,nt1
567 
568  a1=a13(i-1,j-1)
569  a2=a34(i-1,j-1)+a12(i-1,j)
570  a3=a24(i-1,j)
571 
572  g1=g(i-1,j-1)
573  g2=g(i-1,j)
574  g3=g(i-1,j+1)
575 
576  dltk=(dlt(i,j-1)+dlt(i,j))*0.5d0
577  dgdnl=a1*g1+a2*g2+a3*g3
578  dgdn(j)=dgdnl/dltk
579 
580  10 continue
581 
582  dgdn(1)=dgdn(nt1)
583  dgdn(nt)=dgdn(2)
584 
585  do 20 j=2,nt1
586 
587  psb=0.d0
588 
589  do 22 jb=2,nt1
590 
591  psb=psb+binadg(jb,j)*(dgdn(jb)+dgdn(jb+1))*0.5d0
592 
593  22 continue
594 
595  psib(j)=psb
596  psii(i,j)=-psb
597 
598  20 continue
599 
600  psib(1)=psib(nt1)
601  psib(nt)=psib(2)
602 
603  psii(i,1)=-psib(nt1)
604  psii(i,nt)=-psib(2)
605 
606 
607 !!%%%%%%%%%%%%%%%%%%%%%% test for g=r**2
608 
609  do j=2,nt1
610  psib(j)=-r(i,j)**2
611  g(i,j)=-psib(j)
612  enddo
613 
614  psib(1)=psib(nt1)
615  psib(nt)=psib(2)
616 
617  g(i,1)=g(i,nt1)
618  g(i,nt)=g(i,2)
619 
620 !!%%%%%%%%%%%%%%%%%%%%%% test for g=r**2
621 
622 
623  !
624  !!
625  ! !
626  i=nr1 !!!!!!!!!!!!!!!!!!!!! !
627  ! !
628  do 100 j=2,nt1 !!
629  !
630  a7=a24(i,j-1)
631  a8=a34(i,j-1)+a12(i,j)
632  a9=a13(i,j)
633 
634  il=numlin(i,j,nr,nt)
635 
636  right(il)= (a7*psib(j-1)+a8*psib(j)+a9*psib(j+1))
637 
638  100 continue
639 
640 
641 ! do 200 il=1,neq
642 
643 ! right(il)=right(il)!-zpro(il)
644 
645 ! 200 continue
646 
647 
648  return
649  end
650 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
651 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
652 
653  subroutine rightg_test
654 
655 !! riht-hand side for problem L(g)=J
656 
657  include 'double.inc'
658  include 'parevo.inc'
659  parameter(nkp=njlim)
660  include 'dim.inc'
661  include 'compol.inc'
662  include 'compol_add.inc'
663 
664  do il=1,neqp
665 
666  right(il)=0.d0
667 
668  enddo
669 
670  return
671  end
672 
673 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
674 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
675 
676  subroutine f_psib_ext(psex_av)
677 
678  include 'double.inc'
679  include 'parevo.inc'
680  parameter(nkp=njlim)
681  include 'dim.inc'
682  include 'compol.inc'
683  include 'compol_add.inc'
684 
685  dimension psex_b(ntp),rbon(ntp),zbon(ntp)
686  dimension pspl_b(ntp),psi_bon(ntp)
687 
688  i=iplas !!!!!!!!!!!!!!!!!!!!!!
689 
690  do j=1,nt
691  rbon(j)=r(i,j)
692  zbon(j)=z(i,j)
693  psex_b(j)=psie(i,j)
694  pspl_b(j)=psii(i,j)
695  psi_bon(j)=pspl_b(j)+psex_b(j)
696  enddo
697 
698  psex_b(1)=psex_b(nt1)
699  psex_b(nt)=psex_b(2)
700 
701  psex_av = avr_bnd(psex_b)
702 
703  return
704  end
705 
706 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
707 
708  subroutine f_psib_pla(pspl_av)
709 
710  include 'double.inc'
711  include 'parevo.inc'
712  parameter(nkp=njlim)
713  include 'dim.inc'
714  include 'compol.inc'
715  include 'compol_add.inc'
716 
717  common/comaaa/ a12(nrp,ntp),a23(nrp,ntp),a34(nrp,ntp),
718  + a14(nrp,ntp),a13(nrp,ntp),a24(nrp,ntp)
719 
720  common/com_bgr/ bin_adg(ntp,ntp),dg_dn(ntp)
721  dimension pspl_b(ntp)
722 
723  i=iplas !!!!!!!!!!!!!!!!!!!!!!
724 
725  do 10 j=2,nt1
726 
727  pspl_b(j)=psii(i,j)
728 
729  10 continue
730 
731  pspl_b(1)=pspl_b(nt1)
732  pspl_b(nt)=pspl_b(2)
733 
734  pspl_av = avr_bnd(pspl_b)
735 
736  return
737  end
738 
739 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
740 
741  subroutine g_test
742 
743 !! riht-hand side for problem L(g)=J
744 
745  include 'double.inc'
746  include 'parevo.inc'
747  parameter(nkp=njlim)
748  include 'dim.inc'
749  include 'compol.inc'
750  include 'compol_add.inc'
751 
752  dimension dif_g(nrp,ntp)
753 
754  errg_max=0.d0
755 
756  do i=1,nr
757  do j=2,nt1
758 
759  dif_g(i,j)=g(i,j)-r(i,j)**2
760  dif_abs=dabs(dif_g(i,j))
761  errg_max=dmax1(errg_max,dif_abs)
762 
763  psi(i,j)=dif_g(i,j)
764 
765  enddo
766  enddo
767 
768  do i=1,nr
769  psi(i,1)=psi(i,nt1)
770  psi(i,nt)=psi(i,2)
771 
772  enddo
773 
774  call f_wrd
775 
776  return
777  end
778 
779  subroutine put_psib0(psi0_bnd)
780 
781  include 'double.inc'
782  include 'parevo.inc'
783  parameter(nkp=njlim)
784  include 'dim.inc'
785  include 'compol.inc'
786  include 'compol_add.inc'
787 
788  psibon0=psi0_bnd
789 
790  return
791  end
792 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
793 
794 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
795 
796 
797 
798 
f_rigbon
subroutine f_rigbon
Definition: _rig.f:394
f_rightp
subroutine f_rightp
Definition: _rig.f:260
f_wrd
subroutine f_wrd
Definition: _wrd.f:107
rightp_test
subroutine rightp_test
Definition: _rig.f:550
f_rightg
subroutine f_rightg
Definition: _rig.f:12
helena21_mod::dpdpsi
real(r8) function dpdpsi(psi_n)
Definition: helena21_mod.f90:122
numlin
function numlin(i, j, nr, nt)
Definition: com_sub.f:1040
g_test
subroutine g_test
Definition: _rig.f:741
f_psib_ext
subroutine f_psib_ext(psex_av)
Definition: _rig.f:676
avr_bnd
real *8 function avr_bnd(arr)
Definition: com_sub.f:221
rightg_test
subroutine rightg_test
Definition: _rig.f:653
pres_d_psi
subroutine pres_d_psi
Definition: interfac.f:2772
put_psib0
subroutine put_psib0(psi0_bnd)
Definition: _rig.f:779
toksur
subroutine toksur
Definition: _rig.f:447
rigext
subroutine rigext
Definition: _rig.f:344
f_psib_pla
subroutine f_psib_pla(pspl_av)
Definition: _rig.f:708
f_procof
subroutine f_procof(icq, cur_mu)
Definition: _metric.f:836