Changeset 10 for branches/flexpart91_hasod/src_parallel/readageclasses.f90

Timestamp:

Aug 15, 2013, 3:23:48 PM (11 years ago)

Author:

hasod

Message:

ADD: namelist input implemented for all common input files

File:

• branches/flexpart91_hasod/src_parallel/readageclasses.f90 (modified) (6 diffs)

branches/flexpart91_hasod/src_parallel/readageclasses.f90

@@ -28 +28 @@
   !                                                                            *
   !     Author: A. Stohl                                                       *
-  !                                                                            *
   !     20 March 2000                                                          *
+  !     HSO, 14 August 2013
+  !     Added optional namelist input
   !                                                                            *
   !*****************************************************************************
@@ -46 +47 @@
   integer :: i
 
+  ! namelist help variables
+  integer :: readerror
+
+  ! namelist declaration
+  namelist /ageclass/ &
+    nageclass, &
+    lage
+
+  nageclass=-1 ! preset to negative value to identify failed namelist input
 
   ! If age spectra calculation is switched off, set number of age classes
   ! to 1 and maximum age to a large number
   !**********************************************************************
-
 
   if (lagespectra.ne.1) then
@@ -57 +66 @@
     return
   endif
-
 
   ! If age spectra claculation is switched on,
@@ -66 +74 @@
        status='old',err=999)
 
-  do i=1,13
-    read(unitageclasses,*)
-  end do
-  read(unitageclasses,*) nageclass
+  ! try to read in as a namelist
+  read(unitageclasses,ageclass,iostat=readerror)
 
+  if ((nageclass.lt.0).or.(readerror.ne.0)) then
+    rewind(unitageclasses)
+    do i=1,13
+      read(unitageclasses,*)
+    end do
+    read(unitageclasses,*) nageclass
+    read(unitageclasses,*) lage(1)
+    do i=2,nageclass
+      read(unitageclasses,*) lage(i)
+    end do
+  endif
+
+  close(unitageclasses)
 
   if (nageclass.gt.maxageclass) then
@@ -81 +100 @@
   endif
 
-  read(unitageclasses,*) lage(1)
   if (lage(1).le.0) then
     write(*,*) ' #### FLEXPART MODEL ERROR! AGE OF FIRST      #### '
@@ -90 +108 @@
 
   do i=2,nageclass
-    read(unitageclasses,*) lage(i)
     if (lage(i).le.lage(i-1)) then
       write(*,*) ' #### FLEXPART MODEL ERROR! AGE CLASSES     #### '