Changeset 4565 for palm


Ignore:
Timestamp:
Jun 15, 2020 8:30:38 AM (6 months ago)
Author:
oliver.maas
Message:

added new surface temperature forcing method for bc_pt_b = 'dirichlet': surface temperature pt(0) can be linearly increased by pt_surface_heating_rate (in K/h)

Location:
palm/trunk/SOURCE
Files:
4 edited

Legend:

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

    r4564 r4565  
    2424! -----------------
    2525! $Id$
     26! added check for pt_surface_heating_rate
     27!
     28! 4564 2020-06-12 14:03:36Z raasch
    2629! Vertical nesting method of Huq et al. (2019) removed
    2730!
     
    13841387       CALL message( 'check_parameters', 'PA0066', 1, 2, 0, 6, 0 )
    13851388    ENDIF
     1389    IF ( constant_heatflux  .AND.  pt_surface_heating_rate /= 0.0_wp )  THEN
     1390       WRITE ( message_string, * )  'constant_heatflux = .TRUE. is not allo',                      &
     1391               'wed with pt_surface_heating_rate (/=0) = ', pt_surface_heating_rate
     1392       CALL message( 'check_parameters', 'PA0497', 1, 2, 0, 6, 0 )
     1393    ENDIF
    13861394
    13871395!
  • palm/trunk/SOURCE/modules.f90

    r4536 r4565  
    2525! -----------------
    2626! $Id$
     27! added pt_surface_heating_rate
     28!
     29! 4536 2020-05-17 17:24:13Z raasch
    2730! +qsurf_av, ssurf_av
    2831!
     
    942945                                                               !< boundary of total domain
    943946    REAL(wp) ::  pt_surface = 300.0_wp                         !< namelist parameter
     947    REAL(wp) ::  pt_surface_heating_rate = 0.0_wp              !< namelist parameter
    944948    REAL(wp) ::  pt_surface_initial_change = 0.0_wp            !< namelist parameter
    945949    REAL(wp) ::  q_surface = 0.0_wp                            !< namelist parameter
  • palm/trunk/SOURCE/parin.f90

    r4564 r4565  
    2525! -----------------
    2626! $Id$
     27! added pt_surface_heating_rate
     28!
     29! 4564 2020-06-12 14:03:36Z raasch
    2730! Vertical nesting method of Huq et al. (2019) removed
    2831!
     
    195198             nsor, nsor_ini, nudging, nx, ny, nz, ocean_mode, omega,           &
    196199             omega_sor, outflow_source_plane, passive_scalar,                  &
    197              prandtl_number, psolver, pt_damping_factor,        &
    198              pt_damping_width, pt_reference, pt_surface,                       &
     200             prandtl_number, psolver, pt_damping_factor,  pt_damping_width,    &
     201             pt_reference, pt_surface, pt_surface_heating_rate,                &
    199202             pt_surface_initial_change, pt_vertical_gradient,                  &
    200203             pt_vertical_gradient_level, q_surface, q_surface_initial_change,  &
     
    267270             nsor, nsor_ini, nudging, nx, ny, nz, ocean_mode, omega,           &
    268271             omega_sor, outflow_source_plane, passive_scalar,                  &
    269              prandtl_number, psolver, pt_damping_factor,                       &
    270              pt_damping_width, pt_reference, pt_surface,                       &
     272             prandtl_number, psolver, pt_damping_factor, pt_damping_width,     &
     273             pt_surface_heating_rate, pt_reference, pt_surface,                &
    271274             pt_surface_initial_change, pt_vertical_gradient,                  &
    272275             pt_vertical_gradient_level, q_surface, q_surface_initial_change,  &
  • palm/trunk/SOURCE/time_integration.f90

    r4564 r4565  
    2525! -----------------
    2626! $Id$
     27! added new surface temperature forcing method for bc_pt_b = 'dirichlet':
     28! surface temperature pt(0) can be linearly increased by pt_surface_heating_rate (in K/h)
     29!
     30! 4564 2020-06-12 14:03:36Z raasch
    2731! Vertical nesting method of Huq et al. (2019) removed
    2832!
     
    269273               multi_agent_system_end, multi_agent_system_start, nesting_offline, neutral,         &
    270274               nr_timesteps_this_run, nudging, ocean_mode, passive_scalar, pt_reference,           &
    271                pt_slope_offset, random_heatflux, rans_tke_e, run_coupled, salsa,                   &
     275               pt_slope_offset, pt_surface, pt_surface_heating_rate, pt_surface_initial_change,    &
     276               random_heatflux, rans_tke_e, run_coupled, salsa,                                    &
    272277               simulated_time, simulated_time_chr, skip_time_do2d_xy, skip_time_do2d_xz,           &
    273278               skip_time_do2d_yz, skip_time_do3d, skip_time_domask, skip_time_dopr,                &
     
    944949
    945950!
     951!--       Increase temperature pt(0) according to pt_surface_heating_rate (convert from K/h to K/s)
     952          IF ( pt_surface_heating_rate /= 0.0_wp ) THEN
     953             pt(0,:,:) = pt_surface + pt_surface_initial_change +                                  &
     954                         time_since_reference_point * pt_surface_heating_rate / 3600.0_wp
     955          ENDIF
     956
     957!
    946958!--       Impose a turbulent inflow using the recycling method
    947959          IF ( turbulent_inflow )  CALL inflow_turbulence
Note: See TracChangeset for help on using the changeset viewer.