Changeset 31674b5 in flexpart.git for src/concoutput.f90

Timestamp:

Jul 15, 2015, 9:45:11 AM (9 years ago)

Author:

Ignacio Pisso <Ignacio.Pisso@…>

Branches:

svn-petra

Children:

3c1da52

Parents:

480164e

git-author:

Ignacio Pisso <Ignacio.Pisso@…> (07/15/15 09:19:44)

git-committer:

Ignacio Pisso <Ignacio.Pisso@…> (07/15/15 09:45:11)

Message:

from svn branch petra
r36 | pesei | 2015-02-16 19:21:09 +0100 (Mon, 16 Feb 2015) | 1 line
Changed paths:

A /branches/petra/changelog.txt
M /branches/petra/options/COMMAND
M /branches/petra/options/COMMAND.alternative
M /branches/petra/options/COMMAND.reference
M /branches/petra/options/surfdata.t
M /branches/petra/src/FLEXPART.f90
M /branches/petra/src/advance.f90
M /branches/petra/src/com_mod.f90
M /branches/petra/src/concoutput.f90
M /branches/petra/src/concoutput_nest.f90
M /branches/petra/src/concoutput_surf.f90
M /branches/petra/src/par_mod.f90
M /branches/petra/src/readcommand.f90
M /branches/petra/src/readpaths.f90
M /branches/petra/src/readreceptors.f90
M /branches/petra/src/readreleases.f90
M /branches/petra/src/readspecies.f90
M /branches/petra/src/readwind.f90
M /branches/petra/src/readwind_nests.f90
M /branches/petra/src/skplin.f90
M /branches/petra/src/writeheader.f90
M /branches/petra/src/writeheader_nest.f90
M /branches/petra/src/writeheader_txt.f90

Implement switch for incremental deposition, see ticket:113 and many small changes, see changelog.txt

File:

• src/concoutput.f90 (modified) (4 diffs)

src/concoutput.f90

@@ -1 +1 @@
 !**********************************************************************
-! Copyright 1998,1999,2000,2001,2002,2005,2007,2008,2009,2010         *
+! Copyright 1998-2015                                                 *
 ! Andreas Stohl, Petra Seibert, A. Frank, Gerhard Wotawa,             *
 ! Caroline Forster, Sabine Eckhardt, John Burkhart, Harald Sodemann   *
@@ -47 +47 @@
   !                                                                            *
   !     2008 new sparse matrix format                                          *
+  !     PS, 2/2015: option to produce incremental deposition output 
+  !                 access= -> position=
   !                                                                            *
   !*****************************************************************************
@@ -108 +110 @@
   write(adate,'(i8.8)') jjjjmmdd
   write(atime,'(i6.6)') ihmmss
-  open(unitdates,file=path(2)(1:length(2))//'dates', ACCESS='APPEND')
+  open(unitdates,file=path(2)(1:length(2))//'dates', position='append')
   write(unitdates,'(a)') adate//atime
   close(unitdates)  
@@ -588 +590 @@
   !*************************
 
-  do ks=1,nspec
-  do kp=1,maxpointspec_act
-    do i=1,numreceptor
-      creceptor(i,ks)=0.
-    end do
-    do jy=0,numygrid-1
-      do ix=0,numxgrid-1
-        do l=1,nclassunc
-          do nage=1,nageclass
-            do kz=1,numzgrid
-              gridunc(ix,jy,kz,ks,kp,l,nage)=0.
-            end do
-          end do
-        end do
-      end do
-    end do
-  end do
-  end do
-
-
+  creceptor=0.
+  gridunc=0.
+  if (ldep_incr) then ! incremental deposition output
+    wetgridunc=0.
+    drygridunc=0.
+  endif
+  
 end subroutine concoutput