Changeset 7d02c2f in flexpart.git

Timestamp:

Mar 1, 2018, 8:42:11 AM (6 years ago)

Author:

Espen Sollum ATMOS <eso@…>

Branches:

master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie

Children:

91ce5e7

Parents:

01e6052

Message:

Added BW scavenging for nested domains

File:

• src/concoutput_nest.f90 (modified) (1 diff)

src/concoutput_nest.f90

@@ -199 +199 @@
 
   write(anspec,'(i3.3)') ks
-  if ((iout.eq.1).or.(iout.eq.3).or.(iout.eq.5)) then
-    if (ldirect.eq.1) then
-      open(unitoutgrid,file=path(2)(1:length(2))//'grid_conc_nest_' &
+ 
+  if (DRYBKDEP.or.WETBKDEP) then !scavdep output
+      if (DRYBKDEP) &
+      open(unitoutgrid,file=path(2)(1:length(2))//'grid_drydep_nest_'//adate// &
+           atime//'_'//anspec,form='unformatted')
+      if (WETBKDEP) &
+      open(unitoutgrid,file=path(2)(1:length(2))//'grid_wetdep_nest_'//adate// &
+           atime//'_'//anspec,form='unformatted')
+      write(unitoutgrid) itime
+  else
+     if ((iout.eq.1).or.(iout.eq.3).or.(iout.eq.5)) then
+       if (ldirect.eq.1) then
+         open(unitoutgrid,file=path(2)(1:length(2))//'grid_conc_nest_' &
            //adate// &
            atime//'_'//anspec,form='unformatted')
-    else
-      open(unitoutgrid,file=path(2)(1:length(2))//'grid_time_nest_' &
+       else
+         open(unitoutgrid,file=path(2)(1:length(2))//'grid_time_nest_' &
            //adate// &
            atime//'_'//anspec,form='unformatted')
-    endif
-     write(unitoutgrid) itime
-   endif
+       endif
+       write(unitoutgrid) itime
+     endif
+  endif
 
   if ((iout.eq.2).or.(iout.eq.3)) then      ! mixing ratio