Ignore:
Timestamp:
Oct 21, 2014 10:53:05 AM (7 years ago)
Author:
kanani
Message:

New:
---
Subroutine init_plant_canopy added to module plant_canopy_model_mod. (plant_canopy_model)
Alternative method for lad-profile construction added, also, new parameters added.
(header, package_parin, plant_canopy_model, read_var_list, write_var_list)
plant_canopy_model-dependency added to several subroutines. (Makefile)
New package/namelist canopy_par for canopy-related parameters added. (package_parin)

Changed:
---
Code structure of the plant canopy model changed, all canopy-model related code
combined to module plant_canopy_model_mod. (check_parameters, init_3d_model,
modules, timestep)
Module plant_canopy_model_mod added in USE-lists of some subroutines. (check_parameters,
header, init_3d_model, package_parin, read_var_list, user_init_plant_canopy, write_var_list)
Canopy initialization moved to new subroutine init_plant_canopy. (check_parameters,
init_3d_model, plant_canopy_model)
Calculation of canopy timestep-criterion removed, instead, the canopy
drag is now directly limited in the calculation of the canopy tendency terms.
(plant_canopy_model, timestep)
Some parameters renamed. (check_parameters, header, init_plant_canopy,
plant_canopy_model, read_var_list, write_var_list)
Unnecessary 3d-arrays removed. (init_plant_canopy, plant_canopy_model, user_init_plant_canopy)
Parameter checks regarding canopy initialization added. (check_parameters)
All canopy steering parameters moved from namelist inipar to canopy_par. (package_parin, parin)
Some redundant MPI communication removed. (init_plant_canopy)

