Line data Source code
1 : !!****m* ABINIT/m_multibinit_dataset
2 : !! NAME
3 : !! m_multibinit_dataset
4 : !!
5 : !! FUNCTION
6 : !! module with the type for the input of multibinit (should be clean)
7 : !!
8 : !! COPYRIGHT
9 : !! Copyright (C) 2014-2026 ABINIT group (AM)
10 : !! This file is distributed under the terms of the
11 : !! GNU General Public License, see ~abinit/COPYING
12 : !! or http://www.gnu.org/copyleft/gpl.txt .
13 : !!
14 : !! SOURCE
15 :
16 : #if defined HAVE_CONFIG_H
17 : #include "config.h"
18 : #endif
19 :
20 : #include "abi_common.h"
21 :
22 : module m_multibinit_dataset
23 :
24 : use defs_basis
25 : use m_abicore
26 : use m_errors
27 :
28 : use m_xmpi
29 : use m_parser, only : intagm, instrng
30 : use m_fstrings, only : replace, inupper
31 : use m_dtset, only : chkvars
32 : use m_ddb, only : DDB_QTOL
33 : use m_scup_dataset
34 :
35 : implicit none
36 :
37 : private
38 :
39 : public :: multibinit_dtset_type
40 : public :: multibinit_dtset_init
41 : public :: multibinit_dtset_free
42 : public :: outvars_multibinit
43 : public :: invars_multibinit_filenames
44 : public :: invars_multibinit_filenames_from_input_file
45 : public :: invars10
46 : !!***
47 :
48 : !----------------------------------------------------------------------
49 :
50 : !!****t* m_multibinit_dataset/multibinit_dtset_type
51 : !! NAME
52 : !! multibinit_dtset_type
53 : !!
54 : !! FUNCTION
55 : !! The multibinit_dtset_type structured datatype
56 : !! gather all the input variables for the multibinit code.
57 : !!
58 : !! SOURCE
59 :
60 : type multibinit_dtset_type
61 :
62 : ! Integer
63 : integer :: asr
64 : integer :: analyze_anh_pot
65 : integer :: brav
66 : integer :: chneut
67 : integer :: confinement
68 : integer :: conf_power_disp
69 : integer :: conf_power_strain
70 : integer :: dipdip
71 : integer :: eivec
72 : integer :: elphflag
73 : integer :: enunit
74 : integer :: bound_model
75 : integer :: bound_maxCoeff
76 : integer :: bound_SPCoupling
77 : integer :: bound_AnhaStrain
78 : integer :: bound_step
79 : integer :: bound_option
80 : integer :: fit_anhaStrain
81 : integer :: fit_dispterms
82 : integer :: fit_SPCoupling
83 : integer :: fit_SPC_maxS
84 : integer :: fit_generateCoeff
85 : integer :: fit_iatom
86 : integer :: fit_initializeData
87 : integer :: fit_coeff
88 : integer :: fit_option
89 : integer :: fit_ncoeff
90 : integer :: fit_ncoeff_per_iatom
91 : integer :: fit_ncoeff_per_cycle
92 : integer :: fit_nbancoeff
93 : integer :: fit_nfixcoeff
94 : integer :: fit_nimposecoeff
95 : integer :: fit_EFS(3)
96 : integer :: sel_EFS(3)
97 : integer :: opt_EFS(3)
98 : integer :: bound_EFS(3)
99 : integer :: opt_effpot
100 : integer :: opt_ncoeff
101 : integer :: ts_option
102 : integer :: hmctt
103 : integer :: hmcsst
104 : integer :: ifcana
105 : integer :: ifcflag
106 : integer :: ifcout
107 : ! TODO hexu: why integer dtion?
108 : integer :: dtion
109 : integer :: dynamics
110 : integer :: dyn_chksym
111 : integer :: natifc
112 : integer :: natom
113 : integer :: ncoeff
114 : integer :: nctime
115 : integer :: nefield
116 : integer :: ntime
117 : integer :: nnos
118 : integer :: nph1l
119 : integer :: nph2l
120 : integer :: nqshft
121 : integer :: nsphere
122 : integer :: optcell
123 : integer :: prt_model
124 : integer :: dipdip_prt
125 : integer :: coeff_file_rw
126 : integer :: prt_phfrq
127 : integer :: prt_ifc
128 : integer :: randomseed
129 : integer :: strcpling ! Print the 3rd order in xml file
130 : integer :: prtsrlr ! print the short-range/long-range decomposition of phonon freq.
131 : integer :: rfmeth
132 : integer :: restartxf
133 : integer :: symdynmat
134 : integer :: sys_dim
135 : integer :: test_effpot
136 : integer :: test_prt_ph
137 : integer :: dipdip_range(3)
138 : integer :: fit_grid(3)
139 : integer :: fit_rangePower(2)
140 : integer :: bound_rangePower(2)
141 : integer :: bound_cell(3)
142 : integer :: ncell(3)
143 : integer :: ncellmat(3,3)
144 : integer :: ngqpt(9) ! ngqpt(9) instead of ngqpt(3) is needed in wght9.f
145 : integer :: ng2qpt(3)
146 : integer :: kptrlatt(3,3)
147 : integer :: kptrlatt_fine(3,3)
148 : integer :: qrefine(3)
149 : logical :: prt_GF_csv
150 :
151 : ! parameter for lattice
152 : integer :: latt_var_temperature
153 : integer :: latt_temperature_nstep
154 :
155 : ! parameter for hybrid lattice_lwf
156 : integer :: latt_lwf_anharmonic
157 :
158 :
159 : ! parameters for lwf
160 : integer :: lwf_constraint
161 : integer :: lwf_dynamics
162 : integer :: lwf_init_state
163 : integer :: lwf_ntime
164 : integer :: lwf_nctime
165 : !integer :: lwf_self_bound_order
166 : integer :: lwf_temperature_nstep ! var temperature number of steps
167 : integer :: lwf_var_temperature
168 :
169 : ! parameters for spin
170 : ! integer :: spin_calc_traj_obs
171 : integer :: spin_calc_thermo_obs
172 : !integer :: spin_calc_correlation_obs
173 : integer :: spin_dipdip
174 : integer :: spin_dynamics
175 : integer :: spin_init_state
176 : integer :: spin_nctime
177 : integer :: spin_ntime_pre
178 : integer :: spin_ntime
179 : ! integer :: spin_nmatom !TODO hexu: is it needed?
180 : ! integer :: spin_n1l
181 : ! integer :: spin_n2l
182 : integer :: spin_sia_add
183 : integer :: spin_temperature_nstep ! var temperature number of steps
184 : integer :: spin_var_temperature
185 : integer :: spin_write_traj
186 : integer :: efield_type
187 :
188 : ! parameters for spin-lattice coupling
189 : integer :: slc_coupling
190 :
191 : ! Real(dp)
192 : ! Inhomogeneous electric Field variables --Fernando Start
193 : real(dp) :: efield_period
194 : real(dp) :: efield_gmean(3)
195 : real(dp) :: efield_gvel(3)
196 : real(dp) :: efield_sigma
197 : real(dp) :: efield_background(3)
198 : ! Inhomogeneous electric Field variables --Fernando End
199 : real(dp) :: bmass
200 : real(dp) :: conf_power_fact_disp
201 : real(dp) :: conf_power_fact_strain
202 : real(dp) :: delta_df
203 : real(dp) :: dyn_tolsym
204 : real(dp) :: energy_reference
205 : !real(dp) :: efield(3)
206 : real(dp) :: bound_cutoff
207 : real(dp) :: bound_Temp
208 : real(dp) :: fit_cutoff
209 : ! real(dp):: fit_min_bound_coeff
210 : real(dp):: fit_drop_rate
211 : real(dp) :: fit_tolMSDF
212 : real(dp) :: fit_tolMSDS
213 : real(dp) :: fit_tolMSDE
214 : real(dp) :: fit_tolMSDFS
215 : real(dp) :: fit_tolGF
216 : real(dp) :: strprecon
217 : real(dp) :: strfact
218 : real(dp) :: temperature
219 : real(dp) :: rifcsph
220 : real(dp) :: conf
221 : real(dp) :: tolmxf
222 : real(dp) :: acell(3)
223 : real(dp) :: fit_factors(3)
224 : real(dp) :: fit_weight_T
225 : real(dp) :: opt_factors(3)
226 : real(dp) :: bound_factors(3)
227 : real(dp) :: bound_penalty
228 : real(dp) :: strten_reference(6)
229 : real(dp) :: strtarget(6)
230 : real(dp) :: conf_cutoff_strain(6)
231 : real(dp) :: rprim(3,3)
232 :
233 : ! lattice (new) related
234 : real(dp) :: latt_friction ! langevin dynamics friction
235 : real(dp) :: latt_taut ! Berendsen taut
236 : !real(dp) :: latt_taup !
237 : !real(dp) :: latt_compressibility
238 : !integer :: latt_mask(3)
239 : real(dp) :: latt_temperature_start
240 : real(dp) :: latt_temperature_end
241 :
242 : ! lwf related
243 : real(dp) :: lwf_dt
244 : real(dp) :: lwf_mc_avg_amp
245 : real(dp) :: lwf_taut
246 : real(dp) :: lwf_temperature
247 : !real(dp) :: lwf_self_bound_coeff
248 : real(dp) :: lwf_temperature_start ! var temperature start
249 : real(dp) :: lwf_temperature_end ! var temperature end
250 :
251 : ! parameters for spin
252 : real(dp) :: spin_dt
253 : real(dp) :: spin_damping
254 : real(dp) :: spin_sia_k1amp
255 : real(dp) :: spin_temperature
256 : ! TODO hexu: add spin convergence tol. (or remove it)
257 : real(dp) :: spin_temperature_start ! var temperature start
258 : real(dp) :: spin_temperature_end ! var temperature end
259 : !real(dp) :: spin_tolavg !average
260 : !real(dp) :: spin_tolvar !covariance
261 :
262 : real(dp) :: spin_mag_field(3) ! external magnetic field
263 : real(dp) :: spin_projection_qpoint(3) ! qpoint to check if spin configuration is random
264 : real(dp) :: spin_sia_k1dir(3)
265 : real(dp) :: spin_init_qpoint(3) ! qpoint to specify initial spin configuration
266 : real(dp) :: spin_init_rotate_axis(3) ! rotation axis to specify initial spin configuration
267 : real(dp) :: spin_init_orientation(3) ! spin orientation in primitive cell which is then rotated
268 :
269 : ! Integer arrays
270 : integer, allocatable :: atifc(:)
271 : ! atifc(natom)
272 : integer, allocatable :: fit_fixcoeff(:)
273 : ! fit_fixcoeffs(fit_nfixcoeff)
274 : integer, allocatable :: fit_imposecoeff(:)
275 : ! fit_fixcoeffs(fit_nimposecoeff)
276 :
277 : integer, allocatable :: fit_bancoeff(:)
278 : ! fit_bancoeffs(fit_nbancoeff)
279 :
280 : integer, allocatable :: iatfix(:,:)
281 : ! iatfix(3,natom) atom fix contraints for Broyden
282 :
283 : integer, allocatable :: fit_max_nbody_inrange(:)
284 : integer, allocatable :: fit_max_nbody(:)
285 :
286 :
287 : integer, allocatable :: opt_coeff(:)
288 : ! opt_coeff(opt_ncoeff)
289 :
290 : !integer, allocatable :: spin_sublattice(:) ! TODO hexu: difficult to use, better in xml?
291 :
292 : ! Logical array
293 : logical :: fit_on(3)
294 : ! fit_on(1) == TRUE, fit on energy, fit_on(2,3)=TRUE fit on forces stresses, fit_on(1,2,3)=TRUE fit on EFS
295 :
296 : logical :: sel_on(3)
297 : ! sel_on(1) == TRUE, select on energy, sel_on(2,3)=TRUE select on forces stresses, sel_on(1,2,3)=TRUE select on EFS
298 :
299 : logical :: opt_on(3)
300 : ! opt_on(1) == TRUE, optimize on energy, sel_on(2,3)=TRUE optimize on forces stresses, opt_on(1,2,3)=TRUE optimize on EFS
301 :
302 : real(dp), allocatable :: qmass(:)
303 : ! qmass(nnos)
304 :
305 :
306 :
307 : ! Real arrays
308 : real(dp), allocatable :: coefficients(:)
309 : ! coefficients(ncoeff)
310 :
311 : real(dp), allocatable :: conf_cutoff_disp(:)
312 : ! conf_cuttoff(natom)
313 :
314 : real(dp),allocatable :: q1shft(:,:)
315 : !q1shft(3,nqshft) SHIFT for Q point
316 :
317 : real(dp), allocatable :: qnrml1(:)
318 : ! qnrml1(nph1l)
319 :
320 : real(dp), allocatable :: qnrml2(:)
321 : ! qnrml1(nph1l)
322 :
323 : real(dp), allocatable :: qph1l(:,:)
324 : ! qph1l(3,nph1l)
325 :
326 : real(dp), allocatable :: qph2l(:,:)
327 : ! qph2l(3,nph2l)
328 :
329 : real(dp), allocatable :: efield(:, :)
330 : ! efield(3,nefield)
331 :
332 : real(dp), allocatable :: efield_lambda(:,:)
333 : ! efield_lambda(3,nefield)
334 :
335 : real(dp), allocatable :: efield_phase(:)
336 : ! efield_phase(nefield)
337 :
338 :
339 :
340 : !MS all variables for scale-up are put into their one type
341 : type(scup_dtset_type) :: scup_dtset
342 :
343 : ! spin part
344 : !real(dp), allocatable :: gilbert_damping ! if not provided in xml or override is needed.
345 : !real(dp), allocatable :: gyro_ratio(:) ! if not provided in xml
346 :
347 : !real(dp), allocatable :: qspin1l(:,:)
348 : !real(dp), allocatable :: qspin2l(:,:)
349 : ! characters
350 : character(len=fnlen) :: latt_init_hist_fname
351 : character(len=fnlen) :: latt_pot_fname
352 : character(len=fnlen) :: latt_harm_pot_fname
353 : character(len=fnlen) :: latt_anharm_pot_fname
354 : character(len=fnlen) :: latt_training_set_fname
355 : character(len=fnlen) :: latt_test_set_fname
356 : character(len=fnlen) :: latt_ddb_fnames(12)
357 :
358 : character(len=fnlen) :: lwf_init_hist_fname
359 : character(len=fnlen) :: lwf_pot_fname
360 :
361 : character(len=fnlen) :: outdata_prefix
362 :
363 : character(len=fnlen) :: spin_pot_fname
364 : character(len=fnlen) :: spin_init_hist_fname
365 :
366 : character(len=fnlen) :: slc_pot_fname
367 :
368 : end type multibinit_dtset_type
369 : !!***
370 :
371 : contains
372 : !!***
373 :
374 : !!****f* m_multibinit_dataset/multibinit_dtset_init
375 : !!
376 : !! NAME
377 : !! multibinit_dtset_init
378 : !!
379 : !! FUNCTION
380 : !! Init the dtset datatype
381 : !!
382 : !! INPUTS
383 : !! natom=number of atoms, needed for atifc
384 : !!
385 : !! OUTPUT
386 : !! multibinit_dtset <type(multibinit_dtset_type)> = datatype with all the input variables
387 : !!
388 : !! NOTES
389 : !! Should be executed by one processor only.
390 : !!
391 : !! SOURCE
392 :
393 0 : subroutine multibinit_dtset_init(multibinit_dtset,natom)
394 :
395 : !Arguments -------------------------------
396 : !scalars
397 : integer,intent(in) :: natom
398 : type(multibinit_dtset_type),intent(inout) :: multibinit_dtset
399 : !Local variables -------------------------
400 : !scalars
401 : !arrays
402 :
403 : !*********************************************************************
404 :
405 : !copy natom to multibinit_dtset
406 0 : multibinit_dtset%natom=natom
407 :
408 : !=====================================================================
409 : !Scalars
410 : !=====================================================================
411 0 : multibinit_dtset%asr=2
412 0 : multibinit_dtset%analyze_anh_pot=0
413 0 : multibinit_dtset%brav=1
414 0 : multibinit_dtset%bmass=0
415 : ! Inhomogeneous electric Field variables --Fernando Start
416 0 : multibinit_dtset%efield_period=1000000000000000.0
417 0 : multibinit_dtset%efield_gmean(:)=(/0.0,0.0,0.0/)
418 0 : multibinit_dtset%efield_gvel(:)=(/0.0,0.0,0.0/)
419 0 : multibinit_dtset%efield_type=0
420 0 : multibinit_dtset%efield_sigma=0.01
421 0 : multibinit_dtset%efield_background(:)= (/0.0,0.0,0.0/)
422 : ! Inhomogeneous electric Field variables --Fernando End
423 0 : multibinit_dtset%chneut=0
424 0 : multibinit_dtset%confinement=0
425 0 : multibinit_dtset%conf_power_disp=0
426 0 : multibinit_dtset%conf_power_strain=0
427 0 : multibinit_dtset%conf_power_fact_disp=100
428 0 : multibinit_dtset%conf_power_fact_strain=100
429 0 : multibinit_dtset%delta_df= 1d-02
430 0 : multibinit_dtset%dipdip=1
431 0 : multibinit_dtset%dipdip_prt=0
432 0 : multibinit_dtset%dtion=100
433 0 : multibinit_dtset%dynamics=0
434 0 : multibinit_dtset%dyn_chksym=0
435 0 : multibinit_dtset%dyn_tolsym=1d-10
436 0 : multibinit_dtset%eivec=0
437 0 : multibinit_dtset%energy_reference= zero
438 0 : multibinit_dtset%enunit=0
439 0 : multibinit_dtset%fit_anhaStrain=0
440 0 : multibinit_dtset%bound_model=0
441 0 : multibinit_dtset%bound_penalty=1.001d+0
442 0 : multibinit_dtset%bound_anhaStrain=0
443 0 : multibinit_dtset%bound_cutoff=0
444 0 : multibinit_dtset%bound_maxCoeff=4
445 0 : multibinit_dtset%bound_temp=325
446 0 : multibinit_dtset%bound_step=1000
447 0 : multibinit_dtset%bound_SPCoupling=1
448 0 : multibinit_dtset%fit_coeff=0
449 : ! multibinit_dtset%fit_min_bound_coeff=0.0_dp
450 0 : multibinit_dtset%fit_drop_rate=0.0_dp
451 0 : multibinit_dtset%fit_cutoff=0
452 0 : multibinit_dtset%fit_nbancoeff=0
453 0 : multibinit_dtset%fit_ncoeff=0
454 0 : multibinit_dtset%fit_ncoeff_per_iatom=0
455 0 : multibinit_dtset%fit_ncoeff_per_cycle=1
456 0 : multibinit_dtset%fit_weight_T=-0.1_dp
457 :
458 0 : multibinit_dtset%fit_iatom=0
459 0 : multibinit_dtset%ts_option=0
460 0 : multibinit_dtset%fit_nfixcoeff=0
461 0 : multibinit_dtset%fit_nimposecoeff=0
462 0 : multibinit_dtset%fit_EFS = (/ 0, 1, 1 /)
463 0 : multibinit_dtset%sel_EFS = (/ 0, 1, 1 /)
464 0 : multibinit_dtset%opt_EFS = (/ 0, 1, 1 /)
465 0 : multibinit_dtset%bound_EFS = (/ 0, 1, 1 /)
466 0 : multibinit_dtset%fit_option=0
467 0 : multibinit_dtset%fit_SPCoupling=1
468 0 : multibinit_dtset%fit_dispterms=1
469 0 : multibinit_dtset%fit_SPC_maxS=1
470 0 : multibinit_dtset%fit_generateCoeff=1
471 0 : multibinit_dtset%fit_initializeData=0
472 0 : multibinit_dtset%fit_tolMSDE=zero
473 0 : multibinit_dtset%fit_tolMSDS=zero
474 0 : multibinit_dtset%fit_tolMSDF=zero
475 0 : multibinit_dtset%fit_tolMSDFS=zero
476 0 : multibinit_dtset%fit_tolGF=zero
477 0 : multibinit_dtset%ifcana=0
478 0 : multibinit_dtset%ifcflag=1
479 0 : multibinit_dtset%ifcout=-1
480 0 : multibinit_dtset%hmctt=0
481 0 : multibinit_dtset%hmcsst=0
482 0 : multibinit_dtset%prtsrlr=0
483 : ! Langevin friction
484 0 : multibinit_dtset%latt_friction=1d-4
485 : ! Berendsen taut
486 0 : multibinit_dtset%latt_taut=1000.0
487 : !multibinit_dtset%latt_taup=1000.0
488 : !multibinit_dtset%latt_compressibility=0.0
489 :
490 :
491 0 : multibinit_dtset%nefield=1
492 0 : multibinit_dtset%ntime=200
493 0 : multibinit_dtset%nctime=1
494 0 : multibinit_dtset%natifc=natom
495 0 : multibinit_dtset%ncoeff=0
496 0 : multibinit_dtset%nph1l=1
497 0 : multibinit_dtset%nph2l=0
498 0 : multibinit_dtset%nqshft=1
499 0 : multibinit_dtset%nnos=0
500 0 : multibinit_dtset%nsphere=0
501 0 : multibinit_dtset%optcell=0
502 0 : multibinit_dtset%opt_effpot=0
503 : !multibinit_dtset%opt_coeff=0
504 0 : multibinit_dtset%prt_model=0
505 0 : multibinit_dtset%prt_phfrq=0
506 0 : multibinit_dtset%prt_ifc = 0
507 0 : multibinit_dtset%prt_GF_csv = .FALSE.
508 0 : multibinit_dtset%randomseed = 0
509 0 : multibinit_dtset%strcpling = -1
510 0 : multibinit_dtset%qrefine=1
511 0 : multibinit_dtset%restartxf=0
512 0 : multibinit_dtset%rfmeth=1
513 0 : multibinit_dtset%rifcsph=zero
514 0 : multibinit_dtset%strfact=100.0d0
515 0 : multibinit_dtset%strprecon=1.0d0
516 0 : multibinit_dtset%symdynmat=1
517 0 : multibinit_dtset%sys_dim=1
518 0 : multibinit_dtset%temperature=325
519 0 : multibinit_dtset%test_effpot=0
520 0 : multibinit_dtset%test_prt_ph=0
521 0 : multibinit_dtset%tolmxf=2.0d-5
522 0 : multibinit_dtset%latt_temperature_start=0.0
523 0 : multibinit_dtset%latt_temperature_end=0.0
524 0 : multibinit_dtset%latt_temperature_nstep=0
525 0 : multibinit_dtset%latt_var_temperature=0
526 :
527 0 : multibinit_dtset%latt_lwf_anharmonic = 0
528 :
529 0 : multibinit_dtset%lwf_constraint = 0
530 0 : multibinit_dtset%lwf_dynamics = 0
531 0 : multibinit_dtset%lwf_nctime = 1
532 0 : multibinit_dtset%lwf_ntime = 0
533 0 : multibinit_dtset%lwf_init_state = 0
534 0 : multibinit_dtset%lwf_init_hist_fname=""
535 :
536 0 : multibinit_dtset%latt_init_hist_fname=""
537 0 : multibinit_dtset%latt_pot_fname=""
538 0 : multibinit_dtset%latt_harm_pot_fname=""
539 0 : multibinit_dtset%latt_anharm_pot_fname=""
540 0 : multibinit_dtset%latt_training_set_fname=""
541 0 : multibinit_dtset%latt_test_set_fname=""
542 0 : multibinit_dtset%latt_ddb_fnames(12)=""
543 :
544 :
545 0 : multibinit_dtset%lwf_pot_fname=""
546 0 : multibinit_dtset%slc_pot_fname=""
547 0 : multibinit_dtset%lwf_dt=0
548 : !multibinit_dtset%lwf_self_bound_order=0
549 0 : multibinit_dtset%lwf_taut=0.0_dp
550 0 : multibinit_dtset%lwf_temperature=0.0_dp
551 0 : multibinit_dtset%lwf_mc_avg_amp=0.0_dp
552 : !multibinit_dtset%lwf_self_bound_coeff=0.0_dp
553 0 : multibinit_dtset%lwf_temperature_start=0.0
554 0 : multibinit_dtset%lwf_temperature_end= 0.0
555 0 : multibinit_dtset%lwf_temperature_nstep= 0
556 0 : multibinit_dtset%lwf_var_temperature=0
557 :
558 0 : multibinit_dtset%spin_init_hist_fname=""
559 0 : multibinit_dtset%spin_pot_fname=""
560 :
561 : !multibinit_dtset%spin_calc_traj_obs=0
562 0 : multibinit_dtset%spin_calc_thermo_obs=1
563 : !multibinit_dtset%spin_calc_correlation_obs=0
564 0 : multibinit_dtset%spin_dipdip=0
565 0 : multibinit_dtset%spin_dt=100
566 0 : multibinit_dtset%spin_dynamics=0
567 0 : multibinit_dtset%spin_init_hist_fname=""
568 0 : multibinit_dtset%spin_init_state=1
569 0 : multibinit_dtset%spin_ntime_pre=0
570 0 : multibinit_dtset%spin_ntime=10000
571 0 : multibinit_dtset%spin_nctime=100
572 : !multibinit_dtset%spin_nmatom=0
573 : !multibinit_dtset%spin_n1l=1
574 : !multibinit_dtset%spin_n2l=0
575 :
576 :
577 0 : multibinit_dtset%spin_damping=-1.0
578 0 : multibinit_dtset%spin_sia_add=0
579 0 : multibinit_dtset%spin_sia_k1amp=zero
580 0 : multibinit_dtset%spin_temperature=325
581 0 : multibinit_dtset%spin_temperature_start=0.0
582 0 : multibinit_dtset%spin_temperature_end= 0.0
583 0 : multibinit_dtset%spin_temperature_nstep= 0
584 : !multibinit_dtset%spin_tolavg=1d-2 ! TODO hexu: to be decided. should it be a function of temperature?
585 : !multibinit_dtset%spin_tolvar=1d-3 ! TODO hexu: as above.
586 :
587 0 : multibinit_dtset%spin_var_temperature=0
588 0 : multibinit_dtset%spin_write_traj=1
589 :
590 0 : multibinit_dtset%slc_pot_fname=""
591 0 : multibinit_dtset%slc_coupling=0
592 :
593 0 : multibinit_dtset%outdata_prefix=""
594 : !=======================================================================
595 : !Arrays
596 : !=======================================================================
597 0 : multibinit_dtset%acell(:) = one
598 0 : multibinit_dtset%conf_cutoff_strain(1:6) = zero
599 0 : multibinit_dtset%dipdip_range(:)= (/0,0,0/)
600 : !multibinit_dtset%efield(:)= (/0,0,0/)
601 0 : multibinit_dtset%fit_grid(:)= 1
602 0 : multibinit_dtset%fit_rangePower(:)= (/3,4/)
603 0 : multibinit_dtset%bound_rangePower(:)= (/6,6/)
604 0 : multibinit_dtset%bound_cell(:)= (/6,6,6/)
605 0 : multibinit_dtset%ncell(:)= 0
606 0 : multibinit_dtset%ncellmat(:,:)= 0
607 0 : multibinit_dtset%ngqpt(:) = 0
608 0 : multibinit_dtset%ng2qpt(:)= 0
609 0 : multibinit_dtset%strtarget(1:6) = zero
610 0 : multibinit_dtset%qmass(:)= zero
611 0 : multibinit_dtset%rprim(:,:)= zero
612 0 : multibinit_dtset%strten_reference(:)= zero
613 :
614 0 : multibinit_dtset%spin_mag_field(:)=zero
615 0 : multibinit_dtset%spin_projection_qpoint(:)=zero
616 0 : multibinit_dtset%spin_init_qpoint(:)=zero
617 0 : multibinit_dtset%spin_init_rotate_axis(:)=(/1.0, 0.0, 0.0/)
618 0 : multibinit_dtset%spin_init_orientation(:)=(/0.0, 0.0, 1.0/)
619 :
620 :
621 0 : multibinit_dtset%spin_sia_k1dir(:)=(/0.0,0.0,1.0/)
622 :
623 :
624 0 : ABI_MALLOC(multibinit_dtset%atifc,(natom))
625 0 : multibinit_dtset%atifc(:)=0
626 0 : ABI_MALLOC(multibinit_dtset%conf_cutoff_disp,(multibinit_dtset%natom))
627 0 : multibinit_dtset%conf_cutoff_disp(:)=zero
628 0 : ABI_MALLOC(multibinit_dtset%q1shft,(3,multibinit_dtset%nqshft))
629 0 : multibinit_dtset%q1shft(:,:) = zero
630 0 : ABI_MALLOC(multibinit_dtset%iatfix,(3,natom))
631 :
632 : ! multibinit_dtset%latt_mask(:) = 0
633 :
634 0 : end subroutine multibinit_dtset_init
635 : !!***
636 :
637 : !!****f* m_multibinit_dataset/multibinit_dtset_free
638 : !!
639 : !! NAME
640 : !! multibinit_dtset_free
641 : !!
642 : !! FUNCTION
643 : !! deallocate remaining arrays in the multibinit_dtset datastructure
644 : !!
645 : !! INPUTS
646 : !! multibinit_dtset <type(multibinit_dtset_type)> = multibinit_dataset structure
647 : !!
648 : !! OUTPUTS
649 : !! multibinit_dtset <type(multibinit_dtset_type)> = multibinit_dataset structure
650 : !!
651 : !! SOURCE
652 :
653 73 : subroutine multibinit_dtset_free(multibinit_dtset)
654 :
655 : !Arguments ------------------------------------
656 : !scalars
657 : type(multibinit_dtset_type), intent(inout) :: multibinit_dtset
658 :
659 : ! *************************************************************************
660 :
661 73 : ABI_SFREE(multibinit_dtset%atifc)
662 73 : ABI_SFREE(multibinit_dtset%conf_cutoff_disp)
663 73 : ABI_SFREE(multibinit_dtset%fit_fixcoeff)
664 73 : ABI_SFREE(multibinit_dtset%fit_imposecoeff)
665 73 : ABI_SFREE(multibinit_dtset%fit_bancoeff)
666 73 : ABI_SFREE(multibinit_dtset%opt_coeff)
667 :
668 73 : ABI_SFREE(multibinit_dtset%fit_max_nbody)
669 73 : ABI_SFREE(multibinit_dtset%fit_max_nbody_inrange)
670 :
671 73 : ABI_SFREE(multibinit_dtset%qmass)
672 73 : ABI_SFREE(multibinit_dtset%coefficients)
673 73 : ABI_SFREE(multibinit_dtset%qnrml1)
674 73 : ABI_SFREE(multibinit_dtset%qnrml2)
675 73 : ABI_SFREE(multibinit_dtset%qph1l)
676 73 : ABI_SFREE(multibinit_dtset%qph2l)
677 :
678 73 : ABI_SFREE(multibinit_dtset%q1shft)
679 73 : ABI_SFREE(multibinit_dtset%iatfix)
680 :
681 73 : ABI_SFREE(multibinit_dtset%efield)
682 73 : ABI_SFREE(multibinit_dtset%efield_lambda)
683 73 : ABI_SFREE(multibinit_dtset%efield_phase)
684 :
685 :
686 : !ABI_SFREE(multibinit_dtset%gilbert_damping)
687 : !ABI_SFREE(multibinit_dtset%gyro_ratio)
688 : !ABI_SFREE(multibinit_dtset%qph1l_spin)
689 : !ABI_SFREE(multibinit_dtset%qph2l_spin)
690 :
691 : !=====================================================================
692 : !Free Scale-up dataset
693 : !=====================================================================
694 :
695 73 : call scup_dtset_free(multibinit_dtset%scup_dtset)
696 :
697 73 : end subroutine multibinit_dtset_free
698 : !!***
699 :
700 : !----------------------------------------------------------------------
701 : !===============================================================
702 : ! Parsing of input variables for Multibinit
703 : ! read only the latt_inp_ddb_fname, and the outdata_prefix.
704 : ! It contains the reference structure, so that the natom can be decided
705 : !> @ fname: the name of the ddb file
706 : !> @ string: the input file in string
707 : !> @ lenstr: the length of string
708 : !===============================================================
709 73 : subroutine invars_multibinit_filenames( string, lenstr, outdata_prefix, sys_fname )
710 : character(len=fnlen), optional, intent(inout) :: outdata_prefix, sys_fname
711 : character(len=*), intent(inout) :: string
712 : integer, intent(in) :: lenstr
713 : integer :: jdtset,marr,tread
714 : !arrays
715 73 : integer,allocatable :: intarr(:)
716 73 : real(dp),allocatable :: dprarr(:)
717 73 : marr=300
718 73 : ABI_MALLOC(intarr,(marr))
719 73 : ABI_MALLOC(dprarr,(marr))
720 73 : jdtset=1
721 :
722 73 : if(present(outdata_prefix)) then
723 73 : outdata_prefix=""
724 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'outdata_prefix',tread,'KEY',&
725 73 : & key_value=outdata_prefix)
726 73 : if(.not. tread==1) outdata_prefix=""
727 : end if
728 :
729 73 : if(present(sys_fname)) then
730 73 : sys_fname=""
731 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_harm_pot_fname',tread,'KEY',&
732 73 : & key_value=sys_fname)
733 73 : if(.not. tread==1) then
734 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_pot_fname',tread,'KEY',&
735 6 : & key_value=sys_fname)
736 6 : if(.not. tread==1) then
737 1 : sys_fname=""
738 : end if
739 : end if
740 : end if
741 73 : ABI_SFREE(intarr)
742 73 : ABI_SFREE(dprarr)
743 73 : end subroutine invars_multibinit_filenames
744 :
745 :
746 73 : subroutine invars_multibinit_filenames_from_input_file(input_path, outdata_prefix, sys_fname)
747 : character(len=fnlen), optional, intent(inout) :: outdata_prefix, sys_fname, input_path
748 : character(len=strlen) :: string, raw_string
749 : integer :: lenstr, option
750 :
751 : integer :: master, my_rank, comm, nproc, ierr
752 : logical :: iam_master
753 73 : master = 0
754 73 : comm = xmpi_world
755 73 : nproc = xmpi_comm_size(comm)
756 73 : my_rank = xmpi_comm_rank(comm)
757 73 : iam_master = (my_rank == master)
758 :
759 :
760 73 : option=1
761 73 : if (iam_master) then
762 49 : call instrng (input_path,lenstr,option,strlen,string, raw_string)
763 : !To make case-insensitive, map characters to upper case:
764 49 : call inupper(string(1:lenstr))
765 : !Check whether the string only contains valid keywords
766 49 : call chkvars(string)
767 : end if
768 73 : call xmpi_bcast(string, master, comm, ierr)
769 73 : call xmpi_bcast(raw_string, master, comm, ierr)
770 73 : call xmpi_bcast(lenstr, master, comm, ierr)
771 :
772 : ! Save input string in global variable so that we can access it in ntck_open_create
773 73 : ABI_MALLOC_TYPE_SCALAR(character(len=len_trim(string)), INPUT_STRING)
774 73 : INPUT_STRING = string(1:len_trim(string))
775 :
776 73 : call invars_multibinit_filenames( string=string, lenstr=lenstr, sys_fname=sys_fname, outdata_prefix=outdata_prefix)
777 73 : end subroutine invars_multibinit_filenames_from_input_file
778 :
779 : !!****f* m_multibinit_dataset/invars10
780 : !!
781 : !! NAME
782 : !! invars10
783 : !!
784 : !! FUNCTION
785 : !! Open input file for the multibinit code, then reads or echoes the input information.
786 : !!
787 : !! INPUTS
788 : !! lenstr=actual length of string
789 : !! natom=number of atoms, needed for atifc
790 : !! string*(*)=string of characters containing all input variables and data
791 : !!
792 : !! OUTPUT
793 : !! multibinit_dtset <type(multibinit_dtset_type)> = datatype with all the input variables
794 : !!
795 : !! NOTES
796 : !! Should be executed by one processor only.
797 : !!
798 : !! SOURCE
799 :
800 73 : subroutine invars10(multibinit_dtset,lenstr,natom,string)
801 :
802 : !Arguments -------------------------------
803 : !scalars
804 : integer,intent(in) :: lenstr,natom
805 : character(len=*),intent(in) :: string
806 : type(multibinit_dtset_type),intent(inout) :: multibinit_dtset
807 :
808 : !Local variables -------------------------
809 : !Dummy arguments for subroutine 'intagm' to parse input file
810 : !Set routine version number here:
811 : !scalars
812 : integer :: iatifc,ii,iph1,iph2,jdtset,jj,marr,tread,idir,natfix,iatom
813 : integer :: natom_sc
814 : character(len=500) :: message
815 : !arrays
816 : integer,allocatable :: intarr(:)
817 73 : real(dp),allocatable :: dprarr(:),work(:)
818 : ! strings
819 : character(len=fnlen*12) :: lattddb_string
820 : integer :: sidx(13), cnt, i1, i2
821 :
822 73 : real(dp), allocatable :: dptmp(:)
823 :
824 : !*********************************************************************
825 73 : marr=30
826 73 : ABI_MALLOC(intarr,(marr))
827 73 : ABI_MALLOC(dprarr,(marr))
828 :
829 73 : jdtset=1
830 :
831 : !copy natom to multibinit_dtset
832 73 : multibinit_dtset%natom=natom
833 :
834 : !=====================================================================
835 : !start reading in dimensions and non-dependent variables
836 : !=====================================================================
837 :
838 : !A
839 73 : multibinit_dtset%asr=2
840 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'asr',tread,'INT')
841 73 : if(tread==1) multibinit_dtset%asr=intarr(1)
842 73 : if(multibinit_dtset%asr<-2.or.multibinit_dtset%asr>6)then
843 : write(message, '(a,i8,a,a,a,a,a)' )&
844 0 : & 'asr is',multibinit_dtset%asr,', but the only allowed values',ch10,&
845 0 : & 'are 0, 1, 2, 3, 4, 5, 6, -1 or -2 .',ch10,&
846 0 : & 'Action: correct asr in your input file.'
847 : ! Note : negative values are allowed when the acoustic sum rule
848 : ! is to be applied after the analysis of IFCs
849 : ! 3,4 are for rotational invariance (under development)
850 : ! 5 is for hermitian imposition of the ASR
851 0 : ABI_ERROR(message)
852 : end if
853 :
854 73 : multibinit_dtset%analyze_anh_pot=0
855 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'analyze_anh_pot',tread,'INT')
856 73 : if(tread==1) multibinit_dtset%analyze_anh_pot=intarr(1)
857 73 : if(multibinit_dtset%analyze_anh_pot < 0 .or. multibinit_dtset%analyze_anh_pot > 1)then
858 : write(message, '(a,i8,a,a,a,a,a)' )&
859 0 : & 'analyze_anh_pot is',multibinit_dtset%analyze_anh_pot,', but the only allowed values',ch10,&
860 0 : & 'are 0 and 1 .',ch10,&
861 0 : & 'Action: correct analyze_anh_pot in your input file.'
862 0 : ABI_ERROR(message)
863 : end if
864 :
865 : !B
866 73 : multibinit_dtset%brav=1
867 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'brav',tread,'INT')
868 73 : if(tread==1) multibinit_dtset%brav=intarr(1)
869 73 : if(multibinit_dtset%brav/=1)then
870 : write(message, '(a,i8,a,a,a,a,a)' )&
871 0 : & 'brav is',multibinit_dtset%brav,', but the only allowed values',ch10,&
872 0 : & 'are 1 for multibinit (not implemented) .',ch10,&
873 0 : & 'Action: correct brav in your input file.'
874 0 : ABI_ERROR(message)
875 : end if
876 :
877 73 : multibinit_dtset%bmass=0
878 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'bmass',tread,'DPR')
879 73 : if(tread==1) multibinit_dtset%bmass=dprarr(1)
880 73 : if(multibinit_dtset%bmass<0)then
881 : write(message, '(a,f10.2,a,a,a,a,a)' )&
882 0 : & 'bmass is',multibinit_dtset%bmass,', but the only allowed values',ch10,&
883 0 : & 'is superior to 0.',ch10,&
884 0 : & 'Action: correct bmass in your input file.'
885 0 : ABI_ERROR(message)
886 : end if
887 :
888 : ! Inhomogeneous electric Field variables --Fernando Start
889 :
890 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'efield_gmean',tread,'TIM')
891 76 : if(tread==1) multibinit_dtset%efield_gmean(1:3)= dprarr(1:3)
892 292 : do ii=1,3
893 292 : if(multibinit_dtset%efield_gmean(ii)<0)then
894 : write(message, '(a,i0,a,i0,4a,i0,a)' )&
895 0 : & 'efield_gmean(',ii,') is ',multibinit_dtset%efield_gmean(ii),', which is lower',&
896 0 : & ' than 0',&
897 0 : & ch10,'Action: correct efield_gmean(',ii,') in your input file.'
898 0 : ABI_ERROR(message)
899 : end if
900 : end do
901 292 : multibinit_dtset%efield_gvel= (/0.0,0.0,0.0/)
902 73 : if(3>marr)then
903 0 : marr=3
904 0 : ABI_FREE(intarr)
905 0 : ABI_FREE(dprarr)
906 0 : ABI_MALLOC(intarr,(marr))
907 0 : ABI_MALLOC(dprarr,(marr))
908 : end if
909 :
910 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'efield_gvel',tread,'TIM')
911 76 : if(tread==1) multibinit_dtset%efield_gvel(1:3)= dprarr(1:3)
912 292 : do ii=1,3
913 292 : if(multibinit_dtset%efield_gvel(ii)<0)then
914 : write(message, '(a,i0,a,i0,4a,i0,a)' )&
915 0 : & 'efield_gvel(',ii,') is ',multibinit_dtset%efield_gvel(ii),', which is lower',&
916 0 : & ' than 0',&
917 0 : & ch10,'Action: correct efield_gvel(',ii,') in your input file.'
918 0 : ABI_ERROR(message)
919 : end if
920 : end do
921 :
922 73 : multibinit_dtset%efield_type=0
923 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'efield_type',tread,'INT')
924 73 : if(tread==1) multibinit_dtset%efield_type=intarr(1)
925 73 : if(multibinit_dtset%efield_type<0.or.multibinit_dtset%efield_type>6)then
926 : write(message, '(a,i8,a,a,a,a,a)' )&
927 0 : & 'efield_type is',multibinit_dtset%efield_type,', but the only allowed values',ch10,&
928 0 : & 'are 0, 1 or 2, 3, 4, 5 or 6 .',ch10,&
929 0 : & 'Action: correct efield_type in your input file.'
930 0 : ABI_ERROR(message)
931 : end if
932 :
933 73 : if (multibinit_dtset%efield_type/=6 .and. multibinit_dtset%efield_type/=0)then
934 5 : multibinit_dtset%nefield=1
935 68 : else if (multibinit_dtset%efield_type==6)then
936 1 : multibinit_dtset%nefield=2
937 1 : call read_int_var( multibinit_dtset%nefield, 'nefield')
938 1 : if (multibinit_dtset%nefield/=2)then
939 : write(message, '(a,i0,a,a,a)' )&
940 0 : & 'nefield is',multibinit_dtset%nefield, &
941 0 : & ', which can only be 2 for efield_type=6 .',ch10,&
942 0 : & 'Action: correct nefield in your input file.'
943 0 : ABI_ERROR(message)
944 : end if
945 :
946 : end if
947 : ! read efield
948 : call read_dpr_array_var(dptmp, size=multibinit_dtset%nefield*3, &
949 73 : & var_name='efield', type='DPR', default=0.0_dp)
950 219 : ABI_MALLOC(multibinit_dtset%efield, (3,multibinit_dtset%nefield))
951 320 : multibinit_dtset%efield= reshape(dptmp, (/3,multibinit_dtset%nefield/))
952 73 : ABI_FREE(dptmp)
953 :
954 :
955 : ! read efield_phase
956 : call read_dpr_array_var( multibinit_dtset%efield_phase, &
957 73 : & size=multibinit_dtset%nefield, var_name='efield_phase', type='DPR', default=0.0_dp)
958 :
959 : ! read efield_lambda
960 : call read_dpr_array_var(dptmp, size=multibinit_dtset%nefield*3, var_name='efield_lambda', &
961 73 : & type='DPR', default=0.0_dp)
962 219 : ABI_MALLOC(multibinit_dtset%efield_lambda, (3,multibinit_dtset%nefield))
963 320 : multibinit_dtset%efield_lambda= reshape(dptmp, (/3,multibinit_dtset%nefield/))
964 73 : ABI_FREE(dptmp)
965 :
966 :
967 :
968 73 : multibinit_dtset%efield_period=1000000000000000.0
969 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'efield_period',tread,'TIM')
970 73 : if(tread==1) multibinit_dtset%efield_period=dprarr(1)
971 73 : if(multibinit_dtset%efield_period<0)then
972 : write(message, '(a,f10.2,a,a,a,a,a)' )&
973 0 : & 'efield_period is',multibinit_dtset%efield_period,', but the only allowed values',ch10,&
974 0 : & 'is superior to 0.',ch10,&
975 0 : & 'Action: correct efield_period in your input file.'
976 0 : ABI_ERROR(message)
977 : end if
978 :
979 73 : multibinit_dtset%efield_sigma=0.01
980 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'efield_sigma',tread,'TIM')
981 73 : if(tread==1) multibinit_dtset%efield_sigma=dprarr(1)
982 73 : if(multibinit_dtset%efield_period<0)then
983 : write(message, '(a,f10.2,a,a,a,a,a)' )&
984 0 : & 'efield_sigma is',multibinit_dtset%efield_sigma,', but the only allowed values',ch10,&
985 0 : & 'is superior to 0.',ch10,&
986 0 : & 'Action: correct efield_sigma in your input file.'
987 0 : ABI_ERROR(message)
988 : end if
989 292 : multibinit_dtset%efield_background= zero
990 73 : if(3>marr)then
991 0 : marr=3
992 0 : ABI_FREE(intarr)
993 0 : ABI_FREE(dprarr)
994 0 : ABI_MALLOC(intarr,(marr))
995 0 : ABI_MALLOC(dprarr,(marr))
996 : end if
997 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'efield_background',tread,'DPR')
998 88 : if(tread==1) multibinit_dtset%efield_background(1:3)= dprarr(1:3)
999 :
1000 : ! Inhomogeneous electric Field variables --Fernando End
1001 :
1002 292 : multibinit_dtset%bound_EFS=(/0,1,1/)
1003 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'bound_EFS',tread,'INT')
1004 73 : if(tread==1) multibinit_dtset%bound_EFS(1:3)=intarr(1:3)
1005 584 : if(any(multibinit_dtset%bound_EFS<0) .or. any(multibinit_dtset%bound_EFS>1))then
1006 : write(message, '(a,i8,a,a,a)' )&
1007 0 : & 'bound_EFS is',multibinit_dtset%bound_EFS,', but the only allowed values are 0 and 1',ch10,&
1008 0 : & 'Action: correct bound_EFS in your input file.'
1009 0 : ABI_ERROR(message)
1010 : end if
1011 :
1012 292 : multibinit_dtset%bound_factors=(/1,1,1/)
1013 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'bound_factors',tread,'DPR')
1014 76 : if(tread==1) multibinit_dtset%bound_factors(1:3)=dprarr(1:3)
1015 292 : if(any(multibinit_dtset%bound_factors<0))then
1016 : write(message, '(a,i8,a,a,a)' )&
1017 0 : & 'bound_factors is',multibinit_dtset%bound_factors,', but the only allowed values are positive',ch10,&
1018 0 : & 'Action: correct bound_factors in your input file.'
1019 0 : ABI_ERROR(message)
1020 : end if
1021 :
1022 :
1023 : !C
1024 73 : multibinit_dtset%chneut=0
1025 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'chneut',tread,'INT')
1026 73 : if(tread==1) multibinit_dtset%chneut=intarr(1)
1027 73 : if(multibinit_dtset%chneut<0.or.multibinit_dtset%chneut>2)then
1028 : write(message, '(a,i8,a,a,a,a,a)' )&
1029 0 : & 'chneut is',multibinit_dtset%chneut,', but the only allowed values',ch10,&
1030 0 : & 'are 0, 1 or 2 .',ch10,&
1031 0 : & 'Action: correct chneut in your input file.'
1032 0 : ABI_ERROR(message)
1033 : end if
1034 :
1035 73 : multibinit_dtset%confinement=0
1036 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'confinement',tread,'INT')
1037 73 : if(tread==1) multibinit_dtset%confinement=intarr(1)
1038 73 : if(multibinit_dtset%confinement<0.or.multibinit_dtset%confinement>2)then
1039 : write(message, '(a,i8,a,a,a,a,a)' )&
1040 0 : & 'confinement is',multibinit_dtset%confinement,', but the only allowed values',ch10,&
1041 0 : & 'are 0, 1 or 2 .',ch10,&
1042 0 : & 'Action: correct confinement in your input file.'
1043 0 : ABI_ERROR(message)
1044 : end if
1045 :
1046 73 : multibinit_dtset%conf_power_disp=0
1047 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'conf_power_disp',tread,'INT')
1048 73 : if(tread==1) multibinit_dtset%conf_power_disp=intarr(1)
1049 73 : if(multibinit_dtset%conf_power_disp<0)then
1050 : write(message, '(a,i8,a,a,a,a,a)' )&
1051 0 : & 'conf_power_disp is',multibinit_dtset%conf_power_disp,', but the only allowed values',ch10,&
1052 0 : & 'positive .',ch10,&
1053 0 : & 'Action: correct conf_power_disp in your input file.'
1054 0 : ABI_ERROR(message)
1055 : end if
1056 :
1057 73 : multibinit_dtset%conf_power_strain=0
1058 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'conf_power_strain',tread,'INT')
1059 73 : if(tread==1) multibinit_dtset%conf_power_strain=intarr(1)
1060 73 : if(multibinit_dtset%conf_power_strain<0)then
1061 : write(message, '(a,i8,a,a,a,a,a)' )&
1062 0 : & 'conf_power_strain is',multibinit_dtset%conf_power_strain,', but the only allowed values',ch10,&
1063 0 : & 'are positive .',ch10,&
1064 0 : & 'Action: correct conf_power_strain in your input file.'
1065 0 : ABI_ERROR(message)
1066 : end if
1067 :
1068 73 : multibinit_dtset%conf_power_fact_disp=100
1069 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'conf_power_fact_disp',tread,'DPR')
1070 73 : if(tread==1) multibinit_dtset%conf_power_fact_disp=dprarr(1)
1071 :
1072 73 : multibinit_dtset%conf_power_fact_strain=100
1073 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'conf_power_fact_strain',tread,'DPR')
1074 73 : if(tread==1) multibinit_dtset%conf_power_fact_strain=dprarr(1)
1075 :
1076 : !D
1077 :
1078 73 : multibinit_dtset%dipdip=1
1079 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'dipdip',tread,'INT')
1080 73 : if(tread==1) multibinit_dtset%dipdip=intarr(1)
1081 : ! if(multibinit_dtset%dipdip>1.or.multibinit_dtset%dipdip<0)then
1082 : ! write(message, '(a,i8,a,a,a,a,a)' )&
1083 : !& 'dipdip is',multibinit_dtset%dipdip,', but the only allowed values',ch10,&
1084 : !& 'is 1.',ch10,&
1085 : !& 'Action: correct dipdip in your input file.'
1086 : ! ABI_ERROR(message)
1087 : ! end if
1088 :
1089 :
1090 :
1091 73 : multibinit_dtset%dipdip_prt=0
1092 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'dipdip_prt',tread,'INT')
1093 73 : if(tread==1) multibinit_dtset%dipdip_prt=intarr(1)
1094 73 : if(multibinit_dtset%dipdip_prt<0.or.multibinit_dtset%dipdip_prt>1)then
1095 : write(message, '(a,i8,a,a,a,a,a)' )&
1096 0 : & 'dipdip_prt is',multibinit_dtset%prtsrlr,', but the only allowed values',ch10,&
1097 0 : 'are 0 or 1.',ch10,&
1098 0 : & 'Action: correct dipdip_prt in your input file.'
1099 0 : ABI_ERROR(message)
1100 : end if
1101 :
1102 73 : multibinit_dtset%coeff_file_rw=0
1103 : ! call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'coeff_file_rw',tread,'INT')
1104 : ! if(tread==1) multibinit_dtset%coeff_file_rw=intarr(1)
1105 : ! if(multibinit_dtset%coeff_file_rw<0.or.multibinit_dtset%coeff_file_rw>2)then
1106 : ! write(message, '(a,i8,a,a,a,a,a)' )&
1107 : !& 'coeff_file_rw is',multibinit_dtset%prtsrlr,', but the only allowed values',ch10,&
1108 : ! 'are 0 or 1 2.',ch10,&
1109 : !& 'Action: correct coeff_file_rw in your input file.'
1110 : ! ABI_ERROR(message)
1111 : ! end if
1112 :
1113 73 : multibinit_dtset%dtion=100
1114 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'dtion',tread,'INT')
1115 73 : if(tread==1) multibinit_dtset%dtion=intarr(1)
1116 73 : if(multibinit_dtset%dtion<1)then
1117 : write(message, '(a,i8,a,a,a,a,a)' )&
1118 0 : & 'dtion is',multibinit_dtset%dtion,', but the only allowed values',ch10,&
1119 0 : & 'is superior to 1.',ch10,&
1120 0 : & 'Action: correct dtion in your input file.'
1121 0 : ABI_ERROR(message)
1122 : end if
1123 :
1124 :
1125 73 : multibinit_dtset%delta_df= 1d-02
1126 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'delta_df',tread,'DPR')
1127 73 : if(tread==1) multibinit_dtset%delta_df=dprarr(1)
1128 73 : if(multibinit_dtset%delta_df<0)then
1129 : write(message, '(a,es10.2,a,a,a,a,a)' )&
1130 0 : & 'delta_df is',multibinit_dtset%delta_df,', but the only allowed values',ch10,&
1131 0 : & 'are superior to 0 .',ch10,&
1132 0 : & 'Action: correct delta_df in your input file.'
1133 0 : ABI_ERROR(message)
1134 : end if
1135 :
1136 : !E
1137 73 : multibinit_dtset%energy_reference= zero
1138 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'energy_reference',tread,'DPR')
1139 73 : if(tread==1) multibinit_dtset%energy_reference=dprarr(1)
1140 :
1141 73 : multibinit_dtset%enunit=0
1142 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'enunit',tread,'INT')
1143 73 : if(tread==1) multibinit_dtset%enunit=intarr(1)
1144 73 : if(multibinit_dtset%enunit<0.or.multibinit_dtset%enunit>2)then
1145 : write(message, '(a,i0,a,a,a,a,a)' )&
1146 0 : & 'enunit is',multibinit_dtset%enunit,', but the only allowed values',ch10,&
1147 0 : & 'are 0, 1 or 2.',ch10,&
1148 0 : & 'Action: correct enunit in your input file.'
1149 0 : ABI_ERROR(message)
1150 : end if
1151 :
1152 : !F
1153 73 : multibinit_dtset%fit_option=0
1154 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_option',tread,'INT')
1155 73 : if(tread==1) multibinit_dtset%fit_option=intarr(1)
1156 73 : if(multibinit_dtset%fit_option<0.or.multibinit_dtset%fit_option>2)then
1157 : write(message, '(a,i8,a,a,a,a,a)' )&
1158 0 : & 'fit_option is',multibinit_dtset%fit_option,', but the only allowed values',ch10,&
1159 0 : & 'are 0, 1 or 2 for multibinit.',ch10,&
1160 0 : & 'Action: correct fit_option in your input file.'
1161 0 : ABI_ERROR(message)
1162 : end if
1163 :
1164 73 : multibinit_dtset%fit_iatom=0
1165 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_iatom',tread,'INT')
1166 73 : if(tread==1) multibinit_dtset%fit_iatom=intarr(1)
1167 73 : if(multibinit_dtset%fit_iatom<-2)then
1168 : write(message, '(a,i8,a,a,a,a,a)' )&
1169 0 : & 'fit_iatom is',multibinit_dtset%fit_iatom,', but the only allowed values',ch10,&
1170 0 : & 'are larger than -1 for multibinit.',ch10,&
1171 0 : & 'Action: correct fit_iatom in your input file.'
1172 0 : ABI_ERROR(message)
1173 : end if
1174 :
1175 73 : multibinit_dtset%fit_ncoeff=0
1176 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_ncoeff',tread,'INT')
1177 73 : if(tread==1) multibinit_dtset%fit_ncoeff=intarr(1)
1178 73 : if(multibinit_dtset%fit_ncoeff<0)then
1179 : write(message, '(a,i8,a,a,a,a,a)' )&
1180 0 : & 'fit_ncoeff is',multibinit_dtset%fit_ncoeff,', but the only allowed values',ch10,&
1181 0 : & 'are positives for multibinit.',ch10,&
1182 0 : & 'Action: correct fit_ncoeff in your input file.'
1183 0 : ABI_ERROR(message)
1184 : end if
1185 :
1186 73 : multibinit_dtset%fit_ncoeff_per_cycle=1
1187 : !call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_ncoeff_per_cycle',tread,'INT')
1188 : !if(tread==1) multibinit_dtset%fit_ncoeff_per_cycle=intarr(1)
1189 : !if(multibinit_dtset%fit_ncoeff_per_cycle<0)then
1190 : ! write(message, '(a,i8,a,a,a,a,a)' )&
1191 : ! & 'fit_ncoeff_per_cycle is',multibinit_dtset%fit_ncoeff_per_cycle,', but the only allowed values',ch10,&
1192 : ! & 'are positives for multibinit.',ch10,&
1193 : ! & 'Action: correct fit_ncoeff_per_cycle in your input file.'
1194 : ! ABI_ERROR(message)
1195 : !end if
1196 :
1197 73 : multibinit_dtset%fit_ncoeff_per_iatom=0
1198 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_ncoeff_per_iatom',tread,'INT')
1199 73 : if(tread==1) multibinit_dtset%fit_ncoeff_per_iatom=intarr(1)
1200 73 : if(multibinit_dtset%fit_ncoeff_per_iatom<0)then
1201 : write(message, '(a,i8,a,a,a,a,a)' )&
1202 0 : & 'fit_ncoeff_per_iatom is',multibinit_dtset%fit_ncoeff_per_iatom,', but the only allowed values',ch10,&
1203 0 : & 'are positives for multibinit.',ch10,&
1204 0 : & 'Action: correct fit_ncoeff_per_iatom in your input file.'
1205 0 : ABI_ERROR(message)
1206 : end if
1207 :
1208 :
1209 73 : multibinit_dtset%fit_nbancoeff=0
1210 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_nbancoeff',tread,'INT')
1211 73 : if(tread==1) multibinit_dtset%fit_nbancoeff=intarr(1)
1212 73 : if(multibinit_dtset%fit_nbancoeff<0)then
1213 : write(message, '(a,i8,a,a,a,a,a)' )&
1214 0 : & 'fit_nbancoeff is',multibinit_dtset%fit_nbancoeff,', but the only allowed values',ch10,&
1215 0 : & 'are 0 or positive values for multibinit.',ch10,&
1216 0 : & 'Action: correct fit_nbancoeff in your input file.'
1217 0 : ABI_ERROR(message)
1218 : end if
1219 :
1220 73 : multibinit_dtset%fit_nfixcoeff=0
1221 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_nfixcoeff',tread,'INT')
1222 73 : if(tread==1) multibinit_dtset%fit_nfixcoeff=intarr(1)
1223 73 : if(multibinit_dtset%fit_nfixcoeff<-2)then
1224 : write(message, '(a,i8,a,a,a,a,a)' )&
1225 0 : & 'fit_nfixcoeff is',multibinit_dtset%fit_nfixcoeff,', but the only allowed values',ch10,&
1226 0 : & 'are -1 or positives for multibinit.',ch10,&
1227 0 : & 'Action: correct fit_nfixcoeff in your input file.'
1228 0 : ABI_ERROR(message)
1229 : end if
1230 :
1231 73 : multibinit_dtset%fit_nimposecoeff=0
1232 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_nimposecoeff',tread,'INT')
1233 73 : if(tread==1) multibinit_dtset%fit_nimposecoeff=intarr(1)
1234 73 : if(multibinit_dtset%fit_nimposecoeff<-2)then
1235 : write(message, '(a,i8,a,a,a,a,a)' )&
1236 0 : & 'fit_nimposecoeff is',multibinit_dtset%fit_nimposecoeff,', but the only allowed values',ch10,&
1237 0 : & 'are -1 or positives for multibinit.',ch10,&
1238 0 : & 'Action: correct fit_nimposecoeff in your input file.'
1239 0 : ABI_ERROR(message)
1240 : end if
1241 :
1242 292 : multibinit_dtset%fit_EFS=(/0,1,1/)
1243 292 : multibinit_dtset%fit_on = (/ .TRUE.,.TRUE.,.FALSE. /)
1244 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'fit_EFS',tread,'INT')
1245 79 : if(tread==1) multibinit_dtset%fit_EFS(1:3)=intarr(1:3)
1246 73 : if(multibinit_dtset%fit_EFS(1) == 1)multibinit_dtset%fit_on(3) = .TRUE.
1247 73 : if(multibinit_dtset%fit_EFS(2) == 0)multibinit_dtset%fit_on(1) = .FALSE.
1248 73 : if(multibinit_dtset%fit_EFS(3) == 0)multibinit_dtset%fit_on(2) = .FALSE.
1249 584 : if(any(multibinit_dtset%fit_EFS<0) .or. any(multibinit_dtset%fit_EFS>1))then
1250 : write(message, '(a,i8,a,a,a)' )&
1251 0 : & 'fit_EFS is',multibinit_dtset%fit_EFS,', but the only allowed values are 0 and 1',ch10,&
1252 0 : & 'Action: correct fit_EFS in your input file.'
1253 0 : ABI_ERROR(message)
1254 : end if
1255 :
1256 292 : multibinit_dtset%fit_factors=(/1,1,1/)
1257 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'fit_factors',tread,'DPR')
1258 76 : if(tread==1) multibinit_dtset%fit_factors(1:3)=dprarr(1:3)
1259 292 : if(any(multibinit_dtset%fit_factors<0))then
1260 : write(message, '(a,i8,a,a,a)' )&
1261 0 : & 'fit_factors is',multibinit_dtset%fit_factors,', but the only allowed values are positive',ch10,&
1262 0 : & 'Action: correct fit_factors in your input file.'
1263 0 : ABI_ERROR(message)
1264 : end if
1265 292 : multibinit_dtset%sel_EFS=(/0,1,1/)
1266 292 : multibinit_dtset%sel_on = (/ .TRUE.,.TRUE.,.FALSE. /)
1267 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'sel_EFS',tread,'INT')
1268 79 : if(tread==1) multibinit_dtset%sel_EFS(1:3)=intarr(1:3)
1269 73 : if(multibinit_dtset%sel_EFS(1) == 1)multibinit_dtset%sel_on(3) = .TRUE.
1270 73 : if(multibinit_dtset%sel_EFS(2) == 0)multibinit_dtset%sel_on(1) = .FALSE.
1271 73 : if(multibinit_dtset%sel_EFS(3) == 0)multibinit_dtset%sel_on(2) = .FALSE.
1272 584 : if(any(multibinit_dtset%sel_EFS<0) .or. any(multibinit_dtset%sel_EFS>1))then
1273 : write(message, '(a,i8,a,a,a)' )&
1274 0 : & 'sel_EFS is',multibinit_dtset%sel_EFS,', but the only allowed values are 0 and 1',ch10,&
1275 0 : & 'Action: correct sel_EFS in your input file.'
1276 0 : ABI_ERROR(message)
1277 : end if
1278 :
1279 73 : multibinit_dtset%ts_option=0
1280 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'ts_option',tread,'INT')
1281 73 : if(tread==1) multibinit_dtset%ts_option=intarr(1)
1282 73 : if(multibinit_dtset%ts_option<0.or.multibinit_dtset%ts_option>1)then
1283 : write(message, '(a,i8,a,a,a,a,a)' )&
1284 0 : & 'ts_option is',multibinit_dtset%ts_option,', but the only allowed values',ch10,&
1285 0 : & 'are zero and one for multibinit.',ch10,&
1286 0 : & 'Action: correct ts_option in your input file.'
1287 0 : ABI_ERROR(message)
1288 : end if
1289 :
1290 73 : multibinit_dtset%ifcana=0
1291 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'ifcana',tread,'INT')
1292 73 : if(tread==1) multibinit_dtset%ifcana=intarr(1)
1293 73 : if(multibinit_dtset%ifcana<0.or.multibinit_dtset%ifcana>1)then
1294 : write(message, '(a,i0,a,a,a,a,a)' )&
1295 0 : & 'ifcana is',multibinit_dtset%ifcana,', but the only allowed values',ch10,&
1296 0 : & 'are 0 or 1.',ch10,&
1297 0 : & 'Action: correct ifcana in your input file.'
1298 0 : ABI_ERROR(message)
1299 : end if
1300 :
1301 73 : multibinit_dtset%ifcflag=1
1302 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'ifcflag',tread,'INT')
1303 73 : if(tread==1) multibinit_dtset%ifcflag=intarr(1)
1304 73 : if(multibinit_dtset%ifcflag<0.or.multibinit_dtset%ifcflag>1)then
1305 : write(message, '(a,i0,a,a,a,a,a)' )&
1306 0 : & 'ifcflag is',multibinit_dtset%ifcflag,', but the only allowed values',ch10,&
1307 0 : & 'are 0 or 1.',ch10,&
1308 0 : & 'Action: correct ifcflag in your input file.'
1309 0 : ABI_ERROR(message)
1310 : end if
1311 :
1312 73 : multibinit_dtset%prtsrlr=0
1313 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'prtsrlr',tread,'INT')
1314 73 : if(tread==1) multibinit_dtset%prtsrlr=intarr(1)
1315 73 : if(multibinit_dtset%prtsrlr<0.or.multibinit_dtset%prtsrlr>1)then
1316 : write(message, '(a,i8,a,a,a,a,a)' )&
1317 0 : & 'prtsrlr is',multibinit_dtset%prtsrlr,', but the only allowed values',ch10,&
1318 0 : & 'are 0 or 1.',ch10,&
1319 0 : & 'Action: correct prtsrlr in your input file.'
1320 0 : ABI_ERROR(message)
1321 : end if
1322 :
1323 73 : multibinit_dtset%ifcout=2000000 ! or -1 -> max number of ifc
1324 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'ifcout',tread,'INT')
1325 73 : if(tread==1) multibinit_dtset%ifcout=intarr(1)
1326 73 : if(multibinit_dtset%ifcout<-1)then
1327 : write(message, '(a,i0,a,a,a)' )&
1328 0 : & 'ifcout is',multibinit_dtset%ifcout,', which is lower than -1 (default = all ifc) .',ch10,&
1329 0 : & 'Action: correct ifcout in your input file.'
1330 0 : ABI_ERROR(message)
1331 : end if
1332 73 : multibinit_dtset%hmctt=0 ! or -1 -> max number of ifc
1333 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'hmctt',tread,'INT')
1334 73 : if(tread==1) multibinit_dtset%hmctt=intarr(1)
1335 73 : if(multibinit_dtset%hmctt< 0)then
1336 : write(message, '(a,i0,a,a,a)' )&
1337 0 : & 'hmctt is',multibinit_dtset%hmctt,', but only positive values are allowed .',ch10,&
1338 0 : & 'Action: correct hmctt in your input file.'
1339 0 : ABI_ERROR(message)
1340 : end if
1341 :
1342 73 : multibinit_dtset%hmcsst=0 ! or -1 -> max number of ifc
1343 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'hmcsst',tread,'INT')
1344 73 : if(tread==1) multibinit_dtset%hmcsst=intarr(1)
1345 73 : if(multibinit_dtset%hmcsst<0)then
1346 : write(message, '(a,i0,a,a,a)' )&
1347 0 : & 'hmcsst is',multibinit_dtset%hmcsst,', but only positive values are allowed .',ch10,&
1348 0 : & 'Action: correct hmcsst in your input file.'
1349 0 : ABI_ERROR(message)
1350 : end if
1351 :
1352 73 : multibinit_dtset%nctime=1
1353 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'nctime',tread,'INT')
1354 73 : if(tread==1) multibinit_dtset%nctime=intarr(1)
1355 73 : if(multibinit_dtset%nctime<=0)then
1356 : write(message, '(a,i0,a,a,a)' )&
1357 0 : & 'nctime is',multibinit_dtset%ntime,', which is not positive .',ch10,&
1358 0 : & 'Action: correct nctime in your input file.'
1359 0 : ABI_ERROR(message)
1360 : end if
1361 :
1362 :
1363 :
1364 73 : multibinit_dtset%ntime=200
1365 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'ntime',tread,'INT')
1366 73 : if(tread==1) multibinit_dtset%ntime=intarr(1)
1367 73 : if(multibinit_dtset%ntime<0)then
1368 : write(message, '(a,i0,a,a,a)' )&
1369 0 : & 'ntime is',multibinit_dtset%ntime,', which is lower than 0 .',ch10,&
1370 0 : & 'Action: correct ntime in your input file.'
1371 0 : ABI_ERROR(message)
1372 : end if
1373 :
1374 73 : multibinit_dtset%dynamics=0
1375 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'dynamics',tread,'INT')
1376 73 : if(tread==1) multibinit_dtset%dynamics=intarr(1)
1377 : ! >100: the builtin multibinit lattice movers
1378 : if(multibinit_dtset%dynamics/=0.and.multibinit_dtset%dynamics/=6.and.&
1379 : & multibinit_dtset%dynamics/=12.and.multibinit_dtset%dynamics/=13.and.&
1380 : & multibinit_dtset%dynamics/=27.and.&
1381 : & multibinit_dtset%dynamics/=9.and.&
1382 : & multibinit_dtset%dynamics/=7.and.&
1383 : & multibinit_dtset%dynamics/=1.and.&
1384 : & multibinit_dtset%dynamics/=2.and.&
1385 : & multibinit_dtset%dynamics/=22.and.&
1386 : & multibinit_dtset%dynamics/=24.and.multibinit_dtset%dynamics/=25 .and. &
1387 : & multibinit_dtset%dynamics/=101.and.multibinit_dtset%dynamics/=102 .and. &
1388 73 : & multibinit_dtset%dynamics/=103.and.multibinit_dtset%dynamics/=120 &
1389 : ) then
1390 : write(message, '(a,i8,a,a,a,a,a)' )&
1391 0 : & 'dynamics is ',multibinit_dtset%dynamics,', but the only allowed values',ch10,&
1392 0 : & 'are 1,2,6,7,9,12,13, 22,24,25,101,102, 103 or 120 (see ionmov in abinit documentation).',ch10,&
1393 0 : & 'Action: correct dynamics in your input file.'
1394 0 : ABI_ERROR(message)
1395 : end if
1396 :
1397 73 : if(multibinit_dtset%dynamics==120) then
1398 : write(message, '(a,i8,a)' )&
1399 0 : & 'dynamics is ',multibinit_dtset%dynamics,'The atoms will not move. For test only!'
1400 0 : ABI_WARNING(message)
1401 : end if
1402 :
1403 73 : multibinit_dtset%dyn_chksym=0
1404 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'dyn_chksym',tread,'INT')
1405 73 : if(tread==1) multibinit_dtset%dyn_chksym=intarr(1)
1406 73 : if(multibinit_dtset%dyn_chksym<0 .or. multibinit_dtset%dyn_chksym>1)then
1407 : write(message, '(a,i0,a,a,a)' )&
1408 0 : & 'dyn_chksym is',multibinit_dtset%dyn_chksym,', but the only allowed values are 0 and 1.',ch10,&
1409 0 : & 'Action: correct dyn_chksym in your input file.'
1410 0 : ABI_ERROR(message)
1411 : end if
1412 :
1413 73 : multibinit_dtset%dyn_tolsym=1d-10
1414 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'dyn_tolsym',tread,'DPR')
1415 73 : if(tread==1) multibinit_dtset%dyn_tolsym=dprarr(1)
1416 73 : if(multibinit_dtset%dyn_tolsym<0)then
1417 : write(message, '(a,i0,a,a,a)' )&
1418 0 : & 'dyn_tolsym is',multibinit_dtset%dyn_tolsym,', but the only allowed values are positive.',ch10,&
1419 0 : & 'Action: correct dyn_tolsym in your input file.'
1420 0 : ABI_ERROR(message)
1421 : end if
1422 : !L
1423 : ! multibinit_dtset%latt_compressibility=0.0
1424 : ! call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_compressibility',tread,'DPR')
1425 : ! if(tread==1) multibinit_dtset%latt_compressibility=dprarr(1)
1426 :
1427 73 : multibinit_dtset%latt_friction=1e-4
1428 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_friction',tread,'DPR')
1429 73 : if(tread==1) multibinit_dtset%latt_friction=dprarr(1)
1430 73 : multibinit_dtset%latt_lwf_anharmonic=0
1431 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_lwf_anharmonic',tread,'INT')
1432 73 : if(tread==1) multibinit_dtset%latt_lwf_anharmonic=intarr(1)
1433 73 : if( .not. (0 <= multibinit_dtset%latt_lwf_anharmonic .and. multibinit_dtset%latt_lwf_anharmonic< 2) ) then
1434 : write(message, '(a,i8,a,a,a,a,a)' )&
1435 0 : & 'latt_lwf_anharmonic is ',multibinit_dtset%latt_lwf_anharmonic,', but the only allowed values',ch10,&
1436 0 : & 'are 0 and 1',ch10,&
1437 0 : & 'Action: correct latt_lwf_anharmonic in your input file.'
1438 0 : ABI_ERROR(message)
1439 : end if
1440 :
1441 73 : multibinit_dtset%latt_taut=0.0
1442 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_taut',tread,'TIM')
1443 73 : if(tread==1) multibinit_dtset%latt_taut=dprarr(1)
1444 73 : if(multibinit_dtset%latt_taut<0)then
1445 : write(message, '(a,f10.1,a,a,a,a,a)' )&
1446 0 : & 'latt_taut is ',multibinit_dtset%latt_taut,'. The only allowed values',ch10,&
1447 0 : & 'are non-negative values.',ch10,&
1448 0 : & 'Action: correct latt_taut in your input file.'
1449 0 : ABI_ERROR(message)
1450 : end if
1451 :
1452 :
1453 :
1454 :
1455 73 : multibinit_dtset%latt_temperature_start=0.0
1456 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_temperature_start',tread,'DPR')
1457 73 : if(tread==1) multibinit_dtset%latt_temperature_start=dprarr(1)
1458 73 : if(multibinit_dtset%latt_temperature_start<0.0)then
1459 : write(message, '(a,f10.1,a,a,a,a,a)' )&
1460 0 : & 'latt_temperature_start is ',multibinit_dtset%latt_temperature_start,'. The only allowed values',ch10,&
1461 0 : & 'are positives values.',ch10,&
1462 0 : & 'Action: correct latt_semperature_start in your input file.'
1463 0 : ABI_ERROR(message)
1464 : end if
1465 :
1466 73 : multibinit_dtset%latt_temperature_end=0.0
1467 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_temperature_end',tread,'DPR')
1468 73 : if(tread==1) multibinit_dtset%latt_temperature_end=dprarr(1)
1469 73 : if(multibinit_dtset%latt_temperature_end<0)then
1470 : write(message, '(a,f10.1,a,a,a,a,a)' )&
1471 0 : & 'latt_temperature_end is ',multibinit_dtset%latt_temperature_end,'. The only allowed values',ch10,&
1472 0 : & 'are positives values.',ch10,&
1473 0 : & 'Action: correct latt_semperature_end in your input file.'
1474 0 : ABI_ERROR(message)
1475 : end if
1476 :
1477 73 : multibinit_dtset%latt_temperature_nstep=1
1478 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_temperature_nstep',tread,'INT')
1479 73 : if(tread==1) multibinit_dtset%latt_temperature_nstep=intarr(1)
1480 73 : if(multibinit_dtset%latt_temperature_nstep<=0)then
1481 : write(message, '(a,i0,a,a,a,a)' )&
1482 0 : & 'latt_temperature_nstep is',multibinit_dtset%latt_temperature_nstep,', while it should be larger than 0',ch10,&
1483 0 : & 'Action: correct latt_temperature_nstep in your input file.'
1484 0 : ABI_ERROR(message)
1485 : end if
1486 :
1487 73 : multibinit_dtset%latt_var_temperature=0
1488 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_var_temperature',tread,'INT')
1489 73 : if(tread==1) multibinit_dtset%latt_var_temperature=intarr(1)
1490 73 : if(multibinit_dtset%latt_var_temperature/=0.and.multibinit_dtset%latt_var_temperature/=1)then
1491 : write(message, '(a,i0,a,a,a,a,a)' )&
1492 0 : & 'latt_var_temperature is',multibinit_dtset%latt_var_temperature,'. The only allowed values',ch10,&
1493 0 : & 'are 0, or 1.',ch10,&
1494 0 : & 'Action: correct latt_var_temperature in your input file.'
1495 0 : ABI_ERROR(message)
1496 : end if
1497 :
1498 : ! multibinit_dtset%latt_taup=1000
1499 : ! call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_taup',tread,'TIM')
1500 : ! if(tread==1) multibinit_dtset%latt_taup=dprarr(1)
1501 :
1502 73 : multibinit_dtset%lwf_dt= 1d-16
1503 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'lwf_dt',tread,'TIM')
1504 73 : if(tread==1) multibinit_dtset%lwf_dt=dprarr(1)
1505 73 : if(multibinit_dtset%lwf_dt<0)then
1506 : write(message, '(a,es10.2,a,a,a,a,a)' )&
1507 0 : & 'lwf_dt is',multibinit_dtset%lwf_dt,', but the only allowed values',ch10,&
1508 0 : & 'are superior to 0 .',ch10,&
1509 0 : & 'Action: correct lwf_dt in your input file.'
1510 0 : ABI_ERROR(message)
1511 : end if
1512 :
1513 73 : multibinit_dtset%lwf_init_state=0
1514 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'lwf_init_state',tread,'INT')
1515 73 : if(tread==1) multibinit_dtset%lwf_init_state=intarr(1)
1516 73 : if( .not. (multibinit_dtset%lwf_init_state <= 4) ) then
1517 : write(message, '(a,i8,a,a,a,a,a)' )&
1518 0 : & 'lwf_init_state is ',multibinit_dtset%lwf_init_state,', but the only allowed values',ch10,&
1519 0 : & 'are 0, 1,2,3,4 and negative values.',ch10,&
1520 0 : & 'Action: correct lwf_init_state in your input file.'
1521 0 : ABI_ERROR(message)
1522 : end if
1523 :
1524 :
1525 73 : multibinit_dtset%lwf_init_hist_fname=""
1526 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'lwf_init_hist_fname',tread,'KEY',&
1527 73 : & key_value=multibinit_dtset%lwf_init_hist_fname)
1528 73 : if(.not. tread==1) multibinit_dtset%lwf_init_hist_fname="lwf_init_hist.nc"
1529 :
1530 :
1531 73 : multibinit_dtset%spin_init_hist_fname=""
1532 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_init_hist_fname',tread,'KEY',&
1533 73 : & key_value=multibinit_dtset%spin_init_hist_fname)
1534 73 : if(.not. tread==1) multibinit_dtset%spin_init_hist_fname="spin_init_hist.nc"
1535 :
1536 73 : multibinit_dtset%latt_init_hist_fname=""
1537 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_init_hist_fname',tread,'KEY',&
1538 73 : & key_value=multibinit_dtset%latt_init_hist_fname)
1539 73 : if(.not. tread==1) multibinit_dtset%latt_init_hist_fname="latt_init_hist.nc"
1540 :
1541 :
1542 :
1543 :
1544 73 : multibinit_dtset%slc_pot_fname=""
1545 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'slc_pot_fname',tread,'KEY',&
1546 73 : & key_value=multibinit_dtset%slc_pot_fname)
1547 :
1548 :
1549 73 : multibinit_dtset%lwf_constraint=0
1550 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'lwf_constraint',tread,'INT')
1551 73 : if(tread==1) multibinit_dtset%lwf_constraint=intarr(1)
1552 73 : if( .not. (multibinit_dtset%lwf_constraint == 0 .or. multibinit_dtset%lwf_constraint == 1) ) then
1553 : write(message, '(a,i8,a,a,a,a,a)' )&
1554 0 : & 'lwf_constraint is ',multibinit_dtset%lwf_constraint,', but the only allowed values',ch10,&
1555 0 : & 'are 0 or 1',ch10,&
1556 0 : & 'Action: correct lwf_constraint in your input file.'
1557 0 : ABI_ERROR(message)
1558 : end if
1559 :
1560 :
1561 73 : multibinit_dtset%lwf_dynamics=0
1562 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'lwf_dynamics',tread,'INT')
1563 73 : if(tread==1) multibinit_dtset%lwf_dynamics=intarr(1)
1564 73 : if( .not. (multibinit_dtset%lwf_dynamics <= 3) ) then
1565 : write(message, '(a,i8,a,a,a,a,a)' )&
1566 0 : & 'lwf_dynamics is ',multibinit_dtset%lwf_dynamics,', but the only allowed values',ch10,&
1567 0 : & 'are 0, 1, 2, 3 and negative values.',ch10,&
1568 0 : & 'Action: correct lwf_dynamics in your input file.'
1569 0 : ABI_ERROR(message)
1570 : end if
1571 :
1572 73 : multibinit_dtset%lwf_mc_avg_amp=0.05
1573 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'lwf_mc_avg_amp',tread,'DPR')
1574 73 : if(tread==1) multibinit_dtset%lwf_mc_avg_amp=dprarr(1)
1575 73 : if(multibinit_dtset%lwf_mc_avg_amp<0)then
1576 : write(message, '(a,f10.1,a,a,a,a,a)' )&
1577 0 : & 'lwf_mc_avg_amp is ',multibinit_dtset%lwf_mc_avg_amp,'. The only allowed values',ch10,&
1578 0 : & 'are non-negative values.',ch10,&
1579 0 : & 'Action: correct lwf_mc_avg_amp in your input file.'
1580 0 : ABI_ERROR(message)
1581 : end if
1582 :
1583 :
1584 :
1585 :
1586 73 : multibinit_dtset%lwf_nctime=1
1587 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'lwf_nctime',tread,'INT')
1588 73 : if(tread==1) multibinit_dtset%lwf_nctime=intarr(1)
1589 73 : if( .not. (multibinit_dtset%lwf_nctime >= 0) ) then
1590 : write(message, '(a,i8,a,a,a,a,a)' )&
1591 0 : & 'lwf_nctime is ',multibinit_dtset%lwf_nctime,', but the only non-negative values allowed.'
1592 0 : ABI_ERROR(message)
1593 : end if
1594 :
1595 :
1596 73 : multibinit_dtset%lwf_ntime=0
1597 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'lwf_ntime',tread,'INT')
1598 73 : if(tread==1) multibinit_dtset%lwf_ntime=intarr(1)
1599 73 : if( .not. (multibinit_dtset%lwf_ntime >= 0) ) then
1600 : write(message, '(a,i8,a,a,a,a,a)' )&
1601 0 : & 'lwf_ntime is ',multibinit_dtset%lwf_ntime,', but the only non-negative values allowed.'
1602 0 : ABI_ERROR(message)
1603 : end if
1604 :
1605 :
1606 : !multibinit_dtset%lwf_self_bound_order=0
1607 : !call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'lwf_self_bound_order',tread,'INT')
1608 : !if(tread==1) multibinit_dtset%lwf_self_bound_order=intarr(1)
1609 :
1610 :
1611 : !multibinit_dtset%lwf_self_bound_coeff=0.0
1612 : !call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'lwf_self_bound_coeff',tread,'DPR')
1613 : !if(tread==1) multibinit_dtset%lwf_self_bound_coeff=dprarr(1)
1614 :
1615 73 : multibinit_dtset%lwf_taut=0.0
1616 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'lwf_taut',tread,'TIM')
1617 73 : if(tread==1) multibinit_dtset%lwf_taut=dprarr(1)
1618 73 : if(multibinit_dtset%lwf_taut<0)then
1619 : write(message, '(a,f10.1,a,a,a,a,a)' )&
1620 0 : & 'lwf_taut is ',multibinit_dtset%lwf_taut,'. The only allowed values',ch10,&
1621 0 : & 'are non-negative values.',ch10,&
1622 0 : & 'Action: correct lwf_taut in your input file.'
1623 0 : ABI_ERROR(message)
1624 : end if
1625 :
1626 :
1627 73 : multibinit_dtset%lwf_temperature=0.0
1628 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'lwf_temperature',tread,'DPR')
1629 73 : if(tread==1) multibinit_dtset%lwf_temperature=dprarr(1)
1630 73 : if(multibinit_dtset%lwf_temperature<0)then
1631 : write(message, '(a,f10.1,a,a,a,a,a)' )&
1632 0 : & 'lwf_temperature is ',multibinit_dtset%lwf_temperature,'. The only allowed values',ch10,&
1633 0 : & 'are non-negative values.',ch10,&
1634 0 : & 'Action: correct lwf_temperature in your input file.'
1635 0 : ABI_ERROR(message)
1636 : end if
1637 :
1638 73 : multibinit_dtset%lwf_temperature_start=0.0
1639 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'lwf_temperature_start',tread,'DPR')
1640 73 : if(tread==1) multibinit_dtset%lwf_temperature_start=dprarr(1)
1641 73 : if(multibinit_dtset%lwf_temperature_start<0.0)then
1642 : write(message, '(a,f10.1,a,a,a,a,a)' )&
1643 0 : & 'lwf_temperature_start is ',multibinit_dtset%lwf_temperature_start,'. The only allowed values',ch10,&
1644 0 : & 'are positives values.',ch10,&
1645 0 : & 'Action: correct lwf_semperature_start in your input file.'
1646 0 : ABI_ERROR(message)
1647 : end if
1648 :
1649 73 : multibinit_dtset%lwf_temperature_end=0.0
1650 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'lwf_temperature_end',tread,'DPR')
1651 73 : if(tread==1) multibinit_dtset%lwf_temperature_end=dprarr(1)
1652 73 : if(multibinit_dtset%lwf_temperature_end<0)then
1653 : write(message, '(a,f10.1,a,a,a,a,a)' )&
1654 0 : & 'lwf_temperature_end is ',multibinit_dtset%lwf_temperature_end,'. The only allowed values',ch10,&
1655 0 : & 'are positives values.',ch10,&
1656 0 : & 'Action: correct lwf_semperature_end in your input file.'
1657 0 : ABI_ERROR(message)
1658 : end if
1659 :
1660 73 : multibinit_dtset%lwf_temperature_nstep=1
1661 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'lwf_temperature_nstep',tread,'INT')
1662 73 : if(tread==1) multibinit_dtset%lwf_temperature_nstep=intarr(1)
1663 73 : if(multibinit_dtset%lwf_temperature_nstep<=0)then
1664 : write(message, '(a,i0,a,a,a,a)' )&
1665 0 : & 'lwf_temperature_nstep is',multibinit_dtset%lwf_temperature_nstep,', while it should be larger than 0',ch10,&
1666 0 : & 'Action: correct lwf_temperature_nstep in your input file.'
1667 0 : ABI_ERROR(message)
1668 : end if
1669 :
1670 73 : multibinit_dtset%lwf_var_temperature=0
1671 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'lwf_var_temperature',tread,'INT')
1672 73 : if(tread==1) multibinit_dtset%lwf_var_temperature=intarr(1)
1673 73 : if(multibinit_dtset%lwf_var_temperature/=0.and.multibinit_dtset%lwf_var_temperature/=1)then
1674 : write(message, '(a,i0,a,a,a,a,a)' )&
1675 0 : & 'lwf_var_temperature is',multibinit_dtset%lwf_var_temperature,'. The only allowed values',ch10,&
1676 0 : & 'are 0, or 1.',ch10,&
1677 0 : & 'Action: correct lwf_var_temperature in your input file.'
1678 0 : ABI_ERROR(message)
1679 : end if
1680 :
1681 :
1682 :
1683 73 : multibinit_dtset%spin_init_hist_fname=""
1684 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_init_hist_fname',tread,'KEY',&
1685 73 : & key_value=multibinit_dtset%spin_init_hist_fname)
1686 73 : if(.not. tread==1) multibinit_dtset%spin_init_hist_fname="spin_init_hist.nc"
1687 :
1688 73 : multibinit_dtset%latt_init_hist_fname=""
1689 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_init_hist_fname',tread,'KEY',&
1690 73 : & key_value=multibinit_dtset%latt_init_hist_fname)
1691 73 : if(.not. tread==1) multibinit_dtset%latt_init_hist_fname="latt_init_hist.nc"
1692 :
1693 :
1694 :
1695 73 : multibinit_dtset%spin_pot_fname=""
1696 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_pot_fname',tread,'KEY',&
1697 73 : & key_value=multibinit_dtset%spin_pot_fname)
1698 73 : if(.not. tread==1) multibinit_dtset%spin_pot_fname=""
1699 :
1700 73 : multibinit_dtset%latt_pot_fname=""
1701 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_pot_fname',tread,'KEY',&
1702 73 : & key_value=multibinit_dtset%latt_pot_fname)
1703 73 : if(.not. tread==1) multibinit_dtset%latt_pot_fname=""
1704 :
1705 73 : multibinit_dtset%latt_harm_pot_fname=""
1706 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_harm_pot_fname',tread,'KEY',&
1707 73 : & key_value=multibinit_dtset%latt_harm_pot_fname)
1708 73 : if(.not. tread==1) multibinit_dtset%latt_harm_pot_fname=""
1709 :
1710 73 : multibinit_dtset%latt_anharm_pot_fname=""
1711 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_anharm_pot_fname',tread,'KEY',&
1712 73 : & key_value=multibinit_dtset%latt_anharm_pot_fname)
1713 73 : if(.not. tread==1) multibinit_dtset%latt_anharm_pot_fname=""
1714 :
1715 73 : multibinit_dtset%latt_training_set_fname=""
1716 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_training_set_fname',tread,'KEY',&
1717 73 : & key_value=multibinit_dtset%latt_training_set_fname)
1718 73 : if(.not. tread==1) multibinit_dtset%latt_training_set_fname=""
1719 :
1720 73 : multibinit_dtset%latt_test_set_fname=""
1721 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_test_set_fname',tread,'KEY',&
1722 73 : & key_value=multibinit_dtset%latt_test_set_fname)
1723 73 : if(.not. tread==1) multibinit_dtset%latt_test_set_fname=""
1724 :
1725 :
1726 949 : multibinit_dtset%latt_ddb_fnames=""
1727 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_ddb_fnames',tread,'KEY',&
1728 73 : & key_value=lattddb_string)
1729 73 : if(.not. tread==1) then
1730 73 : lattddb_string=""
1731 949 : multibinit_dtset%latt_ddb_fnames(:)=""
1732 : else
1733 0 : sidx(1) = 1; sidx(13) = len(lattddb_string)
1734 0 : cnt = 1
1735 0 : do ii=1,len(lattddb_string)
1736 0 : if (lattddb_string(ii:ii) == ",") then
1737 0 : lattddb_string(ii:ii) = " "
1738 0 : cnt = cnt + 1
1739 0 : sidx(cnt) = ii
1740 0 : ABI_CHECK(cnt <= 12, "Too many commas in latt_ddb_fnames!")
1741 : end if
1742 : end do
1743 :
1744 0 : do ii=1,cnt
1745 0 : i1 = sidx(ii)
1746 0 : i2 = sidx(ii + 1)
1747 0 : cnt = len(adjustl(trim(lattddb_string(i1:i2))))
1748 0 : ABI_CHECK(cnt <= fnlen, "latt_ddb_fnames path too small, increase fnlen")
1749 0 : multibinit_dtset%latt_ddb_fnames(ii) = adjustl(trim(lattddb_string(i1:i2)))
1750 : ! if (len_trim(pp_dirpath) > 0) then
1751 : ! if (len_trim(pp_dirpath) + len_trim(pseudo_paths(ii)) > fnlen) then
1752 : ! MSG_ERROR(sjoin("String of len fnlen:", itoa(fnlen), " too small to contain full pseudo path"))
1753 : ! end if
1754 : ! pseudo_paths(ii) = strcat(pp_dirpath, pseudo_paths(ii))
1755 : !end if
1756 : end do
1757 : end if
1758 :
1759 73 : multibinit_dtset%lwf_pot_fname=""
1760 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'lwf_pot_fname',tread,'KEY',&
1761 73 : & key_value=multibinit_dtset%lwf_pot_fname)
1762 73 : if(.not. tread==1) multibinit_dtset%lwf_pot_fname=""
1763 :
1764 73 : multibinit_dtset%slc_pot_fname=""
1765 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'slc_pot_fname',tread,'KEY',&
1766 73 : & key_value=multibinit_dtset%slc_pot_fname)
1767 73 : if(.not. tread==1) multibinit_dtset%slc_pot_fname=""
1768 :
1769 73 : multibinit_dtset%outdata_prefix=""
1770 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'outdata_prefix',tread,'KEY',&
1771 73 : & key_value=multibinit_dtset%outdata_prefix)
1772 73 : if(.not. tread==1) multibinit_dtset%outdata_prefix="multibinit"
1773 :
1774 : !N
1775 73 : multibinit_dtset%natifc=natom
1776 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'natifc',tread,'INT')
1777 73 : if(tread==1) multibinit_dtset%natifc=intarr(1)
1778 73 : if(multibinit_dtset%natifc<0)then
1779 : write(message, '(a,i0,a,a,a)' )&
1780 0 : & 'natifc is',multibinit_dtset%natifc,', which is lower than 0 .',ch10,&
1781 0 : & 'Action: correct natifc in your input file.'
1782 0 : ABI_ERROR(message)
1783 : end if
1784 :
1785 73 : if(multibinit_dtset%natifc>natom)then
1786 : write(message, '(a,i0,a,a,a,i0,a,a,a)' )&
1787 0 : & 'The number of atom ifc in the input files',multibinit_dtset%natifc,',',ch10,&
1788 0 : & 'is larger than the number of atoms',natom,'.',ch10,&
1789 0 : & 'Action: change natifc in the input file.'
1790 0 : ABI_ERROR(message)
1791 : end if
1792 :
1793 73 : multibinit_dtset%ncoeff=0
1794 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'ncoeff',tread,'INT')
1795 73 : if(tread==1) multibinit_dtset%ncoeff=intarr(1)
1796 73 : if(multibinit_dtset%ncoeff<0)then
1797 : write(message, '(a,i0,a,a,a)' )&
1798 0 : & 'ncoeff is',multibinit_dtset%ncoeff,', which is lower than 0 .',ch10,&
1799 0 : & 'Action: correct ncoeff in your input file.'
1800 0 : ABI_ERROR(message)
1801 : end if
1802 :
1803 292 : multibinit_dtset%ng2qpt(:)=0
1804 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'ng2qpt',tread,'INT')
1805 73 : if(tread==1) multibinit_dtset%ng2qpt(:)=intarr(1:3)
1806 292 : do ii=1,3
1807 292 : if(multibinit_dtset%ng2qpt(ii)<0)then
1808 : write(message, '(a,i0,a,i0,a,a,a,i0,a)' )&
1809 0 : & 'ng2qpt(',ii,') is',multibinit_dtset%ng2qpt(ii),', which is lower than 0 .',ch10,&
1810 0 : & 'Action: correct ng2qpt(',ii,') in your input file.'
1811 0 : ABI_ERROR(message)
1812 : end if
1813 : end do
1814 :
1815 292 : multibinit_dtset%ncell(:)= 1
1816 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'ncell',tread,'INT')
1817 241 : if(tread==1) multibinit_dtset%ncell(1:3)=intarr(1:3)
1818 292 : do ii=1,3
1819 292 : if(multibinit_dtset%ncell(ii)<0.or.multibinit_dtset%ncell(ii)>150)then
1820 : write(message, '(a,i0,a,i0,3a,i0,a)' )&
1821 0 : & 'ncell(',ii,') is ',multibinit_dtset%ncell(ii),', which is lower than 0 of superior than 150.',&
1822 0 : & ch10,'Action: correct ncell(',ii,') in your input file.'
1823 0 : ABI_ERROR(message)
1824 : end if
1825 : end do
1826 :
1827 : ! Set to diagonal ncell. Then if it is specified, overwrite.
1828 : multibinit_dtset%ncellmat(:,:)= reshape([multibinit_dtset%ncell(1), 0, 0,&
1829 : & 0, multibinit_dtset%ncell(2), 0, &
1830 1606 : & 0, 0, multibinit_dtset%ncell(3) ], [3,3])
1831 73 : call intagm(dprarr,intarr,jdtset,marr,9,string(1:lenstr),'ncellmat',tread,'INT')
1832 73 : if(tread==1) then
1833 13 : multibinit_dtset%ncellmat=transpose(reshape(intarr(1:9), [3,3]))
1834 : end if
1835 :
1836 730 : multibinit_dtset%ngqpt(:)= 1
1837 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'ngqpt',tread,'INT')
1838 145 : if(tread==1) multibinit_dtset%ngqpt(1:3)=intarr(1:3)
1839 292 : do ii=1,3
1840 292 : if(multibinit_dtset%ngqpt(ii)<0)then
1841 : write(message, '(a,i0,a,i0,a,a,a,i0,a)' )&
1842 0 : & 'ngqpt(',ii,') is',multibinit_dtset%ngqpt(ii),', which is lower than 0 .',ch10,&
1843 0 : & 'Action: correct ngqpt(',ii,') in your input file.'
1844 0 : ABI_ERROR(message)
1845 : end if
1846 : end do
1847 :
1848 73 : multibinit_dtset%nph1l=1
1849 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'nph1l',tread,'INT')
1850 73 : if(tread==1) multibinit_dtset%nph1l=intarr(1)
1851 73 : if(multibinit_dtset%nph1l<0) then
1852 : write(message, '(a,i0,a,a,a)' )&
1853 0 : & 'nph1l is',multibinit_dtset%nph1l,', which is lower than 0 .',ch10,&
1854 0 : & 'Action: correct nph1l in your input file.'
1855 0 : ABI_ERROR(message)
1856 : end if
1857 :
1858 73 : multibinit_dtset%nph2l=0
1859 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'nph2l',tread,'INT')
1860 73 : if(tread==1) multibinit_dtset%nph2l=intarr(1)
1861 73 : if(multibinit_dtset%nph2l<0)then
1862 : write(message, '(a,i0,a,a,a)' )&
1863 0 : & 'nph2l is',multibinit_dtset%nph2l,', which is lower than 0 .',ch10,&
1864 0 : & 'Action: correct nph2l in your input file.'
1865 0 : ABI_ERROR(message)
1866 : end if
1867 :
1868 73 : multibinit_dtset%nqshft=1
1869 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'nqshft',tread,'INT')
1870 73 : if(tread==1) multibinit_dtset%nqshft=intarr(1)
1871 73 : if(multibinit_dtset%nqshft<0 .or. multibinit_dtset%nqshft==3 .or.&
1872 : & multibinit_dtset%nqshft>=5 )then
1873 : write(message, '(a,i0,a,a,a,a,a)' )&
1874 0 : & 'nqshft is',multibinit_dtset%nqshft,', but the only allowed values',ch10,&
1875 0 : & 'are 1, 2 or 4 .',ch10,&
1876 0 : & 'Action: correct nqshft in your input file.'
1877 0 : ABI_ERROR(message)
1878 : end if
1879 :
1880 73 : multibinit_dtset%nnos=0
1881 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'nnos',tread,'INT')
1882 73 : if(tread==1) multibinit_dtset%nnos=intarr(1)
1883 73 : if(multibinit_dtset%nnos<0)then
1884 : write(message, '(a,i0,a,a,a)' )&
1885 0 : & 'nnos is',multibinit_dtset%nnos,', which is lower than 0',ch10,&
1886 0 : & 'Action: correct nnos in your input file.'
1887 0 : ABI_ERROR(message)
1888 : end if
1889 :
1890 :
1891 73 : multibinit_dtset%nsphere=0
1892 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'nsphere',tread,'INT')
1893 73 : if(tread==1) multibinit_dtset%nsphere=intarr(1)
1894 73 : if(multibinit_dtset%nsphere<0)then
1895 : write(message, '(a,i0,a,a,a)' )&
1896 0 : & 'nsphere is',multibinit_dtset%nsphere,', which is lower than 0',ch10,&
1897 0 : & 'Action: correct nsphere in your input file.'
1898 0 : ABI_ERROR(message)
1899 : end if
1900 :
1901 : !O
1902 :
1903 73 : multibinit_dtset%optcell=0
1904 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'optcell',tread,'INT')
1905 73 : if(tread==1) multibinit_dtset%optcell=intarr(1)
1906 73 : if(multibinit_dtset%optcell<0.or.multibinit_dtset%optcell>9)then
1907 : write(message, '(a,i8,a,a,a,a,a)' )&
1908 0 : & 'optcell is',multibinit_dtset%optcell,', but the only allowed values',ch10,&
1909 0 : & 'are 0, 1, 2, 3, 4, 5, 6 , 7, 8, 9.',ch10,&
1910 0 : & 'Action: correct optcell in your input file.'
1911 0 : ABI_ERROR(message)
1912 : end if
1913 :
1914 :
1915 73 : multibinit_dtset%opt_effpot=0
1916 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'opt_effpot',tread,'INT')
1917 73 : if(tread==1) multibinit_dtset%opt_effpot=intarr(1)
1918 : if(multibinit_dtset%opt_effpot<0 .and. multibinit_dtset%opt_effpot>1)then
1919 : write(message, '(a,i8,a,a,a,a,a)' )&
1920 : & 'opt_effpot is',multibinit_dtset%opt_effpot,', but the only allowed values',ch10,&
1921 : & 'are 0 and 1.',ch10,&
1922 : & 'Action: correct opt_effpot in your input file.'
1923 : ABI_ERROR(message)
1924 : end if
1925 :
1926 :
1927 73 : multibinit_dtset%opt_ncoeff=0
1928 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'opt_ncoeff',tread,'INT')
1929 73 : if(tread==1) multibinit_dtset%opt_ncoeff=intarr(1)
1930 73 : if(multibinit_dtset%opt_ncoeff<0)then
1931 : write(message, '(a,i8,a,a,a,a,a)' )&
1932 0 : & 'opt_ncoeff is',multibinit_dtset%opt_ncoeff,', but the only positive values',ch10,&
1933 0 : & 'are allowed for multibinit.',ch10,&
1934 0 : & 'Action: correct opt_ncoeff in your input file.'
1935 0 : ABI_ERROR(message)
1936 : end if
1937 :
1938 292 : multibinit_dtset%opt_EFS=(/0,1,1/)
1939 292 : multibinit_dtset%opt_on = (/ .TRUE.,.TRUE.,.FALSE. /)
1940 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'opt_EFS',tread,'INT')
1941 73 : if(tread==1) multibinit_dtset%opt_EFS(1:3)=intarr(1:3)
1942 73 : if(multibinit_dtset%opt_EFS(1) == 1)multibinit_dtset%opt_on(3) = .TRUE.
1943 73 : if(multibinit_dtset%opt_EFS(2) == 0)multibinit_dtset%opt_on(1) = .FALSE.
1944 73 : if(multibinit_dtset%opt_EFS(3) == 0)multibinit_dtset%opt_on(2) = .FALSE.
1945 584 : if(any(multibinit_dtset%opt_EFS<0) .or. any(multibinit_dtset%opt_EFS>1))then
1946 : write(message, '(a,i8,a,a,a)' )&
1947 0 : & 'opt_EFS is',multibinit_dtset%opt_EFS,', but the only allowed values are 0 and 1',ch10,&
1948 0 : & 'Action: correct opt_EFS in your input file.'
1949 0 : ABI_ERROR(message)
1950 : endif
1951 :
1952 292 : multibinit_dtset%opt_factors=(/1,1,1/)
1953 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'opt_factors',tread,'DPR')
1954 73 : if(tread==1) multibinit_dtset%opt_factors(1:3)=dprarr(1:3)
1955 292 : if(any(multibinit_dtset%opt_factors<0))then
1956 : write(message, '(a,i8,a,a,a)' )&
1957 0 : & 'opt_factors is',multibinit_dtset%opt_factors,', but the only allowed values are positive',ch10,&
1958 0 : & 'Action: correct opt_factors in your input file.'
1959 0 : ABI_ERROR(message)
1960 : end if
1961 :
1962 : !P
1963 73 : multibinit_dtset%prt_model=0
1964 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'prt_model',tread,'INT')
1965 73 : if(tread==1) multibinit_dtset%prt_model=intarr(1)
1966 73 : if(multibinit_dtset%prt_model<0.or.multibinit_dtset%prt_model>4)then
1967 : write(message, '(a,i8,a,a,a,a,a)' )&
1968 0 : & 'prt_model is',multibinit_dtset%prt_model,', but the only allowed values',ch10,&
1969 0 : & 'are 0, 1 or 2.',ch10,&
1970 0 : & 'Action: correct prt_model in your input file.'
1971 0 : ABI_ERROR(message)
1972 : end if
1973 :
1974 :
1975 73 : multibinit_dtset%prt_phfrq=0
1976 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'prt_phfrq',tread,'INT')
1977 73 : if(tread==1) multibinit_dtset%prt_phfrq=intarr(1)
1978 73 : if(multibinit_dtset%prt_phfrq<0.or.multibinit_dtset%prt_phfrq>2)then
1979 : write(message, '(a,i8,a,a,a,a,a)' )&
1980 0 : & 'prt_phfrq is',multibinit_dtset%prt_phfrq,', but the only allowed values',ch10,&
1981 0 : & 'are 0, 1 or 2.',ch10,&
1982 0 : & 'Action: correct prt_phfrq in your input file.'
1983 0 : ABI_ERROR(message)
1984 : end if
1985 :
1986 73 : multibinit_dtset%fit_initializeData=0
1987 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_initializeData',tread,'INT')
1988 73 : if(tread==1) multibinit_dtset%fit_initializeData=intarr(1)
1989 73 : if(multibinit_dtset%fit_initializeData<0.or.multibinit_dtset%fit_initializeData>1)then
1990 : write(message, '(a,i8,a,a,a,a,a)' )&
1991 0 : & 'fit_initializeData is',multibinit_dtset%fit_initializeData,', but the only allowed values',ch10,&
1992 0 : & 'are 0, 1 or 2.',ch10,&
1993 0 : & 'Action: correct fit_initializeData in your input file.'
1994 0 : ABI_ERROR(message)
1995 : end if
1996 :
1997 :
1998 73 : multibinit_dtset%fit_generateCoeff=1
1999 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_generateCoeff',tread,'INT')
2000 73 : if(tread==1) multibinit_dtset%fit_generateCoeff=intarr(1)
2001 73 : if(multibinit_dtset%fit_generateCoeff<0.or.multibinit_dtset%fit_generateCoeff>1)then
2002 : write(message, '(a,i8,a,a,a,a,a)' )&
2003 0 : & 'fit_generateCoeff is',multibinit_dtset%fit_generateCoeff,', but the only allowed values',ch10,&
2004 0 : & 'are 0, 1 or 2.',ch10,&
2005 0 : & 'Action: correct fit_generateCoeff in your input file.'
2006 0 : ABI_ERROR(message)
2007 : end if
2008 :
2009 : !Default is no output of GF values per processor to csv file
2010 73 : multibinit_dtset%prt_GF_csv = .FALSE.
2011 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'prt_GF_csv',tread,'INT')
2012 73 : if(tread==1)then
2013 1 : if(intarr(1) == 1) multibinit_dtset%prt_GF_csv = .TRUE.
2014 1 : if(intarr(1) == 0) multibinit_dtset%prt_GF_csv = .FALSE.
2015 1 : if(intarr(1) < 0 .or. intarr(1) > 1) then
2016 : write(message, '(a,i0,a,a,a,a,a)' )&
2017 0 : & 'prt_GF_csv is',intarr(1),'. The only allowed values',ch10,&
2018 0 : & 'are 0 (no output) or 1 (print GF values per processor into csv files)',ch10, &
2019 0 : & 'Action: correct prt_GF_csv in your input file.'
2020 0 : ABI_ERROR(message)
2021 : end if
2022 : endif
2023 : !Default is no output of the real space IFC to file
2024 73 : multibinit_dtset%prt_ifc = 0
2025 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'prt_ifc',tread,'INT')
2026 73 : if(tread==1) multibinit_dtset%prt_ifc = intarr(1)
2027 73 : if(multibinit_dtset%prt_ifc < 0 .or. multibinit_dtset%prt_ifc > 1) then
2028 : write(message, '(a,i0,a,a,a,a,a)' )&
2029 0 : & 'prtf_ifc is',multibinit_dtset%prt_ifc,'. The only allowed values',ch10,&
2030 0 : & 'are 0 (no output) or 1 (AI2PS format)',ch10, &
2031 0 : & 'Action: correct prt_ifc in your input file.'
2032 0 : ABI_ERROR(message)
2033 : end if
2034 :
2035 :
2036 73 : multibinit_dtset%randomseed= 0
2037 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'randomseed',tread,'INT')
2038 73 : if(tread==1) multibinit_dtset%randomseed= intarr(1)
2039 :
2040 73 : if(multibinit_dtset%randomseed /=0) then
2041 : write(message, "(a, a, a, a)") &
2042 6 : "The random seed is set to an fixed number, which might lead to wrong result if you're", &
2043 6 : & "performing a series of Monte Carlo/dynamics calculation as the random number will repeat. Only", &
2044 6 : & "do this when you want to repeat previous results. And note that it is only meaningful when", &
2045 12 : & "the number of cpu cores is consistent."
2046 6 : ABI_WARNING(message)
2047 : end if
2048 :
2049 :
2050 :
2051 :
2052 : !Default is no output of the 3rd derivative
2053 73 : multibinit_dtset%strcpling = -1
2054 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'strcpling',tread,'INT')
2055 73 : if(tread==1) multibinit_dtset%strcpling = intarr(1)
2056 73 : if(multibinit_dtset%strcpling < -1 .or. multibinit_dtset%strcpling > 2) then
2057 : write(message, '(a,i0,a,a,a,a,a,a,a)' )&
2058 0 : & 'prtf_3rd is ',multibinit_dtset%strcpling,'. The only allowed values',ch10,&
2059 0 : & 'are 0 (no computation), 1 (only computation)',ch10,&
2060 0 : & 'or 2 (computation and print in xml file)',ch10, &
2061 0 : & 'Action: correct strcpling in your input file.'
2062 0 : ABI_ERROR(message)
2063 : end if
2064 :
2065 : !Q
2066 292 : multibinit_dtset%qrefine=1 ! default is no refinement
2067 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'qrefine',tread,'INT')
2068 73 : if(tread==1) multibinit_dtset%qrefine = intarr(1:3)
2069 292 : do ii=1,3
2070 292 : if(multibinit_dtset%qrefine(ii) < 1) then
2071 : write(message, '(a,3i0,a,a,a,a,a)' )&
2072 0 : & 'qrefine is',multibinit_dtset%qrefine,' The only allowed values',ch10,&
2073 0 : & 'are integers >= 1 giving the refinement of the ngqpt grid',ch10,&
2074 0 : & 'Action: correct qrefine in your input file.'
2075 0 : ABI_ERROR(message)
2076 : end if
2077 : end do
2078 :
2079 : !R
2080 73 : multibinit_dtset%restartxf=0
2081 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'restartxf',tread,'INT')
2082 73 : if(tread==1) multibinit_dtset%restartxf=intarr(1)
2083 73 : if(multibinit_dtset%restartxf < -3 .or. multibinit_dtset%restartxf > 0)then
2084 : write(message, '(a,i8,a,a,a,a,a)' )&
2085 0 : & 'restartxf is',multibinit_dtset%restartxf,', but the only allowed values',ch10,&
2086 0 : & 'is -2 or 0.',ch10,&
2087 0 : & 'Action: correct restartxf in your input file.'
2088 0 : ABI_ERROR(message)
2089 : end if
2090 :
2091 73 : multibinit_dtset%rfmeth=1
2092 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'rfmeth',tread,'INT')
2093 73 : if(tread==1) multibinit_dtset%rfmeth=intarr(1)
2094 73 : if(multibinit_dtset%rfmeth<1.or.multibinit_dtset%rfmeth>2)then
2095 : write(message, '(a,i0,a,a,a,a,a)' )&
2096 0 : & 'rfmeth is',multibinit_dtset%rfmeth,', but the only allowed values',ch10,&
2097 0 : & 'are 1 or 2 . ',ch10,&
2098 0 : & 'Action: correct rfmeth in your input file.'
2099 0 : ABI_ERROR(message)
2100 : end if
2101 :
2102 73 : multibinit_dtset%rifcsph=zero
2103 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'rifcsph',tread,'DPR')
2104 73 : if(tread==1) multibinit_dtset%rifcsph=dprarr(1)
2105 73 : if(multibinit_dtset%rifcsph<-tol12)then
2106 : write(message, '(a,f10.3,a,a,a)' )&
2107 0 : & 'rifcsph is',multibinit_dtset%rifcsph,', which is lower than zero.',ch10,&
2108 0 : & 'Action: correct rifcsph in your input file.'
2109 0 : ABI_ERROR(message)
2110 : end if
2111 : ! R
2112 : ! multibinit_dtset%random_seed=-1
2113 : ! call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'random_seed',tread,'INT')
2114 : ! if(tread==1) multibinit_dtset%random_seed=intarr(1)
2115 :
2116 : !S
2117 :
2118 73 : multibinit_dtset%spin_damping=-1.0
2119 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_damping',tread,'DPR')
2120 73 : if(tread==1) multibinit_dtset%spin_damping=dprarr(1)
2121 :
2122 : ! multibinit_dtset%spin_calc_correlation_obs=0
2123 : ! call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_calc_correlation_obs',tread,'INT')
2124 : ! if(tread==1) multibinit_dtset%spin_calc_correlation_obs=intarr(1)
2125 : ! if(multibinit_dtset%spin_calc_correlation_obs>1.or.multibinit_dtset%spin_calc_correlation_obs<0)then
2126 : ! write(message, '(a,i8,a,a,a,a,a)' )&
2127 : ! & 'spin_calc_correlation_obs is',multibinit_dtset%spin_calc_correlation_obs,', but the only allowed values',ch10,&
2128 : ! & 'is 0 or 1.',ch10,&
2129 : ! & 'Action: correct spin_calc_correlation_obs in your input file.'
2130 : ! ABI_ERROR(message)
2131 : ! end if
2132 : !
2133 73 : multibinit_dtset%spin_calc_thermo_obs=1
2134 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_calc_thermo_obs',tread,'INT')
2135 73 : if(tread==1) multibinit_dtset%spin_calc_thermo_obs=intarr(1)
2136 73 : if(multibinit_dtset%spin_calc_thermo_obs>1.or.multibinit_dtset%spin_calc_thermo_obs<0)then
2137 : write(message, '(a,i8,a,a,a,a,a)' )&
2138 0 : & 'spin_calc_thermo_obs is',multibinit_dtset%spin_calc_thermo_obs,', but the only allowed values',ch10,&
2139 0 : & 'is 0 or 1.',ch10,&
2140 0 : & 'Action: correct spin_calc_thermo_obs in your input file.'
2141 0 : ABI_ERROR(message)
2142 : end if
2143 :
2144 :
2145 : ! multibinit_dtset%spin_calc_traj_obs=0
2146 : ! call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_calc_traj_obs',tread,'INT')
2147 : ! if(tread==1) multibinit_dtset%spin_calc_traj_obs=intarr(1)
2148 : ! if(multibinit_dtset%spin_calc_traj_obs>1.or.multibinit_dtset%spin_calc_traj_obs<0)then
2149 : ! write(message, '(a,i8,a,a,a,a,a)' )&
2150 : ! & 'spin_calc_traj_obs is',multibinit_dtset%spin_calc_traj_obs,', but the only allowed values',ch10,&
2151 : ! & 'is 0 or 1.',ch10,&
2152 : ! & 'Action: correct spin_calc_traj_obs in your input file.'
2153 : ! ABI_ERROR(message)
2154 : ! end if
2155 :
2156 :
2157 73 : multibinit_dtset%spin_dipdip=0
2158 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_dipdip',tread,'INT')
2159 73 : if(tread==1) multibinit_dtset%spin_dipdip=intarr(1)
2160 73 : if(multibinit_dtset%spin_dipdip>1.or.multibinit_dtset%spin_dipdip<0)then
2161 : write(message, '(a,i8,a,a,a,a,a)' )&
2162 0 : & 'spin_dipdip is',multibinit_dtset%spin_dipdip,', but the only allowed values',ch10,&
2163 0 : & 'is 0 or 1.',ch10,&
2164 0 : & 'Action: correct spin_dipdip in your input file.'
2165 0 : ABI_ERROR(message)
2166 : end if
2167 :
2168 :
2169 73 : multibinit_dtset%spin_dt= 1d-16
2170 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_dt',tread,'TIM')
2171 73 : if(tread==1) multibinit_dtset%spin_dt=dprarr(1)
2172 73 : if(multibinit_dtset%spin_dt<0)then
2173 : write(message, '(a,es10.2,a,a,a,a,a)' )&
2174 0 : & 'spin_dt is',multibinit_dtset%spin_dt,', but the only allowed values',ch10,&
2175 0 : & 'are superior to 0 .',ch10,&
2176 0 : & 'Action: correct spin_dt in your input file.'
2177 0 : ABI_ERROR(message)
2178 : end if
2179 :
2180 :
2181 73 : multibinit_dtset%spin_dynamics=0
2182 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_dynamics',tread,'INT')
2183 73 : if(tread==1) multibinit_dtset%spin_dynamics=intarr(1)
2184 73 : if( .not. (multibinit_dtset%spin_dynamics <= 3 .or. multibinit_dtset%spin_dynamics==20) ) then
2185 : write(message, '(a,i8,a,a,a,a,a)' )&
2186 0 : & 'spin_dynamics is ',multibinit_dtset%spin_dynamics,', but the only allowed values',ch10,&
2187 0 : & 'are 0, 1, 2, 3 and 20 and negative values.',ch10,&
2188 0 : & 'Action: correct spin_dynamics in your input file.'
2189 0 : ABI_ERROR(message)
2190 : end if
2191 :
2192 73 : if(multibinit_dtset%spin_dynamics == 20) then
2193 : write(message, '(a,i8,a)' )&
2194 0 : & 'spin_dynamics is ',multibinit_dtset%spin_dynamics,', spins will not move. For test only!!'
2195 0 : ABI_WARNING(message)
2196 : end if
2197 :
2198 :
2199 :
2200 73 : multibinit_dtset%spin_init_state=1
2201 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_init_state',tread,'INT')
2202 73 : if(tread==1) multibinit_dtset%spin_init_state=intarr(1)
2203 73 : if(multibinit_dtset%spin_init_state<1 .or. &
2204 : & multibinit_dtset%spin_init_state>4) then
2205 : write(message, '(a,i8,a,a,a,a,a)' )&
2206 0 : & 'spin_init_state is',multibinit_dtset%spin_init_state,', but the only allowed values',ch10,&
2207 0 : & 'are 1, 2, 3, and 4.',ch10,&
2208 0 : & 'Action: correct spin_init_state in your input file.'
2209 0 : ABI_ERROR(message)
2210 : end if
2211 :
2212 :
2213 :
2214 292 : multibinit_dtset%spin_mag_field= zero
2215 73 : if(3>marr)then
2216 0 : marr=3
2217 0 : ABI_FREE(intarr)
2218 0 : ABI_FREE(dprarr)
2219 0 : ABI_MALLOC(intarr,(marr))
2220 0 : ABI_MALLOC(dprarr,(marr))
2221 : end if
2222 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'spin_mag_field',tread,'BFI')
2223 79 : if(tread==1) multibinit_dtset%spin_mag_field(1:3)= dprarr(1:3)
2224 :
2225 73 : multibinit_dtset%spin_nctime=100
2226 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_nctime',tread,'INT')
2227 73 : if(tread==1) multibinit_dtset%spin_nctime=intarr(1)
2228 73 : if(multibinit_dtset%spin_nctime<=0)then
2229 : write(message, '(a,i0,a,a,a)' )&
2230 0 : & 'spin_nctime is ',multibinit_dtset%spin_nctime,', which is lower than 1 .',ch10,&
2231 0 : & 'Action: correct spin_nctime in your input file.'
2232 0 : ABI_ERROR(message)
2233 : end if
2234 :
2235 73 : multibinit_dtset%spin_ntime_pre=0
2236 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_ntime_pre',tread,'INT')
2237 73 : if(tread==1) multibinit_dtset%spin_ntime_pre=intarr(1)
2238 73 : if(multibinit_dtset%spin_ntime_pre<0)then
2239 : write(message, '(a,i0,a,a,a)' )&
2240 0 : & 'spin_ntime_pre is',multibinit_dtset%spin_ntime_pre,', which is lower than 0 .',ch10,&
2241 0 : & 'Action: correct spin_ntime_pre in your input file.'
2242 0 : ABI_ERROR(message)
2243 : end if
2244 :
2245 :
2246 :
2247 73 : multibinit_dtset%spin_ntime=10000
2248 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_ntime',tread,'INT')
2249 73 : if(tread==1) multibinit_dtset%spin_ntime=intarr(1)
2250 73 : if(multibinit_dtset%spin_ntime<0)then
2251 : write(message, '(a,i0,a,a,a)' )&
2252 0 : & 'spin_ntime is',multibinit_dtset%spin_ntime,', which is lower than 0 .',ch10,&
2253 0 : & 'Action: correct spin_ntime in your input file.'
2254 0 : ABI_ERROR(message)
2255 : end if
2256 :
2257 :
2258 : ! multibinit_dtset%spin_n1l=1
2259 : ! call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_n1l',tread,'INT')
2260 : ! if(tread==1) multibinit_dtset%spin_n1l=intarr(1)
2261 : ! if(multibinit_dtset%spin_n1l<0)then
2262 : ! write(message, '(a,i0,a,a,a)' )&
2263 : ! & 'spin_n1l is',multibinit_dtset%spin_n1l,', which is lower than 0 .',ch10,&
2264 : ! & 'Action: correct spin_n1l in your input file.'
2265 : ! ABI_ERROR(message)
2266 : ! end if
2267 : !
2268 : ! multibinit_dtset%spin_n2l=0
2269 : ! call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_n2l',tread,'INT')
2270 : ! if(tread==1) multibinit_dtset%spin_n2l=intarr(1)
2271 : ! if(multibinit_dtset%spin_n2l<0)then
2272 : ! write(message, '(a,i0,a,a,a)' )&
2273 : ! & 'spin_n2l is',multibinit_dtset%spin_n2l,', which is lower than 0 .',ch10,&
2274 : ! & 'Action: correct spin_n2l in your input file.'
2275 : ! ABI_ERROR(message)
2276 : ! end if
2277 :
2278 292 : multibinit_dtset%spin_init_orientation= [0.0, 0.0, 1.0]
2279 73 : if(3>marr)then
2280 0 : marr=3
2281 0 : ABI_FREE(intarr)
2282 0 : ABI_FREE(dprarr)
2283 0 : ABI_MALLOC(intarr,(marr))
2284 0 : ABI_MALLOC(dprarr,(marr))
2285 : end if
2286 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'spin_init_orientation',tread,'DPR')
2287 79 : if(tread==1) multibinit_dtset%spin_init_orientation(1:3)= dprarr(1:3)
2288 :
2289 292 : multibinit_dtset%spin_projection_qpoint= zero
2290 73 : if(3>marr)then
2291 0 : marr=3
2292 0 : ABI_FREE(intarr)
2293 0 : ABI_FREE(dprarr)
2294 0 : ABI_MALLOC(intarr,(marr))
2295 0 : ABI_MALLOC(dprarr,(marr))
2296 : end if
2297 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'spin_projection_qpoint',tread,'DPR')
2298 79 : if(tread==1) multibinit_dtset%spin_projection_qpoint(1:3)= dprarr(1:3)
2299 :
2300 292 : multibinit_dtset%spin_init_qpoint= zero
2301 73 : if(3>marr)then
2302 0 : marr=3
2303 0 : ABI_FREE(intarr)
2304 0 : ABI_FREE(dprarr)
2305 0 : ABI_MALLOC(intarr,(marr))
2306 0 : ABI_MALLOC(dprarr,(marr))
2307 : end if
2308 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'spin_init_qpoint',tread,'DPR')
2309 79 : if(tread==1) multibinit_dtset%spin_init_qpoint(1:3)= dprarr(1:3)
2310 :
2311 292 : multibinit_dtset%spin_init_rotate_axis= [1.0, 0.0, 0.0]
2312 73 : if(3>marr)then
2313 0 : marr=3
2314 0 : ABI_FREE(intarr)
2315 0 : ABI_FREE(dprarr)
2316 0 : ABI_MALLOC(intarr,(marr))
2317 0 : ABI_MALLOC(dprarr,(marr))
2318 : end if
2319 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'spin_init_rotate_axis',tread,'DPR')
2320 79 : if(tread==1) multibinit_dtset%spin_init_rotate_axis(1:3)= dprarr(1:3)
2321 :
2322 73 : multibinit_dtset%spin_sia_add=0
2323 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_sia_add',tread,'INT')
2324 73 : if(tread==1) multibinit_dtset%spin_sia_add=intarr(1)
2325 73 : if(multibinit_dtset%spin_sia_add <0 .or. multibinit_dtset%spin_sia_add>2 )then
2326 : write(message, '(a,i0,a,a,a)' )&
2327 0 : & 'spin_sia_add is',multibinit_dtset%spin_sia_add,', which is not 0, 1, or 2.',ch10,&
2328 0 : & 'Action: correct spin_sia_add in your input file.'
2329 0 : ABI_ERROR(message)
2330 : end if
2331 :
2332 73 : multibinit_dtset%spin_sia_k1amp=0.0
2333 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_sia_k1amp',tread,'ENE')
2334 73 : if(tread==1) multibinit_dtset%spin_sia_k1amp=dprarr(1)
2335 :
2336 292 : multibinit_dtset%spin_sia_k1dir(:)= [0.0,0.0,1.0]
2337 73 : if(3>marr)then
2338 0 : marr=3
2339 0 : ABI_FREE(intarr)
2340 0 : ABI_FREE(dprarr)
2341 0 : ABI_MALLOC(intarr,(marr))
2342 0 : ABI_MALLOC(dprarr,(marr))
2343 : end if
2344 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'spin_sia_k1dir',tread,'DPR')
2345 73 : if(tread==1) then
2346 0 : dprarr(1:3)=dprarr(1:3)/sqrt(sum(dprarr(1:3)**2))
2347 0 : multibinit_dtset%spin_sia_k1dir(1:3)= dprarr(1:3)
2348 : endif
2349 :
2350 73 : multibinit_dtset%spin_temperature=325
2351 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_temperature',tread,'DPR')
2352 73 : if(tread==1) multibinit_dtset%spin_temperature=dprarr(1)
2353 73 : if(multibinit_dtset%spin_temperature<0)then
2354 : write(message, '(a,f10.1,a,a,a,a,a)' )&
2355 0 : & 'spin_temperature is ',multibinit_dtset%spin_temperature,'. The only allowed values',ch10,&
2356 0 : & 'are non-negative values.',ch10,&
2357 0 : & 'Action: correct spin_temperature in your input file.'
2358 0 : ABI_ERROR(message)
2359 : end if
2360 :
2361 73 : multibinit_dtset%spin_temperature_start=0.0
2362 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_temperature_start',tread,'DPR')
2363 73 : if(tread==1) multibinit_dtset%spin_temperature_start=dprarr(1)
2364 73 : if(multibinit_dtset%spin_temperature_start<0.0)then
2365 : write(message, '(a,f10.1,a,a,a,a,a)' )&
2366 0 : & 'spin_temperature_start is ',multibinit_dtset%spin_temperature_start,'. The only allowed values',ch10,&
2367 0 : & 'are positives values.',ch10,&
2368 0 : & 'Action: correct spin_semperature_start in your input file.'
2369 0 : ABI_ERROR(message)
2370 : end if
2371 :
2372 73 : multibinit_dtset%spin_temperature_end=0.0
2373 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_temperature_end',tread,'DPR')
2374 73 : if(tread==1) multibinit_dtset%spin_temperature_end=dprarr(1)
2375 73 : if(multibinit_dtset%spin_temperature_end<0)then
2376 : write(message, '(a,f10.1,a,a,a,a,a)' )&
2377 0 : & 'spin_temperature_end is ',multibinit_dtset%spin_temperature_end,'. The only allowed values',ch10,&
2378 0 : & 'are positives values.',ch10,&
2379 0 : & 'Action: correct spin_semperature_end in your input file.'
2380 0 : ABI_ERROR(message)
2381 : end if
2382 :
2383 73 : multibinit_dtset%spin_temperature_nstep=1
2384 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_temperature_nstep',tread,'INT')
2385 73 : if(tread==1) multibinit_dtset%spin_temperature_nstep=intarr(1)
2386 73 : if(multibinit_dtset%spin_temperature_nstep<=0)then
2387 : write(message, '(a,i0,a,a,a,a)' )&
2388 0 : & 'spin_temperature_nstep is',multibinit_dtset%spin_temperature_nstep,', while it should be larger than 0',ch10,&
2389 0 : & 'Action: correct spin_temperature_nstep in your input file.'
2390 0 : ABI_ERROR(message)
2391 : end if
2392 :
2393 :
2394 :
2395 : ! multibinit_dtset%spin_tolavg=1d-02
2396 : ! call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_tolavg',tread,'DPR')
2397 : ! if(tread==1) multibinit_dtset%spin_tolavg=dprarr(1)
2398 : ! if(multibinit_dtset%spin_tolavg<=0)then
2399 : ! write(message, '(a,f10.1,a,a,a,a,a)' )&
2400 : ! & 'spin_tolavg is ',multibinit_dtset%spin_tolavg,'. The only allowed values',ch10,&
2401 : ! & 'are positives values.',ch10,&
2402 : ! & 'Action: correct spin_tolavg in your input file.'
2403 : ! ABI_ERROR(message)
2404 : ! end if
2405 : !
2406 : ! multibinit_dtset%spin_tolvar=1d-02
2407 : ! call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_tolvar',tread,'DPR')
2408 : ! if(tread==1) multibinit_dtset%spin_tolvar=dprarr(1)
2409 : ! if(multibinit_dtset%spin_tolvar<=0)then
2410 : ! write(message, '(a,f10.1,a,a,a,a,a)' )&
2411 : ! & 'spin_tolvar is ',multibinit_dtset%spin_tolvar,'. The only allowed values',ch10,&
2412 : ! & 'are positives values.',ch10,&
2413 : ! & 'Action: correct spin_tolvar in your input file.'
2414 : ! ABI_ERROR(message)
2415 : ! end if
2416 :
2417 73 : multibinit_dtset%spin_var_temperature=0
2418 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_var_temperature',tread,'INT')
2419 73 : if(tread==1) multibinit_dtset%spin_var_temperature=intarr(1)
2420 73 : if(multibinit_dtset%spin_var_temperature/=0.and.multibinit_dtset%spin_var_temperature/=1)then
2421 : write(message, '(a,i0,a,a,a,a,a)' )&
2422 0 : & 'spin_var_temperature is',multibinit_dtset%spin_var_temperature,'. The only allowed values',ch10,&
2423 0 : & 'are 0, or 1.',ch10,&
2424 0 : & 'Action: correct spin_var_temperature in your input file.'
2425 0 : ABI_ERROR(message)
2426 : end if
2427 :
2428 73 : multibinit_dtset%spin_write_traj=0
2429 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'spin_write_traj',tread,'INT')
2430 73 : if(tread==1) multibinit_dtset%spin_write_traj=intarr(1)
2431 73 : if(multibinit_dtset%spin_write_traj/=0.and.multibinit_dtset%spin_write_traj/=1)then
2432 : write(message, '(a,i0,a,a,a,a,a)' )&
2433 0 : & 'spin_write_traj is',multibinit_dtset%spin_write_traj,'. The only allowed values',ch10,&
2434 0 : & 'are 0, or 1.',ch10,&
2435 0 : & 'Action: correct spin_write_traj in your input file.'
2436 0 : ABI_ERROR(message)
2437 : end if
2438 :
2439 73 : multibinit_dtset%slc_coupling=0
2440 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'slc_coupling',tread,'INT')
2441 73 : if(tread==1) multibinit_dtset%slc_coupling=intarr(1)
2442 : if(multibinit_dtset%slc_coupling.ne. 1111 .and. &
2443 : & multibinit_dtset%slc_coupling.ne. 1110 .and. &
2444 : & multibinit_dtset%slc_coupling.ne. 1101 .and. &
2445 : & multibinit_dtset%slc_coupling.ne. 1011 .and. &
2446 : & multibinit_dtset%slc_coupling.ne. 111 .and. &
2447 : & multibinit_dtset%slc_coupling.ne. 1100 .and. &
2448 : & multibinit_dtset%slc_coupling.ne. 1010 .and. &
2449 : & multibinit_dtset%slc_coupling.ne. 1001 .and. &
2450 : & multibinit_dtset%slc_coupling.ne. 110 .and. &
2451 : & multibinit_dtset%slc_coupling.ne. 101 .and. &
2452 : & multibinit_dtset%slc_coupling.ne. 11 .and. &
2453 : & multibinit_dtset%slc_coupling.ne. 1000 .and. &
2454 : & multibinit_dtset%slc_coupling.ne. 100 .and. &
2455 : & multibinit_dtset%slc_coupling.ne. 10 .and. &
2456 144 : & multibinit_dtset%slc_coupling.ne. 1 .and. &
2457 : & multibinit_dtset%slc_coupling.ne. 0) then
2458 : write(message, '(a,i8,a,a,a,a,a)' )&
2459 0 : & 'slc_coupling is',multibinit_dtset%slc_coupling,', but the only allowed values',ch10,&
2460 0 : & 'are 1111, 1110, 1101, 1011, 111, 1100, 1010, 1001, 110, 101, 11, 1000, 100, 10, 1, and 0.',ch10,&
2461 0 : & 'Action: correct slc_coupling in your input file.'
2462 0 : ABI_ERROR(message)
2463 : end if
2464 :
2465 :
2466 73 : multibinit_dtset%symdynmat=1
2467 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'symdynmat',tread,'INT')
2468 73 : if(tread==1) multibinit_dtset%symdynmat=intarr(1)
2469 73 : if(multibinit_dtset%symdynmat/=0.and.multibinit_dtset%symdynmat/=1)then
2470 : write(message, '(a,i0,a,a,a,a,a)' )&
2471 0 : & 'symdynmat is',multibinit_dtset%symdynmat,'. The only allowed values',ch10,&
2472 0 : & 'are 0, or 1.',ch10,&
2473 0 : & 'Action: correct symdynmat in your input file.'
2474 0 : ABI_ERROR(message)
2475 : end if
2476 :
2477 73 : multibinit_dtset%sys_dim=1
2478 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'sys_dim',tread,'INT')
2479 73 : if(tread==1) multibinit_dtset%sys_dim=intarr(1)
2480 73 : if(multibinit_dtset%sys_dim<1.or.multibinit_dtset%sys_dim>9)then
2481 : write(message, '(a,i0,5a)' )&
2482 0 : 'sys_dim is ',multibinit_dtset%sys_dim,', but the only allowed values',ch10,&
2483 0 : 'are 1, 2, 3, 4, 5, 6, 7 or 8.',ch10,'Action: correct sys_dim in your input file.'
2484 0 : ABI_ERROR(message)
2485 : end if
2486 :
2487 73 : multibinit_dtset%strfact=100.0d0
2488 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'strfact',tread,'DPR')
2489 73 : if(tread==1) multibinit_dtset%strfact=dprarr(1)
2490 73 : if(multibinit_dtset%strfact<-tol12)then
2491 : write(message, '(a,f10.1,a,a,a,a,a)' )&
2492 0 : & 'strfact is ',multibinit_dtset%strfact,'. The only allowed values',ch10,&
2493 0 : & 'are positives values.',ch10,&
2494 0 : & 'Action: correct strfact in your input file.'
2495 0 : ABI_ERROR(message)
2496 : end if
2497 :
2498 73 : multibinit_dtset%strprecon=1.0d0
2499 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'strprecon',tread,'DPR')
2500 73 : if(tread==1) multibinit_dtset%strprecon=dprarr(1)
2501 73 : if(multibinit_dtset%strprecon<tol8 .or. multibinit_dtset%strprecon >1.0d0)then
2502 : write(message, '(a,f10.1,a,a,a,a,a)' )&
2503 0 : & 'strprecon is ',multibinit_dtset%strprecon,'. The only allowed values',ch10,&
2504 0 : & 'are positives values smaller or equal than 1.',ch10,&
2505 0 : & 'Action: correct strprecon in your input file.'
2506 0 : ABI_ERROR(message)
2507 : end if
2508 :
2509 : !T
2510 :
2511 73 : multibinit_dtset%temperature=325.0d0
2512 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'temperature',tread,'DPR')
2513 73 : if(tread==1) multibinit_dtset%temperature=dprarr(1)
2514 73 : if(multibinit_dtset%temperature<=0)then
2515 : write(message, '(a,f10.1,a,a,a,a,a)' )&
2516 0 : & 'Temperature is ',multibinit_dtset%temperature,'. The only allowed values',ch10,&
2517 0 : & 'are positives values.',ch10,&
2518 0 : & 'Action: correct Temperature in your input file.'
2519 0 : ABI_ERROR(message)
2520 : end if
2521 :
2522 :
2523 :
2524 73 : multibinit_dtset%test_effpot=0
2525 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'test_effpot',tread,'INT')
2526 73 : if(tread==1) multibinit_dtset%test_effpot=intarr(1)
2527 73 : if(multibinit_dtset%test_effpot<0 .or. multibinit_dtset%test_effpot>1)then
2528 : write(message, '(a,i0,a,a,a,a,a)' )&
2529 0 : & 'test_effpot is ',multibinit_dtset%test_effpot,'. The only allowed values',ch10,&
2530 0 : & 'are 0 and 1.',ch10,&
2531 0 : & 'Action: correct test_effpot in your input file.'
2532 0 : ABI_ERROR(message)
2533 : end if
2534 :
2535 73 : multibinit_dtset%test_prt_ph=0
2536 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'test_prt_ph',tread,'INT')
2537 73 : if(tread==1) multibinit_dtset%test_prt_ph=intarr(1)
2538 73 : if(multibinit_dtset%test_prt_ph<0 .or. multibinit_dtset%test_prt_ph>1)then
2539 : write(message, '(a,i0,a,a,a,a,a)' )&
2540 0 : & 'test_prt_ph is ',multibinit_dtset%test_prt_ph,'. The only allowed values',ch10,&
2541 0 : & 'are 0 and 1.',ch10,&
2542 0 : & 'Action: correct test_prt_ph in your input file.'
2543 0 : ABI_ERROR(message)
2544 : end if
2545 :
2546 73 : multibinit_dtset%tolmxf=2.0d-5
2547 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'tolmxf',tread,'DPR')
2548 73 : if(tread==1) multibinit_dtset%tolmxf=dprarr(1)
2549 73 : if(multibinit_dtset%tolmxf<0)then
2550 : write(message, '(a,i0,a,a,a,a,a)' )&
2551 0 : & 'tolmxf is ',multibinit_dtset%tolmxf,'. The only allowed values',ch10,&
2552 0 : & 'are positiv.',ch10,&
2553 0 : & 'Action: correct tolmxf in your input file.'
2554 0 : ABI_ERROR(message)
2555 : end if
2556 :
2557 : !U
2558 :
2559 : !V
2560 :
2561 : !W
2562 :
2563 : !X
2564 :
2565 : !Y
2566 :
2567 : !Z
2568 :
2569 : !=====================================================================
2570 : !end non-dependent variables
2571 : !=====================================================================
2572 :
2573 : !=======================================================================
2574 : !Read in dependent variables (dependent on dimensions above)
2575 : !=======================================================================
2576 :
2577 : !A
2578 292 : multibinit_dtset%acell= one
2579 73 : if(3>marr)then
2580 0 : marr=3
2581 0 : ABI_FREE(intarr)
2582 0 : ABI_FREE(dprarr)
2583 0 : ABI_MALLOC(intarr,(marr))
2584 0 : ABI_MALLOC(dprarr,(marr))
2585 : end if
2586 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'acell',tread,'DPR')
2587 73 : if(tread==1) multibinit_dtset%acell(1:3)= dprarr(1:3)
2588 292 : if(any(multibinit_dtset%acell<=tol10))then
2589 : write(message, '(3a)' )&
2590 0 : & 'There is negative or zero value for cell ',ch10,&
2591 0 : & 'Action: change acell in your input file.'
2592 0 : ABI_ERROR(message)
2593 : end if
2594 :
2595 73 : if(6>marr)then
2596 0 : marr=6
2597 0 : ABI_FREE(intarr)
2598 0 : ABI_FREE(dprarr)
2599 0 : ABI_MALLOC(intarr,(marr))
2600 0 : ABI_MALLOC(dprarr,(marr))
2601 : end if
2602 511 : multibinit_dtset%strtarget(1:6) = zero
2603 73 : call intagm(dprarr,intarr,jdtset,marr,6,string(1:lenstr),'strtarget',tread,'DPR')
2604 91 : if(tread==1) multibinit_dtset%strtarget(1:6)=dprarr(1:6)
2605 :
2606 :
2607 219 : ABI_MALLOC(multibinit_dtset%atifc,(natom))
2608 429 : multibinit_dtset%atifc(:)=0
2609 73 : if(multibinit_dtset%natifc>=1)then
2610 67 : if(multibinit_dtset%natifc>marr)then
2611 0 : marr=multibinit_dtset%natifc
2612 0 : ABI_FREE(intarr)
2613 0 : ABI_FREE(dprarr)
2614 0 : ABI_MALLOC(intarr,(marr))
2615 0 : ABI_MALLOC(dprarr,(marr))
2616 : end if
2617 67 : call intagm(dprarr,intarr,jdtset,marr,multibinit_dtset%natifc,string(1:lenstr),'atifc',tread,'INT')
2618 67 : if(tread==1) then
2619 0 : multibinit_dtset%atifc(1:multibinit_dtset%natifc)= intarr(1:multibinit_dtset%natifc)
2620 : else ! set to the maximum
2621 423 : do iatifc=1,multibinit_dtset%natifc
2622 423 : multibinit_dtset%atifc(iatifc) = iatifc
2623 : end do
2624 : end if
2625 201 : ABI_MALLOC(work,(natom))
2626 423 : work(:)=0
2627 :
2628 423 : do iatifc=1,multibinit_dtset%natifc
2629 356 : if(multibinit_dtset%atifc(iatifc)<=0.or.multibinit_dtset%atifc(iatifc)>natom)then
2630 : write(message, '(a,i0,a,a,a,a,a,i0,a,a,a)' )&
2631 0 : & 'For iatifc=',iatifc,', the number of the atom ifc to be ',ch10,&
2632 0 : & 'analysed is not valid : either negative, ',ch10,&
2633 0 : & 'zero, or larger than natom =',natom,'.',ch10,&
2634 0 : & 'Action: change atifc in your input file.'
2635 0 : ABI_ERROR(message)
2636 : end if
2637 423 : work(multibinit_dtset%atifc(iatifc))=1
2638 : end do
2639 423 : multibinit_dtset%atifc(1:natom)=int(work(:))
2640 67 : ABI_FREE(work)
2641 : end if
2642 :
2643 : !B
2644 :
2645 : !C
2646 219 : ABI_MALLOC(multibinit_dtset%coefficients,(multibinit_dtset%ncoeff))
2647 73 : if (multibinit_dtset%ncoeff/=0)then
2648 8 : if(multibinit_dtset%ncoeff>marr)then
2649 0 : marr=multibinit_dtset%ncoeff
2650 0 : ABI_FREE(intarr)
2651 0 : ABI_FREE(dprarr)
2652 0 : ABI_MALLOC(intarr,(marr))
2653 0 : ABI_MALLOC(dprarr,(marr))
2654 : end if
2655 40 : multibinit_dtset%coefficients(:)=zero
2656 : call intagm(dprarr,intarr,jdtset,marr,multibinit_dtset%ncoeff,&
2657 8 : & string(1:lenstr),'coefficients',tread,'DPR')
2658 8 : if(tread==1)then
2659 40 : do ii=1,multibinit_dtset%ncoeff
2660 40 : multibinit_dtset%coefficients(ii)=dprarr(ii)
2661 : end do
2662 : end if
2663 : end if
2664 :
2665 219 : ABI_MALLOC(multibinit_dtset%conf_cutoff_disp,(multibinit_dtset%natom))
2666 73 : if (multibinit_dtset%natom/=0)then
2667 67 : if(multibinit_dtset%natom>marr)then
2668 0 : marr=multibinit_dtset%natom
2669 0 : ABI_FREE(intarr)
2670 0 : ABI_FREE(dprarr)
2671 0 : ABI_MALLOC(intarr,(marr))
2672 0 : ABI_MALLOC(dprarr,(marr))
2673 : end if
2674 423 : multibinit_dtset%conf_cutoff_disp(:)=zero
2675 : call intagm(dprarr,intarr,jdtset,marr,multibinit_dtset%natom,&
2676 67 : & string(1:lenstr),'conf_cutoff_disp',tread,'DPR')
2677 67 : if(tread==1)then
2678 0 : do ii=1,multibinit_dtset%natom
2679 0 : multibinit_dtset%conf_cutoff_disp(ii)=dprarr(ii)
2680 : end do
2681 : end if
2682 423 : if(any(multibinit_dtset%conf_cutoff_disp<zero))then
2683 : write(message, '(3a)' )&
2684 0 : & 'There is negative value for conf_cutoff_disp ',ch10,&
2685 0 : & 'Action: change acell in your input file.'
2686 0 : ABI_ERROR(message)
2687 : end if
2688 : end if
2689 :
2690 73 : if(6>marr)then
2691 0 : marr=6
2692 0 : ABI_FREE(intarr)
2693 0 : ABI_FREE(dprarr)
2694 0 : ABI_MALLOC(intarr,(marr))
2695 0 : ABI_MALLOC(dprarr,(marr))
2696 : end if
2697 511 : multibinit_dtset%conf_cutoff_strain(1:6) = zero
2698 73 : call intagm(dprarr,intarr,jdtset,marr,6,string(1:lenstr),'conf_cutoff_strain',tread,'DPR')
2699 73 : if(tread==1) multibinit_dtset%conf_cutoff_strain(1:6)=dprarr(1:6)
2700 429 : if(any(multibinit_dtset%conf_cutoff_disp<zero))then
2701 : write(message, '(3a)' )&
2702 0 : & 'There is negative value for conf_cutoff_strain ',ch10,&
2703 0 : & 'Action: change acell in your input file.'
2704 0 : ABI_ERROR(message)
2705 : end if
2706 :
2707 : !D
2708 292 : multibinit_dtset%dipdip_range(:)= (/-375,-375,-375/)
2709 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'dipdip_range',tread,'INT')
2710 73 : if(tread==1) multibinit_dtset%dipdip_range(1:3)=intarr(1:3)
2711 : ! do ii=1,3
2712 : ! if(multibinit_dtset%dipdip_range(ii)<0.or.multibinit_dtset%dipdip_range(ii)>50)then
2713 : ! write(message, '(a,i0,a,i0,4a,i0,a)' )&
2714 : !& 'dipdip_range(',ii,') is ',multibinit_dtset%dipdip_range(ii),', which is lower',&
2715 : !& ' than 0 of superior than 50.',&
2716 : !& ch10,'Action: correct dipdip_range(',ii,') in your input file.'
2717 : ! ABI_ERROR(message)
2718 : ! end if
2719 : ! end do
2720 :
2721 292 : if(any(multibinit_dtset%dipdip_range /=-375)) then
2722 0 : ABI_ERROR("The use of the keyword dipdip_range has been deprecated. Please remove it from the input file!")
2723 : end if
2724 292 : multibinit_dtset%dipdip_range(:)=[0,0,0]
2725 :
2726 :
2727 :
2728 :
2729 : !E
2730 73 : multibinit_dtset%eivec=0
2731 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'eivec',tread,'INT')
2732 73 : if(tread==1) multibinit_dtset%eivec=intarr(1)
2733 73 : if(multibinit_dtset%eivec<0.or.multibinit_dtset%eivec>4)then
2734 : write(message, '(a,i0,a,a,a,a,a)' )&
2735 0 : & 'eivec is',multibinit_dtset%eivec,', but the only allowed values',ch10,&
2736 0 : & 'are 0, 1, 2, 3 or 4.',ch10,&
2737 0 : & 'Action: correct eivec in your input file.'
2738 0 : ABI_ERROR(message)
2739 : end if
2740 :
2741 :
2742 :
2743 : !F
2744 73 : multibinit_dtset%fit_anhaStrain=0
2745 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_anhaStrain',tread,'INT')
2746 73 : if(tread==1) multibinit_dtset%fit_anhaStrain=intarr(1)
2747 : if(multibinit_dtset%fit_anhaStrain<0.and.multibinit_dtset%fit_anhaStrain>1)then
2748 : write(message, '(a,i8,a,a,a,a,a)' )&
2749 : & 'fit_anhaStrain is',multibinit_dtset%fit_anhaStrain,', but the only allowed values',ch10,&
2750 : & 'are 0 or 1 for multibinit.',ch10,&
2751 : & 'Action: correct fit_anhaStrain in your input file.'
2752 : ABI_ERROR(message)
2753 : end if
2754 :
2755 73 : multibinit_dtset%bound_model=0
2756 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'bound_model',tread,'INT')
2757 73 : if(tread==1) multibinit_dtset%bound_model=intarr(1)
2758 : if(multibinit_dtset%bound_model<0.and.multibinit_dtset%bound_model>3)then
2759 : write(message, '(a,i8,a,a,a,a,a)' )&
2760 : & 'bound_model is',multibinit_dtset%bound_model,', but the only allowed values',ch10,&
2761 : & 'are between 0 and 3 for multibinit.',ch10,&
2762 : & 'Action: correct bound_model in your input file.'
2763 : ABI_ERROR(message)
2764 : end if
2765 :
2766 73 : multibinit_dtset%bound_penalty=1.001d+0
2767 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'bound_penalty',tread,'DPR')
2768 73 : if(tread==1) multibinit_dtset%bound_penalty=dprarr(1)
2769 73 : if(multibinit_dtset%bound_penalty<1)then
2770 : write(message, '(a,i8,a,a,a,a,a)' )&
2771 0 : & 'bound_penalty is',multibinit_dtset%bound_penalty,', but the only allowed values',ch10,&
2772 0 : & 'are between larger than 1.',ch10,&
2773 0 : & 'Action: correct bound_penalty in your input file.'
2774 0 : ABI_ERROR(message)
2775 : end if
2776 73 : multibinit_dtset%bound_anhaStrain=0
2777 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_anhaStrain',tread,'INT')
2778 73 : if(tread==1) multibinit_dtset%bound_anhaStrain=intarr(1)
2779 : if(multibinit_dtset%bound_anhaStrain<0.and.multibinit_dtset%bound_anhaStrain>1)then
2780 : write(message, '(a,i8,a,a,a,a,a)' )&
2781 : & 'fit_anhaStrain is',multibinit_dtset%bound_anhaStrain,', but the only allowed values',ch10,&
2782 : & 'are 0 or 1 for multibinit.',ch10,&
2783 : & 'Action: correct fit_anhaStrain in your input file.'
2784 : ABI_ERROR(message)
2785 : end if
2786 :
2787 73 : multibinit_dtset%bound_SPCoupling=1
2788 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'bound_SPCoupling',tread,'INT')
2789 73 : if(tread==1) multibinit_dtset%bound_SPCoupling=intarr(1)
2790 : if(multibinit_dtset%bound_SPCoupling<0.and.multibinit_dtset%bound_SPCoupling>1)then
2791 : write(message, '(a,i8,a,a,a,a,a)' )&
2792 : & 'bound_SPCoupling is',multibinit_dtset%bound_SPCoupling,&
2793 : & ', but the only allowed values',ch10,&
2794 : & 'are 0 or 1 for multibinit.',ch10,&
2795 : & 'Action: correct bound_SPCoupling in your input file.'
2796 : ABI_ERROR(message)
2797 : end if
2798 :
2799 73 : multibinit_dtset%fit_dispterms=1
2800 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_dispterms',tread,'INT')
2801 73 : if(tread==1) multibinit_dtset%fit_dispterms=intarr(1)
2802 : if(multibinit_dtset%fit_dispterms<0.and.multibinit_dtset%fit_dispterms>1)then
2803 : write(message, '(a,i8,a,a,a,a,a)' )&
2804 : & 'fit_dispterms is',multibinit_dtset%fit_dispterms,&
2805 : & ', but the only allowed values',ch10,&
2806 : & 'are 0 or 1 for multibinit.',ch10,&
2807 : & 'Action: correct fit_dispterms in your input file.'
2808 : ABI_ERROR(message)
2809 : end if
2810 :
2811 73 : multibinit_dtset%fit_SPCoupling=1
2812 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_SPCoupling',tread,'INT')
2813 73 : if(tread==1) multibinit_dtset%fit_SPCoupling=intarr(1)
2814 : if(multibinit_dtset%fit_SPCoupling<0.and.multibinit_dtset%fit_SPCoupling>1)then
2815 : write(message, '(a,i8,a,a,a,a,a)' )&
2816 : & 'fit_SPCoupling is',multibinit_dtset%fit_SPCoupling,&
2817 : & ', but the only allowed values',ch10,&
2818 : & 'are 0 or 1 for multibinit.',ch10,&
2819 : & 'Action: correct fit_SPCoupling in your input file.'
2820 : ABI_ERROR(message)
2821 : end if
2822 :
2823 73 : multibinit_dtset%fit_SPC_maxS=1
2824 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_SPC_maxS',tread,'INT')
2825 73 : if(tread==1) multibinit_dtset%fit_SPC_maxS=intarr(1)
2826 73 : if(multibinit_dtset%fit_SPC_maxS<1)then
2827 : write(message, '(a,i8,a,a,a,a,a)' )&
2828 0 : & 'fit_SPC_maxS is',multibinit_dtset%fit_SPC_maxS,&
2829 0 : & ', but the only allowed values',ch10,&
2830 0 : & 'are positiv.',ch10,&
2831 0 : & 'Action: correct fit_SPC_maxS in your input file.'
2832 0 : ABI_ERROR(message)
2833 : end if
2834 :
2835 73 : multibinit_dtset%bound_cutoff=0
2836 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'bound_cutoff',tread,'DPR')
2837 73 : if(tread==1) multibinit_dtset%bound_cutoff=dprarr(1)
2838 73 : if(multibinit_dtset%bound_cutoff<0)then
2839 : write(message, '(a,i8,a,a,a,a,a)' )&
2840 0 : & 'bound_cutoff is',multibinit_dtset%bound_cutoff,', but the only allowed values',ch10,&
2841 0 : & 'are positives for multibinit.',ch10,&
2842 0 : & 'Action: correct bound_cutoff in your input file.'
2843 0 : ABI_ERROR(message)
2844 : end if
2845 :
2846 :
2847 73 : multibinit_dtset%bound_maxCoeff=4
2848 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'bound_maxCoeff',tread,'INT')
2849 73 : if(tread==1) multibinit_dtset%bound_maxCoeff=intarr(1)
2850 : if(multibinit_dtset%bound_maxCoeff<0.and.multibinit_dtset%bound_maxCoeff>1)then
2851 : write(message, '(a,i8,a,a,a,a,a)' )&
2852 : & 'bound_maxCoeff is',multibinit_dtset%bound_maxCoeff,', but the only allowed values',ch10,&
2853 : & 'are 0 or 1 for multibinit.',ch10,&
2854 : & 'Action: correct bound_maxCoeff in your input file.'
2855 : ABI_ERROR(message)
2856 : end if
2857 :
2858 73 : multibinit_dtset%bound_option=1
2859 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'bound_option',tread,'INT')
2860 73 : if(tread==1) multibinit_dtset%bound_option=intarr(1)
2861 : if(multibinit_dtset%bound_option<1.and.multibinit_dtset%bound_option>3)then
2862 : write(message, '(a,i8,a,a,a,a,a)' )&
2863 : & 'bound_option is',multibinit_dtset%bound_option,', but the only allowed values',ch10,&
2864 : & 'are 1,2,3 for multibinit.',ch10,&
2865 : & 'Action: correct bound_option in your input file.'
2866 : ABI_ERROR(message)
2867 : end if
2868 :
2869 :
2870 :
2871 73 : multibinit_dtset%bound_temp=325
2872 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'bound_temp',tread,'DPR')
2873 73 : if(tread==1) multibinit_dtset%bound_temp=dprarr(1)
2874 73 : if(multibinit_dtset%bound_temp<=0)then
2875 : write(message, '(a,f10.1,a,a,a,a,a)' )&
2876 0 : & 'Bound_Temp is ',multibinit_dtset%bound_temp,'. The only allowed values',ch10,&
2877 0 : & 'are positives values.',ch10,&
2878 0 : & 'Action: correct Bound_Temp in your input file.'
2879 0 : ABI_ERROR(message)
2880 : end if
2881 73 : multibinit_dtset%bound_step=1000
2882 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'bound_step',tread,'INT')
2883 73 : if(tread==1) multibinit_dtset%bound_step=intarr(1)
2884 : if(multibinit_dtset%bound_step<0.and.multibinit_dtset%bound_step>1)then
2885 : write(message, '(a,i8,a,a,a,a,a)' )&
2886 : & 'bound_step is',multibinit_dtset%bound_step,', but the only allowed values',ch10,&
2887 : & 'are 0 or 1 for multibinit.',ch10,&
2888 : & 'Action: correct bound_step in your input file.'
2889 : ABI_ERROR(message)
2890 : end if
2891 :
2892 73 : multibinit_dtset%fit_coeff=0
2893 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_coeff',tread,'INT')
2894 73 : if(tread==1) multibinit_dtset%fit_coeff=intarr(1)
2895 : if(multibinit_dtset%fit_coeff<0.and.multibinit_dtset%fit_coeff>1)then
2896 : write(message, '(a,i8,a,a,a,a,a)' )&
2897 : & 'fit_coeff is',multibinit_dtset%fit_coeff,', but the only allowed values',ch10,&
2898 : & 'are 0 or 1 for multibinit.',ch10,&
2899 : & 'Action: correct fit_coeff in your input file.'
2900 : ABI_ERROR(message)
2901 : end if
2902 :
2903 :
2904 73 : multibinit_dtset%fit_cutoff=0
2905 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_cutoff',tread,'DPR')
2906 73 : if(tread==1) multibinit_dtset%fit_cutoff=dprarr(1)
2907 73 : if(multibinit_dtset%fit_cutoff<0)then
2908 : write(message, '(a,i8,a,a,a,a,a)' )&
2909 0 : & 'fit_cutoff is',multibinit_dtset%fit_cutoff,', but the only allowed values',ch10,&
2910 0 : & 'are positives for multibinit.',ch10,&
2911 0 : & 'Action: correct fit_cutoff in your input file.'
2912 0 : ABI_ERROR(message)
2913 : end if
2914 :
2915 : !multibinit_dtset%fit_drop_rate=0.0_dp
2916 : !call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_drop_rate',tread,'DPR')
2917 : !if(tread==1) multibinit_dtset%fit_drop_rate=dprarr(1)
2918 : !if(multibinit_dtset%fit_drop_rate<0.or.multibinit_dtset%fit_drop_rate>1.0)then
2919 : ! write(message, '(a,i0,a,a,a,a)' )&
2920 : ! & 'fit_drop_rate is ',multibinit_dtset%fit_drop_rate,', which is lower',&
2921 : ! & ' than 0 of superior than 1.',&
2922 : ! & ch10,'Action: correct fit_drop_rate in your input file.'
2923 : ! ABI_ERROR(message)
2924 : !end if
2925 :
2926 :
2927 :
2928 :
2929 292 : multibinit_dtset%fit_grid(:)= 1
2930 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'fit_grid',tread,'INT')
2931 73 : if(tread==1) multibinit_dtset%fit_grid(1:3)=intarr(1:3)
2932 292 : do ii=1,3
2933 292 : if(multibinit_dtset%fit_grid(ii)<0.or.multibinit_dtset%fit_grid(ii)>20)then
2934 : write(message, '(a,i0,a,i0,a,a,a,i0,a)' )&
2935 0 : & 'fit_grid(',ii,') is ',multibinit_dtset%fit_grid(ii),', which is lower',&
2936 0 : & ' than 0 of superior than 20.',&
2937 0 : & ch10,'Action: correct fit_grid(',ii,') in your input file.'
2938 0 : ABI_ERROR(message)
2939 : end if
2940 : end do
2941 :
2942 : ! multibinit_dtset%fit_min_bound_coeff=0.0_dp
2943 : ! call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_min_bound_coeff',tread,'DPR')
2944 : ! if(tread==1) multibinit_dtset%fit_min_bound_coeff=dprarr(1)
2945 : !
2946 :
2947 :
2948 :
2949 219 : multibinit_dtset%fit_rangePower(:)= (/3,4/)
2950 73 : call intagm(dprarr,intarr,jdtset,marr,2,string(1:lenstr),'fit_rangePower',tread,'INT')
2951 109 : if(tread==1) multibinit_dtset%fit_rangePower(1:2)=intarr(1:2)
2952 219 : do ii=1,2
2953 219 : if(multibinit_dtset%fit_rangePower(ii)<0.or.multibinit_dtset%fit_rangePower(ii)>20)then
2954 : write(message, '(a,i0,a,i0,a,a,a,i0,a)' )&
2955 0 : & 'fit_rangePower(',ii,') is ',multibinit_dtset%fit_rangePower(ii),', which is lower',&
2956 0 : & ' than 0 of superior than 20.',&
2957 0 : & ch10,'Action: correct fit_rangePower(',ii,') in your input file.'
2958 0 : ABI_ERROR(message)
2959 : end if
2960 : end do
2961 :
2962 219 : multibinit_dtset%bound_rangePower(:)= (/6,6/)
2963 73 : call intagm(dprarr,intarr,jdtset,marr,2,string(1:lenstr),'bound_rangePower',tread,'INT')
2964 75 : if(tread==1) multibinit_dtset%bound_rangePower(1:2)=intarr(1:2)
2965 219 : do ii=1,2
2966 219 : if(multibinit_dtset%bound_rangePower(ii)<=0.or.multibinit_dtset%bound_rangePower(ii)>20)then
2967 : write(message, '(a,i0,a,i0,a,a,a,i0,a)' )&
2968 0 : & 'bound_rangePower(',ii,') is ',multibinit_dtset%bound_rangePower(ii),', which is lower',&
2969 0 : & ' than 0 of superior than 20.',&
2970 0 : & ch10,'Action: correct bound_rangePower(',ii,') in your input file.'
2971 0 : ABI_ERROR(message)
2972 : end if
2973 : end do
2974 :
2975 292 : multibinit_dtset%bound_cell(:)= (/6,6,6/)
2976 73 : call intagm(dprarr,intarr,jdtset,marr,3,string(1:lenstr),'bound_cell',tread,'INT')
2977 73 : if(tread==1) multibinit_dtset%bound_cell(1:3)=intarr(1:3)
2978 292 : do ii=1,3
2979 292 : if(multibinit_dtset%bound_cell(ii)<=0.or.multibinit_dtset%bound_cell(ii)>20)then
2980 : write(message, '(a,i0,a,i0,4a,i0,a)' )&
2981 0 : & 'bound_cell(',ii,') is ',multibinit_dtset%bound_cell(ii),', which is lower',&
2982 0 : & ' than 0 of superior than 20.',&
2983 0 : & ch10,'Action: correct bound_cell(',ii,') in your input file.'
2984 0 : ABI_ERROR(message)
2985 : end if
2986 : end do
2987 :
2988 73 : multibinit_dtset%fit_tolMSDE=0
2989 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_tolMSDE',tread,'DPR')
2990 73 : if(tread==1) multibinit_dtset%fit_tolMSDE=dprarr(1)
2991 73 : if(multibinit_dtset%fit_tolMSDE<0)then
2992 : write(message, '(a,i8,a,a,a,a,a)' )&
2993 0 : & 'fit_tolMSDE is',multibinit_dtset%fit_tolMSDE,', but the only allowed values',ch10,&
2994 0 : & 'are positives for multibinit.',ch10,&
2995 0 : & 'Action: correct fit_tolMSDE in your input file.'
2996 0 : ABI_ERROR(message)
2997 : end if
2998 :
2999 73 : multibinit_dtset%fit_tolMSDF=0
3000 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_tolMSDF',tread,'DPR')
3001 73 : if(tread==1) multibinit_dtset%fit_tolMSDF=dprarr(1)
3002 73 : if(multibinit_dtset%fit_tolMSDF<0)then
3003 : write(message, '(a,i8,a,a,a,a,a)' )&
3004 0 : & 'fit_tolMSDF is',multibinit_dtset%fit_tolMSDF,', but the only allowed values',ch10,&
3005 0 : & 'are positives for multibinit.',ch10,&
3006 0 : & 'Action: correct fit_tolMSDF in your input file.'
3007 0 : ABI_ERROR(message)
3008 : end if
3009 :
3010 73 : multibinit_dtset%fit_tolMSDS=0
3011 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_tolMSDS',tread,'DPR')
3012 73 : if(tread==1) multibinit_dtset%fit_tolMSDS=dprarr(1)
3013 73 : if(multibinit_dtset%fit_tolMSDS<0)then
3014 : write(message, '(a,i8,a,a,a,a,a)' )&
3015 0 : & 'fit_tolMSDS is',multibinit_dtset%fit_tolMSDS,', but the only allowed values',ch10,&
3016 0 : & 'are positives for multibinit.',ch10,&
3017 0 : & 'Action: correct fit_tolMSDS in your input file.'
3018 0 : ABI_ERROR(message)
3019 : end if
3020 :
3021 73 : multibinit_dtset%fit_tolMSDFS=0
3022 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_tolMSDFS',tread,'DPR')
3023 73 : if(tread==1) multibinit_dtset%fit_tolMSDFS=dprarr(1)
3024 73 : if(multibinit_dtset%fit_tolMSDFS<0)then
3025 : write(message, '(a,i8,a,a,a,a,a)' )&
3026 0 : & 'fit_tolMSDFS is',multibinit_dtset%fit_tolMSDFS,', but the only allowed values',ch10,&
3027 0 : & 'are positives for multibinit.',ch10,&
3028 0 : & 'Action: correct fit_tolMSDFS in your input file.'
3029 0 : ABI_ERROR(message)
3030 : end if
3031 :
3032 73 : multibinit_dtset%fit_tolGF=0
3033 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_tolGF',tread,'DPR')
3034 73 : if(tread==1) multibinit_dtset%fit_tolGF=dprarr(1)
3035 73 : if(multibinit_dtset%fit_tolGF<0)then
3036 : write(message, '(a,i8,a,a,a,a,a)' )&
3037 0 : & 'fit_tolGF is',multibinit_dtset%fit_tolGF,', but the only allowed values',ch10,&
3038 0 : & 'are positives for multibinit.',ch10,&
3039 0 : & 'Action: correct fit_tolGF in your input file.'
3040 0 : ABI_ERROR(message)
3041 : end if
3042 :
3043 73 : multibinit_dtset%fit_weight_T=-0.1_dp
3044 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'fit_weight_T',tread,'DPR')
3045 73 : if(tread==1) multibinit_dtset%fit_weight_T=dprarr(1)
3046 :
3047 :
3048 : !G
3049 :
3050 : !H
3051 :
3052 : !I
3053 :
3054 :
3055 : !J
3056 :
3057 : !K
3058 :
3059 : !L
3060 : ! multibinit_dtset%latt_mask(:)=0
3061 : ! call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'latt_mask',tread,'INT')
3062 : ! if(tread==1) then
3063 : ! do ii=1, 3
3064 : ! multibinit_dtset%latt_mask(ii)=intarr(ii)
3065 : ! if(multibinit_dtset%latt_mask(ii) <0 .or. multibinit_dtset%latt_mask(ii) >1)then
3066 : ! write(message, '(a)' )&
3067 : ! & ' latt_mask element should be 0 or 1.'
3068 : ! ABI_ERROR(message)
3069 : ! end if
3070 : ! end do
3071 : ! end if
3072 :
3073 :
3074 : !M
3075 : !N
3076 :
3077 219 : ABI_MALLOC(multibinit_dtset%fit_bancoeff,(multibinit_dtset%fit_nbancoeff))
3078 73 : if (multibinit_dtset%fit_nbancoeff >0)then
3079 0 : if(multibinit_dtset%fit_nbancoeff>marr)then
3080 0 : marr=multibinit_dtset%fit_nbancoeff
3081 0 : ABI_FREE(intarr)
3082 0 : ABI_MALLOC(intarr,(marr))
3083 : end if
3084 0 : multibinit_dtset%fit_bancoeff(:)=0
3085 : call intagm(dprarr,intarr,jdtset,marr,multibinit_dtset%fit_nbancoeff,&
3086 0 : & string(1:lenstr),'fit_bancoeff',tread,'INT')
3087 0 : if(tread==1)then
3088 0 : do ii=1,multibinit_dtset%fit_nbancoeff
3089 0 : multibinit_dtset%fit_bancoeff(ii)=intarr(ii)
3090 : end do
3091 : end if
3092 : end if
3093 :
3094 219 : ABI_MALLOC(multibinit_dtset%fit_fixcoeff,(multibinit_dtset%fit_nfixcoeff))
3095 73 : if (multibinit_dtset%fit_nfixcoeff >0)then
3096 7 : if(multibinit_dtset%fit_nfixcoeff>marr)then
3097 0 : marr=multibinit_dtset%fit_nfixcoeff
3098 0 : ABI_FREE(intarr)
3099 0 : ABI_MALLOC(intarr,(marr))
3100 : end if
3101 21 : multibinit_dtset%fit_fixcoeff(:)=0
3102 : call intagm(dprarr,intarr,jdtset,marr,multibinit_dtset%fit_nfixcoeff,&
3103 7 : & string(1:lenstr),'fit_fixcoeff',tread,'INT')
3104 7 : if(tread==1)then
3105 21 : do ii=1,multibinit_dtset%fit_nfixcoeff
3106 21 : multibinit_dtset%fit_fixcoeff(ii)=intarr(ii)
3107 : end do
3108 : end if
3109 : end if
3110 :
3111 219 : ABI_MALLOC(multibinit_dtset%fit_imposecoeff,(multibinit_dtset%fit_nimposecoeff))
3112 73 : if (multibinit_dtset%fit_nimposecoeff >0)then
3113 0 : if(multibinit_dtset%fit_nimposecoeff>marr)then
3114 0 : marr=multibinit_dtset%fit_nimposecoeff
3115 0 : ABI_FREE(intarr)
3116 0 : ABI_MALLOC(intarr,(marr))
3117 : end if
3118 0 : multibinit_dtset%fit_imposecoeff(:)=0
3119 : call intagm(dprarr,intarr,jdtset,marr,multibinit_dtset%fit_nimposecoeff,&
3120 0 : & string(1:lenstr),'fit_imposecoeff',tread,'INT')
3121 0 : if(tread==1)then
3122 0 : do ii=1,multibinit_dtset%fit_nimposecoeff
3123 0 : multibinit_dtset%fit_imposecoeff(ii)=intarr(ii)
3124 : end do
3125 : end if
3126 : end if
3127 :
3128 : !O
3129 : block
3130 : integer :: minv, maxv, n
3131 73 : minv =multibinit_dtset%fit_rangepower(1)
3132 73 : maxv =multibinit_dtset%fit_rangepower(2)
3133 73 : n=maxv-minv+1
3134 : call read_int_array_var(int_array_var=multibinit_dtset%fit_max_nbody_inrange, &
3135 : & size=n, &
3136 73 : & var_name='fit_max_nbody', type='INT', default=999)
3137 219 : ABI_MALLOC(multibinit_dtset%fit_max_nbody, (maxv))
3138 348 : multibinit_dtset%fit_max_nbody(:)=999
3139 202 : multibinit_dtset%fit_max_nbody(minv:maxv) = multibinit_dtset%fit_max_nbody_inrange
3140 :
3141 :
3142 : end block
3143 :
3144 :
3145 219 : ABI_MALLOC(multibinit_dtset%opt_coeff,(multibinit_dtset%opt_ncoeff))
3146 73 : if (multibinit_dtset%opt_ncoeff >0)then
3147 1 : if(multibinit_dtset%opt_ncoeff>marr)then
3148 0 : marr=multibinit_dtset%opt_ncoeff
3149 0 : ABI_FREE(intarr)
3150 0 : ABI_MALLOC(intarr,(marr))
3151 : end if
3152 3 : multibinit_dtset%opt_coeff(:)=0
3153 : call intagm(dprarr,intarr,jdtset,marr,multibinit_dtset%opt_ncoeff,&
3154 1 : & string(1:lenstr),'opt_coeff',tread,'INT')
3155 1 : if(tread==1)then
3156 3 : do ii=1,multibinit_dtset%opt_ncoeff
3157 3 : multibinit_dtset%opt_coeff(ii)=intarr(ii)
3158 : end do
3159 : end if
3160 : end if
3161 :
3162 : !P
3163 :
3164 : !Q
3165 219 : ABI_MALLOC(multibinit_dtset%qmass,(multibinit_dtset%nnos))
3166 73 : multibinit_dtset%qmass(:)= zero
3167 73 : if(multibinit_dtset%nnos>=1)then
3168 0 : if(multibinit_dtset%nnos>marr)then
3169 0 : marr=multibinit_dtset%nnos
3170 0 : ABI_FREE(intarr)
3171 0 : ABI_FREE(dprarr)
3172 0 : ABI_MALLOC(intarr,(marr))
3173 0 : ABI_MALLOC(dprarr,(marr))
3174 : end if
3175 0 : call intagm(dprarr,intarr,jdtset,marr,multibinit_dtset%nnos,string(1:lenstr),'qmass',tread,'DPR')
3176 0 : if(tread==1) multibinit_dtset%qmass(:)=dprarr(1:multibinit_dtset%nnos)
3177 : end if
3178 :
3179 73 : if (multibinit_dtset%nqshft/=0)then
3180 73 : if(3*multibinit_dtset%nqshft>marr)then
3181 0 : marr=3*multibinit_dtset%nqshft
3182 0 : ABI_FREE(intarr)
3183 0 : ABI_FREE(dprarr)
3184 0 : ABI_MALLOC(intarr,(marr))
3185 0 : ABI_MALLOC(dprarr,(marr))
3186 : end if
3187 219 : ABI_MALLOC(multibinit_dtset%q1shft,(3,multibinit_dtset%nqshft))
3188 365 : multibinit_dtset%q1shft(:,:)=zero
3189 73 : call intagm(dprarr,intarr,jdtset,marr,3*multibinit_dtset%nqshft, string(1:lenstr),'q1shft',tread,'DPR')
3190 73 : if(tread==1) multibinit_dtset%q1shft(1:3,1:multibinit_dtset%nqshft)=&
3191 0 : & reshape(dprarr(1:3*multibinit_dtset%nqshft),(/3,multibinit_dtset%nqshft/))
3192 : end if
3193 :
3194 219 : ABI_MALLOC(multibinit_dtset%qph1l,(3,multibinit_dtset%nph1l))
3195 219 : ABI_MALLOC(multibinit_dtset%qnrml1,(multibinit_dtset%nph1l))
3196 73 : if (multibinit_dtset%nph1l/=0)then
3197 73 : if(4*multibinit_dtset%nph1l>marr)then
3198 8 : marr=4*multibinit_dtset%nph1l
3199 8 : ABI_FREE(intarr)
3200 8 : ABI_FREE(dprarr)
3201 24 : ABI_MALLOC(intarr,(marr))
3202 24 : ABI_MALLOC(dprarr,(marr))
3203 : end if
3204 589 : multibinit_dtset%qph1l(:,:)=zero
3205 202 : multibinit_dtset%qnrml1(:)=zero
3206 73 : call intagm(dprarr,intarr,jdtset,marr,4*multibinit_dtset%nph1l,string(1:lenstr),'qph1l',tread,'DPR')
3207 73 : if(tread==1)then
3208 72 : do iph1=1,multibinit_dtset%nph1l
3209 256 : do ii=1,3
3210 256 : multibinit_dtset%qph1l(ii,iph1)=dprarr(ii+(iph1-1)*4)
3211 : end do
3212 64 : multibinit_dtset%qnrml1(iph1)=dprarr(4+(iph1-1)*4)
3213 72 : if(abs(multibinit_dtset%qnrml1(iph1))<DDB_QTOL)then
3214 : write(message, '(a,a,a,a,a)' )&
3215 0 : & 'The first list of wavevectors ','should not have non-analytical data.',ch10,&
3216 0 : & 'Action: correct the first list',' of wavevectors in the input file.'
3217 0 : ABI_ERROR(message)
3218 : end if
3219 : end do
3220 : end if
3221 : end if
3222 :
3223 219 : ABI_MALLOC(multibinit_dtset%qph2l,(3,multibinit_dtset%nph2l))
3224 219 : ABI_MALLOC(multibinit_dtset%qnrml2,(multibinit_dtset%nph2l))
3225 73 : if (multibinit_dtset%nph2l/=0)then
3226 0 : if(4*multibinit_dtset%nph2l>marr)then
3227 0 : marr=4*multibinit_dtset%nph2l
3228 0 : ABI_FREE(intarr)
3229 0 : ABI_FREE(dprarr)
3230 0 : ABI_MALLOC(intarr,(marr))
3231 0 : ABI_MALLOC(dprarr,(marr))
3232 : end if
3233 0 : multibinit_dtset%qph2l(:,:)=zero
3234 0 : multibinit_dtset%qnrml2(:)=zero
3235 0 : call intagm(dprarr,intarr,jdtset,marr,4*multibinit_dtset%nph2l,string(1:lenstr),'qph2l',tread,'DPR')
3236 0 : if(tread==1)then
3237 0 : do iph2=1,multibinit_dtset%nph2l
3238 0 : do ii=1,3
3239 0 : multibinit_dtset%qph2l(ii,iph2)=dprarr(ii+(iph2-1)*4)
3240 : end do
3241 0 : multibinit_dtset%qnrml2(iph2)=dprarr(4+(iph2-1)*4)
3242 0 : if(abs(multibinit_dtset%qnrml2(iph2))>DDB_QTOL)then
3243 : write(message, '(a,a,a,a,a)' )&
3244 0 : & 'The second list of wavevectors',' should have only non-analytical data.',ch10,&
3245 0 : & 'Action: correct the second list','of wavevectors in the input file.'
3246 0 : ABI_ERROR(message)
3247 : end if
3248 : end do
3249 : end if
3250 : end if
3251 :
3252 : !R
3253 73 : if(9>marr)then
3254 0 : marr=9
3255 0 : ABI_FREE(intarr)
3256 0 : ABI_FREE(dprarr)
3257 0 : ABI_MALLOC(intarr,(marr))
3258 0 : ABI_MALLOC(dprarr,(marr))
3259 : end if
3260 949 : multibinit_dtset%rprim(:,:)= zero
3261 73 : call intagm(dprarr,intarr,jdtset,marr,9,string(1:lenstr),'rprim',tread,'DPR')
3262 73 : if(tread==1) then
3263 0 : multibinit_dtset%rprim(1:3,1:3)= reshape(dprarr(1:9),(/3,3/))
3264 : ! check new rprimd
3265 : if(all(abs(multibinit_dtset%rprim(1,:))<tol16).or.&
3266 0 : & all(abs(multibinit_dtset%rprim(2,:))<tol16).or.all(abs(multibinit_dtset%rprim(3,:))<tol16)) then
3267 : write(message, '(3a)' )&
3268 0 : & ' There is a problem with rprim',ch10,&
3269 0 : & 'Action: correct rprim'
3270 0 : ABI_BUG(message)
3271 : end if
3272 : end if
3273 : !S
3274 :
3275 73 : if(6>marr)then
3276 0 : marr=6
3277 0 : ABI_FREE(intarr)
3278 0 : ABI_FREE(dprarr)
3279 0 : ABI_MALLOC(intarr,(marr))
3280 0 : ABI_MALLOC(dprarr,(marr))
3281 : end if
3282 511 : multibinit_dtset%strten_reference(:)= zero
3283 73 : call intagm(dprarr,intarr,jdtset,marr,6,string(1:lenstr),'strten_reference',tread,'DPR')
3284 73 : if(tread==1) multibinit_dtset%strten_reference(1:6)= dprarr(1:6)
3285 :
3286 : !T
3287 :
3288 : !U
3289 :
3290 : !V
3291 :
3292 : !W
3293 :
3294 : !X
3295 :
3296 : !Y
3297 :
3298 : !Z
3299 :
3300 : !=======================================================================
3301 : ! Read Geometric constraint variables for relaxation
3302 : !=======================================================================
3303 :
3304 73 : natom_sc = natom*multibinit_dtset%ncell(1)*multibinit_dtset%ncell(2)*multibinit_dtset%ncell(3)
3305 219 : ABI_MALLOC(multibinit_dtset%iatfix,(3,natom_sc))
3306 20865 : multibinit_dtset%iatfix(:,:)=0
3307 :
3308 365 : do idir=0,3
3309 :
3310 292 : if(idir==0)then
3311 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'natfix',tread,'INT')
3312 219 : else if(idir==1)then
3313 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'natfixx',tread,'INT')
3314 146 : else if(idir==2)then
3315 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'natfixy',tread,'INT')
3316 73 : else if(idir==3)then
3317 73 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),'natfixz',tread,'INT')
3318 : end if
3319 :
3320 : ! Use natfix also for natfixx,natfixy,natfixz
3321 292 : natfix=0
3322 292 : if(tread==1)then
3323 4 : natfix=intarr(1)
3324 4 : ABI_FREE(intarr)
3325 4 : ABI_FREE(dprarr)
3326 12 : ABI_MALLOC(intarr,(natfix))
3327 12 : ABI_MALLOC(dprarr,(natfix))
3328 4 : marr = natfix
3329 : endif
3330 :
3331 : ! Checks the validity of natfix
3332 292 : if (natfix<0 .or. natfix>natom_sc) then
3333 : write(message, '(a,a,a,i0,a,i4,a,a,a)' )&
3334 0 : & 'The input variables natfix, natfixx, natfixy and natfixz must be',ch10,&
3335 0 : & 'between 0 and natom of the supercell (= ',natom_sc,'), while one of them is ',natfix,'.',ch10,&
3336 0 : & 'Action: correct that occurence in your input file.'
3337 0 : ABI_ERROR(message)
3338 : end if
3339 :
3340 : ! Read iatfix
3341 292 : if(idir==0)then
3342 73 : call intagm(dprarr,intarr,jdtset,marr,natfix,string(1:lenstr),'iatfix',tread,'INT')
3343 219 : else if(idir==1)then
3344 73 : call intagm(dprarr,intarr,jdtset,marr,natfix,string(1:lenstr),'iatfixx',tread,'INT')
3345 146 : else if(idir==2)then
3346 73 : call intagm(dprarr,intarr,jdtset,marr,natfix,string(1:lenstr),'iatfixy',tread,'INT')
3347 73 : else if(idir==3)then
3348 73 : call intagm(dprarr,intarr,jdtset,marr,natfix,string(1:lenstr),'iatfixz',tread,'INT')
3349 : end if
3350 :
3351 : ! If some iatfix was read, natfix must vanish
3352 292 : if (natfix==0 .and. tread==1)then
3353 : write(message, '(a,i1,5a)' )&
3354 0 : & 'For direction ',idir,' the corresponding natfix is zero,',ch10,&
3355 0 : & 'while iatfix specifies some atoms to be fixed.',ch10,&
3356 0 : & 'Action: either specify a non-zero natfix(x,y,z) or suppress iatfix(x,y,z).'
3357 0 : ABI_ERROR(message)
3358 : end if
3359 :
3360 : ! If natfix is non-zero, iatfix must be defined
3361 292 : if (natfix>0 .and. tread==0)then
3362 : write(message, '(a,i1,3a,i0,3a)' )&
3363 0 : & 'For direction ',idir,' no iatfix has been specified,',ch10,&
3364 0 : & 'while natfix specifies that some atoms to be fixed, natfix= ',natfix,'.',ch10,&
3365 0 : & 'Action: either set natfix(x,y,z) to zero or define iatfix(x,y,z).'
3366 0 : ABI_ERROR(message)
3367 : end if
3368 :
3369 365 : if(tread==1)then
3370 11 : do ii=1,natfix
3371 : ! Checks the validity of the input iatfix
3372 7 : if (intarr(ii)<1 .or. intarr(ii)>natom_sc) then
3373 : write(message, '(a,a,a,i0,a,a,a)' )&
3374 0 : & 'The input variables iatfix, iatfixx, iatfixy and iatfixz must be',ch10,&
3375 0 : & 'between 1 and natom of the supercell, while one of them is ',intarr(ii),'.',ch10,&
3376 0 : & 'Action: correct that occurence in your input file.'
3377 0 : ABI_ERROR(message)
3378 : end if
3379 : ! Finally set the value of the internal iatfix array
3380 151 : do iatom=1,natom_sc
3381 147 : if(intarr(ii)==iatom)then
3382 10 : if(idir==0)multibinit_dtset%iatfix(1:3,iatom)=1
3383 7 : if(idir/=0)multibinit_dtset%iatfix(idir,iatom)=1
3384 : end if
3385 : end do
3386 : end do
3387 : end if
3388 : end do
3389 :
3390 : !=======================================================================
3391 : ! Read SCALE UP variables
3392 : !=======================================================================
3393 :
3394 73 : call invars10scup(multibinit_dtset%scup_dtset,lenstr,string)
3395 :
3396 : !=======================================================================
3397 : !Finished reading in variables - deallocate
3398 : !=======================================================================
3399 :
3400 73 : ABI_FREE(dprarr)
3401 73 : ABI_FREE(intarr)
3402 :
3403 : !=======================================================================
3404 : !Check consistency of input variables:
3405 : !=======================================================================
3406 :
3407 73 : if(multibinit_dtset%prtsrlr/=0 .and. multibinit_dtset%ifcflag/=1) then
3408 : write(message, '(3a)' )&
3409 0 : & 'ifcflag must be 1 for the SR/LR decomposition of the phonon frequencies',ch10,&
3410 0 : & 'Action: correct ifcflag in your input file.'
3411 0 : ABI_ERROR(message)
3412 : end if
3413 :
3414 : !FIXME: add check that if freeze_displ /= 0 then you need to be doing ifc and phonon interpolation
3415 :
3416 730 : if (multibinit_dtset%ifcflag > 0 .and. sum(abs(multibinit_dtset%ngqpt)) == 0) then
3417 : write(message, '(3a)' )&
3418 0 : & 'if you want interatomic force constant output, multibinit needs ngqpt input variable ',ch10,&
3419 0 : & 'Action: set ngqpt in your input file.'
3420 0 : ABI_ERROR(message)
3421 : end if
3422 :
3423 : !check that q-grid refinement is a divisor of ngqpt in each direction
3424 584 : if(any(multibinit_dtset%qrefine(:) > 1) .and. &
3425 : & any(abs(dmod(dble(multibinit_dtset%ngqpt(1:3))/dble(multibinit_dtset%qrefine(:)),one)) > tol10)) then
3426 : write(message, '(a,3i0,a,a,a,3i8,a,a)' )&
3427 0 : & 'qrefine is',multibinit_dtset%qrefine,' The only allowed values',ch10,&
3428 0 : & 'are integers which are divisors of the ngqpt grid', multibinit_dtset%ngqpt,ch10,&
3429 0 : & 'Action: correct qrefine in your input file.'
3430 0 : ABI_ERROR(message)
3431 : end if
3432 :
3433 : ! check new rprimd
3434 1022 : if(all(multibinit_dtset%acell(:) > one).and.all(abs(multibinit_dtset%rprim(:,:))<tol16))then
3435 : write(message, '(3a)' )&
3436 0 : & ' acell is defined but there is no rprim',ch10,&
3437 0 : & 'Action: add rprim input'
3438 0 : ABI_BUG(message)
3439 : end if
3440 :
3441 :
3442 : !check the fit_bancoeff and fit_fixcoeff and opt_coeff, fit_imposecoeff
3443 73 : do ii=1,multibinit_dtset%fit_nbancoeff
3444 73 : do jj=ii+1,multibinit_dtset%fit_nbancoeff
3445 0 : if (multibinit_dtset%fit_bancoeff(ii) == multibinit_dtset%fit_bancoeff(jj))then
3446 : write(message, '(a,I0,a,I0,2a)' )&
3447 0 : & ' There is two similar numbers for fit_bancoeff: ',multibinit_dtset%fit_bancoeff(ii),&
3448 0 : & ' and ', multibinit_dtset%fit_bancoeff(jj),ch10,&
3449 0 : & 'Action: change fit_bancoeff'
3450 0 : ABI_BUG(message)
3451 : end if
3452 : end do
3453 : end do
3454 :
3455 87 : do ii=1,multibinit_dtset%fit_nfixcoeff
3456 94 : do jj=ii+1,multibinit_dtset%fit_nfixcoeff
3457 21 : if (multibinit_dtset%fit_fixcoeff(ii) == multibinit_dtset%fit_fixcoeff(jj))then
3458 : write(message, '(a,I0,a,I0,2a)' )&
3459 0 : & ' There is two similar numbers for fit_fixcoeff: ',multibinit_dtset%fit_fixcoeff(ii),&
3460 0 : & ' and ', multibinit_dtset%fit_fixcoeff(jj),ch10,&
3461 0 : & 'Action: change fit_fixcoeff'
3462 0 : ABI_BUG(message)
3463 : end if
3464 : end do
3465 : end do
3466 :
3467 :
3468 75 : do ii=1,multibinit_dtset%opt_ncoeff
3469 76 : do jj=ii+1,multibinit_dtset%opt_ncoeff
3470 3 : if (multibinit_dtset%opt_coeff(ii) == multibinit_dtset%opt_coeff(jj))then
3471 : write(message, '(a,I0,a,I0,2a)' )&
3472 0 : & ' There is two similar numbers for opt_coeff: ',multibinit_dtset%opt_coeff(ii),&
3473 0 : & ' and ', multibinit_dtset%opt_coeff(jj),ch10,&
3474 0 : & 'Action: change opt_coeff'
3475 0 : ABI_BUG(message)
3476 : end if
3477 : end do
3478 : end do
3479 :
3480 73 : do ii=1,multibinit_dtset%fit_nimposecoeff
3481 73 : do jj=ii+1,multibinit_dtset%fit_nimposecoeff
3482 0 : if (multibinit_dtset%fit_imposecoeff(ii) == multibinit_dtset%fit_imposecoeff(jj))then
3483 : write(message, '(a,I0,a,I0,2a)' )&
3484 0 : & ' There is two similar numbers for fit_imposecoeff: ',multibinit_dtset%fit_imposecoeff(ii),&
3485 0 : & ' and ', multibinit_dtset%fit_imposecoeff(jj),ch10,&
3486 0 : & 'Action: change fit_imposecoeff'
3487 0 : ABI_BUG(message)
3488 : end if
3489 : end do
3490 : end do
3491 :
3492 292 : multibinit_dtset%dipdip_range=0
3493 292 : do ii=1,3
3494 292 : if(multibinit_dtset%dipdip_range(ii) < multibinit_dtset%ncell(ii)) then
3495 219 : write(message,'(4a,3I3,3a,3I3,6a)') ch10,&
3496 219 : & ' --- !WARNING',ch10,&
3497 219 : & ' The range of dipdip_range (',multibinit_dtset%dipdip_range(:),')',ch10,&
3498 219 : & ' But the range of the cell for the simulation is',&
3499 219 : & multibinit_dtset%ncell(:),')',ch10,&
3500 219 : & ' dipdip_range is set to ncell.',ch10,&
3501 438 : & ' ---',ch10
3502 219 : multibinit_dtset%dipdip_range(ii) = multibinit_dtset%ncell(ii)
3503 219 : call wrtout(std_out,message,'COLL')
3504 : end if
3505 : ! if(multibinit_dtset%dipdip_range(ii) < multibinit_dtset%bound_cell(ii)) then
3506 : ! write(message,'(4a,3I3,3a,3I3,6a)') ch10,&
3507 : !& ' --- !WARNING',ch10,&
3508 : !& ' The range of dipdip_range (',multibinit_dtset%dipdip_range(:),')',ch10,&
3509 : !& ' But the range of the cell for the simulation is',&
3510 : !& multibinit_dtset%ncell(:),')',ch10,&
3511 : !& ' dipdip_range is set to bound_cell.',ch10,&
3512 : !& ' ---',ch10
3513 : ! multibinit_dtset%dipdip_range(:) = multibinit_dtset%ncell(:)
3514 : ! call wrtout(std_out,message,'COLL')
3515 : ! exit
3516 : ! end if
3517 : end do
3518 :
3519 : !Check if only one tolerance is specify
3520 73 : if(abs(multibinit_dtset%fit_tolMSDF) >zero .and. abs(multibinit_dtset%fit_tolMSDS) >zero) then
3521 : write(message, '(3a)' ) &
3522 0 : & ' There is two tolerance flags for the fit: fit_tolMSDF and fit_tolMSDS',ch10,&
3523 0 : & 'Action: Put only one tolerance flag'
3524 0 : ABI_BUG(message)
3525 : end if
3526 73 : if(abs(multibinit_dtset%fit_tolMSDF) >zero .and. abs(multibinit_dtset%fit_tolMSDE) >zero)then
3527 : write(message, '(3a)' ) &
3528 0 : & ' There is two tolerance flags for the fit: fit_tolMSDF and fit_tolMSDE',ch10,&
3529 0 : & 'Action: Put only one tolerance flag'
3530 0 : ABI_BUG(message)
3531 : end if
3532 73 : if(abs(multibinit_dtset%fit_tolMSDF) >zero .and. abs(multibinit_dtset%fit_tolMSDFS) >zero)then
3533 : write(message, '(3a)' ) &
3534 0 : & ' There is two tolerance flags for the fit: fit_tolMSDF and fit_tolMSDFS',ch10,&
3535 0 : & 'Action: Put only one tolerance flag'
3536 0 : ABI_BUG(message)
3537 : end if
3538 73 : if(abs(multibinit_dtset%fit_tolMSDS) >zero .and. abs(multibinit_dtset%fit_tolMSDE) >zero)then
3539 : write(message, '(3a)' ) &
3540 0 : & ' There is two tolerance flags for the fit: fit_tolMSDS and fit_tolMSDE',ch10,&
3541 0 : & 'Action: Put only one tolerance flag'
3542 0 : ABI_BUG(message)
3543 : end if
3544 73 : if(abs(multibinit_dtset%fit_tolMSDS) >zero .and. abs(multibinit_dtset%fit_tolMSDFS) >zero)then
3545 : write(message, '(3a)' ) &
3546 0 : & ' There is two tolerance flags for the fit: fit_tolMSDS and fit_tolMSDFS',ch10,&
3547 0 : & 'Action: Put only one tolerance flag'
3548 0 : ABI_BUG(message)
3549 : end if
3550 73 : if(abs(multibinit_dtset%fit_tolMSDE) >zero .and. abs(multibinit_dtset%fit_tolMSDFS) >zero)then
3551 : write(message, '(3a)' ) &
3552 0 : & ' There is two tolerance flags for the fit: fit_tolMSDE and fit_tolMSDFS',ch10,&
3553 0 : & 'Action: Put only one tolerance flag'
3554 0 : ABI_BUG(message)
3555 : end if
3556 :
3557 : contains
3558 1 : subroutine read_int_var(int_var, var_name)
3559 : ! Read an integer variable from the input
3560 : integer, intent(inout) :: int_var
3561 : character(*), intent(in) :: var_name
3562 1 : call intagm(dprarr,intarr,jdtset,marr,1,string(1:lenstr),var_name,tread,'INT')
3563 1 : if(tread==1) int_var=intarr(1)
3564 1 : end subroutine read_int_var
3565 :
3566 73 : subroutine read_int_array_var(int_array_var, size, var_name, type, default)
3567 : ! Read an integer array variable from the input
3568 : integer, allocatable, intent(inout) :: int_array_var(:)
3569 : integer, intent(in) :: size, default
3570 : character(*), intent(in) :: var_name, type
3571 219 : ABI_MALLOC(int_array_var,(size))
3572 73 : if (size >0)then
3573 73 : if(size>marr)then
3574 0 : marr=size
3575 0 : ABI_FREE(intarr)
3576 0 : ABI_MALLOC(intarr,(marr))
3577 : end if
3578 202 : int_array_var(:)=default
3579 : call intagm(dprarr,intarr,jdtset,marr,size,&
3580 73 : & string(1:lenstr),var_name,tread, type)
3581 73 : if(tread==1)then
3582 2 : do ii=1,size
3583 2 : int_array_var(ii)=intarr(ii)
3584 : end do
3585 : end if
3586 : end if
3587 73 : end subroutine read_int_array_var
3588 :
3589 219 : subroutine read_dpr_array_var(dpr_array_var, size, var_name, type, default)
3590 : ! Read a double precision real array variable from the input
3591 : real(dp), allocatable, intent(inout) :: dpr_array_var(:)
3592 : integer, intent(in) :: size
3593 : real(dp), intent(in) :: default
3594 : character(*), intent(in) :: var_name, type
3595 657 : ABI_MALLOC(dpr_array_var,(size))
3596 219 : if (size >0)then
3597 18 : if(size>marr)then
3598 0 : marr=size
3599 0 : ABI_FREE(intarr)
3600 0 : ABI_FREE(dprarr)
3601 0 : ABI_MALLOC(intarr,(marr))
3602 0 : ABI_MALLOC(dprarr,(marr))
3603 : end if
3604 67 : dpr_array_var(:)=default
3605 : call intagm(dprarr,intarr,jdtset,marr,size,&
3606 18 : & string(1:lenstr),var_name,tread, type)
3607 18 : if(tread==1)then
3608 51 : do ii=1,size
3609 51 : dpr_array_var(ii)=dprarr(ii)
3610 : end do
3611 : end if
3612 : end if
3613 219 : end subroutine read_dpr_array_var
3614 :
3615 : end subroutine invars10
3616 : !!***
3617 :
3618 :
3619 :
3620 : !----------------------------------------------------------------------
3621 :
3622 : !!****f* m_multibinit_dataset/outvars_multibinit
3623 : !!
3624 : !! NAME
3625 : !! outvars_multibinit
3626 : !!
3627 : !! FUNCTION
3628 : !! Takes as an input the input dtset for multibinit and echoes it to
3629 : !! the output
3630 : !!
3631 : !! INPUTS
3632 : !! multibinit_dtset <type(multibinit_dtset_type)> datatype with all the input variables
3633 : !! nunit=unit number for input or output
3634 : !!
3635 : !! OUTPUT
3636 : !! (only writing)
3637 : !!
3638 : !! NOTES
3639 : !! Should be executed by one processor only.
3640 : !!
3641 : !! SOURCE
3642 :
3643 98 : subroutine outvars_multibinit (multibinit_dtset,nunit)
3644 :
3645 : !Arguments -------------------------------
3646 : !scalars
3647 : integer,intent(in) :: nunit
3648 : type(multibinit_dtset_type),intent(in) :: multibinit_dtset
3649 :
3650 : !Local variables -------------------------
3651 : !Set routine version number here:
3652 : !scalars
3653 : integer :: ii,iph1,iph2,iqshft,natfix
3654 : integer :: n
3655 :
3656 : !*********************************************************************
3657 :
3658 : !Write the heading
3659 7938 : write(nunit,'(a,80a,a)') ch10,('=',ii=1,80),ch10
3660 : write(nunit, '(a,a)' )&
3661 98 : & ' -outvars_multibinit: echo values of input variables ----------------------',ch10
3662 :
3663 : !The flags
3664 98 : if(multibinit_dtset%ifcflag/=0)then
3665 98 : write(nunit,'(a)')' Flags : '
3666 98 : if(multibinit_dtset%ifcflag/=0)write(nunit,'(3x,a9,3i10)')' ifcflag',multibinit_dtset%ifcflag
3667 98 : if(multibinit_dtset%prt_model/=0)write(nunit,'(3x,a9,3i10)')'prt_model',multibinit_dtset%prt_model
3668 98 : if(multibinit_dtset%prt_phfrq/=0)write(nunit,'(3x,a9,3i10)')'prt_phfrq',multibinit_dtset%prt_phfrq
3669 98 : if(multibinit_dtset%strcpling/=0)write(nunit,'(3x,a9,3i10)')' strcpling',multibinit_dtset%strcpling
3670 98 : if(multibinit_dtset%strcpling==2)write(nunit,'(3x,a9,3es8.2)')'delta_df',multibinit_dtset%delta_df
3671 : end if
3672 :
3673 98 : if(multibinit_dtset%dynamics/=0)then
3674 66 : write(nunit,'(a)')' Molecular Dynamics :'
3675 66 : write(nunit,'(3x,a9,3I10.1)')' dynamics',multibinit_dtset%dynamics
3676 66 : write(nunit,'(3x,a9,3I10.1)')' dyn_chksym',multibinit_dtset%dyn_chksym
3677 66 : write(nunit,'(3x,a9,3F10.1)')' temp',multibinit_dtset%temperature
3678 66 : write(nunit,'(3x,a9,3I10.1)')' ntime',multibinit_dtset%ntime
3679 66 : if (multibinit_dtset%nctime /=1)then
3680 10 : write(nunit,'(3x,a9,3I10.1)')' nctime',multibinit_dtset%nctime
3681 : end if
3682 66 : write(nunit,'(3x,a9,3i10)') ' ncell',multibinit_dtset%ncell
3683 66 : write(nunit,'(3x,a9,3i10)') ' dtion',multibinit_dtset%dtion
3684 66 : if (multibinit_dtset%restartxf/=0) then
3685 18 : write(nunit,'(3x,a9,3i10)') 'restartxf',multibinit_dtset%restartxf
3686 : end if
3687 66 : if(multibinit_dtset%dynamics==13)then
3688 32 : write(nunit,'(3x,a9,3i10)')' optcell',multibinit_dtset%optcell
3689 32 : write(nunit,'(3x,a9,3F12.1)')' bmass',multibinit_dtset%bmass
3690 32 : write(nunit,'(3x,a9,3I10)')' nnos',multibinit_dtset%nnos
3691 32 : write(nunit,'(3x,a12)',advance='no')' qmass '
3692 32 : write(nunit,'(3x,15F12.10)') (multibinit_dtset%qmass(ii),ii=1,multibinit_dtset%nnos)
3693 : end if
3694 :
3695 : ! TODO: print the value of the efield here if it is present.
3696 :
3697 30994 : if(any(multibinit_dtset%iatfix /= 0))then
3698 2 : natfix = 0
3699 42 : do ii=1,size(multibinit_dtset%iatfix,2)
3700 134 : if(any(multibinit_dtset%iatfix(:,ii) /= 0))then
3701 14 : natfix = natfix + 1
3702 : endif
3703 : enddo
3704 2 : write(nunit,'(3x,a9,3I10)')' natfix', natfix
3705 : endif
3706 :
3707 : if(multibinit_dtset%dynamics==101)then
3708 : end if
3709 :
3710 66 : if(multibinit_dtset%dynamics==102)then
3711 2 : write(nunit,'(a15,ES15.5)')'latt_friction',multibinit_dtset%latt_friction
3712 : end if
3713 :
3714 66 : if(multibinit_dtset%dynamics==103)then
3715 2 : write(nunit,'(a15,ES15.5)')' latt_taut',multibinit_dtset%latt_taut
3716 : end if
3717 :
3718 66 : if(multibinit_dtset%dynamics==104)then
3719 0 : write(nunit,'(a15,ES15.5)')' latt_taut',multibinit_dtset%latt_taut
3720 : ! write(nunit,'(a15,ES15.5)')' latt_taup',multibinit_dtset%latt_taup
3721 : ! write(nunit,'(a15,ES15.5)')'compressibility',multibinit_dtset%latt_compressibility
3722 : end if
3723 :
3724 66 : if(multibinit_dtset%dynamics==105)then
3725 0 : write(nunit,'(a15,ES15.5)')' latt_taut',multibinit_dtset%latt_taut
3726 : ! write(nunit,'(a15,ES15.5)')' latt_taup',multibinit_dtset%latt_taup
3727 : ! write(nunit,'(a15,ES15.5)')'compressibility',multibinit_dtset%latt_compressibility
3728 : ! write(nunit,'(a15,ES15.5)')' latt_mask',(multibinit_dtset%latt_mask(ii), ii=1, 3)
3729 : end if
3730 : end if
3731 :
3732 98 : if(multibinit_dtset%lwf_dynamics/=0) then
3733 2 : write(nunit,'(a)')' LWF Dynamics :'
3734 2 : write(nunit,'(12x,a16,I12.1)')'lwf_dynamics',multibinit_dtset%lwf_dynamics
3735 2 : write(nunit, '(13x, a15, I12.1)') 'lwf_init_state', multibinit_dtset%lwf_init_state
3736 2 : write(nunit,'(10x, a18, 5x, F10.5)')'lwf_temperature',multibinit_dtset%lwf_temperature
3737 2 : write(nunit,'(13x,a15,ES15.5, a8)') 'lwf_dt',multibinit_dtset%lwf_dt*Time_Sec , ' second' !TODO: use a.u.
3738 2 : write(nunit,'(13x,a15,I10.1)')'lwf_ntime',multibinit_dtset%lwf_ntime
3739 2 : write(nunit,'(13x,a15,I10.1)')'lwf_nctime',multibinit_dtset%lwf_nctime
3740 : !write(nunit,'(8x,a20,I10.1)')'lwf_self_bound_order',multibinit_dtset%lwf_self_bound_order
3741 : !write(nunit,'(8x,a20,F10.5)')'lwf_self_bound_coeff',multibinit_dtset%lwf_self_bound_coeff
3742 2 : write(nunit, '(6x, a22, I12.1)') 'lwf_var_temperature', multibinit_dtset%lwf_var_temperature
3743 2 : write(nunit, '(6x, a22, 5x, F10.5)') 'lwf_temperature_start', multibinit_dtset%lwf_temperature_start
3744 2 : write(nunit, '(6x, a22, 5x, F10.5)') 'lwf_temperature_end', multibinit_dtset%lwf_temperature_end
3745 2 : write(nunit, '(5x, a23, I12.1)') 'lwf_temperature_nstep', multibinit_dtset%lwf_temperature_nstep
3746 :
3747 2 : if(multibinit_dtset%lwf_dynamics==1) then
3748 0 : write(nunit,'(20x,a8,ES15.5, a8)') 'lwf_taut',multibinit_dtset%lwf_dt*Time_Sec , ' second' !TODO: use a.u.
3749 2 : else if(multibinit_dtset%lwf_dynamics==3) then
3750 2 : write(nunit,'(10x, a18, 5x, F10.5)')'lwf_mc_avg_amp',multibinit_dtset%lwf_mc_avg_amp
3751 : end if
3752 : end if
3753 :
3754 98 : if(multibinit_dtset%spin_dynamics/=0) then
3755 4 : write(nunit,'(a)')' Spin Dynamics :'
3756 :
3757 : !write(nunit,'(3x,a25,I12.1)')'spin_calc_correlation_obs',multibinit_dtset%spin_calc_correlation_obs
3758 4 : write(nunit,'(3x,a25,I12.1)')'spin_calc_thermo_obs',multibinit_dtset%spin_calc_thermo_obs
3759 : !write(nunit,'(3x,a25,I12.1)')'spin_calc_traj_obs',multibinit_dtset%spin_calc_traj_obs
3760 4 : write(nunit,'(12x,a16,I12.1)')'spin_dynamics',multibinit_dtset%spin_dynamics
3761 4 : write(nunit,'(10x, a18, 5x, F10.5)')'spin_temperature',multibinit_dtset%spin_temperature
3762 4 : write(nunit,'(10x, a18, 5x, F10.5)')'spin_damping',multibinit_dtset%spin_damping
3763 4 : write(nunit,'(9x,a19,I10.1)')'spin_ntime_pre',multibinit_dtset%spin_ntime_pre
3764 4 : write(nunit,'(13x,a15,I10.1)')'spin_ntime',multibinit_dtset%spin_ntime
3765 4 : write(nunit,'(13x,a15,3I10)') 'ncell',multibinit_dtset%ncell !TODO hexu: duplicate but dynamics can be 0.
3766 4 : write(nunit,'(13x,a15,ES15.5, a8)') 'spin_dt',multibinit_dtset%spin_dt*Time_Sec , ' second' !TODO: use a.u.
3767 : !write(nunit,'(3x,a14,3es10.5)') ' spin_tolavg',multibinit_dtset%spin_tolavg
3768 : !write(nunit,'(3x,a14,3es10.5)') ' spin_tolvar',multibinit_dtset%spin_tolvar
3769 16 : write(nunit,'(13x,a15,3ES15.5,a8)') 'spin_mag_field', (multibinit_dtset%spin_mag_field(ii)/Bfield_Tesla,ii=1,3), ' Tesla'
3770 4 : write(nunit, '(13x, a15, I12.1)') 'spin_sia_add', multibinit_dtset%spin_sia_add
3771 4 : write(nunit, '(13x, a15, ES15.5)') 'spin_sia_k1amp', multibinit_dtset%spin_sia_k1amp
3772 16 : write(nunit, '(13x, a15, 3ES15.5)') 'spin_sia_k1dir', (multibinit_dtset%spin_sia_k1dir(ii), ii=1,3)
3773 16 : write(nunit,'(5x,a23, 3es15.5)') 'spin_projection_qpoint', (multibinit_dtset%spin_projection_qpoint(ii),ii=1,3)
3774 4 : write(nunit, '(13x, a15, I12.1)') 'spin_init_state', multibinit_dtset%spin_init_state
3775 4 : if(multibinit_dtset%spin_init_state==3) then
3776 0 : write(nunit,'(6x,a22, 3es15.5)') 'spin_init_orientation', (multibinit_dtset%spin_init_orientation(ii),ii=1,3)
3777 0 : write(nunit,'(11x,a17, 3es15.5)') 'spin_init_qpoint', (multibinit_dtset%spin_init_qpoint(ii),ii=1,3)
3778 0 : write(nunit,'(6x,a22,3es15.5)') 'spin_init_rotate_axis', (multibinit_dtset%spin_init_rotate_axis(ii),ii=1,3)
3779 : endif
3780 4 : write(nunit, '(6x, a22, I12.1)') 'spin_var_temperature', multibinit_dtset%spin_var_temperature
3781 4 : write(nunit, '(6x, a22, 5x, F10.5)') 'spin_temperature_start', multibinit_dtset%spin_temperature_start
3782 4 : write(nunit, '(6x, a22, 5x, F10.5)') 'spin_temperature_end', multibinit_dtset%spin_temperature_end
3783 4 : write(nunit, '(5x, a23, I12.1)') 'spin_temperature_nstep', multibinit_dtset%spin_temperature_nstep
3784 4 : write(nunit, '(13x, a15, I12.1)') 'spin_write_traj', multibinit_dtset%spin_write_traj
3785 : end if
3786 :
3787 98 : if(multibinit_dtset%slc_coupling/=0) then
3788 4 : write(nunit,'(6x, a22, I12.4)') 'spin-lattice coupling', multibinit_dtset%slc_coupling
3789 : endif
3790 :
3791 98 : if(multibinit_dtset%confinement==1)then
3792 0 : write(nunit,'(a)')' Confinement information :'
3793 0 : write(nunit,'(1x,a22,I5.1)')' conf_power_disp',multibinit_dtset%conf_power_disp
3794 0 : write(nunit,'(1x,a22,I5.1)')' conf_power_strain',multibinit_dtset%conf_power_strain
3795 0 : write(nunit,'(1x,a22,3es16.8)')' conf_power_fact_disp',multibinit_dtset%conf_power_fact_disp
3796 0 : write(nunit,'(1x,a22,3es16.8)')'conf_power_fact_strain',multibinit_dtset%conf_power_fact_strain
3797 0 : write(nunit,'(1x,a22)')' conf_cutoff_disp'
3798 0 : write(nunit,'(19x,3es16.8)') (multibinit_dtset%conf_cutoff_disp(ii),ii=1,multibinit_dtset%natom)
3799 0 : write(nunit,'(1x,a22)')' conf_cutoff_strain'
3800 0 : write(nunit,'(19x,3es16.8)') (multibinit_dtset%conf_cutoff_strain(ii),ii=1,6)
3801 : end if
3802 :
3803 98 : if(multibinit_dtset%fit_coeff/=0)then
3804 20 : write(nunit,'(a)')' Fit the coefficients :'
3805 20 : write(nunit,'(1x,a17,I3.1)')' fit_coeff',multibinit_dtset%fit_coeff
3806 20 : write(nunit,'(1x,a17,I3.1)')'fit_generateCoeff',multibinit_dtset%fit_generateCoeff
3807 20 : if(multibinit_dtset%fit_initializeData==0)then
3808 20 : write(nunit,'(1x,a17,I3.1)')'fit_initializeData',multibinit_dtset%fit_initializeData
3809 : end if
3810 20 : if(multibinit_dtset%fit_tolMSDE > 0)then
3811 0 : write(nunit,'(1x,a17,es16.8)')' fit_tolMSDE',multibinit_dtset%fit_tolMSDE
3812 : end if
3813 20 : if(multibinit_dtset%fit_tolMSDF > 0)then
3814 0 : write(nunit,'(1x,a17,es16.8)')' fit_tolMSDF',multibinit_dtset%fit_tolMSDF
3815 : end if
3816 20 : if(multibinit_dtset%fit_tolMSDS > 0)then
3817 0 : write(nunit,'(1x,a17,es16.8)')' fit_tolMSDS',multibinit_dtset%fit_tolMSDS
3818 : end if
3819 20 : if(multibinit_dtset%fit_tolMSDFS > 0)then
3820 0 : write(nunit,'(1x,a17,es16.8)')' fit_tolMSDFS',multibinit_dtset%fit_tolMSDFS
3821 : end if
3822 20 : if(multibinit_dtset%fit_weight_T > 0.0_dp)then
3823 2 : write(nunit,'(1x,a17,es16.8)')' fit_weight_T',multibinit_dtset%fit_weight_T
3824 : end if
3825 :
3826 20 : write(nunit,'(1x,a17,es16.8)')' fit_cutoff',multibinit_dtset%fit_cutoff
3827 20 : write(nunit,'(1x,a17,es16.8)')' fit_droprate',multibinit_dtset%fit_drop_rate
3828 : !write(nunit,'(1x,a17,es16.8)')' fit_min_bound_coeff',multibinit_dtset%fit_min_bound_coeff
3829 20 : write(nunit,'(1x,a17,I3.1)')' fit_option',multibinit_dtset%fit_option
3830 20 : write(nunit,'(1x,a17,2x,I0)')' fit_iatom',multibinit_dtset%fit_iatom
3831 20 : write(nunit,'(1x,a17,2x,I0)')' fit_ncoeff',multibinit_dtset%fit_ncoeff
3832 20 : write(nunit,'(1x,a17,3i3)') ' fit_grid',multibinit_dtset%fit_grid
3833 20 : write(nunit,'(1x,a17,I3.1)')' ts_option',multibinit_dtset%ts_option
3834 20 : write(nunit,'(1x,a17,2i3)') ' fit_rangePower',multibinit_dtset%fit_rangePower
3835 20 : write(nunit,'(1x,a17,I3)') ' fit_dispterms',multibinit_dtset%fit_dispterms
3836 20 : write(nunit,'(1x,a17,I3)') ' fit_anhaStrain',multibinit_dtset%fit_anhaStrain
3837 20 : write(nunit,'(1x,a17,I3)') ' fit_SPCoupling',multibinit_dtset%fit_SPCoupling
3838 20 : write(nunit,'(1x,a17,I3)') ' fit_SPC_maxS',multibinit_dtset%fit_SPC_maxS
3839 :
3840 :
3841 20 : n=multibinit_dtset%fit_rangePower(2)-multibinit_dtset%fit_rangePower(1)+1
3842 20 : if(n > 0) then
3843 20 : write(nunit,'(1x,a17)',advance='no') ' fit_max_nbody'
3844 42 : write(nunit,'(4x,9i7)') (multibinit_dtset%fit_max_nbody_inrange(ii),ii=1,n)
3845 : end if
3846 :
3847 :
3848 20 : if(multibinit_dtset%fit_nbancoeff /= 0) then
3849 0 : write(nunit,'(1x,a17,I3)') ' fit_nbancoeff',multibinit_dtset%fit_nbancoeff
3850 0 : write(nunit,'(1x,a17)',advance='no')' fit_bancoeff'
3851 0 : write(nunit,'(4x,9i7)') (multibinit_dtset%fit_bancoeff(ii),ii=1,multibinit_dtset%fit_nbancoeff)
3852 : end if
3853 20 : if(multibinit_dtset%fit_nfixcoeff /= 0) then
3854 6 : write(nunit,'(1x,a17,I3)') ' fit_nfixcoeff',multibinit_dtset%fit_nfixcoeff
3855 6 : write(nunit,'(1x,a17)',advance='no')' fit_fixcoeff'
3856 18 : write(nunit,'(4x,9i7)') (multibinit_dtset%fit_fixcoeff(ii),ii=1,multibinit_dtset%fit_nfixcoeff)
3857 : end if
3858 20 : if(multibinit_dtset%fit_nimposecoeff /= 0) then
3859 0 : write(nunit,'(1x,a17,I3)') ' fit_nimposecoeff',multibinit_dtset%fit_nimposecoeff
3860 0 : write(nunit,'(1x,a17)',advance='no')' fit_imposecoeff'
3861 0 : write(nunit,'(4x,9i7)') (multibinit_dtset%fit_imposecoeff(ii),ii=1,multibinit_dtset%fit_nimposecoeff)
3862 : end if
3863 : end if
3864 :
3865 98 : if(multibinit_dtset%opt_effpot/=0)then
3866 2 : write(nunit,'(a)')'Optimize Effective Potential :'
3867 2 : if(multibinit_dtset%opt_ncoeff /= 0) then
3868 2 : write(nunit,'(1x,a17,I3)') ' opt_ncoeff',multibinit_dtset%opt_ncoeff
3869 2 : write(nunit,'(1x,a17)',advance='no')' opt_coeff'
3870 6 : write(nunit,'(4x,9i7)') (multibinit_dtset%opt_coeff(ii),ii=1,multibinit_dtset%opt_ncoeff)
3871 : end if
3872 : end if
3873 :
3874 98 : if(multibinit_dtset%bound_model /=0)then
3875 2 : write(nunit,'(a)')' Bound the coefficients :'
3876 2 : write(nunit,'(1x,a16,I3.1)') ' bound_model',multibinit_dtset%bound_model
3877 2 : write(nunit,'(1x,a16,es12.4)') ' bound_penalty',multibinit_dtset%bound_penalty
3878 2 : write(nunit,'(1x,a16,I3)') 'bound_anhaStrain',multibinit_dtset%bound_anhaStrain
3879 2 : write(nunit,'(1x,a16,I3)') 'bound_SPCoupling',multibinit_dtset%bound_SPCoupling
3880 2 : write(nunit,'(1x,a16,es16.8)')' bound_cutoff',multibinit_dtset%bound_cutoff
3881 2 : write(nunit,'(1x,a16,1x,3I3)') ' bound_cell',multibinit_dtset%bound_cell
3882 2 : write(nunit,'(1x,a16,1x,I3)') ' bound_maxCoeff',multibinit_dtset%bound_maxCoeff
3883 2 : write(nunit,'(1x,a16,1x,I3)') ' bound_option',multibinit_dtset%bound_option
3884 2 : write(nunit,'(1x,a16,es16.8)') ' bound_temp',multibinit_dtset%bound_temp
3885 2 : write(nunit,'(1x,a16,I7)') ' bound_step',multibinit_dtset%bound_step
3886 2 : write(nunit,'(1x,a16,2I3.1)')'bound_rangePower',multibinit_dtset%bound_rangePower
3887 : end if
3888 :
3889 :
3890 : !MS Variables for SCALE-UP
3891 : !Replace by function from m_scup_dataset.F90
3892 98 : if(multibinit_dtset%scup_dtset%scup_elec_model)then
3893 0 : call outvars_scup(multibinit_dtset%scup_dtset,nunit)
3894 : end if
3895 :
3896 : !Write the general information
3897 : if( multibinit_dtset%rfmeth/=1 .or. &
3898 : & multibinit_dtset%enunit/=0 .or. &
3899 : & multibinit_dtset%eivec/=0 .or. &
3900 98 : & multibinit_dtset%asr/=0 .or. &
3901 : & multibinit_dtset%chneut/=0)then
3902 98 : write(nunit,'(a)')' Miscellaneous information :'
3903 98 : if(multibinit_dtset%rfmeth/=1)write(nunit,'(3x,a9,3i10)')' rfmeth',multibinit_dtset%rfmeth
3904 98 : if(multibinit_dtset%enunit/=0)write(nunit,'(3x,a9,3i10)')' enunit',multibinit_dtset%enunit
3905 98 : if(multibinit_dtset%eivec/=0) write(nunit,'(3x,a9,3i10)')' eivec',multibinit_dtset%eivec
3906 98 : if(multibinit_dtset%asr/=0) write(nunit,'(3x,a9,3i10)')' asr',multibinit_dtset%asr
3907 98 : if(multibinit_dtset%chneut/=0)write(nunit,'(3x,a9,3i10)')' chneut',multibinit_dtset%chneut
3908 : end if
3909 :
3910 :
3911 : !For interatomic force constant information
3912 98 : if(multibinit_dtset%ifcflag/=0)then
3913 98 : write(nunit,'(a)')' Interatomic Force Constants Inputs :'
3914 98 : write(nunit,'(3x,a9,3i10)')' dipdip',multibinit_dtset%dipdip
3915 98 : if(multibinit_dtset%dipdip /= 0)then
3916 66 : write(nunit,'(a12,3i10)') 'dipdip_range',multibinit_dtset%dipdip_range
3917 : end if
3918 98 : if(multibinit_dtset%dipdip_prt/=0)then
3919 42 : write(nunit,'(a12,3i10)') 'dipdip_prt',multibinit_dtset%dipdip_prt
3920 : end if
3921 98 : if(multibinit_dtset%nsphere/=0)write(nunit,'(3x,a9,3i10)')' nsphere',multibinit_dtset%nsphere
3922 98 : if(abs(multibinit_dtset%rifcsph)>tol10)write(nunit,'(3x,a9,E16.6)')' nsphere',multibinit_dtset%rifcsph
3923 98 : write(nunit,'(3x,a9,3i10)')' ifcana',multibinit_dtset%ifcana
3924 98 : write(nunit,'(3x,a9,3i10)')' ifcout',multibinit_dtset%ifcout
3925 98 : if(multibinit_dtset%natifc>=1)then
3926 86 : write(nunit,'(3x,a9,3i10)')' natifc',multibinit_dtset%natifc
3927 86 : write(nunit,'(3x,a12)',advance='no')' atifc '
3928 582 : write(nunit,'(3x,15i4)') (multibinit_dtset%atifc(ii)*ii,ii=1,multibinit_dtset%natifc)
3929 :
3930 : end if
3931 98 : write(nunit,'(a)')' Description of grid 1 :'
3932 98 : write(nunit,'(3x,a9,3i10)')' brav',multibinit_dtset%brav
3933 98 : write(nunit,'(3x,a9,3i10)')' ngqpt',multibinit_dtset%ngqpt(1:3)
3934 98 : write(nunit,'(3x,a9,3i10)')' nqshft',multibinit_dtset%nqshft
3935 98 : if (multibinit_dtset%nqshft/=0)then
3936 98 : write(nunit,'(3x,a9)')' q1shft'
3937 196 : do iqshft=1,multibinit_dtset%nqshft
3938 490 : write(nunit,'(19x,4es16.8)') (multibinit_dtset%q1shft(ii,iqshft),ii=1,3)
3939 : end do
3940 : end if
3941 392 : if (any(multibinit_dtset%qrefine(:) > 1)) then
3942 0 : write(nunit,'(3x,a9,3i10)')' qrefine', multibinit_dtset%qrefine
3943 : end if
3944 : end if
3945 :
3946 :
3947 : !List of vector 1 (reduced coordinates)
3948 98 : if(multibinit_dtset%nph1l/=0)then
3949 98 : write(nunit,'(a)')' First list of wavevector (reduced coord.) :'
3950 98 : write(nunit,'(3x,a9,3i10)')' nph1l',multibinit_dtset%nph1l
3951 98 : write(nunit,'(3x,a9)')' qph1l'
3952 252 : do iph1=1,multibinit_dtset%nph1l
3953 : write(nunit,'(19x,3es16.8,2x,es11.3)') &
3954 714 : & (multibinit_dtset%qph1l(ii,iph1),ii=1,3),multibinit_dtset%qnrml1(iph1)
3955 : end do
3956 : end if
3957 :
3958 : !List of vector 2 (cartesian coordinates)
3959 98 : if(multibinit_dtset%nph2l/=0)then
3960 0 : write(nunit,'(a)')' Second list of wavevector (cart. coord.) :'
3961 0 : write(nunit,'(3x,a9,3i10)')' nph2l',multibinit_dtset%nph2l
3962 0 : write(nunit,'(3x,a9)')' qph2l'
3963 0 : do iph2=1,multibinit_dtset%nph2l
3964 : write(nunit,'(19x,3es16.8,2x,es11.3)') &
3965 0 : & (multibinit_dtset%qph2l(ii,iph2),ii=1,3),multibinit_dtset%qnrml2(iph2)
3966 : end do
3967 : end if
3968 :
3969 7938 : write(nunit,'(a,80a,a)') ch10,('=',ii=1,80),ch10
3970 :
3971 98 : end subroutine outvars_multibinit
3972 : !!***
3973 :
3974 :
3975 :
3976 0 : end module m_multibinit_dataset
3977 : !!***
|
