Changeset 4138764 in flexpart.git for src/timemanager.f90

Timestamp:

Apr 7, 2021, 8:45:57 AM (3 years ago)

Author:

Sabine <sabine.eckhardt@…>

Branches:

dev

Children:

1228ef7

Parents:

03adec6 (diff), 759df5f (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent.

Message:

Merge remote-tracking branch 'refs/remotes/origin/dev' into dev

File:

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

src/timemanager.f90

@@ +1 @@
+! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
+! SPDX-License-Identifier: GPL-3.0-or-later
+
 subroutine timemanager(metdata_format)
 
@@ -87 +90 @@
   implicit none
 
+  integer(selected_int_kind(16)) :: idummy,idummy2
   integer :: metdata_format
   integer :: j,ks,kp,l,n,itime=0,nstop,nstop1
 ! integer :: ksp
   integer :: loutnext,loutstart,loutend
-  integer :: ix,jy,ldeltat,itage,nage,idummy
+  integer :: ix,jy,ldeltat,itage,nage
   integer :: i_nan=0,ii_nan,total_nan_intl=0  !added by mc to check instability in CBL scheme 
   real :: outnum,weight,prob_rec(maxspec),prob(maxspec),decfact,wetscav
@@ -376 +380 @@
               call concoutput(itime,outnum,gridtotalunc,wetgridtotalunc,drygridtotalunc)
             endif
-          else  
+          else ! surf only  
             if (verbosity.eq.1) then
-             print*,'call concoutput_surf '
-             call system_clock(count_clock)
-             write(*,*) 'system clock',count_clock - count_clock0   
+              print*,'call concoutput_surf '
+              call system_clock(count_clock)
+              write(*,*) 'system clock',count_clock - count_clock0   
             endif
             if (lnetcdfout.eq.1) then
@@ -395 +399 @@
                 endif
               else
-              call concoutput_surf(itime,outnum,gridtotalunc,wetgridtotalunc,drygridtotalunc)
+                call concoutput_surf(itime,outnum,gridtotalunc,wetgridtotalunc,drygridtotalunc)
               endif
               if (verbosity.eq.1) then
@@ -413 +417 @@
                   call concoutput_inversion_nest(itime,outnum)
                 else 
-                call concoutput_surf_nest(itime,outnum)
-              endif
+                  call concoutput_surf_nest(itime,outnum)
+                endif
               endif
             else
@@ -581 +585 @@
        do ks=1,nspec
          if  ((xscav_frac1(j,ks).lt.0)) then
-            call get_wetscav(itime,lsynctime,loutnext,j,ks,grfraction,idummy,idummy,wetscav)
+            call get_wetscav(itime,lsynctime,loutnext,j,ks,grfraction,idummy,idummy2,wetscav)
             if (wetscav.gt.0) then
                 xscav_frac1(j,ks)=wetscav* &