-
Notifications
You must be signed in to change notification settings - Fork 7
Expand file tree
/
Copy pathdust_data_mod.F90
More file actions
executable file
·97 lines (80 loc) · 3.88 KB
/
dust_data_mod.F90
File metadata and controls
executable file
·97 lines (80 loc) · 3.88 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
!>\file dust_data_mod.F90
!! This file contains the data for the dust flux schemes.
module dust_data_mod
use mpas_kind_types
implicit none
integer, parameter :: ndust = 5
integer, parameter :: ndcls = 3
integer, parameter :: ndsrc = 1
integer, parameter :: maxstypes = 100
integer, parameter :: nsalt = 9
real(RKIND), parameter :: dyn_visc = 1.5E-5
! -- dust parameters
! never used: integer, dimension(ndust), parameter :: ipoint = (/ 3, 2, 2, 2, 2 /)
real(RKIND), dimension(ndust), parameter :: den_dust = (/ 2500., 2650., 2650., 2650., 2650. /)
real(RKIND), dimension(ndust), parameter :: reff_dust = (/ 0.73D-6, 1.4D-6, 2.4D-6, 4.5D-6, 8.0D-6 /)
real(RKIND), dimension(ndust), parameter :: frac_s = (/ 0.1, 0.25, 0.25, 0.25, 0.25 /)
real(RKIND), dimension(ndust), parameter :: lo_dust = (/ 0.1D-6, 1.0D-6, 1.8D-6, 3.0D-6, 6.0D-6 /)
real(RKIND), dimension(ndust), parameter :: up_dust = (/ 1.0D-6, 1.8D-6, 3.0D-6, 6.0D-6,10.0D-6 /)
! never used: real(RKIND), dimension(ndust, 12) :: ch_dust = 0.8e-09_RKIND
! -- default dust parameters
! -- AFWA & GOCART
! -----------+----------+-----------+
! Parameter | FIM-Chem | HRRR-Chem |
! -----------+----------+-----------+
! alpha | 1.0 | 0.5 |
! gamma | 1.6 | 1.0 |
! -----------+----------+-----------+
! Never used:
! real(RKIND), parameter :: afwa_alpha = 0.2
! real(RKIND), parameter :: afwa_gamma = 1.3
! real(RKIND), parameter :: gocart_alpha = 0.3
! real(RKIND), parameter :: gocart_gamma = 1.3
! -- FENGSHA
! Never used:
! real(RKIND), parameter :: fengsha_alpha = 0.3
! real(RKIND), parameter :: fengsha_gamma = 1.3
! -- FENGSHA threshold velocities based on Dale A. Gillette's data
integer, parameter :: fengsha_maxstypes = 13
real(RKIND), dimension(fengsha_maxstypes), parameter :: dust_uthres = &
(/ 0.065, & ! Sand - 1
0.18, & ! Loamy Sand - 2
0.27, & ! Sandy Loam - 3
0.30, & ! Silt Loam - 4
0.35, & ! Silt - 5
0.38, & ! Loam - 6
0.35, & ! Sandy Clay Loam - 7
0.41, & ! Silty Clay Loam - 8
0.41, & ! Clay Loam - 9
0.45, & ! Sandy Clay - 10
0.50, & ! Silty Clay - 11
0.45, & ! Clay - 12
9999.0 /) ! Other - 13
! -- FENGSHA uses precalculated drag partition
integer, parameter :: dust_calcdrag = 1
! -- FENGSHA dust moisture parameterization 1:fecan - 2:shao
integer :: dust_moist_opt = 1
! real(RKIND) :: dust_alpha = 10.0
! real(RKIND) :: dust_gamma = 1.3
! real(RKIND) :: dust_moist_correction = 0.5
! real(RKIND) :: dust_drylimit_factor = 0.25
! -- sea salt parameters
integer, dimension(nsalt), parameter :: spoint = (/ 1, 2, 2, 2, 2, 2, 3, 3, 3 /) ! 1 Clay, 2 Silt, 3 Sand
real(RKIND), dimension(nsalt), parameter :: reff_salt = &
(/ 0.71D-6, 1.37D-6, 2.63D-6, 5.00D-6, 9.50D-6, 18.1D-6, 34.5D-6, 65.5D-6, 125.D-6 /)
real(RKIND), dimension(nsalt), parameter :: den_salt = &
(/ 2500., 2650., 2650., 2650., 2650., 2650., 2650., 2650., 2650. /)
real(RKIND), dimension(nsalt), parameter :: frac_salt = &
(/ 1., 0.2, 0.2, 0.2, 0.2, 0.2, 0.333, 0.333, 0.333 /)
! -- soil vegatation parameters
integer, parameter :: max_soiltyp = 30
real(RKIND), dimension(max_soiltyp), parameter :: &
maxsmc = (/ 0.421, 0.464, 0.468, 0.434, 0.406, 0.465, &
0.404, 0.439, 0.421, 0.000, 0.000, 0.000, &
0.000, 0.000, 0.000, 0.000, 0.000, 0.000, &
0.000, 0.000, 0.000, 0.000, 0.000, 0.000, &
0.000, 0.000, 0.000, 0.000, 0.000, 0.000 /)
! -- other soil parameters
! never used: real(RKIND), dimension(maxstypes) :: porosity
public
end module dust_data_mod