Ignore:
Timestamp:
Mar 21, 2014 11:00:16 AM (8 years ago)
Author:
raasch
Message:

Changed:


-s real64 removed (.mrun.config.hlrnIII)
-r8 removed (.mrun.config.imuk)
deleted: .mrun.config.imuk_ice2_netcdf4 .mrun.config.imuk_hlrn

REAL constants defined as wp-kind in modules

"baroclinicity" renamed "baroclinity", "ocean version" replaced by
"ocean mode"

code parts concerning old output formats "iso2d" and "avs" removed.
netCDF is the only remaining output format.

Errors:


bugfix: duplicate error message 56 removed

File:
1 edited

Legend:

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

    r1321 r1327  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! REAL constants defined as wp-kind
     23! -avs_output, data_output_format, do3d_compress, iso2d_output, netcdf_output
    2324!
    2425! Former revisions:
     
    347348               
    348349
    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. parametrization
    352                  a_term = 9.65,     & !: coef. for terminal velocity (m s-1)
    353                  a_vent = 0.78,     & !: coef. for ventilation effect
    354                  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. parametrization
    357                  b_term = 9.8,      & !: coef. for terminal velocity (m s-1)
    358                  b_vent = 0.308,    & !: coef. for ventilation effect
    359                  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)
    360361                 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. parametrization
    364                  c_const = 0.93,    & !: const. in Taylor-microscale Reynolds number
    365                  c_evap = 0.7,      & !: constant in evaporation
    366                  c_sedimentation = 2.0, & !: Courant number of sedimentation process
    367                  c_term = 600.0,    & !: coef. for terminal velocity (m-1)
    368                  cof(6) = (/ 76.18009172947146,      & !: coefficients in the
    369                              -86.50532032941677,     & !: numerical
    370                              24.01409824083091,      & !: calculation of the
    371                              -1.231739572450155,     & !: gamma function
    372                              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)
    376377                effective_coll_efficiency, & !:
    377                 eps_ros = 1.0E-4,  & !: accuracy of Rosenbrock method
    378                 eps_sb = 1.0E-20,  & !: threshold in two-moments scheme
    379                 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)
    387388                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 concentration
    392                 prec_time_const = 0.001, & !: coef. in Kessler scheme
     389                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
    393394                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 scheme
    396                 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 number
    399                 schmidt_p_1d3,     & !: schmidt**( 1.0 / 3.0 )
    400                 sigma_gc = 1.3,    & !: log-normal geometric standard deviation
    401                 stp = 2.5066282746310005, & !: parameter in gamma function
    402                 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 NaCl
    404                 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)
    409410
    410411    REAL(wp), DIMENSION(:), ALLOCATABLE     ::  hyrho, pt_d_t, t_d_pt 
     
    511512    CHARACTER (LEN=1000) ::  message_string = ' '
    512513
    513     CHARACTER (LEN=7),  DIMENSION(100) ::  do3d_comp_prec = ' '
    514     CHARACTER (LEN=10), DIMENSION(10)  ::  data_output_format = ' '
    515514    CHARACTER (LEN=11), DIMENSION(100) ::  data_output = ' ',    &
    516515                                           data_output_user = ' ', doav = ' '
     
    581580                    mask_i_global, mask_j_global, mask_k_global
    582581
    583     LOGICAL ::  avs_output = .FALSE., &
    584                 bc_lr_cyc =.TRUE., bc_lr_dirrad = .FALSE., &
     582    LOGICAL ::  bc_lr_cyc =.TRUE., bc_lr_dirrad = .FALSE., &
    585583                bc_lr_raddir = .FALSE., bc_ns_cyc = .TRUE., &
    586584                bc_ns_dirrad = .FALSE., bc_ns_raddir = .FALSE.,&
     
    595593                dissipation_control = .FALSE., disturbance_created = .FALSE., &
    596594                do2d_at_begin = .FALSE., do3d_at_begin = .FALSE., &
    597                 do3d_compress = .FALSE., do_sum = .FALSE., &
     595                do_sum = .FALSE., &
    598596                dp_external = .FALSE., dp_smooth = .FALSE., &
    599597                drizzle = .FALSE., dt_fixed = .FALSE., &
     
    604602                inflow_l = .FALSE., inflow_n = .FALSE., &
    605603                inflow_r = .FALSE., inflow_s = .FALSE., &
    606                 iso2d_output = .FALSE., large_scale_forcing = .FALSE., &
     604                large_scale_forcing = .FALSE., &
    607605                large_scale_subsidence = .FALSE., lsf_surf = .TRUE., &
    608606                lsf_vert = .TRUE., lptnudge = .FALSE., lqnudge = .FALSE., &
    609607                lunudge = .FALSE., lvnudge = .FALSE., lwnudge = .FALSE., &
    610608                masking_method = .FALSE., mg_switch_to_pe0 = .FALSE., &
    611                 netcdf_output = .FALSE., neutral = .FALSE., nudging = .FALSE., &
     609                neutral = .FALSE., nudging = .FALSE., &
    612610                ocean = .FALSE., on_device = .FALSE., &
    613611                outflow_l = .FALSE., outflow_n = .FALSE., outflow_r = .FALSE., &
     
    632630                data_output_yz(0:1) = .FALSE.
    633631
    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.0
    704 
    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, &
    709707                 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
    734732
    735733    REAL(wp), DIMENSION(:), ALLOCATABLE ::  dp_smooth_factor
    736734
    737735    REAL(wp), DIMENSION(max_masks,mask_xyz_dimension) :: &
    738             mask_x = -1.0, mask_y = -1.0, mask_z = -1.0
     736            mask_x = -1.0_wp, mask_y = -1.0_wp, mask_z = -1.0_wp
    739737    REAL(wp), DIMENSION(max_masks,3) ::                  &
    740             mask_x_loop = -1.0, mask_y_loop = -1.0, mask_z_loop = -1.0
     738            mask_x_loop = -1.0_wp, mask_y_loop = -1.0_wp, mask_z_loop = -1.0_wp
    741739   
    742740!
     
    787785                use_seperate_pe_for_dvrp_output = .FALSE.
    788786
    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.0
    793 
    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 /)
    799797
    800798#if defined( __decalpha )
    801799    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 /) )
    805803
    806804    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 /) )
    811809
    812810    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 ) /),         &
    816814                                            (/ 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,                       &
    819817                                   interval_values_dvrp_prt,                    &
    820818                                   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 ) /),     &
    823821                                                   (/ 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
    827825#else
    828826    REAL(wp), DIMENSION(2,10)     ::  slicer_range_limits_dvrp
     
    832830    REAL(sp), DIMENSION(2,100) ::  interval_values_dvrp,                       &
    833831                                   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 /
    854852#endif
    855853
     
    10241022    LOGICAL ::  run_control_header_1d = .FALSE., stop_dt_1d = .FALSE.
    10251023
    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, &
    10311029                    vsws1d, z01d, z0h1d
    10321030
     
    12331231    LOGICAL, DIMENSION(:), ALLOCATABLE ::  particle_mask, tail_mask
    12341232
    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_global
     1233    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
    12451243
    12461244    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.9
     1245                    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
    12511249
    12521250    REAL(wp), DIMENSION(:), ALLOCATABLE     ::  log_z_z0
     
    13931391
    13941392    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, &
    14001398                 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
    14061404
    14071405    SAVE
     
    14541452                     plot_spectra_level(100) = 999999
    14551453
    1456     REAL(wp)    ::  time_to_start_sp = 0.0
     1454    REAL(wp)    ::  time_to_start_sp = 0.0_wp
    14571455
    14581456    SAVE
Note: See TracChangeset for help on using the changeset viewer.