88
88
character*12 control_boundry
89
89
external control_boundry
91
logical brillioun_print
91
logical brillioun_print,control_Mulliken
92
92
integer brillioun_nbrillioun,brillioun_nbrillq
93
93
real*8 brillioun_weight_brdcst,brillioun_ks_brdcst
94
94
real*8 brillioun_k_brdcst,brillioun_weight
95
external brillioun_print
95
external brillioun_print,control_Mulliken
96
96
external brillioun_nbrillioun,brillioun_nbrillq
97
97
external brillioun_weight_brdcst,brillioun_ks_brdcst
98
98
external brillioun_k_brdcst,brillioun_weight
148
148
c character*255 cpsp_comment,comment
149
149
c external cpsp_comment
151
integer ion_nconstraints,ion_ndof
152
external ion_nconstraints,ion_ndof
344
347
write(luout,1130) "spin orbit"
350
call v_bwexc_print(luout,control_gga())
347
IF (control_gga().eq.0) THEN
348
write(luout,1131) 'Vosko et al parameterization'
349
ELSE IF (control_gga().eq.10) THEN
351
> 'PBE96 (White and Bird) parameterization'
352
ELSE IF (control_gga().eq.11) THEN
354
> 'BLYP (White and Bird) parameterization'
355
ELSE IF (control_gga().eq.12) THEN
357
> 'revPBE (White and Bird) parameterization'
358
ELSE IF (control_gga().eq.13) THEN
360
> 'PBEsol (White and Bird) parameterization'
362
ELSE IF (control_gga().eq.110) THEN
364
> 'PBE0 (White and Bird) parameterization'
365
ELSE IF (control_gga().eq.112) THEN
367
> 'revPBE0 (White and Bird) parameterization'
368
ELSE IF (control_gga().eq.113) THEN
370
> 'BNL (White and Bird) parameterization'
371
ELSE IF (control_gga().eq.114) THEN
373
> 'HSE (White and Bird) parameterization'
374
ELSE IF (control_gga().eq.200) THEN
375
write(luout,1131) 'Hartree-Fock'
377
write(luout,1131) 'unknown parameterization'
378
call errquit('bad exchange_correlation',0, INPUT_ERR)
380
352
call band_print_HFX(luout)
381
353
write(luout,1140)
382
354
do ia = 1,ion_nkatm()
415
387
write(luout,1200) cx,cy,cz
416
388
write(luout,1210) gx,gy,gz
389
write(luout,1211) ion_nconstraints(),ion_ndof()
418
391
write(luout,1220) en(1),en(cpsi_ispin()),
419
392
> ' (fourier space)'
740
717
> ' ) - atomic mass= ',F7.3,' - fixed')
741
718
1200 FORMAT(5X,' G.C. ',' (',3F11.5,' )')
742
719
1210 FORMAT(5X,' C.O.M.',' (',3F11.5,' )')
720
1211 FORMAT(5X,' number of constraints = ', I6,' ( DOF = ',I6,' )' )
743
721
1220 FORMAT(/' number of electrons: spin up=',F8.2,
744
722
> ' spin down=',F8.2,A)
745
723
1221 FORMAT( ' number of orbitals: spin up=',I8,