~ubuntu-branches/ubuntu/utopic/nwchem/utopic

« back to all changes in this revision

Viewing changes to src/nwdft/rt_tddft/input/rt_tddft_input_print.F

  • Committer: Package Import Robot
  • Author(s): Michael Banck, Daniel Leidert, Andreas Tille, Michael Banck
  • Date: 2013-07-04 12:14:55 UTC
  • mfrom: (1.1.2)
  • Revision ID: package-import@ubuntu.com-20130704121455-5tvsx2qabor3nrui
Tags: 6.3-1
* New upstream release.
* Fixes anisotropic properties (Closes: #696361).
* New features include:
  + Multi-reference coupled cluster (MRCC) approaches
  + Hybrid DFT calculations with short-range HF 
  + New density-functionals including Minnesota (M08, M11) and HSE hybrid
    functionals
  + X-ray absorption spectroscopy (XAS) with TDDFT
  + Analytical gradients for the COSMO solvation model
  + Transition densities from TDDFT 
  + DFT+U and Electron-Transfer (ET) methods for plane wave calculations
  + Exploitation of space group symmetry in plane wave geometry optimizations
  + Local density of states (LDOS) collective variable added to Metadynamics
  + Various new XC functionals added for plane wave calculations, including
    hybrid and range-corrected ones
  + Electric field gradients with relativistic corrections 
  + Nudged Elastic Band optimization method
  + Updated basis sets and ECPs 

[ Daniel Leidert ]
* debian/watch: Fixed.

[ Andreas Tille ]
* debian/upstream: References

[ Michael Banck ]
* debian/upstream (Name): New field.
* debian/patches/02_makefile_flags.patch: Refreshed.
* debian/patches/06_statfs_kfreebsd.patch: Likewise.
* debian/patches/07_ga_target_force_linux.patch: Likewise.
* debian/patches/05_avoid_inline_assembler.patch: Removed, no longer needed.
* debian/patches/09_backported_6.1.1_fixes.patch: Likewise.
* debian/control (Build-Depends): Added gfortran-4.7 and gcc-4.7.
* debian/patches/10_force_gcc-4.7.patch: New patch, explicitly sets
  gfortran-4.7 and gcc-4.7, fixes test suite hang with gcc-4.8 (Closes:
  #701328, #713262).
* debian/testsuite: Added tests for COSMO analytical gradients and MRCC.
* debian/rules (MRCC_METHODS): New variable, required to enable MRCC methods.

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
C
 
2
C     rt_tddft_input_print.F
 
3
C     
 
4
C     Parses input deck for rt-tddft print options.
 
5
C     
 
6
C
 
7
      subroutine rt_tddft_input_print (rtdb)
 
8
      implicit none
 
9
 
 
10
#include "rt_tddft.fh"
 
11
#include "errquit.fh"
 
12
#include "inp.fh"
 
13
#include "rtdb.fh"
 
14
#include "mafdecls.fh"
 
15
#include "stdio.fh"
 
16
 
 
17
      
 
18
C     == Inputs ==
 
19
      integer, intent(in) :: rtdb
 
20
 
 
21
 
 
22
C     == Parameters ==
 
23
      character(*), parameter :: pname = "rt_tddft_input_print: "
 
24
 
 
25
 
 
26
C     == Variables ==
 
27
      logical got_opt
 
28
      character*255 curr_popt
 
29
      logical print_dipole
 
30
      logical print_quad
 
31
      logical print_field
 
32
      logical print_moocc
 
33
      logical print_energy
 
34
      logical print_cputime
 
35
      logical print_charge
 
36
      logical print_convergence
 
37
      logical print_s2
 
38
      
 
39
 
 
40
C     (zero print options if print directive found)
 
41
      print_dipole = .false.
 
42
      print_quad = .false.
 
43
      print_field = .false. 
 
44
      print_moocc = .false.
 
45
      print_energy = .false. 
 
46
      print_cputime = .false.
 
47
      print_charge = .false.
 
48
      print_convergence = .false.
 
49
      print_s2 = .false.
 
50
            
 
51
 
 
52
C     (scan over all words found after "print" directive)
 
53
      got_opt = .true.
 
54
      
 
55
      do while (got_opt)
 
56
         got_opt = inp_a (curr_popt)
 
57
         
 
58
         if ( trim(curr_popt) .eq. "") exit
 
59
         
 
60
C     Now that have an option, toggle print flag if it matches.  If no
 
61
C     match, print an error.
 
62
         if (curr_popt .eq. "dipole") then 
 
63
            print_dipole = .true.
 
64
         elseif (curr_popt .eq. "quadrupole") then 
 
65
            print_quad = .true.
 
66
         elseif (curr_popt .eq. "field") then 
 
67
            print_field = .true.
 
68
         elseif (curr_popt .eq. "moocc") then 
 
69
            print_moocc = .true.
 
70
         elseif (curr_popt .eq. "energy") then 
 
71
            print_energy = .true.
 
72
         elseif (curr_popt .eq. "cputime") then 
 
73
            print_cputime = .true.
 
74
         elseif (curr_popt .eq. "charge") then 
 
75
            print_charge = .true.
 
76
         elseif (curr_popt .eq. "convergence") then 
 
77
            print_convergence = .true.
 
78
         elseif (curr_popt .eq. "s2") then 
 
79
            print_s2 = .true.
 
80
         elseif (curr_popt .eq. "*") then 
 
81
            print_dipole = .true.
 
82
            print_quad = .true.
 
83
            print_field = .true. 
 
84
            print_moocc = .true.
 
85
            print_energy = .true. 
 
86
            print_cputime = .true.
 
87
            print_charge = .true.
 
88
            print_convergence = .true.
 
89
            print_s2 = .true.
 
90
            
 
91
         else
 
92
            call errquit (pname//
 
93
     $           "unknown print option: "//trim(curr_popt)//
 
94
     $           ".  Allowed options are: dipole, field, moocc,"//
 
95
     $           " energy, cputime, charge, convergence, s2, *.",
 
96
     $           0, 0)
 
97
         endif
 
98
      enddo
 
99
 
 
100
C
 
101
C     Save all params to rtdb, overwriting the defaults.
 
102
C     
 
103
      if (.not.rtdb_put(rtdb,'rt_tddft:print_dipole',
 
104
     $     mt_log,1,print_dipole))
 
105
     $     call errquit(pname//'Write failed to rtdb',0,RTDB_ERR)
 
106
 
 
107
      if (.not.rtdb_put(rtdb,'rt_tddft:print_quad',
 
108
     $     mt_log,1,print_quad))
 
109
     $     call errquit(pname//'Write failed to rtdb',0,RTDB_ERR)
 
110
 
 
111
      if (.not.rtdb_put(rtdb,'rt_tddft:print_field',
 
112
     $     mt_log,1,print_field))
 
113
     $     call errquit(pname//'Write failed to rtdb',0,RTDB_ERR)
 
114
 
 
115
      if (.not.rtdb_put(rtdb,'rt_tddft:print_moocc',
 
116
     $     mt_log,1,print_moocc))
 
117
     $     call errquit(pname//'Write failed to rtdb',0,RTDB_ERR)
 
118
 
 
119
      if (.not.rtdb_put(rtdb,'rt_tddft:print_energy',
 
120
     $     mt_log,1,print_energy))
 
121
     $     call errquit(pname//'Write failed to rtdb',0,RTDB_ERR)
 
122
 
 
123
      if (.not.rtdb_put(rtdb,'rt_tddft:print_cputime',
 
124
     $     mt_log,1,print_cputime))
 
125
     $     call errquit(pname//'Write failed to rtdb',0,RTDB_ERR)
 
126
 
 
127
      if (.not.rtdb_put(rtdb,'rt_tddft:print_charge',
 
128
     $     mt_log,1,print_charge))
 
129
     $     call errquit(pname//'Write failed to rtdb',0,RTDB_ERR)
 
130
 
 
131
      if (.not.rtdb_put(rtdb,'rt_tddft:print_convergence',
 
132
     $     mt_log,1,print_convergence))
 
133
     $     call errquit(pname//'Write failed to rtdb',0,RTDB_ERR)
 
134
 
 
135
      if (.not.rtdb_put(rtdb,'rt_tddft:print_s2',
 
136
     $     mt_log,1,print_s2))
 
137
     $     call errquit(pname//'Write failed to rtdb',0,RTDB_ERR)
 
138
      
 
139
      end subroutine
 
140