Changeset 4380 for palm/trunk/SOURCE


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

salsa_mod: minor bug in dry deposition and formatting

File:
1 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
Note: See TracChangeset for help on using the changeset viewer.