Changeset 1115 for palm/trunk/SOURCE/modules.f90
- Timestamp:
- Mar 26, 2013 6:16:16 PM (12 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/modules.f90 (modified) (19 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/modules.f90
r1114 r1115 20 20 ! Current revisions: 21 21 ! ------------------ 22 ! 22 ! unused variables removed 23 23 ! 24 24 ! Former revisions: … … 408 408 c_u_m, c_u_m_l, c_v_m, c_v_m_l, c_w_m, c_w_m_l, ddzu, ddzu_pres, & 409 409 dd2zu, dzu, ddzw, dzw, hyp, inflow_damping_factor, lad, l_grid, & 410 n r_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, zw410 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 413 413 414 414 REAL, DIMENSION(:,:), ALLOCATABLE :: & … … 430 430 tend_nr, tend_q, tend_qr, tric, u_m_l, u_m_n, u_m_r, u_m_s, v_m_l, & 431 431 v_m_n, v_m_r, v_m_s, w_m_l, w_m_n, w_m_r, w_m_s 432 433 432 434 433 #if defined( __nopointer ) 435 434 REAL, DIMENSION(:,:,:), ALLOCATABLE, TARGET :: & 436 e, e_p, nr, nr_p, p, prho, pt, pt_p, q, q_p, q l, ql_c, ql_v, ql_vp,&437 q r, qr_p, rho, sa, sa_p, te_m, tnr_m, tpt_m, tq_m, tqr_m, tsa_m, &438 t u_m, tv_m, tw_m, u, u_p, v, v_p, vpt, w, w_p435 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 439 438 #else 440 439 REAL, DIMENSION(:,:,:), ALLOCATABLE, TARGET :: & 441 440 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, q l_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_3441 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 444 443 445 444 REAL, DIMENSION(:,:,:), POINTER :: & 446 e, e_p, nr, nr_p, prho, pt, pt_p, q, q_p, q l, 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_p445 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 449 448 #endif 450 449 … … 474 473 475 474 REAL, DIMENSION(:,:,:), ALLOCATABLE, TARGET :: & 476 e_av, lpt_av, nr_av, p_av, pc_av, pr_av, prr_av, pt_av, q_av, q l_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_av475 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 479 478 480 479 END MODULE averaging … … 525 524 eps_ros = 1.0E-4, & ! accuracy of Rosenbrock method 526 525 eps_sb = 1.0E-20, & ! threshold in two-moments scheme 527 k_cc = 9.44E09, & ! const. rain-rainkernel (m3 kg-2 s-1)526 k_cc = 9.44E09, & ! const. cloud-cloud kernel (m3 kg-2 s-1) 528 527 k_cr0 = 4.33, & ! const. cloud-rain kernel (m3 kg-1 s-1) 529 528 k_rr = 7.12, & ! const. rain-rain kernel (m3 kg-1 s-1) 530 529 k_br = 1000., & ! const. in breakup parametrization (m-1) 530 k_st = 1.2E8, & ! const. in drizzle parametrization (m-1 s-1) 531 531 kappa_rr = 60.7, & ! const. in collision kernel (kg-1/3) 532 532 kin_vis_air = 1.4086E-5, & ! kin. viscosity of air (m2 s-1) … … 536 536 molecular_weight_of_solute = 0.05844, & ! mol. m. NaCl (kg mol-1) 537 537 molecular_weight_of_water = 0.01801528, & ! mol. m. H2O (kg mol-1) 538 nc = 70.0E6, & ! cloud droplet concentration538 nc_const = 70.0E6, & ! cloud droplet concentration 539 539 prec_time_const = 0.001, & !coef. in Kessler scheme 540 540 pirho_l, dpirho_l, & ! pi * rho_l / 6.0; 6.0 / ( pi * rho_l ) … … 545 545 schmidt = 0.71, & ! Schmidt number 546 546 schmidt_p_1d3, & ! schmidt**( 1.0 / 3.0 ) 547 sigma_gc = 1.3, & ! log-normal geometric standard deviation 547 548 stp = 2.5066282746310005, & ! parameter in gamma function 548 549 thermal_conductivity_l = 2.43E-2, & ! therm. cond. air (J m-1 s-1 K-1) … … 554 555 w_precipitation = 9.65 ! maximum terminal velocity (m s-1) 555 556 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 558 558 559 559 REAL, DIMENSION(:,:), ALLOCATABLE :: precipitation_amount, & … … 630 630 scalar_advec = 'ws-scheme' 631 631 CHARACTER (LEN=20) :: bc_e_b = 'neumann', bc_lr = 'cyclic', & 632 bc_nr_b = 'neumann', bc_nr_t = 'neumann', &633 632 bc_ns = 'cyclic', bc_p_b = 'neumann', & 634 633 bc_p_t = 'dirichlet', bc_pt_b = 'dirichlet', & 635 634 bc_pt_t = 'initial_gradient', & 636 635 bc_q_b = 'dirichlet', bc_q_t = 'neumann', & 637 bc_qr_b = 'neumann', bc_qr_t = 'neumann',&638 636 bc_s_b = 'dirichlet', bc_s_t = 'neumann', & 639 637 bc_sa_t = 'neumann', & … … 681 679 dp_level_ind_b = 0, dvrp_filecount = 0, & 682 680 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, & 685 683 ibc_sa_t, ibc_uv_b, ibc_uv_t, icloud_scheme, & 686 684 inflow_disturbance_begin = -1, inflow_disturbance_end = -1, & … … 705 703 mask_size(max_masks,3) = -1, mask_size_l(max_masks,3) = -1, & 706 704 mask_start_l(max_masks,3) = -1, & 707 nr_vertical_gradient_level_ind(10) = -9999, &708 705 pt_vertical_gradient_level_ind(10) = -9999, & 709 706 q_vertical_gradient_level_ind(10) = -9999, & 710 qr_vertical_gradient_level_ind(10) = -9999, &711 707 sa_vertical_gradient_level_ind(10) = -9999, & 712 708 section(100,3), section_xy(100) = -9999, & … … 738 734 constant_top_momentumflux = .FALSE., & 739 735 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., & 742 737 data_output_2d_on_each_pe = .TRUE., & 743 738 dissipation_control = .FALSE., disturbance_created = .FALSE., & … … 779 774 averaging_interval = 0.0, averaging_interval_pr = 9999999.9, & 780 775 averaging_interval_sp = 9999999.9, bc_pt_t_val, bc_q_t_val, & 781 b c_qr_t_val, bc_nr_t_val, bottom_salinityflux = 0.0, &776 bottom_salinityflux = 0.0, & 782 777 building_height = 50.0, building_length_x = 50.0, & 783 778 building_length_y = 50.0, building_wall_left = 9999999.9, & … … 799 794 dt_do2d_xy = 9999999.9, dt_do2d_xz = 9999999.9, & 800 795 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, & 802 797 dt_run_control = 60.0, dt_3d = -1.0, dz = -1.0, & 803 798 dz_max = 9999999.9, dz_stretch_factor = 1.08, & … … 810 805 maximum_cpu_time_allowed = 0.0, & 811 806 molecular_viscosity = 1.461E-5, & 812 nr_surface = 0.0, nr_surface_initial_change = 0.0, &813 807 old_dt = 1.0E-10, omega = 7.29212E-5, omega_sor = 1.8, & 814 808 particle_maximum_age = 9999999.9, & … … 819 813 pt_surface = 300.0, pt_surface_initial_change = 0.0, & 820 814 q_surface = 0.0, q_surface_initial_change = 0.0, & 821 qr_surface = 0.0, qr_surface_initial_change = 0.0, &822 815 rayleigh_damping_factor = -1.0, rayleigh_damping_height = -1.0, & 823 816 recycling_width = 9999999.9, residual_limit = 1.0E-4, & … … 832 825 surface_heatflux = 9999999.9, surface_pressure = 1013.25, & 833 826 surface_scalarflux = 9999999.9, surface_waterflux = 9999999.9, & 834 surface_waterflux_nr = 0.0, surface_waterflux_qr = 0.0, &835 827 s_surface = 0.0, s_surface_initial_change = 0.0, & 836 828 termination_time_needed = -1.0, time_coupling = 0.0, & … … 854 846 lad_vertical_gradient_level(10) = -9999999.9, & 855 847 mask_scale(3), & 856 nr_vertical_gradient(10) = 0.0, &857 nr_vertical_gradient_level(10) = -1.0, &858 848 pt_vertical_gradient(10) = 0.0, & 859 849 pt_vertical_gradient_level(10) = -9999999.9, & 860 850 q_vertical_gradient(10) = 0.0, & 861 851 q_vertical_gradient_level(10) = -1.0, & 862 qr_vertical_gradient(10) = 0.0, &863 qr_vertical_gradient_level(10) = -1.0, &864 852 s_vertical_gradient(10) = 0.0, & 865 853 s_vertical_gradient_level(10) = -1.0, & … … 877 865 volume_flow(1:2) = 0.0, volume_flow_area(1:2) = 0.0, & 878 866 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, & 880 868 wall_qflux(0:4) = 0.0, wall_qrflux(0:4) = 0.0, & 881 869 wall_salinityflux(0:4) = 0.0, wall_scalarflux(0:4) = 0.0, &
Note: See TracChangeset
for help on using the changeset viewer.
