Changes in src/timemanager.f90 [2eefa58:c7d1052] in flexpart.git

File:

• src/timemanager.f90 (modified) (5 diffs)

src/timemanager.f90

@@ -408 +408 @@
 #endif
             else
-              if (linversionout.eq.1) then
-                call concoutput_inversion(itime,outnum,gridtotalunc,wetgridtotalunc,drygridtotalunc)
-                if (verbosity.eq.1) then
-                  print*,'called concoutput_inversion'
-                  call system_clock(count_clock)
-                  write(*,*) 'system clock',count_clock - count_clock0 
-                endif
-              else
               call concoutput_surf(itime,outnum,gridtotalunc,wetgridtotalunc,drygridtotalunc)
-              endif
               if (verbosity.eq.1) then
                 print*,'called concoutput_surf '
@@ -431 +422 @@
                 call concoutput_nest(itime,outnum)
               else 
-                if(linversionout.eq.1) then
-                  call concoutput_inversion_nest(itime,outnum)
-                else 
                 call concoutput_surf_nest(itime,outnum)
-              endif
               endif
             else
@@ -464 +451 @@
 45      format(i13,' Seconds simulated: ',i13, ' Particles:    Uncertainty: ',3f7.3)
 46      format(' Simulated ',f7.1,' hours (',i13,' s), ',i13, ' particles')
-        if (ipout.ge.1) then
-          if (mod(itime,ipoutfac*loutstep).eq.0) call partoutput(itime) ! dump particle positions
-          if (ipout.eq.3) call partoutput_average(itime) ! dump particle positions
-        endif
+        if (ipout.ge.1) call partoutput(itime)    ! dump particle positions
         loutnext=loutnext+loutstep
         loutstart=loutnext-loutaver/2
@@ -625 +609 @@
 !        write (*,*) 'advance: ',prob(1),xmass1(j,1),ztra1(j)
 
-  ! Calculate average position for particle dump output
-  !****************************************************
-
-        if (ipout.eq.3) call partpos_average(itime,j)
-
-
   ! Calculate the gross fluxes across layer interfaces
   !***************************************************
@@ -752 +730 @@
   if (ipout.eq.2) call partoutput(itime)     ! dump particle positions
 
-  if (linit_cond.ge.1) then
-    if(linversionout.eq.1) then
-      call initial_cond_output_inversion(itime)   ! dump initial cond. field
-    else
-      call initial_cond_output(itime)   ! dump initial cond. fielf
-    endif
-  endif
+  if (linit_cond.ge.1) call initial_cond_output(itime)   ! dump initial cond. field
 
   !close(104)