Ignore:
Timestamp:
Mar 26, 2013 6:16:16 PM (9 years ago)
Author:
hoffmann
Message:

optimization of two-moments cloud physics

File:
1 edited

Legend:

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

    r1114 r1115  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! unused variables removed
    2323!
    2424! Former revisions:
     
    408408          c_u_m, c_u_m_l, c_v_m, c_v_m_l, c_w_m, c_w_m_l, ddzu, ddzu_pres,     &
    409409          dd2zu, dzu, ddzw, dzw, hyp, inflow_damping_factor, lad, l_grid,      &
    410           nr_init, ptdf_x, ptdf_y, pt_init, q_init, qr_init, rdf, rdf_sc,      &
    411           sa_init, ug, u_init, u_nzb_p1_for_vfc, vg, v_init, v_nzb_p1_for_vfc, &
    412           w_subs, zu, zw
     410          nc_1d, nr_1d, ptdf_x, ptdf_y, pt_1d, pt_init, q_1d, q_init, qc_1d,   &
     411          qr_1d, rdf, rdf_sc, sa_init, ug, u_init, u_nzb_p1_for_vfc, vg,      &
     412          v_init, v_nzb_p1_for_vfc, w_subs, zu, zw
    413413
    414414    REAL, DIMENSION(:,:), ALLOCATABLE ::                                       &
     
    430430          tend_nr, tend_q, tend_qr, tric, u_m_l, u_m_n, u_m_r, u_m_s, v_m_l,   &
    431431          v_m_n, v_m_r, v_m_s, w_m_l, w_m_n, w_m_r, w_m_s
    432            
    433432
    434433#if defined( __nopointer )
    435434    REAL, DIMENSION(:,:,:), ALLOCATABLE, TARGET ::                             &
    436           e, e_p, nr, nr_p, p, prho, pt, pt_p, q, q_p, ql, ql_c, ql_v, ql_vp,  &
    437           qr, qr_p, rho, sa, sa_p, te_m, tnr_m, tpt_m, tq_m, tqr_m, tsa_m,     &
    438           tu_m, tv_m, tw_m, u, u_p, v, v_p, vpt, w, w_p
     435          e, e_p, nr, nr_p, p, prho, pt, pt_p, q, q_p, qc, ql, ql_c, ql_v,     &
     436          ql_vp, qr, qr_p, rho, sa, sa_p, te_m, tnr_m, tpt_m, tq_m, tqr_m,     &
     437          tsa_m, tu_m, tv_m, tw_m, u, u_p, v, v_p, vpt, w, w_p
    439438#else
    440439    REAL, DIMENSION(:,:,:), ALLOCATABLE, TARGET ::                             &
    441440          e_1, e_2, e_3, p, prho_1, nr_1, nr_2, nr_3, pt_1, pt_2, pt_3, q_1,   &
    442           q_2, q_3, ql_v, ql_vp, ql_1, ql_2, qr_1, qr_2, qr_3, rho_1, sa_1,    &
    443           sa_2, sa_3, u_1, u_2, u_3, v_1, v_2, v_3, vpt_1, w_1, w_2, w_3
     441          q_2, q_3, qc_1, ql_v, ql_vp, ql_1, ql_2, qr_1, qr_2, qr_3, rho_1,    &
     442          sa_1, sa_2, sa_3, u_1, u_2, u_3, v_1, v_2, v_3, vpt_1, w_1, w_2, w_3
    444443
    445444    REAL, DIMENSION(:,:,:), POINTER ::                                         &
    446           e, e_p, nr, nr_p, prho, pt, pt_p, q, q_p, ql, ql_c, qr, qr_p, rho,   &
    447           sa, sa_p, te_m, tnr_m, tpt_m, tq_m, tqr_m, tsa_m, tu_m, tv_m, tw_m,  &
    448           u, u_p, v, v_p, vpt, w, w_p
     445          e, e_p, nr, nr_p, prho, pt, pt_p, q, q_p, qc, ql, ql_c, qr, qr_p,    &
     446          rho, sa, sa_p, te_m, tnr_m, tpt_m, tq_m, tqr_m, tsa_m, tu_m, tv_m,   &
     447          tw_m, u, u_p, v, v_p, vpt, w, w_p
    449448#endif
    450449
     
    474473
    475474    REAL, DIMENSION(:,:,:), ALLOCATABLE, TARGET :: &
    476           e_av, lpt_av, nr_av, p_av, pc_av, pr_av, prr_av, pt_av, q_av, ql_av, &
    477           ql_c_av, ql_v_av, ql_vp_av, qr_av, qv_av, rho_av, s_av, sa_av, u_av, &
    478           v_av, vpt_av, w_av
     475          e_av, lpt_av, nr_av, p_av, pc_av, pr_av, prr_av, pt_av, q_av, qc_av, &
     476          ql_av, ql_c_av, ql_v_av, ql_vp_av, qr_av, qv_av, rho_av, s_av, sa_av,&
     477          u_av, v_av, vpt_av, w_av
    479478 
    480479 END MODULE averaging
     
    525524             eps_ros = 1.0E-4,  & ! accuracy of Rosenbrock method
    526525             eps_sb = 1.0E-20,  & ! threshold in two-moments scheme
    527              k_cc = 9.44E09,    & ! const. rain-rain kernel (m3 kg-2 s-1)
     526             k_cc = 9.44E09,    & ! const. cloud-cloud kernel (m3 kg-2 s-1)
    528527             k_cr0 = 4.33,      & ! const. cloud-rain kernel (m3 kg-1 s-1)
    529528             k_rr = 7.12,       & ! const. rain-rain kernel (m3 kg-1 s-1)
    530529             k_br = 1000.,      & ! const. in breakup parametrization (m-1)
     530             k_st = 1.2E8,      & ! const. in drizzle parametrization (m-1 s-1)
    531531             kappa_rr = 60.7,   & ! const. in collision kernel (kg-1/3)
    532532             kin_vis_air = 1.4086E-5, & ! kin. viscosity of air (m2 s-1)
     
    536536             molecular_weight_of_solute = 0.05844, & ! mol. m. NaCl (kg mol-1)
    537537             molecular_weight_of_water = 0.01801528, & ! mol. m. H2O (kg mol-1)
    538              nc = 70.0E6,       & ! cloud droplet concentration
     538             nc_const = 70.0E6, & ! cloud droplet concentration
    539539             prec_time_const = 0.001, & !coef. in Kessler scheme
    540540             pirho_l, dpirho_l, & ! pi * rho_l / 6.0; 6.0 / ( pi * rho_l )
     
    545545             schmidt = 0.71,    & ! Schmidt number
    546546             schmidt_p_1d3,     & ! schmidt**( 1.0 / 3.0 )
     547             sigma_gc = 1.3,    & ! log-normal geometric standard deviation
    547548             stp = 2.5066282746310005, & ! parameter in gamma function
    548549             thermal_conductivity_l = 2.43E-2, & ! therm. cond. air (J m-1 s-1 K-1)
     
    554555             w_precipitation = 9.65      ! maximum terminal velocity (m s-1)
    555556
    556     REAL, DIMENSION(:), ALLOCATABLE     ::  dr, hyrho, lambda_r, mu_r, pt_d_t, &
    557                                             sed_nr, sed_q, sed_qr, t_d_pt, xr 
     557    REAL, DIMENSION(:), ALLOCATABLE     ::  hyrho, pt_d_t, t_d_pt 
    558558
    559559    REAL, DIMENSION(:,:), ALLOCATABLE   ::  precipitation_amount, &
     
    630630                             scalar_advec = 'ws-scheme'
    631631    CHARACTER (LEN=20)   ::  bc_e_b = 'neumann', bc_lr = 'cyclic', &
    632                              bc_nr_b = 'neumann', bc_nr_t = 'neumann', &
    633632                             bc_ns = 'cyclic', bc_p_b = 'neumann', &
    634633                             bc_p_t = 'dirichlet', bc_pt_b = 'dirichlet', &
    635634                             bc_pt_t = 'initial_gradient', &
    636635                             bc_q_b = 'dirichlet', bc_q_t = 'neumann', &
    637                              bc_qr_b = 'neumann', bc_qr_t = 'neumann',&
    638636                             bc_s_b = 'dirichlet', bc_s_t = 'neumann', &
    639637                             bc_sa_t = 'neumann', &
     
    681679                dp_level_ind_b = 0, dvrp_filecount = 0, &
    682680                dz_stretch_level_index, gamma_mg, gathered_size, &
    683                 grid_level, ibc_e_b, ibc_nr_b, ibc_nr_t, ibc_p_b, ibc_p_t, &
    684                 ibc_pt_b, ibc_pt_t, ibc_q_b, ibc_q_t, ibc_qr_b, ibc_qr_t, &
     681                grid_level, ibc_e_b, ibc_p_b, ibc_p_t, &
     682                ibc_pt_b, ibc_pt_t, ibc_q_b, ibc_q_t, &
    685683                ibc_sa_t, ibc_uv_b, ibc_uv_t, icloud_scheme, &
    686684                inflow_disturbance_begin = -1, inflow_disturbance_end = -1, &
     
    705703                mask_size(max_masks,3) = -1, mask_size_l(max_masks,3) = -1, &
    706704                mask_start_l(max_masks,3) = -1, &
    707                 nr_vertical_gradient_level_ind(10) = -9999, &
    708705                pt_vertical_gradient_level_ind(10) = -9999, &
    709706                q_vertical_gradient_level_ind(10) = -9999, &
    710                 qr_vertical_gradient_level_ind(10) = -9999, &
    711707                sa_vertical_gradient_level_ind(10) = -9999, &
    712708                section(100,3), section_xy(100) = -9999, &
     
    738734                constant_top_momentumflux = .FALSE., &
    739735                constant_top_salinityflux = .TRUE., &
    740                 constant_waterflux = .TRUE., constant_waterflux_nr = .TRUE., &
    741                 constant_waterflux_qr = .TRUE., create_disturbances = .TRUE., &
     736                constant_waterflux = .TRUE., create_disturbances = .TRUE., &
    742737                data_output_2d_on_each_pe = .TRUE., &
    743738                dissipation_control = .FALSE., disturbance_created = .FALSE., &
     
    779774             averaging_interval = 0.0, averaging_interval_pr = 9999999.9, &
    780775             averaging_interval_sp = 9999999.9, bc_pt_t_val, bc_q_t_val, &
    781              bc_qr_t_val, bc_nr_t_val, bottom_salinityflux = 0.0, &
     776             bottom_salinityflux = 0.0, &
    782777             building_height = 50.0, building_length_x = 50.0, &
    783778             building_length_y = 50.0, building_wall_left = 9999999.9, &
     
    799794             dt_do2d_xy = 9999999.9, dt_do2d_xz = 9999999.9, &
    800795             dt_do2d_yz = 9999999.9, dt_do3d = 9999999.9, dt_dvrp = 9999999.9, &
    801              dt_max = 20.0, dt_restart = 9999999.9, &
     796             dt_max = 20.0, dt_micro = -1.0, dt_restart = 9999999.9, &
    802797             dt_run_control = 60.0, dt_3d = -1.0, dz = -1.0, &
    803798             dz_max = 9999999.9, dz_stretch_factor = 1.08, &
     
    810805             maximum_cpu_time_allowed = 0.0,  &
    811806             molecular_viscosity = 1.461E-5, &
    812              nr_surface = 0.0, nr_surface_initial_change = 0.0, &
    813807             old_dt = 1.0E-10, omega = 7.29212E-5, omega_sor = 1.8, &
    814808             particle_maximum_age = 9999999.9, &
     
    819813             pt_surface = 300.0, pt_surface_initial_change = 0.0, &
    820814             q_surface = 0.0, q_surface_initial_change = 0.0, &
    821              qr_surface = 0.0, qr_surface_initial_change = 0.0, &
    822815             rayleigh_damping_factor = -1.0, rayleigh_damping_height = -1.0, &
    823816             recycling_width = 9999999.9, residual_limit = 1.0E-4, &
     
    832825             surface_heatflux = 9999999.9, surface_pressure = 1013.25, &
    833826             surface_scalarflux = 9999999.9, surface_waterflux = 9999999.9, &
    834              surface_waterflux_nr = 0.0, surface_waterflux_qr = 0.0, &
    835827             s_surface = 0.0, s_surface_initial_change = 0.0, &
    836828             termination_time_needed = -1.0, time_coupling = 0.0, &
     
    854846             lad_vertical_gradient_level(10) = -9999999.9, &
    855847             mask_scale(3), &
    856              nr_vertical_gradient(10) = 0.0, &
    857              nr_vertical_gradient_level(10) = -1.0, &
    858848             pt_vertical_gradient(10) = 0.0, &
    859849             pt_vertical_gradient_level(10) = -9999999.9, &
    860850             q_vertical_gradient(10) = 0.0, &
    861851             q_vertical_gradient_level(10) = -1.0, &
    862              qr_vertical_gradient(10) = 0.0, &
    863              qr_vertical_gradient_level(10) = -1.0, &
    864852             s_vertical_gradient(10) = 0.0, &
    865853             s_vertical_gradient_level(10) = -1.0, &
     
    877865             volume_flow(1:2) = 0.0, volume_flow_area(1:2) = 0.0, &
    878866             volume_flow_initial(1:2) = 0.0, wall_heatflux(0:4) = 0.0, &
    879              wall_humidityflux(0:4) = 0.0, wall_nrflux(0:4) = 0.0, & 
     867             wall_humidityflux(0:4) = 0.0, wall_nrflux(0:4) = 0.0, &
    880868             wall_qflux(0:4) = 0.0, wall_qrflux(0:4) = 0.0, &
    881869             wall_salinityflux(0:4) = 0.0, wall_scalarflux(0:4) = 0.0, &
Note: See TracChangeset for help on using the changeset viewer.