Changeset 660bcb7 in flexpart.git for src/readspecies.f90

Timestamp:

Nov 24, 2014, 1:36:35 PM (9 years ago)

Author:

Anne Fouilloux <annefou@…>

Branches:

NetCDF

Children:

3073eaf

Parents:

4bf4a69

Message:

NETCDF outputs from svn trunk (hasod: ADD: netcdf module file)
I did not take changes in advance.f90; it will be added later.
I also changed OPENs where status was set to new and access is set to APPEND (had problems on abel.uio.no with intel compilers 2013.sp1.3)
It should contains technical changes only and results should be identical.

File:

• src/readspecies.f90 (modified) (2 diffs)

src/readspecies.f90

@@ -239 +239 @@
 20   continue
 
+
+  ! Read in daily and day-of-week variation of emissions, if available
+  !*******************************************************************
+  ! HSO: This is not yet implemented as namelist parameters
+  ! default values set to 1
+
+  do j=1,24           ! initialize everything to no variation
+    area_hour(i,j)=1.
+    point_hour(i,j)=1.
+  end do
+  do j=1,7
+    area_dow(i,j)=1.
+    point_dow(i,j)=1.
+  end do
+
   if (readerror.ne.0) then ! text format input
-
-    ! Read in daily and day-of-week variation of emissions, if available
-    !*******************************************************************
-    ! HSO: This is not yet implemented as namelist parameters
-
-    do j=1,24           ! initialize everything to no variation
-      area_hour(i,j)=1.
-      point_hour(i,j)=1.
-    end do
-    do j=1,7
-      area_dow(i,j)=1.
-      point_dow(i,j)=1.
-    end do
 
     read(unitspecies,*,end=22)
@@ -269 +271 @@
   ! namelist output if requested
   if (nmlout.eqv..true.) then
-    open(unitspecies,file=path(2)(1:length(2))//'SPECIES_'//aspecnumb//'.namelist',access='append',status='new',err=1000)
+    open(unitspecies,file=path(2)(1:length(2))//'SPECIES_'//aspecnumb//'.namelist',access='append',status='unknown',err=1000)
     write(unitspecies,nml=species_params)
     close(unitspecies)