Changes between Version 1 and Version 2 of doc/app/examples/rrtmg


Ignore:
Timestamp:
May 4, 2015 1:01:50 PM (10 years ago)
Author:
maronga
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • doc/app/examples/rrtmg

    v1 v2  
    1 Coming soon.
     1=== Run with land surface model and coupled to the RRTMG radiation model ===
     2
     3When using the land surface model (LSM) in combination with the RRTMG radiation scheme, at least two additional NetCDF files are required in the INPUT directory with suffixes _rlw (longwave radiation) and _rsw (shortwave radiation). They have to be placed in the INPUT directory on the host where the simulation will run. The layout of the data files must be equal to the standard files provided by RRTMG (see [source:palm/trunk/lib/rrtmg/data here]). A sounding file for trace gas data can be added optionally, having the suffic _rsnd.
     4
     5In addition, the PALM configuration file .mrun.config has to be extended as follows
     6{{{
     7...
     8#----------------------------------------------------------------------------
     9# List of input-files
     10#----------------------------------------------------------------------------
     11...
     12rrtmg_lw.nc              in:locopt     d3#:d3f  $base_data/$fname/INPUT    _rlw nc
     13rrtmg_sw.nc              in:locopt     d3#:d3f  $base_data/$fname/INPUT    _rsw nc
     14RAD_SND_DATA             in:locopt     d3#:d3f  $base_data/$fname/INPUT    _rsnd nc
     15...
     16}}}
     17
     18With the following parameter file example_lsm_rrtmg_p3d a simulation with LSM and RRTMG can be performed.
     19{{{
     20 &inipar nx = 63, ny = 63, nz = 96,
     21         dx = 50.0, dy = 50.0, dz = 25.0,
     22         dz_stretch_level = 1800.0,
     23
     24         fft_method = 'fftw',
     25
     26         initializing_actions = 'set_constant_profiles',
     27
     28         ug_surface = 2.5, vg_surface = 0.0,
     29         pt_surface = 292.5,
     30         pt_vertical_gradient       =   0.0,  12.0,   0.5,
     31         pt_vertical_gradient_level =   0.0, 800.0, 825.0,
     32
     33         humidity = .TRUE.,
     34         cloud_physics = .TRUE.,
     35         precipitation = .TRUE.,
     36
     37         q_surface = 0.005,
     38
     39         surface_pressure = 1000.0,
     40
     41         phi = 51.96,
     42
     43         prandtl_layer = .TRUE.,
     44         bc_pt_b = 'dirichlet',
     45         bc_q_b = 'dirichlet',
     46         bc_e_b = 'neumann', /
     47
     48
     49 &d3par  end_time = 18000.0,
     50
     51         create_disturbances = .T.,
     52         dt_disturb = 150.0, disturbance_energy_limit = 0.01,
     53
     54         data_output_2d_on_each_pe = .T.,
     55
     56         dt_run_control = 1.0,
     57
     58         data_output = 'ql', 'qv', 'vpt', 'pt', 'shf*_xy', 'qsws*_xy',
     59                       'lwp*_xy', 'prr', 'qc', 'qr',
     60                       't_soil', 't_soil_av', 'pt_av',
     61                       'm_soil', 'm_soil_av',
     62                       'c_liq*_xy', 'c_liq*_xy_av',
     63                       'c_soil*_xy', 'c_soil*_xy_av',
     64                       'c_veg*_xy', 'c_veg*_xy_av',
     65                       'lai*_xy', 'lai*_xy_av',
     66                       'ghf_eb*_xy', 'ghf_eb*_xy_av',
     67                       'shf_eb*_xy', 'shf_eb*_xy_av',
     68                       'qsws_eb*_xy', 'qsws_eb*_xy_av',
     69                       'qsws_liq_eb*_xy', 'qsws_liq_eb*_xy_av',
     70                       'qsws_soil_eb*_xy', 'qsws_soil_eb*_xy_av',
     71                       'qsws_veg_eb*_xy', 'qsws_veg_eb*_xy_av',
     72                       'rad_net*_xy', 'rad_net*_xy_av',
     73                       'rad_sw_in', 'rad_sw_out', 'rad_lw_in', 'rad_lw_out',
     74
     75        data_output_pr = '#pt', '#qv', 'wpt', '#t_soil', '#m_soil', 'u', 'v', 'rad_sw_in', 'rad_sw_out', 'rad_lw_in', 'rad_lw_out',
     76         
     77         dt_data_output = 600.0,
     78
     79         dt_data_output_av = 600.0,
     80         averaging_interval = 10.0,
     81         dt_dopr = 600.0,
     82/
     83
     84 &lsm_par veg_type = 2, soil_type = 7,
     85
     86          conserve_water_content = .T.,
     87          dewfall = .T.,
     88
     89          c_surface = 20000.0,
     90          vegetation_coverage = 0.9,
     91          z0_eb = 0.1,
     92          z0h_eb = 0.02,
     93          lambda_surface_stable = 5.0,
     94          lambda_surface_unstable = 5.0,
     95          min_canopy_resistance = 110.0,
     96          min_soil_resistance = 50.0,
     97          leaf_area_index = 2.0,
     98
     99          zs = 0.07, 0.34, 1.47, 2.86,
     100          soil_moisture = 0.3, 0.3, 0.3, 0.3,
     101          root_fraction = 0.25, 0.25, 0.25, 0.25,
     102          soil_temperature = 290.0, 287.0, 285.0, 283.0, 283.0,
     103/
     104
     105 &radiation_par radiation_scheme = 'rrtmg',
     106                lambda = 0.0,
     107                albedo_type = 5,
     108                constant_albedo = .F.,
     109                day_init = 287,
     110                time_utc_init = 28800.0,
     111                dt_radiation = 60.0,
     112/
     113}}}