''' 
This is a config file for Aeolus2. Here, all parameters and user settings can be defined. This file is then read
by Aeolus2_1.py main script. 
#Structure:
#1) Path specs, where to read/write data, how to name output files
#2) User settings: submodule activations
#3) Scaling and physical parameters
#4) Discretization Parameters
#5) Integration Parameters
'''


### Imports
#import numpy as np

### path specification ###
#input_folder          = "initial0g"               # location of the input files (used for restart of the model)
#start_file            = "output_1"            # output_00022
#output_folder         = 'initial0g'           

#n_iterations          = 2000000               # total number of iterations

#    --------------------- *** User setting & parameterization *** -----------
#smooth_run            = 1                     # 1: latitudinal fine resolution (768*384) with 1 min numerical time step; 
#fast_run              = 0                     # 1: time step=step min, some second oreder nonliear terms have been modified, resolution: (384*192)
#super_fast            = 0                     # 1: time step=step min, some second oreder nonliear terms have been modified, resolution: (192*96)   

moist_convection      = 1
#topography            = 1                     # 1: inclusion of topography
#restart_run           = 1                     # 1: restart from previous run except TQG_adjustment
Radiative_Transfer    = 1                     # 1: RRTMG_LW/SW: rapid radiative transfer model that utilizes the correlated-k approach to calculate longwave fluxes
unparallel_TQG_adj    = 0                     # 1: provides fully nonlinear thermo-quasi-geostrophic balanced state of initial state or any other initialization

summer_solstice       = 0                     # 1: just insolation of summer solstice will be imposed
winter_solstice       = 0                     # 1: just insolation of winter solstice will be imposed  

CO2ppmv               = 280. 
T_s0                  = 250.0277       # 262.1009#280.
L_d_eq                = 4262500.       # Equatorial Rossby deformation radius= sqrt(c/beta), where c=sqrt(BgH), 

#if topography==1:
delta1  = 0.60#0.35#0.57 # H1/H0
#else:
#   delta1  = 0.5#0.35

# Physical parameters
B1                    = 1.0 # Mean value of g*theta/theta_o at the lower (first) layer
B2                    = 1.1907 #Jan. 1980 Mean value of (theta2)/(theta1)+\delta(b2-b1) at the middle (second) layer
B3                    = 1.2 # Mean value of g*theta/theta_o at the upper (third) layer, it does not exist in two-layer configuration
#g = g0*np.array([[B1, B1, B1],[B1, B2, B2],[B1,B2,B3]])