Ignore:
Timestamp:
May 31, 2019 3:19:05 PM (2 years ago)
Author:
monakurppa
Message:

remove salsa_util_mod.f90 and correct some bugs in salsa and salsa test case

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/TESTS/cases/urban_environment_salsa/INPUT/urban_environment_salsa_p3d

    r3924 r4012  
    1919!-- Numerics
    2020!------------------------------------------------------------------------------
    21     psolver              = 'multigrid',
    22     momentum_advec       = 'ws-scheme',  ! default advection scheme
    23     scalar_advec         = 'ws-scheme',
    24     conserve_volume_flow = .T.,
     21    fft_method           = 'fftw',
    2522
    2623!
     
    2825!------------------------------------------------------------------------------     
    2926    humidity             = .T.,
    30     q_surface            = 0.008,
     27    q_surface            = 0.001,
    3128
    3229!
    3330!-- Initialization
    3431!------------------------------------------------------------------------------
    35     initializing_actions = 'set_constant_profiles',
    36 
    37     pt_surface           = 293.15,
     32    initializing_actions = 'set_constant_profiles',   !'set_constant_profiles', or 'inifor',
     33                                                      ! or both (met. from inifor, chem constant profiles)
     34
     35    pt_surface           = 277.15,
    3836    ug_surface           = 1.0,
    3937    vg_surface           = 0.0,
    4038
    41     roughness_length     = 0.05,
    42     neutral              = .T.,
    43    
    44     day_of_year_init     =  65,          ! March 06, to avoid extreme temperatures
    45     time_utc_init        = 71700.0,      !(19:55, UTC in seconds)
    46 
    47 !
    48 !-- Boundary conditions
    49 !------------------------------------------------------------------------------
    50     bc_lr                = 'cyclic',
    51     bc_ns                = 'cyclic',
    52     bc_uv_t              = 'dirichlet',  ! channel flow boundary condition
     39    day_of_year_init     = 65,             ! march 06, to avoid extreme temperatures
     40    time_utc_init        = 36000.0,        ! (0.0 s = 00 UTC midnight),
     41
     42!
     43!-- Wall/soil spinup
     44!------------------------------------------------------------------------------   
     45    spinup_time         = 1200.0,
     46    spinup_pt_mean      = 284.15,
     47    spinup_pt_amplitude = 10.0,             ! give a dirunal cycle from 6 - 16 °C
     48    dt_spinup           = 120.0,
     49    data_output_during_spinup = .F.,
    5350
    5451!
     
    6057!-- Physics
    6158!------------------------------------------------------------------------------
    62     longitude            = 13.4,         ! Berlin longitude
    63     latitude             = 52.5,         ! Berlin latitude
    64      
     59    longitude            = 13.4,
     60    latitude             = 52.5,
     61
    6562 / ! end of inipar namelist
    6663
     
    7976
    8077!
    81 !-- General output settings
    82 !------------------------------------------------------------------------------
    83     netcdf_data_format    = 4,           ! use NetCDF4
    84 
    85 !
    8678!-- Run-control/timeseries output settings
    8779!------------------------------------------------------------------------------
    8880    dt_run_control        =  0.0,
    89     dt_dots               =  0.0,
    90 
     81    dt_dots               =  10.0,
     82
     83!
     84!-- Profile output settings
     85!------------------------------------------------------------------------------
     86    skip_time_dopr        =  0.0,
     87    dt_dopr               = 60.0,
     88    averaging_interval_pr = 60.0,
     89    dt_averaging_input_pr =  0.0,
     90
     91    data_output_pr        = '#u',  'u*2',  'wu',  'w*u*',  'w"u"',
     92                            '#v',  'v*2',  'wv',  'w*v*',  'w"v"',
     93!                            'w',   'w*2',
     94                            '#theta', 'theta*2', 'wtheta', 'w*theta*', 'w"theta"',
     95                            '#q',  'q*2',
    9196!
    9297!-- 2D/3D output settings
    9398!------------------------------------------------------------------------------
    9499    do3d_at_begin      = .T.,
    95 
    96     dt_data_output     = 10.0,
    97 
     100    do2d_at_begin      = .T.,
     101
     102    dt_data_output     = 10.0,
     103    dt_data_output_av  = 10.0,
     104    averaging_interval = 10.0,
     105    dt_averaging_input =  0.0,
     106
     107    section_xy         = 0,
    98108
    99109    data_output = 'u', 'v', 'w',
    100                   'Ntot', 'LDSA', 'PM2.5', 's_OC', 's_NH', 's_NO',
    101                   'g_H2SO4','g_HNO3','g_NH3','g_OCNV','g_OCSV',
     110                  'Ntot', 'PM2.5', 'LDSA',
     111                  's_OC','g_OCNV','g_OCSV',
    102112                  'N_bin3', 'm_bin4',
    103113
     114
    104115 / ! end of d3par namelist
    105116
     117
     118!-------------------------------------------------------------------------------
     119!-- RADIATION MODEL PARAMETER NAMELIST
     120!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/radpar
     121!-------------------------------------------------------------------------------
     122  &radiation_parameters
     123
     124     radiation_scheme =  'clear-sky',    !'clear-sky' or 'rrtmg',
     125     albedo_type = 5,
     126     constant_albedo = .F.,
     127
     128     dt_radiation = 60.0,
     129
     130     surface_reflections = .T.,
     131     nrefsteps = 3,
     132
     133     rad_angular_discretization = .F.,
     134     raytrace_mpi_rma = .F.,
     135
     136  /
     137
     138
     139!-------------------------------------------------------------------------------
     140!-- LAND SURFACE MODEL PARAMETER NAMELIST
     141!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/lsmpar
     142!-------------------------------------------------------------------------------
     143 &land_surface_parameters
     144
     145    constant_roughness     = .T.,
     146    aero_resist_kray       = .T.,
     147
     148    vegetation_type        = 2,
     149    soil_type              = 3,
     150    conserve_water_content = .T.,
     151
     152    root_fraction          = 0.49, 0.26, 0.25, 0.0, 0.0, 0.0, 0.0, 0.0,
     153    soil_temperature       = 277.15, 277.0, 277.0, 277.0, 277.0, 277.0, 277.0, 277.0,
     154    soil_moisture          = 0.30,   0.30,  0.30,  0.30,  0.30, 0.30,  0.30,  0.30,
     155    deep_soil_temperature  = 277.0,
     156   
     157 / ! end of lsm_par namelist
     158
     159
     160!-------------------------------------------------------------------------------
     161!-- URBAN SURFACE MODEL PARAMETER NAMELIST
     162!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/usmpar
     163!-------------------------------------------------------------------------------
     164 &urban_surface_parameters
     165
     166    usm_material_model     = .T.,
     167    usm_wall_mod = .T.,
     168
     169 / ! end of urban_surface_par namelist
     170 
    106171
    107172!-------------------------------------------------------------------------------
     
    112177 
    113178    canopy_mode       = 'read_from_file_3d',
    114     canopy_drag_coeff = 0.5,
     179    canopy_drag_coeff = 0.3,
    115180     
    116181 / ! end of canopy_par namelist
Note: See TracChangeset for help on using the changeset viewer.