Changeset 4380


Ignore:
Timestamp:
Jan 17, 2020 11:39:51 PM (17 months ago)
Author:
monakurppa
Message:

salsa_mod: minor bug in dry deposition and formatting

Location:
palm/trunk
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/salsa_mod.f90

    r4364 r4380  
    2626! -----------------
    2727! $Id$
     28! - Error in saving the surface fluxes in an array that is applied in the
     29!   deposition scheme
     30! - Corrections in the header: aerosol bin diameters and lower bin limits not
     31!   printed correctly
     32!
     33! 4364 2020-01-08 02:12:31Z monakurppa
    2834! Set time coordinate in the input data relative to origin_time rather than to
    2935! 00:00:00 UTC
     
    12811287    IF ( lspartition )  WRITE( io, 11 )
    12821288    IF ( nldepo )       WRITE( io, 12 ) nldepo_pcm, nldepo_surf
    1283     WRITE( io, 13 )  reglim, nbin, bin_low_limits
     1289    WRITE( io, 13 )  reglim, nbin, ( aero(:)%vlolim / api6 )**0.33333333_wp
     1290    WRITE( io, 25 )  aero(:)%dmid
    12841291    IF ( init_aerosol_type == 0 )  WRITE( io, 14 ) nsect
    12851292    WRITE( io, 15 ) ncc, listspec, mass_fracs_a, mass_fracs_b
     
    131013173   FORMAT  (/'    Timestep: dt_salsa = ', F6.2, '  s')
    131113184   FORMAT  (/'    Array shape (z,y,x,bins):'/                                                     &
    1312               '       aerosol_number:  ', 4(I3))
    1313 5   FORMAT  (/'       aerosol_mass:    ', 4(I3),/                                                  &
     1319              '       aerosol_number:  ', 4(I5))
     13205   FORMAT  (/'       aerosol_mass:    ', 4(I5),/                                                  &
    13141321              '       (advect_particle_water = ', L1, ')')
    1315 6   FORMAT   ('       salsa_gas: ', 4(I3),/                                                        &
     13226   FORMAT   ('       salsa_gas: ', 4(I5),/                                                        &
    13161323              '       (salsa_gases_from_chem = ', L1, ')')
    131713247   FORMAT  (/'    Aerosol dynamic processes included: ')
     
    1323133013  FORMAT  (/'    Aerosol bin subrange limits (in metres): ',  3(ES10.2E3), /                     &
    13241331              '    Number of size bins for each aerosol subrange: ', 2I3,/                         &
    1325               '    Aerosol bin limits (in metres): ', 9(ES10.2E3))
     1332              '    Aerosol bin lower limits (in metres): ', 12(ES10.2E3))
     133325  FORMAT  (/'    Bin geometric mean diameters (in metres): ', 12(ES10.2E3))
    1326133414  FORMAT   ('    Initial number concentration in bins at the lowest level (#/m**3):', 9(ES10.2E3))
    1327133515  FORMAT  (/'    Number of chemical components used: ', I1,/                                     &
     
    93319339!--          Save number fluxes in the end
    93329340             surface%answs(m,ib) = surface%answs(m,ib) + source_array(j,i,ib) * rho_air_zw(k-1)
    9333              aerosol_number(ib)%source(j,i) = aerosol_number(ib)%source(j,i) + surface%answs(m,ib)
     9341             aerosol_number(ib)%source(j,i) = surface%answs(m,ib)
    93349342
    93359343          ENDDO  ! ib
     
    93779385       surface%amsws(surf_num,ic) = surface%amsws(surf_num,ic) + mass_frac * nsource *             &
    93789386                                    aero(ib)%core * prho * rho_air_zw(k-1)
    9379        aerosol_mass(ic)%source(j,i) = aerosol_mass(ic)%source(j,i) + surface%amsws(surf_num,ic)
     9387       aerosol_mass(ic)%source(j,i) = surface%amsws(surf_num,ic)
    93809388
    93819389    END SUBROUTINE set_mass_flux
  • palm/trunk/TESTS/cases/urban_environment_salsa/MONITORING/urban_environment_salsa_rc

    r4364 r4380  
    11
    22 ******************************    --------------------------------------------
    3  * PALM 6.0  Rev: 4363M       *    atmosphere - run without 1D - prerun
     3 * PALM 6.0  Rev: 4379M       *    atmosphere - run without 1D - prerun
    44 ******************************    --------------------------------------------
    55
    6  Date:               2020-01-07    Run:       urban_environment_salsa           
    7  Time:                 21:00:12    Run-No.:   00
     6 Date:               2020-01-17    Run:       urban_environment_salsa           
     7 Time:                 18:25:39    Run-No.:   00
    88 Run on host:           default
    99 Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
     
    259259
    260260    Array shape (z,y,x,bins):
    261        aerosol_number:   60 19 19  8
    262 
    263        aerosol_mass:     60 19 19 32
     261       aerosol_number:     60   19   19    8
     262
     263       aerosol_mass:       60   19   19  32
    264264       (advect_particle_water = T)
    265        salsa_gas:  60 19 19  5
     265       salsa_gas:    60   19   19    5
    266266       (salsa_gases_from_chem = F)
    267267
     
    272272    Aerosol bin subrange limits (in metres):  3.00E-009 1.00E-008 2.50E-006
    273273    Number of size bins for each aerosol subrange:   1  7
    274     Aerosol bin limits (in metres):  0.00E+000 0.00E+000 0.00E+000 0.00E+000 0.00E+000 0.00E+000 0.00E+000 0.00E+000
     274    Aerosol bin lower limits (in metres):  3.00E-009 1.00E-008 2.20E-008 4.84E-008 1.07E-007 2.35E-007 5.16E-007 1.14E-006
     275
     276    Bin geometric mean diameters (in metres):  5.48E-009 1.48E-008 3.26E-008 7.18E-008 1.58E-007 3.48E-007 7.66E-007 1.69E-006
    275277
    276278    Number of chemical components used: 3
Note: See TracChangeset for help on using the changeset viewer.