Index | Rundeck: E6TomaF40 | Created: Fri May 17 03:30:39 EDT 2024 |
TRDUST_MOD | File: TRDUST_COM.f |
Author : Jan Perlwitz, Reha Cakmur, Ina Tegen
Version: 3.0
Subroutines:
Global Variables:
- by6 : used from constant
- im : used from resolution
- jm : used from resolution
- lm : used from resolution
- ntm_dust : used from tracer_com
- ntm_clay : used from tracer_com
- trname : used from oldtracer_mod
- max_len_name : used from oldtracer_mod
- by8 :
real(kind=8), parameter
- 0.25d0/2d0
Initial Value= 0.25D0/2D0
- i :
integer
- by4 :
real(kind=8), parameter
- 1D0/4D0
Initial Value= 1D0/4D0
- dust_names :
character(len=max_len_name), allocatable, dimension(:), dimension(:)
- names of soil dust aerosol tracers
- daridsoils :
real(kind=8), parameter
- median particle diameter by volume of arid soils
Initial Value= 3.4d0
- sigmaaridsoils :
real(kind=8), parameter
- particle diameter standard deviation by volume of
Initial Value= 3.0d0
- cv :
real(kind=8), parameter
- normalization constant for emitted volume size distribution
Initial Value= 12.63
- lambda :
real(kind=8), parameter
- crack propagation length for emitted volume size distribution
Initial Value= 12
- ndustbins :
integer, parameter
- number of size classes for soil dust aerosol tracers
Initial Value= 6
Used by: | :local_dust_emission | GHY_TRACERS:ghy_tracers_set_cell | SOCPBL |
- dustbounds :
real(kind=8), parameter, dimension(ndustbins+1)
- particle diameter bounds of soil dust aerosol bins [1 um]
Initial Value= (/ 0.1d0, 2.0d0, 4.0d0, 8.0d0, 16.0d0, 32.0d0, 64.0d0 /)
- nsubclays :
integer, parameter
- number of sub bins in clay size class of soil dust aerosols
Initial Value= 4
Used by: | :hetcdust | :init_ijlts_diag | :init_ijts_diag | :init_rad | :sulfdust | RAD_COM |
- subclaybounds :
real(kind=8), parameter, dimension(nsubclays+1)
- particle diameter bounds of clay sub bins [1 um]
Initial Value= (/ 0.1d0, 0.2d0, 0.5d0, 1.d0, 2.d0 /)
- ndustbinsfull :
integer, parameter
- number of soil dust bins for radiation
Initial Value= nSubClays+nDustBins-1
Used by: | :hetcdust | :sulfdust |
- dustboundsradia :
real(kind=8), parameter, dimension(ndustbinsfull+1)
- particle diameter bounds of soil dust aerosol bins
Initial Value= (/ subClayBounds, dustBounds( 3:nDustBins +1 ) /)
- subclayweights :
real(kind=8), dimension(nsubclays, ntm_clay)
- weights for masses in the sub bins of the clay size
Initial Value= reshape
Used by: | :hetcdust | :init_rad | :sulfdust |
- imdust :
integer
- : 0: scheme using PDF of wind speed (default)
Initial Value= 0
Used by: | :dust_emission_constraints | :hetcdust | :init_ijts_diag | :init_jls_diag | :init_rad | :local_dust_emission | :sulfdust | GHY_TRACERS |
- cwicub :
real(kind=8), parameter
- uplift factor [kg*s**2/m**5] for all size classes
Initial Value= 52.D-9
Used by: | :local_dust_emission |
- fclwicub :
real(kind=8)
- fraction [1] of uplifted clay for scheme using cubes of
Initial Value= By6
Used by: | :local_dust_emission |
- fsiwicub :
real(kind=8)
- fractions [1] of uplifted silt for scheme using cubes of
Initial Value= By8
Used by: | :local_dust_emission |
- hbaij :
real(kind=8), allocatable, dimension(:, :)
- accumulated precipitation - evaporation balance [kg/m^2]
Used by: | :tracer_ic | :tracerio | PBL_DRV:dust_emission_prep | PBL_DRV:save_dust_emission_vars |
- ricntd :
real(kind=8), allocatable, dimension(:, :)
- no. of hours with negative precipitation - evaporation balance [1]
Used by: | :tracer_ic | :tracerio | PBL_DRV:dust_emission_prep | PBL_DRV:save_dust_emission_vars |
- dryhr :
real(kind=8), allocatable, dimension(:, :)
- number of hours with evaporation-precipitation greater Zero
Used by: | GHY_TRACERS:ghy_tracers_set_cell |
- frclay :
real(kind=8), allocatable, dimension(:, :)
- fraction of clay
Used by: | GHY_TRACERS:ghy_tracers_set_cell |
- frsilt :
real(kind=8), allocatable, dimension(:, :)
- fraction of silt
Used by: | GHY_TRACERS:ghy_tracers_set_cell |
- vtrsh :
real(kind=8), allocatable, dimension(:, :)
- threshold wind speed above which dust emis. is allowed [m/s]
Used by: | GHY_TRACERS:ghy_tracers_set_cell |
- prelay :
real(kind=8), allocatable, dimension(:, :, :)
- distributed array with some prec info needed for simple wet
Used by: | :condse |
- cwipdf :
real(kind=8), parameter
- uplift factor [kg*s**2/m**5] for all size classes of soil dust
Initial Value= 12.068996D-9
Used by: | :local_dust_emission |
- scaledustemission :
real(kind=8)
- Initial Value
= 1.d0
Used by: | :local_dust_emission |
- fracclaypdfscheme :
real(kind=8)
- Initial Value
= 1.d0
Used by: | :local_dust_emission |
- fracsiltpdfscheme :
real(kind=8)
- Initial Value
= 1.d0
Used by: | :local_dust_emission |
- vegetationers :
integer
- Initial Value
= 1
- wsubtke_dust_on :
integer
- Initial Value
= 1
Used by: | :dust_emission_constraints |
- wsubwd_dust_on :
integer
- Initial Value
= 1
Used by: | :dust_emission_constraints |
- wsubwm_dust_on :
integer
- Initial Value
= 1
Used by: | :dust_emission_constraints |
- ers_data :
real(kind=8), allocatable, dimension(:, :, :)
- field of ERS data
Used by: | GHY_TRACERS:ghy_tracers_set_cell |
- prefdustsources :
integer
- Initial Value
= 1
- dustsourcefunction :
real(kind=8), allocatable, dimension(:, :)
- distribution of preferred dust sources
Used by: | GHY_TRACERS:ghy_tracers_set_cell |
- x3 :
real(kind=8), dimension(Lkm)
- indices of lock up table for emission
Used by: | :dust_emission_constraints |
- x2 :
real(kind=8), dimension(Ljm)
- indices of lock up table for emission
Used by: | :dust_emission_constraints |
- x1 :
real(kind=8), dimension(Lim)
- indices of lock up table for emission
Used by: | :dust_emission_constraints |
- lim :
integer, parameter
- Initial Value
= 234
Used by: | :dust_emission_constraints |
- ljm :
integer, parameter
- Initial Value
= 234
Used by: | :dust_emission_constraints |
- lkm :
integer, parameter
- Initial Value
= 22
Used by: | :dust_emission_constraints |
- table :
real(kind=8), dimension( Lim, Ljm, Lkm )
- array of lock up table for emission local to each grid box
Used by: | :dust_emission_constraints |
- wsubtke_com :
real(kind=8), allocatable, dimension(:, :)
- distributed array of subscale turbulent term
- wsubwd_com :
real(kind=8), allocatable, dimension(:, :)
- distributed array of subscale dry convective term
- wsubwm_com :
real(kind=8), allocatable, dimension(:, :)
- distributed array of subscale moist convective term
- naerocomdust :
integer, parameter
- Number of AEROCOM dust classes
Initial Value= 4
- d_dust :
real(kind=8), allocatable, dimension(:, :, :, :)
- prescribed daily dust emissions [kg] (e.g. AEROCOM)
Used by: | GHY_TRACERS:ghy_tracers_set_cell |
- nminerals :
integer, parameter
- number of different minerals
Initial Value= 8
- mineralnames :
character(len=4), parameter, dimension(nminerals)
- names of Minerals
Initial Value= (/'Illi', 'Kaol', 'Smec', 'Calc', 'Quar', 'Feld', 'Feox', 'Gyps' /)
- dustbinnamesin :
character(len=5), parameter, dimension(ndustbins)
- dust bin names in mineral fraction input file
Initial Value= (/ 'Clay ', 'Silt1', 'Silt2', 'Silt3', 'Silt4', 'Silt5' /)
- mineralfractions :
real(kind=8), allocatable, dimension(:, :, :)
- distribution of mineral fractions in soils for each
Used by: | GHY_TRACERS:ghy_tracers_set_cell |
- calceffectiveradius :
integer
- flag whether to calculate or prescribe
Initial Value= 0
- radiusclay :
real(kind=8), parameter, dimension(nsubclays)
- dry radius of clay dust/minerals [um]
Initial Value= (/ 0.132d0, 0.23d0, 0.416d0, 0.766d0 /)
- radiussil1 :
real(kind=8), parameter
- dry radius of silt1 dust/minerals [um]
Initial Value= 1.386d0
- radiussil2 :
real(kind=8), parameter
- dry radius of silt2 dust/minerals [um]
Initial Value= 2.773d0
- radiussil3 :
real(kind=8), parameter
- dry radius of silt3 dust/minerals [um]
Initial Value= 5.545d0
- radiussil4 :
real(kind=8), parameter
- dry radius of silt4 dust/minerals [um]
Initial Value= 11.090d0
- radiussil5 :
real(kind=8), parameter
- dry radius of silt5 dust/minerals [um]
Initial Value= 22.0d0
- radiusminerals :
real(kind=8), parameter, dimension(ndustbinsfull)
- dry radius of all dust/minerals [um]
Initial Value= (/ radiusClay, radiusSil1, radiusSil2, radiusSil3, radiusSil4, radiusSil5 /)
Used by: | :hetcdust | :sulfdust |
- dryeffradminerals :
real(kind=8), parameter, dimension(ndustbinsfull)
- dry effective radius of minerals for radiation [um]
Initial Value= (/ radiusMinerals /)
Used by: | :init_rad |
- ndustemij :
integer, parameter
- index of dust emission in ijts_isrc
Initial Value= 1
Used by: | :init_ijts_diag | GHY_TRACERS |
- ndustem2ij :
integer, parameter
- index of dust emission according to cubic scheme
Initial Value= 2
Used by: | :init_ijts_diag | GHY_TRACERS |
- ndustturbij :
integer, parameter
- index of dust dry turbulent deposition in ijts_isrc
Initial Value= 3
- ndustev1ij :
integer, parameter
- index of number of dust events below threshold wind
Initial Value= 1
Used by: | :init_ijts_diag | GHY_TRACERS |
- ndustev2ij :
integer, parameter
- index of number of dust events above threshold wind
Initial Value= 2
Used by: | :init_ijts_diag | GHY_TRACERS |
- ndustwthij :
integer, parameter
- index of threshold velocity in ijts_spec
Initial Value= 3
Used by: | :init_ijts_diag | GHY_TRACERS |
- to_conc_soildust :
integer
- : For printout of 3D soil dust aerosol concentration
Initial Value= 0
- ndustemjl :
integer, parameter
- index of dust emission in jls_source
Initial Value= 1
Used by: | :init_jls_diag | GHY_TRACERS |
- ndustem2jl :
integer, parameter
- index of dust emission according to cubic scheme
Initial Value= 2
Used by: | :init_jls_diag | GHY_TRACERS |
- ndustturbjl :
integer, parameter
- index of dust dry turbulent deposition in jls_source
Initial Value= 3
- ndustev1jl :
integer, parameter
- index of number of dust events below threshold wind
Initial Value= 1
Used by: | :init_jls_diag | GHY_TRACERS |
- ndustev2jl :
integer, parameter
- index of number of dust events above threshold wind
Initial Value= 2
Used by: | :init_jls_diag | GHY_TRACERS |
- ndustwthjl :
integer, parameter
- index of threshold velocity in ijts_spec
Initial Value= 3
Used by: | :init_jls_diag | GHY_TRACERS |
- dustdiagsubdd :
- structured type for specific subdaily dust
- dustdiagsubdd_acc :
type(dustdiagsubdd)
- structured variable to accumulate data for