sam/input.f at master · tslund/sam · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
      subroutine input( fname, write_params, ierr )

      include 'sam.h'

      character( *) fname
      character( 8) default_flag
      logical       write_params

      call set_labels( labels, values, required, fixed, L_params,
     &                 n_inputs_r, n_inputs_i, n_inputs,
     &                 n_dynpar_r, n_dynpar_i, n_dynpar, n_params, ierr)
      if( ierr .ne. 0 ) return

      call read_params(fname, labels, values, found, n_inputs, ierr)
      if( ierr .ne. 0 ) return

c   *** Assign input data to the /params/ common block

      call assign_params( loc(xL), loc(Nx), values,
     &                    n_inputs_r, n_inputs_i )

      call check_inputs( ierr )
      if( ierr .ne. 0 ) return

c   *** Write the input parameters to file params.out

      if( write_params ) then

         open(newunit=i_unit,file='params.out')

         do i=1, n_inputs
            default_flag = ' '
            if( .not. found(i) ) default_flag = 'DEFAULT'
            write(i_unit,*) labels(i), values(i), default_flag
         end do

         close(i_unit)

      end if

      return
      end