-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathstats_lh_sfc_module.F90
137 lines (101 loc) · 4.74 KB
/
stats_lh_sfc_module.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
!-----------------------------------------------------------------------
! $Id$
!===============================================================================
module stats_lh_sfc_module
implicit none
private ! Set Default Scope
public :: stats_init_lh_sfc
! Constant parameters
integer, parameter, public :: nvarmax_lh_sfc = 10 ! Maximum variables allowed
contains
!-----------------------------------------------------------------------
subroutine stats_init_lh_sfc( vars_lh_sfc, & ! intent(in)
l_error, & ! intent(inout)
stats_metadata, stats_lh_sfc ) ! intent(inout)
! Description:
! Initializes array indices for stats_lh_sfc
! References:
! None
!-----------------------------------------------------------------------
use constants_clubb, only: &
fstderr ! Constant(s)
use stats_type_utilities, only: &
stat_assign ! Procedure
use stats_type, only: &
stats ! Type
use stats_variables, only: &
stats_metadata_type
implicit none
! External
intrinsic :: trim
!--------------------- Input Variable ---------------------
character(len= * ), dimension(nvarmax_lh_sfc), intent(in) :: vars_lh_sfc
!--------------------- InOut Variables ---------------------
type (stats_metadata_type), intent(inout) :: &
stats_metadata
type (stats), target, intent(inout) :: &
stats_lh_sfc
logical, intent(inout) :: l_error
!--------------------- Local Varables ---------------------
integer :: i, k
!--------------------- Begin Code ---------------------
! Default initialization for array indices for stats_sfc is zero (see module
! stats_variables)
! Assign pointers for statistics variables stats_sfc
k = 1
do i = 1, stats_lh_sfc%num_output_fields
select case ( trim( vars_lh_sfc(i) ) )
case ( 'lh_morr_snow_rate' )
stats_metadata%ilh_morr_snow_rate = k
call stat_assign( var_index=stats_metadata%ilh_morr_snow_rate, var_name="lh_morr_snow_rate", & !intent(in)
var_description="Snow+Ice+Graupel fallout rate from Morrison scheme [mm/day]", & ! In
var_units="mm/day", l_silhs=.true., & ! intent(in)
grid_kind=stats_lh_sfc ) ! intent(inout)
k = k + 1
case ( 'lh_vwp' )
stats_metadata%ilh_vwp = k
call stat_assign( var_index=stats_metadata%ilh_vwp, var_name="lh_vwp", & ! intent(in)
var_description="Vapor water path [kg/m^2]", var_units="kg/m^2", & ! intent(in)
l_silhs=.true., & ! intent(in)
grid_kind=stats_lh_sfc ) ! intent(inout)
k = k + 1
case ( 'lh_lwp' )
stats_metadata%ilh_lwp = k
call stat_assign( var_index=stats_metadata%ilh_lwp, var_name="lh_lwp", & ! intent(in)
var_description="Liquid water path [kg/m^2]", var_units="kg/m^2", & ! intent(in)
l_silhs=.true., & ! intent(in)
grid_kind=stats_lh_sfc ) ! intent(inout)
k = k + 1
case ( 'k_lh_start' )
stats_metadata%ik_lh_start = k
call stat_assign( var_index=stats_metadata%ik_lh_start, var_name="k_lh_start", & ! intent(in)
var_description="Index of height level for SILHS sampling preferentially within &
&cloud [integer]", var_units="integer", l_silhs=.true., &
grid_kind=stats_lh_sfc ) ! intent(inout)
k = k + 1
case ( 'lh_sample_weights_sum' )
stats_metadata%ilh_sample_weights_sum = k
call stat_assign( var_index=stats_metadata%ilh_sample_weights_sum, & ! intent(in)
var_name="lh_sample_weights_sum", & ! intent(in)
var_description="Sum of the sample point weights [-]", var_units="-", & ! intent(in)
l_silhs=.true., & ! intent(in)
grid_kind=stats_lh_sfc ) ! intent(inout)
k = k + 1
case ( 'lh_sample_weights_avg' )
stats_metadata%ilh_sample_weights_avg = k
call stat_assign( var_index=stats_metadata%ilh_sample_weights_avg, & ! intent(in)
var_name="lh_sample_weights_avg", & ! intent(in)
var_description="Average of the sample point weights [-]", & ! intent(in)
var_units="-", & ! intent(in)
l_silhs=.true., & ! intent(in)
grid_kind=stats_lh_sfc ) ! intent(inout)
k = k + 1
case default
write(fstderr,*) 'Error: unrecognized variable in vars_lh_sfc: ', &
trim( vars_lh_sfc(i) )
l_error = .true. ! This will stop the run.
end select
end do ! i = 1, stats_lh_sfc%num_output_fields
return
end subroutine stats_init_lh_sfc
end module stats_lh_sfc_module