Skip to content
Snippets Groups Projects
Select Git revision
  • 500cc4cf3393c07784616b4a14dbda56eb4001c5
  • master default protected
  • dev-lkugler
  • teaching-2024
  • old_config_2023-05 protected
  • v2025.2
  • v2024.6
  • v2024.2.20
8 results

cfg.py

Blame
  • cfg.py 2.06 KiB
    import numpy as np
    from config import clusters  # from . = problem in archivedir
    cluster = clusters.vsc  # change cluster configuration here
    
    
    class ExperimentConfiguration(object):
        def __init__(self):
            pass
    
    
    exp = ExperimentConfiguration()
    exp.expname = "exp_v1.16_Pwbub-1_40mem"
    exp.model_dx = 2000
    exp.n_ens = 40
    exp.n_nodes = 10
    
    n_obs = 121  #961 900: 10km resoltn # radar: n_obs for each observation height level
    
    vis = dict(plotname='VIS 0.6µm',  plotunits='[1]',
               kind='MSG_4_SEVIRI_BDRF', sat_channel=1, n_obs=n_obs, 
               error_generate=0.03, error_assimilate=0.06,
               cov_loc_radius_km=32)
    
    wv73 = dict(plotname='Brightness temperature WV 7.3µm',  plotunits='[K]',
                kind='MSG_4_SEVIRI_TB', sat_channel=6, n_obs=n_obs, 
                error_generate=1., error_assimilate=False,
                cov_loc_radius_km=32)
    
    ir108 = dict(plotname='Brightness temperature IR 10.8µm', plotunits='[K]',
                 kind='MSG_4_SEVIRI_TB', sat_channel=9, n_obs=n_obs, 
                 error_generate=5., error_assimilate=10.,
                 cov_loc_radius_km=32)
    
    radar = dict(plotname='Radar reflectivity', plotunits='[dBz]',
                 kind='RADAR_REFLECTIVITY', n_obs=n_obs, 
                 error_generate=2.5, error_assimilate=5.,
                 heights=np.arange(1000, 15001, 1000),
                 cov_loc_radius_km=32, cov_loc_vert_km=4)
    
    t2m = dict(plotname='SYNOP Temperature', plotunits='[K]',
               kind='SYNOP_TEMPERATURE', n_obs=n_obs, 
               error_generate=0.1, error_assimilate=1.,
               cov_loc_radius_km=20, cov_loc_vert_km=3)
    
    psfc = dict(plotname='SYNOP Pressure', plotunits='[dBz]',
                kind='SYNOP_SURFACE_PRESSURE', n_obs=n_obs, 
                error_generate=50., error_assimilate=100.,
                cov_loc_radius_km=32, cov_loc_vert_km=5)
    
    
    exp.observations = [] #wv73, vis] # 108, wv73, vis]
    exp.update_vars = ['T', 'QVAPOR', 'QCLOUD', 'QICE','CLDFRA']
    #exp.update_vars = ['U', 'V', 'T', 'PH', 'MU', 'QVAPOR', 'QCLOUD', 'QICE', 'TSK', 'CLDFRA']
    
    # directory paths depend on the name of the experiment
    cluster.expname = exp.expname