Changeset b4d29ce in flexpart.git for src/readcommand.f90

Timestamp:

Jul 3, 2014, 2:55:50 PM (10 years ago)

Author:

Harald Sodemann <harald.sodemann@…>

Branches:

master, 10.4.1_pesei, FPv9.3.1, FPv9.3.1b_testing, FPv9.3.2, GFS_025, NetCDF, bugfixes+enhancements, deposition, dev, fp9.3.1-20161214-nc4, grib2nc4_repair, inputlist, laptop, release-10, release-10.4.1, scaling-bug, svn-petra, svn-trunk, univie

Children:

4b093cc, 75dfd42, 326b31b

Parents:

87910af

Message:

• Implemented optional namelist input for COMMAND, RELEASES, SPECIES, AGECLASSES,OUTGRID,OUTGRID_NEST,RECEPTORS
• Implemented com_mod switch nmlout to write input files as namelist to the output directory (.true. by default)
• Proposed updated startup and runtime output (may change back to previous info if desired)

git-svn-id: ​http://flexpart.flexpart.eu:8088/svn/FlexPart90/trunk@27 ef8cc7e1-21b7-489e-abab-c1baa636049d

File:

• src/readcommand.f90 (modified) (8 diffs)

src/readcommand.f90

@@ -29 +29 @@
   !                                                                            *
   !     18 May 1996                                                            *
+  !     HSO, 1 July 2014                                                       *
+  !     Added optional namelist input                                          *
   !                                                                            *
   !*****************************************************************************
@@ -80 +82 @@
   character(len=50) :: line
   logical :: old
-  logical :: nml_COMMAND=.true. , nmlout=.true.  !.false.
   integer :: readerror
 
@@ -111 +112 @@
 
   ! Presetting namelist command
-  ldirect=1
+  ldirect=0
   ibdate=20000101
   ibtime=0
@@ -142 +143 @@
   ! Namelist input first: try to read as namelist file
   !**************************************************************************
-  open(unitcommand,file=path(1)(1:length(1))//'COMMAND',status='old', &
-         form='formatted',iostat=readerror)
-   ! If fail, check if file does not exist
-   if (readerror.ne.0) then
-     print*,'***ERROR: file COMMAND not found in ' 
-     print*, path(1)(1:length(1))//'COMMAND'
-     print*, 'Check your pathnames file.'
-     stop
-   endif
-
-   ! print error code
-   !write(*,*) 'readcommand > readerror open=' , readerror
-   !probe first line  
-   read (unitcommand,901) line
-   !write(*,*) 'index(line,COMMAND) =', index(line,'COMMAND') 
-
+  open(unitcommand,file=path(1)(1:length(1))//'COMMAND',status='old',form='formatted',err=999)
+
+  ! try namelist input (default)
+  read(unitcommand,command,iostat=readerror)
+  close(unitcommand)
+
+  ! distinguish namelist from fixed text input
+  if ((readerror.ne.0).or.(ldirect.eq.0)) then ! parse as text file format
  
-   !default is namelist input 
-   ! distinguish namelist from fixed text input
-   if (index(line,'COMMAND') .eq. 0) then
-   nml_COMMAND = .false. 
-   !write(*,*) 'COMMAND file does not contain the string COMMAND in the first line'      
-   endif
-   !write(*,*) 'readcommand > read as namelist? ' , nml_COMMAND
-   rewind(unitcommand)
-   read(unitcommand,command,iostat=readerror)
-
-  close(unitcommand)
-
-  !write(*,*) 'readcommand > readerror read=' , readerror
-  ! If error in namelist format, try to open with old input code
-  ! if (readerror.ne.0) then
-  ! IP 21/5/2014 the previous line cause the old long format 
-  ! to be confused with namelist input
-  
-  ! use text input 
-  if (nml_COMMAND .eqv. .false.) then 
-
-    open(unitcommand,file=path(1)(1:length(1))//'COMMAND',status='old', &
-         err=999)
+    open(unitcommand,file=path(1)(1:length(1))//'COMMAND',status='old', err=999)
 
     ! Check the format of the COMMAND file (either in free format,
@@ -193 +164 @@
     if (index(line,'LDIRECT') .eq. 0) then
       old = .false.
-    !write(*,*) 'readcommand old short'
+      write(*,*) 'COMMAND in old short format, please update to namelist format'
     else
       old = .true.
-    !write(*,*) 'readcommand old long'
+      write(*,*) 'COMMAND in old long format, please update to namelist format'
     endif
     rewind(unitcommand)
@@ -206 +177 @@
     call skplin(7,unitcommand)
     if (old) call skplin(1,unitcommand)
-
     read(unitcommand,*) ldirect
     if (old) call skplin(3,unitcommand)
@@ -265 +235 @@
     write(unitcommand,nml=command)
     close(unitcommand)
-     ! open(unitheader,file=path(2)(1:length(2))//'header_nml',status='new',err=999)
-     ! write(unitheader,NML=COMMAND)
-     !close(unitheader) 
   endif
 
@@ -616 +583 @@
 
 1000   write(*,*) ' #### FLEXPART MODEL ERROR! FILE "COMMAND"    #### '
-       write(*,*) ' #### CANNOT BE OPENED IN THE DIRECTORY       #### '
-        write(*,'(a)') path(2)(1:length(1))
-        stop
+  write(*,*) ' #### CANNOT BE OPENED IN THE DIRECTORY       #### '
+  write(*,'(a)') path(2)(1:length(1))
+  stop
 end subroutine readcommand