source: LMDZ5/branches/Cold_pool_death/libf/phylmd/Cold_pool_death/derivs_param.F90 @ 3899

Last change on this file since 3899 was 3899, checked in by lguez, 3 years ago

Old work on cold pool death.

File size: 2.0 KB
Line 
1module derivs_param
2
3  implicit none
4
5  real eps_cp_m ! entrainment coefficient at the border of the cold pool
6  !$omp THREADPRIVATE(eps_cp)
7
8  real beta_m ! entrainment at the top of the mixed layer
9  !$omp THREADPRIVATE(beta)
10
11  logical av_thvu_constant_m
12  !$omp THREADPRIVATE(av_thvu_constant)
13
14  real k_star_m
15  !$omp THREADPRIVATE(k_star)
16
17  real theta_vs_m ! surface virtual potential temperature, in K
18  !$omp THREADPRIVATE(theta_vs_m)
19
20  real rate_uns_m ! mass rate in one unsaturated downdraft, in kg s-1
21  !$omp THREADPRIVATE(rate_uns)
22
23  real av_rho_m ! average mass density in the cold pool, in kg m-3
24  !$omp THREADPRIVATE(av_rho_m)
25
26  real theta_v_uns_m
27  ! virtual potential temperature of the unsaturated downdraft, in K
28  !$omp THREADPRIVATE(theta_v_uns_m)
29
30  real C_dw_m ! drag coefficient under cold pools, in m s-1
31  !$omp THREADPRIVATE(C_dw_m)
32
33contains
34
35  subroutine set_module_var(eps_cp, beta, av_thvu_constant, k_star, theta_vs, &
36       rate_uns, av_rho, theta_v_uns, c_dw)
37
38    real, intent(in):: eps_cp
39    ! entrainment coefficient at the border of the cold pool
40
41    real, intent(in):: beta ! entrainment at the top of the mixed layer
42    logical, intent(in):: av_thvu_constant
43    real, intent(in):: k_star
44
45    real, intent(in):: theta_vs
46    ! virtual potential temperature at the surface, in K
47
48    real, intent(in):: rate_uns
49    ! mass rate in one unsaturated downdraft, in kg s-1
50
51    real, intent(in):: av_rho ! average mass density in the cold pool, in kg m-3
52
53    real, intent(in):: theta_v_uns
54    ! virtual potential temperature of the unsaturated downdraft, in K
55   
56    real, intent(in):: C_dw ! drag coefficient under cold pools, in m s-1
57
58    !----------------------------------------------
59
60    eps_cp_m = eps_cp
61    beta_m = beta
62    av_thvu_constant_m = av_thvu_constant
63    k_star_m = k_star
64    theta_vs_m = theta_vs
65    rate_uns_m = rate_uns
66    av_rho_m = av_rho
67    theta_v_uns_m = theta_v_uns
68    c_dw_m = c_dw
69
70  end subroutine set_module_var
71
72end module derivs_param
Note: See TracBrowser for help on using the repository browser.