V3FIT
•All Classes Namespaces Files Functions Variables Enumerations Macros Pages
vsetup.f
1  SUBROUTINE vsetup (iseq_count)
2  USE vmec_main
3  USE vacmod
4  USE realspace
5  USE mgrid_mod, ONLY: nbcoil_max, nlim_max, nextcur, mgrid_mode
6  USE gmres_mod, ONLY: nfcn
7 
8  IMPLICIT NONE
9 !-----------------------------------------------
10 ! D u m m y A r g u m e n t s
11 !C-----------------------------------------------
12  INTEGER, INTENT(IN) :: iseq_count
13 ! REAL(dp) :: vseton, vsetoff
14 C-----------------------------------------------
15 !
16 ! Reset default initial values
17 !
18 ! m=1 constraint (=t: apply correct, polar constraint; =f, apply approx. constraint)
19  lconm1 = .true.
20 ! lconm1 = .false.
21 ! IF (lrfp) lconm1 = .false. !SPH102109 converges better
22 
23 ! 2d preconditioner
24  nfcn = 0
25 
26  lrecon = .true.
27  loldout = .false.
28  ledge_dump = .false.
29 
30  z00 = zero
31  mgrid_mode = 'S' !Assume scaled mode; read in from mgrid in free-bdy mode
32  nextcur = 0
33 !
34 ! SOME DEFAULT VALUES WHICH MAY BE READ IN
35 !
36 
37 !
38 ! Reconstruction stuff
39 !
40  delbsq = one
41  lpofr = .true.
42  imse = -1
43  itse = 0
44  isnodes = 0
45  ipnodes = 0
46  iopt_raxis = 1
47  imatch_phiedge = 1
48  nflxs = 0
49  nbfld = 0
50  mseangle_offset = zero
51  mseangle_offsetm = zero
52  pres_offset = zero
53  sigma_current = 1.e30_dp
54  sigma_delphid = 1.e30_dp
55  tensi = one
56  tensp = one
57  tensi2 = zero
58  fpolyi = one
59  presfac = one
60  phidiam = 1.e30_dp
61 
62  mseprof = one
63  indxflx = 0
64  indxbfld = 0
65  sigma_stark = 1.1*cbig
66  sigma_thom = 1.1*cbig
67  sigma_flux = 1.1*cbig
68  sigma_b = 1.1*cbig
69 
70 !
71 ! FREE-BOUNDARY STUFF, ONLY INITIALIZED FIRST TIME
72 !
73  IF (iseq_count .eq. 0) THEN
74  nbcoil_max = 0
75  nlim_max = 0
76  END IF
77 
78  END SUBROUTINE vsetup
Generated on Thu Mar 5 2020 15:49:27 for V3FIT by   doxygen 1.8.17