Changeset 4380 for palm/trunk/SOURCE/salsa_mod.f90

Timestamp:

Jan 17, 2020 11:39:51 PM (5 years ago)

Author:

monakurppa

Message:

salsa_mod: minor bug in dry deposition and formatting

File:

• palm/trunk/SOURCE/salsa_mod.f90 (modified) (6 diffs)

palm/trunk/SOURCE/salsa_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! - Error in saving the surface fluxes in an array that is applied in the
+!   deposition scheme
+! - Corrections in the header: aerosol bin diameters and lower bin limits not
+!   printed correctly
+! 
+! 4364 2020-01-08 02:12:31Z monakurppa
 ! Set time coordinate in the input data relative to origin_time rather than to
 ! 00:00:00 UTC
@@ -1281 +1287 @@
     IF ( lspartition )  WRITE( io, 11 )
     IF ( nldepo )       WRITE( io, 12 ) nldepo_pcm, nldepo_surf
-    WRITE( io, 13 )  reglim, nbin, bin_low_limits
+    WRITE( io, 13 )  reglim, nbin, ( aero(:)%vlolim / api6 )**0.33333333_wp
+    WRITE( io, 25 )  aero(:)%dmid
     IF ( init_aerosol_type == 0 )  WRITE( io, 14 ) nsect
     WRITE( io, 15 ) ncc, listspec, mass_fracs_a, mass_fracs_b
@@ -1310 +1317 @@
 3   FORMAT  (/'    Timestep: dt_salsa = ', F6.2, '  s')
 4   FORMAT  (/'    Array shape (z,y,x,bins):'/                                                     &
-              '       aerosol_number:  ', 4(I3)) 
-5   FORMAT  (/'       aerosol_mass:    ', 4(I3),/                                                  &
+              '       aerosol_number:  ', 4(I5)) 
+5   FORMAT  (/'       aerosol_mass:    ', 4(I5),/                                                  &
               '       (advect_particle_water = ', L1, ')')
-6   FORMAT   ('       salsa_gas: ', 4(I3),/                                                        &
+6   FORMAT   ('       salsa_gas: ', 4(I5),/                                                        &
               '       (salsa_gases_from_chem = ', L1, ')')
 7   FORMAT  (/'    Aerosol dynamic processes included: ')
@@ -1323 +1330 @@
 13  FORMAT  (/'    Aerosol bin subrange limits (in metres): ',  3(ES10.2E3), /                     &
               '    Number of size bins for each aerosol subrange: ', 2I3,/                         &
-              '    Aerosol bin limits (in metres): ', 9(ES10.2E3))
+              '    Aerosol bin lower limits (in metres): ', 12(ES10.2E3))
+25  FORMAT  (/'    Bin geometric mean diameters (in metres): ', 12(ES10.2E3))
 14  FORMAT   ('    Initial number concentration in bins at the lowest level (#/m**3):', 9(ES10.2E3))
 15  FORMAT  (/'    Number of chemical components used: ', I1,/                                     &
@@ -9331 +9339 @@
 !--          Save number fluxes in the end
              surface%answs(m,ib) = surface%answs(m,ib) + source_array(j,i,ib) * rho_air_zw(k-1)
-             aerosol_number(ib)%source(j,i) = aerosol_number(ib)%source(j,i) + surface%answs(m,ib)
+             aerosol_number(ib)%source(j,i) = surface%answs(m,ib)
 
           ENDDO  ! ib
@@ -9377 +9385 @@
        surface%amsws(surf_num,ic) = surface%amsws(surf_num,ic) + mass_frac * nsource *             &
                                     aero(ib)%core * prho * rho_air_zw(k-1)
-       aerosol_mass(ic)%source(j,i) = aerosol_mass(ic)%source(j,i) + surface%amsws(surf_num,ic)
+       aerosol_mass(ic)%source(j,i) = surface%amsws(surf_num,ic)
 
     END SUBROUTINE set_mass_flux