A Modular Arbitrary-Order Ocean-Atmosphere Model -- Stochastic implementation
stoch_params.f90
Go to the documentation of this file.
1 
2 ! stoch_params.f90
3 !
4 !> The stochastic models parameters module.
5 !
6 !> @copyright
7 !> 2018 Jonathan Demaeyer.
8 !> See LICENSE.txt for license information.
9 !
10 !---------------------------------------------------------------------------
11 !
12 !> @remark
13 !>
14 !>
15 !
16 !---------------------------------------------------------------------------
17 
18 
20  USE params, only: dt
21  IMPLICIT NONE
22 
23  PUBLIC
24 
25  REAL(KIND=8) :: mnuti !< Multiplicative noise update time interval
26 
27  REAL(KIND=8) :: t_trans_stoch !< Transient time period of the stochastic model evolution
28  REAL(KIND=8) :: q_ar !< Atmospheric resolved component noise amplitude
29  REAL(KIND=8) :: q_au !< Atmospheric unresolved component noise amplitude
30  REAL(KIND=8) :: q_or !< Oceanic resolved component noise amplitude
31  REAL(KIND=8) :: q_ou !< Oceanic unresolved component noise amplitude
32  REAL(KIND=8) :: dtn !< Square root of the timestep
33  REAL(KIND=8) :: eps_pert !< Perturbation parameter for the coupling
34  REAL(KIND=8) :: tdelta !< Time separation parameter
35 
36  REAL(KIND=8) :: muti !< Memory update time interval
37  REAL(KIND=8) :: meml !< Time over which the memory kernel is integrated
38  REAL(KIND=8) :: t_trans_mem !< Transient time period to initialize the memory term
39  CHARACTER(len=4) :: x_int_mode !< Integration mode for the resolved component
40  REAL(KIND=8) :: dts !< Intrisic resolved dynamics time step
41 
42  INTEGER :: mems !< Number of steps in the memory kernel integral
43  REAL(KIND=8) :: dtsn !< Square root of the intrisic resolved dynamics time step
44 
45  REAL(KIND=8) :: maxint !< Upper integration limit of the correlations
46  CHARACTER(LEN=4) :: load_mode !< Loading mode for the correlations
47  CHARACTER(LEN=4) :: int_corr_mode !< Correlation integration mode
48 
49  CHARACTER(len=4) :: mode !< Stochastic mode parameter
50 
51 
52 
53 
54 CONTAINS
55 
56  !> Stochastic parameters initialization routine
57  SUBROUTINE init_stoch_params
58 
59  namelist /mtvparams/ mnuti
60  namelist /stparams/ q_ar,q_au,q_or,q_ou,eps_pert,tdelta,t_trans_stoch
61  namelist /wlparams/ muti,meml,x_int_mode,dts,t_trans_mem
62  namelist /corr_init_mode/ load_mode,int_corr_mode,maxint
63  namelist /stoch_int_params/ mode
64 
65 
66  OPEN(8, file="stoch_params.nml", status='OLD', recl=80, delim='APOSTROPHE')
67  READ(8,nml=mtvparams)
68  READ(8,nml=wlparams)
69  READ(8,nml=stparams)
70  READ(8,nml=stoch_int_params)
71  READ(8,nml=corr_init_mode)
72  CLOSE(8)
73 
74  dtn=sqrt(dt)
75  dtsn=sqrt(dts)
76  mems=ceiling(meml/muti)
77 
80 
81  END SUBROUTINE init_stoch_params
82 END MODULE stoch_params
The stochastic models parameters module.
real(kind=8) q_au
Atmospheric unresolved component noise amplitude.
character(len=4) load_mode
Loading mode for the correlations.
subroutine init_stoch_params
Stochastic parameters initialization routine.
real(kind=8) muti
Memory update time interval.
real(kind=8) q_or
Oceanic resolved component noise amplitude.
character(len=4) x_int_mode
Integration mode for the resolved component.
real(kind=8) t_trans_stoch
Transient time period of the stochastic model evolution.
character(len=4) int_corr_mode
Correlation integration mode.
real(kind=8) dtn
Square root of the timestep.
real(kind=8) tdelta
Time separation parameter.
real(kind=8) q_ou
Oceanic unresolved component noise amplitude.
real(kind=8) dtsn
Square root of the intrisic resolved dynamics time step.
real(kind=8) meml
Time over which the memory kernel is integrated.
The model parameters module.
Definition: params.f90:18
real(kind=8) maxint
Upper integration limit of the correlations.
real(kind=8) eps_pert
Perturbation parameter for the coupling.
real(kind=8) dt
Integration time step.
Definition: params.f90:77
character(len=4) mode
Stochastic mode parameter.
real(kind=8) mnuti
Multiplicative noise update time interval.
real(kind=8) q_ar
Atmospheric resolved component noise amplitude.
real(kind=8) t_trans_mem
Transient time period to initialize the memory term.
real(kind=8) dts
Intrisic resolved dynamics time step.
integer mems
Number of steps in the memory kernel integral.