V3FIT
initialize_coilsin.f
1  SUBROUTINE initialize_coilsin
2  USE stel_kinds
3  USE coilsnamin
4 
5  IMPLICIT NONE
6 !-----------------------------------------------
7 ! L o c a l V a r i a b l e s
8 !-----------------------------------------------
9  INTEGER :: i
10 !-----------------------------------------------
11  bcoil_file = "bcoilxyz.dat"
12  lvf = .false.
13  lvfc = .true.
14  lcc_vf = .true.
15  lvfvar = .false.
16  lvfr = .false.
17  lvfz = .false.
18  ltfc = .false.
19  ltfcv = .false.
20  lbcoil = .false.
21  lbcoil_cur = .true.
22  lbnorm = .false.
23  lmodcur = .false.
24  lsadcur = .false.
25  lsadshape = .true.
26  lmodular = .false.
27  lsaddle = .true.
28  lsmod = .false.
29  lspline = .false.
30  lctrlpt = .false.
31  lsplbkp = .false.
32  lncsx = .false.
33  lqos = .false.
34  lsymm = .true.
35  laccess = .false.
36  laxis = .false.
37  lpolcur = .false.
38  lp_bg = .true.
39  lrestart = .false.
40  ls_cur = .true.
41 
42  nwdim = 400
43  nvf_fix = 2 !(=2 for QPS)
44  n_access = 0
45  nf_phi = 0
46  nf_rho = 0
47  nmod_coils_per_period = 0
48  nsad_coils_per_period = 0
49  nsad_u = 0
50  nsad_v = 0
51  numsurf = 0
52  numsurf_sad = 0
53  nstep = 10
54  niter_opt = 1
55  nopt_alg = 0
56  nopt_wsurf = -1
57  nfils = 1
58  nvar_vc = 1
59  nvar_uc = 1
60 
61  i_pol = 0
62  i_tfc = 0
63  rhoc = 0
64  rhos = 0
65  phic = 0
66  phis = 0
67  curmod = 0
68  dcc_wgt = 0
69  dcc_exp = 0
70  dcc_tgt = 0
71  rc_wgt = 0
72  rc_exp = 0
73  rc_tgt = 0
74  lmod_wgt = 0
75  lmod_tgt = 0
76  lsad_wgt = 0
77  lsad_tgt = 0
78  sad_v_c = 0
79  sad_v_s = 0
80  sad_u_c = 0
81  sad_u_s = 0
82  sad_v0 = 0
83  sad_u0 = 0
84  cursad = 0
85  cc_vf = 0
86  r_ext = 0
87  ymin_wgt = 0
88  ymin_tgt = 0
89  rmax_wgt = 0
90  rmax_tgt = 0
91  rs_wgt = 0
92  rs_exp = 0
93  rs_tgt = 0
94  dac_wgt = 0
95  dac_exp = 0
96  dac_tgt = 0
97  cvf_wgt = 0
98  cvf_tgt = 0
99  rvf_wgt = 0
100  rvf_tgt = 0
101  rc_vf = 0
102  zc_vf = 0
103  rcfc_vf = 0
104  rcfs_vf = 0
105  cs_wgt = 0
106  cs_tgt = 0
107  csc_wgt = 0
108  csc_tgt = 0
109  scd_wgt = 0
110  scd_tgt = 0
111  dpc_wgt = 0
112  vacfld_wgt = 0
113  mc_bg = 0
114  bcoil_cur = 0
115  dscxp_wgt = 0
116  dscxp_exp = 0
117  dscxp_tgt = 0
118  delt = 0.032_dp
119  deln = 0.055_dp ! NCSX value
120  bkp_wgt = 10000
121  bkp_tgt = 0.02_dp
122  mxb_wgt = 0
123  csad_scl = 1
124  nmid = 0
125  nsmid = 0
126  n_access = 0
127  num_vf = 0
128  nrvf_c = 0
129  dsc_wgt = 0
130  sc_dmin_wgt = 0
131 
132  DO i = 1, ncdim
133  nsad_group(i) = i
134  END DO
135 
136  END SUBROUTINE initialize_coilsin