← Back to branch summary

~jose-soler/siesta/unfolding

« back to all changes in this revision

Viewing changes to Pseudo/atom/compat_params.f

  • Committer: Alberto Garcia
  • Date: 2016-01-25 16:00:16 UTC
  • mto: (483.3.1 rel-4.0)
  • mto: This revision was merged to the branch mainline in revision 485.
  • Revision ID: albertog@icmab.es-20160125160016-c1qivg1zw06kgyfd
Prepare GPL release

* Include proper headers

* Add Docs/Contributors.txt and NOTICE.txt files.

* Update READMEs and LICENSE files in several directories.

* Remove Pseudo/atom, Util/test-xml

* Remove DOM files from Src/xmlparser

Show diffs side-by-side

added added

removed removed

Lines of Context:
Pseudo/atom/compat_params.f
1
 
c
2
 
      subroutine compat_params(str)
3
 
c
4
 
c     Some internal parameters tend to change over time... This routine
5
 
c     assigns values to them based on a set of directives, specified
6
 
c     at the top of the input file as:
7
 
c     
8
 
c     %define COMPAT_UCB
9
 
c     %define NEW_CC
10
 
c  
11
 
c
12
 
 
13
 
c       
14
 
c     COMPAT_UCB :  Revert to the standard circa 1990 UCB values. Note
15
 
c                that these correspond to the first released version
16
 
c                of Jose Luis Martins code, not to the old Froyen
17
 
c                version (that would be 'froyen' --to be implemented)
18
 
c     (The default is: to use a denser grid up to larger radii.
19
 
c                Use a larger value for the ps's ecuts.
20
 
c                Use the Soler-Balbas XC package)
21
 
c
22
 
c     NEW_CC     : New core-correction scheme
23
 
c     OLD_CC     : Old core-correction scheme  (see wrapup.f)
24
 
c
25
 
c     The default is to use the new CC scheme only for GGA calculations.
26
 
c
27
 
c
28
 
c     For compatibility with an interim scheme using strings in the
29
 
c     input file, this routine accepts an argument "str". If not
30
 
c     empty, the user is warned that those strings now carry no
31
 
c     weight.
32
 
c
33
 
 
34
 
      character*(*) str
35
 
 
36
 
      include 'compat.h'
37
 
 
38
 
      logical leqi, defined
39
 
      external leqi, defined
40
 
c
41
 
      if (str .ne. " ") then
42
 
         write(6,'(a,a20)')
43
 
     $        '** WARNING: Compatibility string obsolete: ',
44
 
     $        str
45
 
         stop 'COMPAT'
46
 
      endif
47
 
 
48
 
      if (defined('COMPAT_UCB')) then
49
 
 
50
 
         write(6,'(a)') '*** UCB compatibility mode ***'
51
 
         aa_def = 6.d0
52
 
         bb_def = 40.d0
53
 
         rmax_def = 80.d0
54
 
         ecuts = 1.2d-4
55
 
         use_excorr = .true.
56
 
c
57
 
      else
58
 
 
59
 
         aa_def = 6.d0
60
 
         bb_def = 80.d0
61
 
         rmax_def = 120.d0 
62
 
         ecuts = 1.d-3
63
 
         use_excorr = .false.
64
 
 
65
 
      endif
66
 
c
67
 
c     Flag to use the old excorr subroutine.
68
 
c
69
 
      if (defined('USE_OLD_EXCORR')) use_excorr = .true.
70
 
c
71
 
c     Avoid cutting off ionic unscreened pseudopotentials
72
 
c
73
 
      if (defined('NO_PS_CUTOFFS')) ecuts = 0.d0
74
 
c
75
 
      use_old_cc = defined('OLD_CC')
76
 
      use_new_cc = defined('NEW_CC')
77
 
 
78
 
      end
79