~ubuntu-branches/ubuntu/saucy/nwchem/saucy

« back to all changes in this revision

Viewing changes to src/smd/smd-subgroups/smd_res.F

  • Committer: Package Import Robot
  • Author(s): Michael Banck, Michael Banck, Daniel Leidert
  • Date: 2012-02-09 20:02:41 UTC
  • mfrom: (1.1.1)
  • Revision ID: package-import@ubuntu.com-20120209200241-jgk03qfsphal4ug2
Tags: 6.1-1
* New upstream release.

[ Michael Banck ]
* debian/patches/02_makefile_flags.patch: Updated.
* debian/patches/02_makefile_flags.patch: Use internal blas and lapack code.
* debian/patches/02_makefile_flags.patch: Define GCC4 for LINUX and LINUX64
  (Closes: #632611 and LP: #791308).
* debian/control (Build-Depends): Added openssh-client.
* debian/rules (USE_SCALAPACK, SCALAPACK): Removed variables (Closes:
  #654658).
* debian/rules (LIBDIR, USE_MPIF4, ARMCI_NETWORK): New variables.
* debian/TODO: New file.
* debian/control (Build-Depends): Removed libblas-dev, liblapack-dev and
  libscalapack-mpi-dev.
* debian/patches/04_show_testsuite_diff_output.patch: New patch, shows the
  diff output for failed tests.
* debian/patches/series: Adjusted.
* debian/testsuite: Optionally run all tests if "all" is passed as option.
* debian/rules: Run debian/testsuite with "all" if DEB_BUILD_OPTIONS
  contains "checkall".

[ Daniel Leidert ]
* debian/control: Used wrap-and-sort. Added Vcs-Svn and Vcs-Browser fields.
  (Priority): Moved to extra according to policy section 2.5.
  (Standards-Version): Bumped to 3.9.2.
  (Description): Fixed a typo.
* debian/watch: Added.
* debian/patches/03_hurd-i386_define_path_max.patch: Added.
  - Define MAX_PATH if not defines to fix FTBFS on hurd.
* debian/patches/series: Adjusted.

Show diffs side-by-side

added added

removed removed

Lines of Context:
1
 
      subroutine smd_res_init(namespace,rtdb)
2
 
      implicit none
3
 
#include "errquit.fh"
4
 
#include "inp.fh"
5
 
#include "mafdecls.fh"
6
 
#include "rtdb.fh"
7
 
#include "util.fh"
8
 
#include "global.fh"
9
 
c     
10
 
      character*(*) namespace
11
 
      integer rtdb
12
 
c
13
 
      character*32 pname
14
 
      character*80 tag
15
 
      character*255 filename
16
 
      integer nr,nt
17
 
      integer i_s,i_t,i_ind
18
 
      logical result
19
 
c
20
 
      pname = "smd_res_init"
21
 
c
22
 
      write(*,*) "in "//pname
23
 
c
24
 
      if(.not.rtdb_cget(rtdb,'smd:coordfile',1,filename))
25
 
     >  call errquit(
26
 
     >       pname//'error reading coordinate file',0, RTDB_ERR)
27
 
      write(*,*) "found coordinates",filename
28
 
c
29
 
c     get number of residues
30
 
c     ----------------------
31
 
      call smd_pdb_nres(filename,nr)
32
 
      write(*,*) "number of res",nr
33
 
c
34
 
c     get number of atoms
35
 
c     -------------------
36
 
      call smd_atom_ntot(namespace,nt)
37
 
      write(*,*) "number of atoms",nt
38
 
c
39
 
c     create residue data structures
40
 
c     ------------------------------
41
 
      call smd_data_create(namespace,"res:name",16*nr,MT_BYTE)
42
 
      call smd_data_create(namespace,"res:size",nr,MT_INT)
43
 
      call smd_data_create(namespace,"res:",nr,MT_INT)
44
 
 
45
 
      tag = "atom:res"
46
 
      call smd_data_get_index(namespace,tag,i_ind,result)
47
 
      if(.not. result) 
48
 
     >  call errquit(
49
 
     >       pname//'error getting index for '//tag,0, RTDB_ERR)
50
 
 
51
 
      tag = "res:size"
52
 
      call smd_data_get_index(namespace,tag,i_s,result)
53
 
      if(.not. result) 
54
 
     >  call errquit(
55
 
     >       pname//'error getting index for '//tag,0, RTDB_ERR)
56
 
 
57
 
      tag = "res:name"
58
 
      call smd_data_get_index(namespace,tag,i_t,result)
59
 
      if(.not. result) 
60
 
     >  call errquit(
61
 
     >       pname//'error getting index for '//tag,0, RTDB_ERR)
62
 
 
63
 
 
64
 
      call smd_pdb_read_res(filename,nt,nr,
65
 
     +                         byte_mb(i_t),
66
 
     +                         int_mb(i_ind),
67
 
     +                         int_mb(i_s))
68
 
 
69
 
      return
70
 
      end