1 | # Stochastic forcing option, for ARW only |
---|
2 | |
---|
3 | # 3D arrays |
---|
4 | state real VERTSTRUCC ikj dyn_em 1 - rhd "VERTSTRUCC " "vertical structure for stoch. forcing " "" |
---|
5 | state real VERTSTRUCS ikj dyn_em 1 - rhd "VERTSTRUCS " "vertical structure for stoch. forcing " "" |
---|
6 | |
---|
7 | state real ru_tendf_stoch ikjb dyn_em 1 X rhdf=(p2c:dt) "ru_tendf_stoch" "stochastic forcing, U " "" |
---|
8 | state real rv_tendf_stoch ikjb dyn_em 1 Y rhdf=(p2c:dt) "rv_tendf_stoch" "stochastic forcing, V " "" |
---|
9 | state real rt_tendf_stoch ikjb dyn_em 1 - rhdf=(p2c:dt) "rt_tendf_stoch" "stochastic forcing, T " "" |
---|
10 | |
---|
11 | # 2d arrays |
---|
12 | state real SPSTREAMFORCC ij misc 1 - r "SPSTREAMFORCC" "real spect. coeff. of stoch. streamfunction perturb." "" |
---|
13 | state real SPSTREAMFORCS ij misc 1 - r "SPSTREAMFORCS" "imag. spect. coeff. of stoch. streamfunction perturb." "" |
---|
14 | state real SPTFORCC ij misc 1 - r "SPTFORCC" "real spect. coeff. of stoch. temperature perturb." "" "" |
---|
15 | state real SPTFORCS ij misc 1 - r "SPTFORCS" "imag. spect. coeff. of stoch. temperature perturb." "" "" |
---|
16 | state real SPSTREAM_AMP ij misc 1 - r "SPSTREAM_AMP" "amplitude of stoch. streamfunction perturb." "" "" |
---|
17 | state real SPT_AMP ij misc 1 - r "SPT_AMP" "amplitude of stoch. temperature perturb." "" "" |
---|
18 | |
---|
19 | # 1d arrays for FFT transpose |
---|
20 | state real RU_REAL ikj dyn_em 1 XYZ |
---|
21 | state real RU_IMAG ikj dyn_em 1 XYZ |
---|
22 | state real RU_REAL_xxx ikjx dyn_em 1 XYZ |
---|
23 | state real RU_REAL_yyy ikjy dyn_em 1 XYZ |
---|
24 | state real RU_IMAG_xxx ikjx dyn_em 1 XYZ |
---|
25 | state real RU_IMAG_yyy ikjy dyn_em 1 XYZ |
---|
26 | state real RV_REAL ikj dyn_em 1 XYZ |
---|
27 | state real RV_IMAG ikj dyn_em 1 XYZ |
---|
28 | state real RV_REAL_xxx ikjx dyn_em 1 XYZ |
---|
29 | state real RV_REAL_yyy ikjy dyn_em 1 XYZ |
---|
30 | state real RV_IMAG_xxx ikjx dyn_em 1 XYZ |
---|
31 | state real RV_IMAG_yyy ikjy dyn_em 1 XYZ |
---|
32 | state real RT_REAL ikj dyn_em 1 XYZ |
---|
33 | state real RT_IMAG ikj dyn_em 1 XYZ |
---|
34 | state real RT_REAL_xxx ikjx dyn_em 1 XYZ |
---|
35 | state real RT_REAL_yyy ikjy dyn_em 1 XYZ |
---|
36 | state real RT_IMAG_xxx ikjx dyn_em 1 XYZ |
---|
37 | state real RT_IMAG_yyy ikjy dyn_em 1 XYZ |
---|
38 | |
---|
39 | xpose XPOSE_STOCH_BACK_U_REAL dyn_em RU_REAL,RU_REAL_xxx,RU_REAL_yyy |
---|
40 | xpose XPOSE_STOCH_BACK_U_IMAG dyn_em RU_IMAG,RU_IMAG_xxx,RU_IMAG_yyy |
---|
41 | xpose XPOSE_STOCH_BACK_V_REAL dyn_em RV_REAL,RV_REAL_xxx,RV_REAL_yyy |
---|
42 | xpose XPOSE_STOCH_BACK_V_IMAG dyn_em RV_IMAG,RV_IMAG_xxx,RV_IMAG_yyy |
---|
43 | xpose XPOSE_STOCH_BACK_T_REAL dyn_em RT_REAL,RT_REAL_xxx,RT_REAL_yyy |
---|
44 | xpose XPOSE_STOCH_BACK_T_IMAG dyn_em RT_IMAG,RT_IMAG_xxx,RT_IMAG_yyy |
---|
45 | |
---|
46 | # Variables |
---|
47 | state integer SEED1 - misc 1 - rh "SEED1" "RANDOM SEED NUMBER 1" "" |
---|
48 | state integer SEED2 - misc 1 - rh "SEED2" "RANDOM SEED NUMBER 2" "" |
---|
49 | state logical did_stoch - misc 1 - r "DID_STOCH" "Logical to tell us that we already did the initialization for dom 1" "" |
---|
50 | |
---|
51 | # Namelist parameters |
---|
52 | rconfig integer stoch_force_opt namelist,physics max_domains 0 - "stochastic forcing option: 0=none, 1=backscatter" |
---|
53 | rconfig integer stoch_vertstruc_opt namelist,physics max_domains 0 - "vertical structure for stochastic forcing: 0=constant, 1=random phase" |
---|
54 | rconfig integer nens namelist,physics 1 1 - "ensemble member" "" "" |
---|
55 | rconfig real tot_backscat_psi namelist,physics max_domains 115200 - "total backscattered energy 240*240*2" "m2 s-3" "" |
---|
56 | rconfig real tot_backscat_t namelist,physics max_domains 2.0E-06 - "total backscattered temperature" "Cp m2 s-3" "" |
---|
57 | |
---|
58 | # Package declarations |
---|
59 | package no_stoch_force stoch_force_opt==0 - - |
---|
60 | package stoch_backscatter stoch_force_opt==1 - state:ru_tendf_stoch,rv_tendf_stoch,rt_tendf_stoch,SPSTREAMFORCC,SPSTREAMFORCS,SPTFORCC,SPTFORCS,SPSTREAM_AMP,SPT_AMP,VERTSTRUCC,VERTSTRUCS,RU_IMAG,RU_REAL_xxx,RU_REAL_yyy,RU_IMAG_xxx,RU_IMAG_yyy,RV_IMAG,RV_REAL_xxx,RV_REAL_yyy,RV_IMAG_xxx,RV_IMAG_yyy,RT_IMAG,RT_REAL_xxx,RT_REAL_yyy,RT_IMAG_xxx,RT_IMAG_yyy,RU_REAL,RV_REAL,RT_REAL |
---|