Changeset 1327 for palm/trunk/SOURCE/modules.f90
- Timestamp:
- Mar 21, 2014 11:00:16 AM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/modules.f90
r1321 r1327 20 20 ! Current revisions: 21 21 ! ------------------ 22 ! 22 ! REAL constants defined as wp-kind 23 ! -avs_output, data_output_format, do3d_compress, iso2d_output, netcdf_output 23 24 ! 24 25 ! Former revisions: … … 347 348 348 349 349 REAL(wp) :: a_1 = 8.69E-4 , & !: coef. in turb. parametrization (cm-2 s3)350 a_2 = -7.38E-5 , & !: coef. in turb. parametrization (cm-2 s3)351 a_3 = -1.40E-2 , & !: coef. in turb. parametrization352 a_term = 9.65 , & !: coef. for terminal velocity (m s-1)353 a_vent = 0.78 , & !: coef. for ventilation effect354 b_1 = 11.45E-6 , & !: coef. in turb. parametrization (m)355 b_2 = 9.68E-6 , & !: coef. in turb. parametrization (m)356 b_3 = 0.62 , & !: coef. in turb. parametrization357 b_term = 9.8 , & !: coef. for terminal velocity (m s-1)358 b_vent = 0.308 , & !: coef. for ventilation effect359 beta_cc = 3.09E-4 , & !: coef. in turb. parametrization (cm-2 s3)350 REAL(wp) :: a_1 = 8.69E-4_wp, & !: coef. in turb. parametrization (cm-2 s3) 351 a_2 = -7.38E-5_wp, & !: coef. in turb. parametrization (cm-2 s3) 352 a_3 = -1.40E-2_wp, & !: coef. in turb. parametrization 353 a_term = 9.65_wp, & !: coef. for terminal velocity (m s-1) 354 a_vent = 0.78_wp, & !: coef. for ventilation effect 355 b_1 = 11.45E-6_wp, & !: coef. in turb. parametrization (m) 356 b_2 = 9.68E-6_wp, & !: coef. in turb. parametrization (m) 357 b_3 = 0.62_wp, & !: coef. in turb. parametrization 358 b_term = 9.8_wp, & !: coef. for terminal velocity (m s-1) 359 b_vent = 0.308_wp, & !: coef. for ventilation effect 360 beta_cc = 3.09E-4_wp, & !: coef. in turb. parametrization (cm-2 s3) 360 361 bfactor, & 361 c_1 = 4.82E-6 , & !: coef. in turb. parametrization (m)362 c_2 = 4.8E-6 , & !: coef. in turb. parametrization (m)363 c_3 = 0.76 , & !: coef. in turb. parametrization364 c_const = 0.93 , & !: const. in Taylor-microscale Reynolds number365 c_evap = 0.7 , & !: constant in evaporation366 c_sedimentation = 2.0 , & !: Courant number of sedimentation process367 c_term = 600.0 , & !: coef. for terminal velocity (m-1)368 cof(6) = (/ 76.18009172947146 , & !: coefficients in the369 -86.50532032941677 , & !: numerical370 24.01409824083091 , & !: calculation of the371 -1.231739572450155 , & !: gamma function372 0.1208650973866179E-2 , &373 -0.5395239384953E-5 /), &374 cp = 1005.0 , & !: heat capacity of dry air (J kg-1 K-1)375 diff_coeff_l = 0.23E-4 , & !: diffusivity of water vapor (m2 s-1)362 c_1 = 4.82E-6_wp, & !: coef. in turb. parametrization (m) 363 c_2 = 4.8E-6_wp, & !: coef. in turb. parametrization (m) 364 c_3 = 0.76_wp, & !: coef. in turb. parametrization 365 c_const = 0.93_wp, & !: const. in Taylor-microscale Reynolds number 366 c_evap = 0.7_wp, & !: constant in evaporation 367 c_sedimentation = 2.0_wp, & !: Courant number of sedimentation process 368 c_term = 600.0_wp, & !: coef. for terminal velocity (m-1) 369 cof(6) = (/ 76.18009172947146_wp, & !: coefficients in the 370 -86.50532032941677_wp, & !: numerical 371 24.01409824083091_wp, & !: calculation of the 372 -1.231739572450155_wp, & !: gamma function 373 0.1208650973866179E-2_wp, & 374 -0.5395239384953E-5_wp /), & 375 cp = 1005.0_wp, & !: heat capacity of dry air (J kg-1 K-1) 376 diff_coeff_l = 0.23E-4_wp, & !: diffusivity of water vapor (m2 s-1) 376 377 effective_coll_efficiency, & !: 377 eps_ros = 1.0E-4 , & !: accuracy of Rosenbrock method378 eps_sb = 1.0E-20 , & !: threshold in two-moments scheme379 k_cc = 9.44E09 , & !: const. cloud-cloud kernel (m3 kg-2 s-1)380 k_cr0 = 4.33 , & !: const. cloud-rain kernel (m3 kg-1 s-1)381 k_rr = 7.12 , & !: const. rain-rain kernel (m3 kg-1 s-1)382 k_br = 1000. , & !: const. in breakup parametrization (m-1)383 k_st = 1.2E8 , & !: const. in drizzle parametrization (m-1 s-1)384 kappa_rr = 60.7 , & !: const. in collision kernel (kg-1/3)385 kin_vis_air = 1.4086E-5 , & !: kin. viscosity of air (m2 s-1)386 l_v = 2.5E+06 , & !: latent heat of vaporization (J kg-1)378 eps_ros = 1.0E-4_wp, & !: accuracy of Rosenbrock method 379 eps_sb = 1.0E-20_wp, & !: threshold in two-moments scheme 380 k_cc = 9.44E09_wp, & !: const. cloud-cloud kernel (m3 kg-2 s-1) 381 k_cr0 = 4.33_wp, & !: const. cloud-rain kernel (m3 kg-1 s-1) 382 k_rr = 7.12_wp, & !: const. rain-rain kernel (m3 kg-1 s-1) 383 k_br = 1000._wp, & !: const. in breakup parametrization (m-1) 384 k_st = 1.2E8_wp, & !: const. in drizzle parametrization (m-1 s-1) 385 kappa_rr = 60.7_wp, & !: const. in collision kernel (kg-1/3) 386 kin_vis_air = 1.4086E-5_wp, & !: kin. viscosity of air (m2 s-1) 387 l_v = 2.5E+06_wp, & !: latent heat of vaporization (J kg-1) 387 388 l_d_cp, l_d_r, l_d_rv, & !: l_v / cp, l_v / r_d, l_v / r_v 388 mass_of_solute = 1.0E-17 , & !: soluted NaCl (kg)389 molecular_weight_of_solute = 0.05844 , & !: mol. m. NaCl (kg mol-1)390 molecular_weight_of_water = 0.01801528 , & !: mol. m. H2O (kg mol-1)391 nc_const = 70.0E6 , & !: cloud droplet concentration392 prec_time_const = 0.001 , & !: coef. in Kessler scheme389 mass_of_solute = 1.0E-17_wp, & !: soluted NaCl (kg) 390 molecular_weight_of_solute = 0.05844_wp, & !: mol. m. NaCl (kg mol-1) 391 molecular_weight_of_water = 0.01801528_wp, & !: mol. m. H2O (kg mol-1) 392 nc_const = 70.0E6_wp, & !: cloud droplet concentration 393 prec_time_const = 0.001_wp, & !: coef. in Kessler scheme 393 394 pirho_l, dpirho_l, & !: pi * rho_l / 6.0; 6.0 / ( pi * rho_l ) 394 rho_l = 1.0E3 , & !: density of water (kg m-3)395 ql_crit = 0.0005 , & !: coef. in Kessler scheme396 r_d = 287.0 , & !: sp. gas const. dry air (J kg-1 K-1)397 r_v = 461.51 , & !: sp. gas const. water vapor (J kg-1 K-1)398 schmidt = 0.71 , & !: Schmidt number399 schmidt_p_1d3, & !: schmidt**( 1.0 / 3.0 )400 sigma_gc = 1.3 , & !: log-normal geometric standard deviation401 stp = 2.5066282746310005 , & !: parameter in gamma function402 thermal_conductivity_l = 2.43E-2 , & !: therm. cond. air (J m-1 s-1 K-1)403 vanthoff = 2.0 , & !: van't Hoff factor for NaCl404 x0 = 2.6E-10 , & !: separating drop mass (kg)405 xrmin = 2.6E-10 , & !: minimum rain drop size (kg)406 xrmax = 5.0E-6 , & !: maximum rain drop site (kg)407 dt_precipitation = 100.0 , & !: timestep precipitation (s)408 w_precipitation = 9.65 !: maximum terminal velocity (m s-1)395 rho_l = 1.0E3_wp, & !: density of water (kg m-3) 396 ql_crit = 0.0005_wp, & !: coef. in Kessler scheme 397 r_d = 287.0_wp, & !: sp. gas const. dry air (J kg-1 K-1) 398 r_v = 461.51_wp, & !: sp. gas const. water vapor (J kg-1 K-1) 399 schmidt = 0.71_wp, & !: Schmidt number 400 schmidt_p_1d3, & !: schmidt**( 1.0 / 3.0 ) 401 sigma_gc = 1.3_wp, & !: log-normal geometric standard deviation 402 stp = 2.5066282746310005_wp, & !: parameter in gamma function 403 thermal_conductivity_l = 2.43E-2_wp, & !: therm. cond. air (J m-1 s-1 K-1) 404 vanthoff = 2.0_wp, & !: van't Hoff factor for NaCl 405 x0 = 2.6E-10_wp, & !: separating drop mass (kg) 406 xrmin = 2.6E-10_wp, & !: minimum rain drop size (kg) 407 xrmax = 5.0E-6_wp, & !: maximum rain drop site (kg) 408 dt_precipitation = 100.0_wp, & !: timestep precipitation (s) 409 w_precipitation = 9.65_wp !: maximum terminal velocity (m s-1) 409 410 410 411 REAL(wp), DIMENSION(:), ALLOCATABLE :: hyrho, pt_d_t, t_d_pt … … 511 512 CHARACTER (LEN=1000) :: message_string = ' ' 512 513 513 CHARACTER (LEN=7), DIMENSION(100) :: do3d_comp_prec = ' '514 CHARACTER (LEN=10), DIMENSION(10) :: data_output_format = ' '515 514 CHARACTER (LEN=11), DIMENSION(100) :: data_output = ' ', & 516 515 data_output_user = ' ', doav = ' ' … … 581 580 mask_i_global, mask_j_global, mask_k_global 582 581 583 LOGICAL :: avs_output = .FALSE., & 584 bc_lr_cyc =.TRUE., bc_lr_dirrad = .FALSE., & 582 LOGICAL :: bc_lr_cyc =.TRUE., bc_lr_dirrad = .FALSE., & 585 583 bc_lr_raddir = .FALSE., bc_ns_cyc = .TRUE., & 586 584 bc_ns_dirrad = .FALSE., bc_ns_raddir = .FALSE.,& … … 595 593 dissipation_control = .FALSE., disturbance_created = .FALSE., & 596 594 do2d_at_begin = .FALSE., do3d_at_begin = .FALSE., & 597 do 3d_compress = .FALSE., do_sum = .FALSE., &595 do_sum = .FALSE., & 598 596 dp_external = .FALSE., dp_smooth = .FALSE., & 599 597 drizzle = .FALSE., dt_fixed = .FALSE., & … … 604 602 inflow_l = .FALSE., inflow_n = .FALSE., & 605 603 inflow_r = .FALSE., inflow_s = .FALSE., & 606 iso2d_output = .FALSE.,large_scale_forcing = .FALSE., &604 large_scale_forcing = .FALSE., & 607 605 large_scale_subsidence = .FALSE., lsf_surf = .TRUE., & 608 606 lsf_vert = .TRUE., lptnudge = .FALSE., lqnudge = .FALSE., & 609 607 lunudge = .FALSE., lvnudge = .FALSE., lwnudge = .FALSE., & 610 608 masking_method = .FALSE., mg_switch_to_pe0 = .FALSE., & 611 ne tcdf_output = .FALSE., neutral = .FALSE., nudging = .FALSE., &609 neutral = .FALSE., nudging = .FALSE., & 612 610 ocean = .FALSE., on_device = .FALSE., & 613 611 outflow_l = .FALSE., outflow_n = .FALSE., outflow_r = .FALSE., & … … 632 630 data_output_yz(0:1) = .FALSE. 633 631 634 REAL(wp) :: advected_distance_x = 0.0 , advected_distance_y = 0.0, &635 alpha_surface = 0.0 , atmos_ocean_sign = 1.0, &636 averaging_interval = 0.0 , averaging_interval_pr = 9999999.9, &637 averaging_interval_sp = 9999999.9 , bc_pt_t_val, bc_q_t_val, &638 bottom_salinityflux = 0.0 , &639 building_height = 50.0 , building_length_x = 50.0, &640 building_length_y = 50.0 , building_wall_left = 9999999.9, &641 building_wall_south = 9999999.9 , canyon_height = 50.0, &642 canyon_width_x = 9999999.9 , canyon_width_y = 9999999.9, &643 canyon_wall_left = 9999999.9 , canyon_wall_south = 9999999.9, &644 cthf = 0.0 , cfl_factor = -1.0, cos_alpha_surface, &645 coupling_start_time = 0.0 , disturbance_amplitude = 0.25, &646 disturbance_energy_limit = 0.01 , &647 disturbance_level_b = -9999999.9 , &648 disturbance_level_t = -9999999.9 , &649 dp_level_b = 0.0 , drag_coefficient = 0.0, &650 dt = -1.0 , dt_averaging_input = 0.0, &651 dt_averaging_input_pr = 9999999.9 , dt_coupling = 9999999.9, &652 dt_data_output = 9999999.9 , &653 dt_data_output_av = 9999999.9 , dt_disturb = 9999999.9, &654 dt_dopr = 9999999.9 , dt_dopr_listing = 9999999.9, &655 dt_dopts = 9999999.9 , dt_dosp = 9999999.9, dt_dots = 9999999.9, &656 dt_do2d_xy = 9999999.9 , dt_do2d_xz = 9999999.9, &657 dt_do2d_yz = 9999999.9 , dt_do3d = 9999999.9, dt_dvrp = 9999999.9, &658 dt_max = 20.0 , dt_micro = -1.0, dt_restart = 9999999.9, &659 dt_run_control = 60.0 , dt_3d = -1.0, dz = -1.0, &660 dz_max = 9999999.9 , dz_stretch_factor = 1.08, &661 dz_stretch_level = 100000.0 , e_init = 0.0, e_min = 0.0, &662 end_time = 0.0 , &663 f = 0.0 , fs = 0.0, g = 9.81, inflow_damping_height = 9999999.9, &664 inflow_damping_width = 9999999.9 , kappa = 0.4, km_constant = -1.0,&665 lad_surface = 0.0 , leaf_surface_concentration = 0.0, &666 mask_scale_x = 1.0 , mask_scale_y = 1.0, mask_scale_z = 1.0, &667 maximum_cpu_time_allowed = 0.0 , &668 molecular_viscosity = 1.461E-5 , &669 old_dt = 1.0E-10 , omega = 7.29212E-5, omega_sor = 1.8, &670 particle_maximum_age = 9999999.9 , &671 phi = 55.0 , prandtl_number = 1.0, &672 precipitation_amount_interval = 9999999.9 , prho_reference, &673 pt_damping_factor = 0.0 , pt_damping_width = 0.0, &674 pt_reference = 9999999.9 , pt_slope_offset = 0.0, &675 pt_surface = 300.0 , pt_surface_initial_change = 0.0, &676 q_surface = 0.0 , q_surface_initial_change = 0.0, &677 rayleigh_damping_factor = -1.0 , rayleigh_damping_height = -1.0, &678 recycling_width = 9999999.9 , residual_limit = 1.0E-4, &679 restart_time = 9999999.9 , rho_reference, rho_surface, &680 rif_max = 1.0 , rif_min = -5.0, roughness_length = 0.1, &681 sa_surface = 35.0 , scalar_exchange_coefficient = 0.0, &682 simulated_time = 0.0 , simulated_time_at_begin, sin_alpha_surface, &683 skip_time_data_output = 0.0 , skip_time_data_output_av = 9999999.9,&684 skip_time_dopr = 9999999.9 , skip_time_dosp = 9999999.9, &685 skip_time_do2d_xy = 9999999.9 , skip_time_do2d_xz = 9999999.9, &686 skip_time_do2d_yz = 9999999.9 , skip_time_do3d = 9999999.9, &687 surface_heatflux = 9999999.9 , surface_pressure = 1013.25, &688 surface_scalarflux = 9999999.9 , surface_waterflux = 9999999.9, &689 s_surface = 0.0 , s_surface_initial_change = 0.0, &690 termination_time_needed = -1.0 , time_coupling = 0.0, &691 time_disturb = 0.0 , time_dopr = 0.0, time_dopr_av = 0.0, &692 time_dopr_listing = 0.0 , time_dopts = 0.0, time_dosp = 0.0, &693 time_dosp_av = 0.0 , time_dots = 0.0, time_do2d_xy = 0.0, &694 time_do2d_xz = 0.0 , time_do2d_yz = 0.0, time_do3d = 0.0, &695 time_do_av = 0.0 , time_do_sla = 0.0, time_dvrp = 0.0, &696 time_restart = 9999999.9 , time_run_control = 0.0,&697 time_since_reference_point, top_heatflux = 9999999.9 , &698 top_momentumflux_u = 9999999.9 , &699 top_momentumflux_v = 9999999.9 , top_salinityflux = 9999999.9, &700 ug_surface = 0.0 , u_bulk = 0.0, u_gtrans = 0.0, &701 vg_surface = 0.0 , vpt_reference = 9999999.9, &702 v_bulk = 0.0 , v_gtrans = 0.0, wall_adjustment_factor = 1.8, &703 z_max_do2d = -1.0 , z0h_factor = 1.0704 705 REAL(wp) :: do2d_xy_last_time(0:1) = -1.0 , do2d_xz_last_time(0:1) = -1.0, &706 do2d_yz_last_time(0:1) = -1.0 , dpdxy(1:2) = 0.0, &707 dt_domask(max_masks) = 9999999.9 , lad_vertical_gradient(10) = 0.0,&708 lad_vertical_gradient_level(10) = -9999999.9 , &632 REAL(wp) :: advected_distance_x = 0.0_wp, advected_distance_y = 0.0_wp, & 633 alpha_surface = 0.0_wp, atmos_ocean_sign = 1.0_wp, & 634 averaging_interval = 0.0_wp, averaging_interval_pr = 9999999.9_wp, & 635 averaging_interval_sp = 9999999.9_wp, bc_pt_t_val, bc_q_t_val, & 636 bottom_salinityflux = 0.0_wp, & 637 building_height = 50.0_wp, building_length_x = 50.0_wp, & 638 building_length_y = 50.0_wp, building_wall_left = 9999999.9_wp, & 639 building_wall_south = 9999999.9_wp, canyon_height = 50.0_wp, & 640 canyon_width_x = 9999999.9_wp, canyon_width_y = 9999999.9_wp, & 641 canyon_wall_left = 9999999.9_wp, canyon_wall_south = 9999999.9_wp, & 642 cthf = 0.0_wp, cfl_factor = -1.0_wp, cos_alpha_surface, & 643 coupling_start_time = 0.0_wp, disturbance_amplitude = 0.25_wp, & 644 disturbance_energy_limit = 0.01_wp, & 645 disturbance_level_b = -9999999.9_wp, & 646 disturbance_level_t = -9999999.9_wp, & 647 dp_level_b = 0.0_wp, drag_coefficient = 0.0_wp, & 648 dt = -1.0_wp, dt_averaging_input = 0.0_wp, & 649 dt_averaging_input_pr = 9999999.9_wp, dt_coupling = 9999999.9_wp, & 650 dt_data_output = 9999999.9_wp, & 651 dt_data_output_av = 9999999.9_wp, dt_disturb = 9999999.9_wp, & 652 dt_dopr = 9999999.9_wp, dt_dopr_listing = 9999999.9_wp, & 653 dt_dopts = 9999999.9_wp, dt_dosp = 9999999.9_wp, dt_dots = 9999999.9_wp, & 654 dt_do2d_xy = 9999999.9_wp, dt_do2d_xz = 9999999.9_wp, & 655 dt_do2d_yz = 9999999.9_wp, dt_do3d = 9999999.9_wp, dt_dvrp = 9999999.9_wp, & 656 dt_max = 20.0_wp, dt_micro = -1.0_wp, dt_restart = 9999999.9_wp, & 657 dt_run_control = 60.0_wp, dt_3d = -1.0_wp, dz = -1.0_wp, & 658 dz_max = 9999999.9_wp, dz_stretch_factor = 1.08_wp, & 659 dz_stretch_level = 100000.0_wp, e_init = 0.0_wp, e_min = 0.0_wp, & 660 end_time = 0.0_wp, & 661 f = 0.0_wp, fs = 0.0_wp, g = 9.81_wp, inflow_damping_height = 9999999.9_wp, & 662 inflow_damping_width = 9999999.9_wp, kappa = 0.4_wp, km_constant = -1.0_wp,& 663 lad_surface = 0.0_wp, leaf_surface_concentration = 0.0_wp, & 664 mask_scale_x = 1.0_wp, mask_scale_y = 1.0_wp, mask_scale_z = 1.0_wp, & 665 maximum_cpu_time_allowed = 0.0_wp, & 666 molecular_viscosity = 1.461E-5_wp, & 667 old_dt = 1.0E-10_wp, omega = 7.29212E-5_wp, omega_sor = 1.8_wp, & 668 particle_maximum_age = 9999999.9_wp, & 669 phi = 55.0_wp, prandtl_number = 1.0_wp, & 670 precipitation_amount_interval = 9999999.9_wp, prho_reference, & 671 pt_damping_factor = 0.0_wp, pt_damping_width = 0.0_wp, & 672 pt_reference = 9999999.9_wp, pt_slope_offset = 0.0_wp, & 673 pt_surface = 300.0_wp, pt_surface_initial_change = 0.0_wp, & 674 q_surface = 0.0_wp, q_surface_initial_change = 0.0_wp, & 675 rayleigh_damping_factor = -1.0_wp, rayleigh_damping_height = -1.0_wp, & 676 recycling_width = 9999999.9_wp, residual_limit = 1.0E-4_wp, & 677 restart_time = 9999999.9_wp, rho_reference, rho_surface, & 678 rif_max = 1.0_wp, rif_min = -5.0_wp, roughness_length = 0.1_wp, & 679 sa_surface = 35.0_wp, scalar_exchange_coefficient = 0.0_wp, & 680 simulated_time = 0.0_wp, simulated_time_at_begin, sin_alpha_surface, & 681 skip_time_data_output = 0.0_wp, skip_time_data_output_av = 9999999.9_wp,& 682 skip_time_dopr = 9999999.9_wp, skip_time_dosp = 9999999.9_wp, & 683 skip_time_do2d_xy = 9999999.9_wp, skip_time_do2d_xz = 9999999.9_wp, & 684 skip_time_do2d_yz = 9999999.9_wp, skip_time_do3d = 9999999.9_wp, & 685 surface_heatflux = 9999999.9_wp, surface_pressure = 1013.25_wp, & 686 surface_scalarflux = 9999999.9_wp, surface_waterflux = 9999999.9_wp, & 687 s_surface = 0.0_wp, s_surface_initial_change = 0.0_wp, & 688 termination_time_needed = -1.0_wp, time_coupling = 0.0_wp, & 689 time_disturb = 0.0_wp, time_dopr = 0.0_wp, time_dopr_av = 0.0_wp, & 690 time_dopr_listing = 0.0_wp, time_dopts = 0.0_wp, time_dosp = 0.0_wp, & 691 time_dosp_av = 0.0_wp, time_dots = 0.0_wp, time_do2d_xy = 0.0_wp, & 692 time_do2d_xz = 0.0_wp, time_do2d_yz = 0.0_wp, time_do3d = 0.0_wp, & 693 time_do_av = 0.0_wp, time_do_sla = 0.0_wp, time_dvrp = 0.0_wp, & 694 time_restart = 9999999.9_wp, time_run_control = 0.0_wp,& 695 time_since_reference_point, top_heatflux = 9999999.9_wp, & 696 top_momentumflux_u = 9999999.9_wp, & 697 top_momentumflux_v = 9999999.9_wp, top_salinityflux = 9999999.9_wp, & 698 ug_surface = 0.0_wp, u_bulk = 0.0_wp, u_gtrans = 0.0_wp, & 699 vg_surface = 0.0_wp, vpt_reference = 9999999.9_wp, & 700 v_bulk = 0.0_wp, v_gtrans = 0.0_wp, wall_adjustment_factor = 1.8_wp, & 701 z_max_do2d = -1.0_wp, z0h_factor = 1.0_wp 702 703 REAL(wp) :: do2d_xy_last_time(0:1) = -1.0_wp, do2d_xz_last_time(0:1) = -1.0_wp, & 704 do2d_yz_last_time(0:1) = -1.0_wp, dpdxy(1:2) = 0.0_wp, & 705 dt_domask(max_masks) = 9999999.9_wp, lad_vertical_gradient(10) = 0.0_wp,& 706 lad_vertical_gradient_level(10) = -9999999.9_wp, & 709 707 mask_scale(3), & 710 pt_vertical_gradient(10) = 0.0 , &711 pt_vertical_gradient_level(10) = -9999999.9 , &712 q_vertical_gradient(10) = 0.0 , &713 q_vertical_gradient_level(10) = -1.0 , &714 s_vertical_gradient(10) = 0.0 , &715 s_vertical_gradient_level(10) = -1.0 , &716 sa_vertical_gradient(10) = 0.0 , &717 sa_vertical_gradient_level(10) = -9999999.9 , &718 skip_time_domask(max_masks) = 9999999.9 , threshold(20) = 0.0, &719 time_domask(max_masks) = 0.0 , &720 tsc(10) = (/ 1.0 , 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0/), &721 u_profile(100) = 9999999.9 , uv_heights(100) = 9999999.9, &722 v_profile(100) = 9999999.9 , &723 ug_vertical_gradient(10) = 0.0 , &724 ug_vertical_gradient_level(10) = -9999999.9 , &725 vg_vertical_gradient(10) = 0.0 , &726 vg_vertical_gradient_level(10) = -9999999.9 , &727 volume_flow(1:2) = 0.0 , volume_flow_area(1:2) = 0.0, &728 volume_flow_initial(1:2) = 0.0 , wall_heatflux(0:4) = 0.0, &729 wall_humidityflux(0:4) = 0.0 , wall_nrflux(0:4) = 0.0, &730 wall_qflux(0:4) = 0.0 , wall_qrflux(0:4) = 0.0, &731 wall_salinityflux(0:4) = 0.0 , wall_scalarflux(0:4) = 0.0, &732 subs_vertical_gradient(10) = 0.0 , &733 subs_vertical_gradient_level(10) = -9999999.9 708 pt_vertical_gradient(10) = 0.0_wp, & 709 pt_vertical_gradient_level(10) = -9999999.9_wp, & 710 q_vertical_gradient(10) = 0.0_wp, & 711 q_vertical_gradient_level(10) = -1.0_wp, & 712 s_vertical_gradient(10) = 0.0_wp, & 713 s_vertical_gradient_level(10) = -1.0_wp, & 714 sa_vertical_gradient(10) = 0.0_wp, & 715 sa_vertical_gradient_level(10) = -9999999.9_wp, & 716 skip_time_domask(max_masks) = 9999999.9_wp, threshold(20) = 0.0_wp, & 717 time_domask(max_masks) = 0.0_wp, & 718 tsc(10) = (/ 1.0_wp, 1.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp /), & 719 u_profile(100) = 9999999.9_wp, uv_heights(100) = 9999999.9_wp, & 720 v_profile(100) = 9999999.9_wp, & 721 ug_vertical_gradient(10) = 0.0_wp, & 722 ug_vertical_gradient_level(10) = -9999999.9_wp, & 723 vg_vertical_gradient(10) = 0.0_wp, & 724 vg_vertical_gradient_level(10) = -9999999.9_wp, & 725 volume_flow(1:2) = 0.0_wp, volume_flow_area(1:2) = 0.0_wp, & 726 volume_flow_initial(1:2) = 0.0_wp, wall_heatflux(0:4) = 0.0_wp, & 727 wall_humidityflux(0:4) = 0.0_wp, wall_nrflux(0:4) = 0.0_wp, & 728 wall_qflux(0:4) = 0.0_wp, wall_qrflux(0:4) = 0.0_wp, & 729 wall_salinityflux(0:4) = 0.0_wp, wall_scalarflux(0:4) = 0.0_wp, & 730 subs_vertical_gradient(10) = 0.0_wp, & 731 subs_vertical_gradient_level(10) = -9999999.9_wp 734 732 735 733 REAL(wp), DIMENSION(:), ALLOCATABLE :: dp_smooth_factor 736 734 737 735 REAL(wp), DIMENSION(max_masks,mask_xyz_dimension) :: & 738 mask_x = -1.0 , mask_y = -1.0, mask_z = -1.0736 mask_x = -1.0_wp, mask_y = -1.0_wp, mask_z = -1.0_wp 739 737 REAL(wp), DIMENSION(max_masks,3) :: & 740 mask_x_loop = -1.0 , mask_y_loop = -1.0, mask_z_loop = -1.0738 mask_x_loop = -1.0_wp, mask_y_loop = -1.0_wp, mask_z_loop = -1.0_wp 741 739 742 740 ! … … 787 785 use_seperate_pe_for_dvrp_output = .FALSE. 788 786 789 REAL(wp) :: clip_dvrp_l = 9999999.9 , clip_dvrp_n = 9999999.9, &790 clip_dvrp_r = 9999999.9 , clip_dvrp_s = 9999999.9, &791 superelevation = 1.0 , superelevation_x = 1.0, &792 superelevation_y = 1.0 , vc_alpha = 38.0793 794 REAL(wp), DIMENSION(2) :: color_interval = (/ 0.0 , 1.0/), &795 dvrpsize_interval = (/ 0.0 , 1.0/)796 797 REAL(wp), DIMENSION(3) :: groundplate_color = (/ 0.0 , 0.6, 0.0/), &798 topography_color = (/ 0.8 , 0.7, 0.6/)787 REAL(wp) :: clip_dvrp_l = 9999999.9_wp, clip_dvrp_n = 9999999.9_wp, & 788 clip_dvrp_r = 9999999.9_wp, clip_dvrp_s = 9999999.9_wp, & 789 superelevation = 1.0_wp, superelevation_x = 1.0_wp, & 790 superelevation_y = 1.0_wp, vc_alpha = 38.0_wp 791 792 REAL(wp), DIMENSION(2) :: color_interval = (/ 0.0_wp, 1.0_wp /), & 793 dvrpsize_interval = (/ 0.0_wp, 1.0_wp /) 794 795 REAL(wp), DIMENSION(3) :: groundplate_color = (/ 0.0_wp, 0.6_wp, 0.0_wp /), & 796 topography_color = (/ 0.8_wp, 0.7_wp, 0.6_wp /) 799 797 800 798 #if defined( __decalpha ) 801 799 REAL(wp), DIMENSION(2,10) :: slicer_range_limits_dvrp = RESHAPE( (/ & 802 -1.0 , 1.0, -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, &803 -1.0 , 1.0, -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, &804 -1.0 , 1.0, -1.0, 1.0/), (/ 2, 10 /) )800 -1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, & 801 -1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, & 802 -1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp /), (/ 2, 10 /) ) 805 803 806 804 REAL(wp), DIMENSION(3,10) :: isosurface_color = RESHAPE( (/ & 807 0.9 , 0.9, 0.9, 0.8, 0.1, 0.1, 0.1, 0.1, 0.8, &808 0.1 , 0.8, 0.1, 0.6, 0.1, 0.1, 0.1, 0.1, 0.6, &809 0.1 , 0.6, 0.1, 0.4, 0.1, 0.1, 0.1, 0.1, 0.4, &810 0.1 , 0.4, 0.1/), (/ 3, 10 /) )805 0.9_wp, 0.9_wp, 0.9_wp, 0.8_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.8_wp, & 806 0.1_wp, 0.8_wp, 0.1_wp, 0.6_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.6_wp, & 807 0.1_wp, 0.6_wp, 0.1_wp, 0.4_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.4_wp, & 808 0.1_wp, 0.4_wp, 0.1_wp /), (/ 3, 10 /) ) 811 809 812 810 REAL(sp), DIMENSION(2,100) :: interval_values_dvrp, interval_h_dvrp = & 813 RESHAPE( (/ 270.0 , 225.0, 225.0, 180.0, &814 70.0 , 25.0, 25.0, -25.0, &815 ( 0.0 , i9 = 1, 192 ) /), &811 RESHAPE( (/ 270.0_sp, 225.0_sp, 225.0_sp, 180.0_sp, & 812 70.0_sp, 25.0_sp, 25.0_sp, -25.0_sp, & 813 ( 0.0_sp, i9 = 1, 192 ) /), & 816 814 (/ 2, 100 /) ), & 817 interval_l_dvrp = 0.5 , interval_s_dvrp = 1.0,&818 interval_a_dvrp = 0.0 , &815 interval_l_dvrp = 0.5_sp, interval_s_dvrp = 1.0_sp,& 816 interval_a_dvrp = 0.0_sp, & 819 817 interval_values_dvrp_prt, & 820 818 interval_h_dvrp_prt = RESHAPE( & 821 (/ 270.0 , 225.0, 225.0, 180.0, 70.0, 25.0, &822 25.0 , -25.0, ( 0.0, i9 = 1, 192 ) /), &819 (/ 270.0_sp, 225.0_sp, 225.0_sp, 180.0_sp, 70.0_sp, 25.0_sp, & 820 25.0_sp, -25.0_sp, ( 0.0_sp, i9 = 1, 192 ) /), & 823 821 (/ 2, 100 /) ), & 824 interval_l_dvrp_prt = 0.5 , &825 interval_s_dvrp_prt = 1.0 , &826 interval_a_dvrp_prt = 0.0 822 interval_l_dvrp_prt = 0.5_sp, & 823 interval_s_dvrp_prt = 1.0_sp, & 824 interval_a_dvrp_prt = 0.0_sp 827 825 #else 828 826 REAL(wp), DIMENSION(2,10) :: slicer_range_limits_dvrp … … 832 830 REAL(sp), DIMENSION(2,100) :: interval_values_dvrp, & 833 831 interval_values_dvrp_prt, interval_h_dvrp, & 834 interval_h_dvrp_prt, interval_l_dvrp = 0.5 , &835 interval_l_dvrp_prt = 0.5 , interval_s_dvrp = 1.0, &836 interval_s_dvrp_prt = 1.0 , interval_a_dvrp = 0.0, &837 interval_a_dvrp_prt = 0.0 838 839 DATA slicer_range_limits_dvrp / -1.0 , 1.0, -1.0, 1.0, -1.0, 1.0, &840 -1.0 , 1.0, -1.0, 1.0, -1.0, 1.0, &841 -1.0 , 1.0, -1.0, 1.0, -1.0, 1.0, &842 -1.0 , 1.0/843 844 DATA isosurface_color / 0.9 , 0.9, 0.9, 0.8, 0.1, 0.1, 0.1, 0.1, 0.8, &845 0.1 , 0.8, 0.1, 0.6, 0.1, 0.1, 0.1, 0.1, 0.6, &846 0.1 , 0.6, 0.1, 0.4, 0.1, 0.1, 0.1, 0.1, 0.4, &847 0.1 , 0.4, 0.1/848 849 DATA interval_h_dvrp / 270.0 , 225.0, 225.0, 180.0, 70.0, 25.0, &850 25.0 , -25.0, 192 * 0.0/851 852 DATA interval_h_dvrp_prt / 270.0 , 225.0, 225.0, 180.0, 70.0, 25.0, &853 25.0 , -25.0, 192 * 0.0/832 interval_h_dvrp_prt, interval_l_dvrp = 0.5_sp, & 833 interval_l_dvrp_prt = 0.5_sp, interval_s_dvrp = 1.0_sp, & 834 interval_s_dvrp_prt = 1.0_sp, interval_a_dvrp = 0.0_sp, & 835 interval_a_dvrp_prt = 0.0_sp 836 837 DATA slicer_range_limits_dvrp / -1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, & 838 -1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, & 839 -1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, & 840 -1.0_wp, 1.0_wp / 841 842 DATA isosurface_color / 0.9_wp, 0.9_wp, 0.9_wp, 0.8_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.8_wp, & 843 0.1_wp, 0.8_wp, 0.1_wp, 0.6_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.6_wp, & 844 0.1_wp, 0.6_wp, 0.1_wp, 0.4_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.4_wp, & 845 0.1_wp, 0.4_wp, 0.1_wp / 846 847 DATA interval_h_dvrp / 270.0_wp, 225.0_wp, 225.0_wp, 180.0_wp, 70.0_wp, 25.0_wp, & 848 25.0_wp, -25.0_wp, 192 * 0.0_wp / 849 850 DATA interval_h_dvrp_prt / 270.0_wp, 225.0_wp, 225.0_wp, 180.0_wp, 70.0_wp, 25.0_wp, & 851 25.0_wp, -25.0_wp, 192 * 0.0_wp / 854 852 #endif 855 853 … … 1024 1022 LOGICAL :: run_control_header_1d = .FALSE., stop_dt_1d = .FALSE. 1025 1023 1026 REAL(wp) :: damp_level_1d = -1.0 , dt_1d = 60.0, dt_max_1d = 300.0, &1027 dt_pr_1d = 9999999.9 , dt_run_control_1d = 60.0, &1028 end_time_1d = 864000.0 , old_dt_1d = 1.0E-10, &1029 qs1d, simulated_time_1d = 0.0 , time_pr_1d = 0.0, &1030 time_run_control_1d = 0.0 , ts1d, us1d, usws1d, &1024 REAL(wp) :: damp_level_1d = -1.0_wp, dt_1d = 60.0_wp, dt_max_1d = 300.0_wp, & 1025 dt_pr_1d = 9999999.9_wp, dt_run_control_1d = 60.0_wp, & 1026 end_time_1d = 864000.0_wp, old_dt_1d = 1.0E-10_wp, & 1027 qs1d, simulated_time_1d = 0.0_wp, time_pr_1d = 0.0_wp, & 1028 time_run_control_1d = 0.0_wp, ts1d, us1d, usws1d, & 1031 1029 vsws1d, z01d, z0h1d 1032 1030 … … 1233 1231 LOGICAL, DIMENSION(:), ALLOCATABLE :: particle_mask, tail_mask 1234 1232 1235 REAL(wp) :: c_0 = 3.0 , dt_min_part = 0.0002, dt_prel = 9999999.9, &1236 dt_sort_particles = 0.0 , dt_write_particle_data = 9999999.9, &1237 dvrp_psize = 9999999.9 , end_time_prel = 9999999.9, &1238 initial_weighting_factor = 1.0 , &1239 maximum_tailpoint_age = 100000.0 , &1240 minimum_tailpoint_distance = 0.0 , &1241 particle_advection_start = 0.0 , sgs_wfu_part = 0.3333333, &1242 sgs_wfv_part = 0.3333333 , sgs_wfw_part = 0.3333333, &1243 time_prel = 0.0 , time_sort_particles = 0.0, &1244 time_write_particle_data = 0.0 , z0_av_global1233 REAL(wp) :: c_0 = 3.0_wp, dt_min_part = 0.0002_wp, dt_prel = 9999999.9_wp, & 1234 dt_sort_particles = 0.0_wp, dt_write_particle_data = 9999999.9_wp, & 1235 dvrp_psize = 9999999.9_wp, end_time_prel = 9999999.9_wp, & 1236 initial_weighting_factor = 1.0_wp, & 1237 maximum_tailpoint_age = 100000.0_wp, & 1238 minimum_tailpoint_distance = 0.0_wp, & 1239 particle_advection_start = 0.0_wp, sgs_wfu_part = 0.3333333_wp, & 1240 sgs_wfv_part = 0.3333333_wp, sgs_wfw_part = 0.3333333_wp, & 1241 time_prel = 0.0_wp, time_sort_particles = 0.0_wp, & 1242 time_write_particle_data = 0.0_wp, z0_av_global 1245 1243 1246 1244 REAL(wp), DIMENSION(max_number_of_particle_groups) :: & 1247 density_ratio = 9999999.9 , pdx = 9999999.9, pdy = 9999999.9, &1248 pdz = 9999999.9 , psb = 9999999.9, psl = 9999999.9, &1249 psn = 9999999.9 , psr = 9999999.9, pss = 9999999.9, &1250 pst = 9999999.9 , radius = 9999999.91245 density_ratio = 9999999.9_wp, pdx = 9999999.9_wp, pdy = 9999999.9_wp, & 1246 pdz = 9999999.9_wp, psb = 9999999.9_wp, psl = 9999999.9_wp, & 1247 psn = 9999999.9_wp, psr = 9999999.9_wp, pss = 9999999.9_wp, & 1248 pst = 9999999.9_wp, radius = 9999999.9_wp 1251 1249 1252 1250 REAL(wp), DIMENSION(:), ALLOCATABLE :: log_z_z0 … … 1393 1391 1394 1392 REAL(wp) :: cross_ts_uymax(20) = & 1395 (/ 999.999 , 999.999, 999.999, 999.999, 999.999, &1396 999.999 , 999.999, 999.999, 999.999, 999.999, &1397 999.999 , 999.999, 999.999, 999.999, 999.999, &1398 999.999 , 999.999, 999.999, 999.999, 999.999/),&1399 cross_ts_uymax_computed(20) = 999.999 , &1393 (/ 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, & 1394 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, & 1395 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, & 1396 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp /),& 1397 cross_ts_uymax_computed(20) = 999.999_wp, & 1400 1398 cross_ts_uymin(20) = & 1401 (/ 999.999 , 999.999, 999.999, -5.000, 999.999, &1402 999.999 , 0.000, 999.999, 999.999, 999.999, &1403 999.999 , 999.999, 999.999, 999.999, 999.999, &1404 999.999 , 999.999, 999.999, 999.999, 999.999/),&1405 cross_ts_uymin_computed(20) = 999.999 1399 (/ 999.999_wp, 999.999_wp, 999.999_wp, -5.000_wp, 999.999_wp, & 1400 999.999_wp, 0.000_wp, 999.999_wp, 999.999_wp, 999.999_wp, & 1401 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, & 1402 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp /),& 1403 cross_ts_uymin_computed(20) = 999.999_wp 1406 1404 1407 1405 SAVE … … 1454 1452 plot_spectra_level(100) = 999999 1455 1453 1456 REAL(wp) :: time_to_start_sp = 0.0 1454 REAL(wp) :: time_to_start_sp = 0.0_wp 1457 1455 1458 1456 SAVE
Note: See TracChangeset
for help on using the changeset viewer.