-
Notifications
You must be signed in to change notification settings - Fork 12
/
parameters_model.F90
154 lines (112 loc) · 4.53 KB
/
parameters_model.F90
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
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
!-------------------------------------------------------------------------------
! $Id$
!===============================================================================
module parameters_model
! Description:
! Contains model parameters that are determined at run time rather than
! compile time.
!
! References:
! None
!-------------------------------------------------------------------------------
use clubb_precision, only: &
core_rknd
implicit none
private ! Default scope
integer, parameter :: &
sp = selected_real_kind(6) ! 32-bit floating point number
! Maximum magnitude of PDF parameter 'mixt_frac'.
real( kind = core_rknd ), public :: mixt_frac_max_mag
!$omp threadprivate(mixt_frac_max_mag)
! Model parameters and constraints setup in the namelists
real( kind = core_rknd ), public :: &
T0 = 300._core_rknd, & ! Reference temperature (usually 300) [K]
ts_nudge = 0._core_rknd ! Timescale of u/v nudging [s]
#ifdef GFDL
real( kind = core_rknd ), public :: & ! h1g, 2010-06-15
cloud_frac_min ! minimum cloud fraction for droplet #
!$omp threadprivate( cloud_frac_min )
#endif
!$omp threadprivate(T0, ts_nudge)
real( kind = core_rknd), public :: &
rtm_min = epsilon( rtm_min ), & ! Value below which rtm will be nudged [kg/kg]
rtm_nudge_max_altitude = 10000._core_rknd ! Highest altitude at which to nudge rtm [m]
!$omp threadprivate( rtm_min, rtm_nudge_max_altitude )
integer, public :: &
sclr_dim = 0, & ! Number of passive scalars
edsclr_dim = 0, & ! Number of eddy-diff. passive scalars
hydromet_dim = 0 ! Number of hydrometeor species
!$omp threadprivate(sclr_dim, edsclr_dim, hydromet_dim)
real( kind = core_rknd ), dimension(:), allocatable, public :: &
sclr_tol ! Threshold(s) on the passive scalars [units vary]
!$omp threadprivate(sclr_tol)
real( kind = sp ), public :: PosInf
!$omp threadprivate(PosInf)
public :: setup_parameters_model
contains
!-------------------------------------------------------------------------------
subroutine setup_parameters_model &
( T0_in, ts_nudge_in, &
hydromet_dim_in, &
sclr_dim_in, sclr_tol_in, edsclr_dim_in &
#ifdef GFDL
, cloud_frac_min_in & ! hlg, 2010-6-15
#endif
)
! Description:
! Sets parameters to their initial values
!
! References:
! None
!-------------------------------------------------------------------------------
use parameters_tunable, only: &
Skw_max_mag
use clubb_precision, only: &
core_rknd ! Variable(s)
implicit none
! External
intrinsic :: sqrt, allocated, transfer
! Constants
integer(kind=4), parameter :: nanbits = 2139095040
! Input Variables
real( kind = core_rknd ), intent(in) :: &
T0_in, & ! Ref. temperature [K]
ts_nudge_in ! Timescale for u/v nudging [s]
#ifdef GFDL
real( kind = core_rknd ), intent(in) :: cloud_frac_min_in ! h1g, 2010-06-15
#endif
integer, intent(in) :: &
hydromet_dim_in, & ! Number of hydrometeor species
sclr_dim_in, & ! Number of passive scalars
edsclr_dim_in ! Number of eddy-diff. passive scalars
real( kind = core_rknd ), intent(in), dimension(sclr_dim_in) :: &
sclr_tol_in ! Threshold on passive scalars
! --- Begin Code ---
! Formula from subroutine pdf_closure, where sigma_sqd_w = 0.4 and Skw =
! Skw_max_mag in this formula. Note that this is constant, but can't appear
! with a Fortran parameter attribute, so we define it here.
mixt_frac_max_mag = 1.0_core_rknd &
- ( 0.5_core_rknd * ( 1.0_core_rknd - Skw_max_mag / &
sqrt( 4.0_core_rknd * ( 1.0_core_rknd - 0.4_core_rknd )**3 &
+ Skw_max_mag**2 ) ) ) ! Known magic number
T0 = T0_in
ts_nudge = ts_nudge_in
hydromet_dim = hydromet_dim_in
sclr_dim = sclr_dim_in
edsclr_dim = edsclr_dim_in
! In a tuning run, this array has the potential to be allocated already
if ( .not. allocated( sclr_tol ) ) then
allocate( sclr_tol(1:sclr_dim) )
else
deallocate( sclr_tol )
allocate( sclr_tol(1:sclr_dim) )
end if
sclr_tol(1:sclr_dim) = sclr_tol_in(1:sclr_dim)
PosInf = transfer( nanbits, PosInf )
#ifdef GFDL
cloud_frac_min = cloud_frac_min_in ! h1g, 2010-06-15
#endif
return
end subroutine setup_parameters_model
!-------------------------------------------------------------------------------
end module parameters_model