Bugfix:
---
Missing KIND-attribute for REAL constant added. (check_parameters)
DO-WHILE-loop for lad-profile output restricted. (header)
Removed double-listing of use_upstream_for_tke in ONLY-list of module
control_parameters. (prognostic_equations)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/modules.f90

    r1469 r1484  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! Changes due to new module structure of the plant canopy model:
     23!   canopy-model related parameters/variables moved to module
     24!   plant_canopy_model_mod
    2325!
    2426! Former revisions:
     
    3840! 1429 2014-07-15 12:53:45Z knoop
    3941! +ensemble_member_nr to prepare the random_generator for ensemble runs
    40 ! 
     42!
    4143! 1382 2014-04-30 12:15:41Z boeske
    4244! Renamed variables which store large scale forcing tendencies
     
    294296    REAL(wp), DIMENSION(:), ALLOCATABLE ::                                     &
    295297          c_u_m, c_u_m_l, c_v_m, c_v_m_l, c_w_m, c_w_m_l, ddzu, ddzu_pres,     &
    296           dd2zu, dzu, ddzw, dzw, hyp, inflow_damping_factor, lad, l_grid,      &
     298          dd2zu, dzu, ddzw, dzw, hyp, inflow_damping_factor, l_grid,           &
    297299          nc_1d, nr_1d, ptdf_x, ptdf_y, p_surf, pt_surf, pt_1d, pt_init,       &
    298300          qsws_surf, q_1d, q_init, q_surf, qc_1d, qr_1d, rdf, rdf_sc,          &
     
    313315
    314316    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::                                 &
    315           canopy_heat_flux, cdc, d, de_dx, de_dy, de_dz, diss, diss_l_e,       &
     317          d, de_dx, de_dy, de_dz, diss, diss_l_e,                              &
    316318          diss_l_nr, diss_l_pt, diss_l_q, diss_l_qr, diss_l_sa, diss_l_u,      &
    317319          diss_l_v, diss_l_w, flux_l_e, flux_l_nr, flux_l_pt, flux_l_q,        &
    318           flux_l_qr, flux_l_sa, flux_l_u, flux_l_v, flux_l_w, kh, km, lad_s,   &
    319           lad_u, lad_v, lad_w, lai, l_wall, p_loc, sec, sls, tend, tric,       &
     320          flux_l_qr, flux_l_sa, flux_l_u, flux_l_v, flux_l_w, kh, km,          &
     321          l_wall, p_loc, tend, tric,                                           &
    320322          u_m_l, u_m_n, u_m_r, u_m_s, v_m_l, v_m_n, v_m_r, v_m_s, w_m_l,       &
    321323          w_m_n, w_m_r, w_m_s
     
    533535                             bc_sa_t = 'neumann', &
    534536                             bc_uv_b = 'dirichlet', bc_uv_t = 'dirichlet', &
    535                              canopy_mode = 'block', &
    536537                             cloud_scheme = 'seifert_beheng', &
    537538                             coupling_mode = 'uncoupled', &
     
    586587                     nr_timesteps_this_run = 0, &
    587588                     nsor = 20, nsor_ini = 100, n_sor, normalizing_region = 0, &
    588                      nz_do3d = -9999, pch_index = 0, prt_time_count = 0, &
     589                     nz_do3d = -9999, prt_time_count = 0, &
    589590                     recycling_plane, runnr = 0, &
    590591                     skip_do_avs = 0, subdomain_size, terminate_coupled = 0, &
     
    596597                     do3d_no(0:1) = 0, do3d_time_count(0:1), &
    597598                     domask_no(max_masks,0:1) = 0, domask_time_count(max_masks,0:1),&
    598                      lad_vertical_gradient_level_ind(10) = -9999, &
    599599                     mask_size(max_masks,3) = -1, mask_size_l(max_masks,3) = -1, &
    600600                     mask_start_l(max_masks,3) = -1, &
     
    652652                outflow_l = .FALSE., outflow_n = .FALSE., outflow_r = .FALSE., &
    653653                outflow_s = .FALSE., passive_scalar = .FALSE., &
    654                 plant_canopy = .FALSE., prandtl_layer = .TRUE., &
     654                prandtl_layer = .TRUE., &
    655655                precipitation = .FALSE., &
    656656                radiation = .FALSE., random_heatflux = .FALSE., &
     
    682682                 canyon_width_x = 9999999.9_wp, canyon_width_y = 9999999.9_wp, &
    683683                 canyon_wall_left = 9999999.9_wp, canyon_wall_south = 9999999.9_wp, &
    684                  cthf = 0.0_wp, cfl_factor = -1.0_wp, cos_alpha_surface, &
     684                 cfl_factor = -1.0_wp, cos_alpha_surface, &
    685685                 coupling_start_time = 0.0_wp, disturbance_amplitude = 0.25_wp, &
    686686                 disturbance_energy_limit = 0.01_wp, &
    687687                 disturbance_level_b = -9999999.9_wp, &
    688688                 disturbance_level_t = -9999999.9_wp, &
    689                  dp_level_b = 0.0_wp, drag_coefficient = 0.0_wp, &
     689                 dp_level_b = 0.0_wp, &
    690690                 dt = -1.0_wp, dt_averaging_input = 0.0_wp, &
    691691                 dt_averaging_input_pr = 9999999.9_wp, dt_coupling = 9999999.9_wp, &
     
    703703                 f = 0.0_wp, fs = 0.0_wp, g = 9.81_wp, inflow_damping_height = 9999999.9_wp, &
    704704                 inflow_damping_width = 9999999.9_wp, kappa = 0.4_wp, km_constant = -1.0_wp,&
    705                  lad_surface = 0.0_wp, leaf_surface_concentration = 0.0_wp, &
    706705                 mask_scale_x = 1.0_wp, mask_scale_y = 1.0_wp, mask_scale_z = 1.0_wp, &
    707706                 maximum_cpu_time_allowed = 0.0_wp,  &
     
    719718                 restart_time = 9999999.9_wp, rho_reference, rho_surface, &
    720719                 rif_max = 1.0_wp, rif_min = -5.0_wp, roughness_length = 0.1_wp, &
    721                  sa_surface = 35.0_wp, scalar_exchange_coefficient = 0.0_wp, &
     720                 sa_surface = 35.0_wp, &
    722721                 simulated_time = 0.0_wp, simulated_time_at_begin, sin_alpha_surface, &
    723722                 skip_time_data_output = 0.0_wp, skip_time_data_output_av = 9999999.9_wp,&
     
    745744    REAL(wp) ::  do2d_xy_last_time(0:1) = -1.0_wp, do2d_xz_last_time(0:1) = -1.0_wp, &
    746745                 do2d_yz_last_time(0:1) = -1.0_wp, dpdxy(1:2) = 0.0_wp, &
    747                  dt_domask(max_masks) = 9999999.9_wp, lad_vertical_gradient(10) = 0.0_wp,&
    748                  lad_vertical_gradient_level(10) = -9999999.9_wp, &
     746                 dt_domask(max_masks) = 9999999.9_wp, &
    749747                 mask_scale(3), &
    750748                 pt_vertical_gradient(10) = 0.0_wp, &
Note: See TracChangeset for help on using the changeset viewer.