Changeset 4573 for palm/trunk/SOURCE


Ignore:
Timestamp:
Jun 24, 2020 1:08:47 PM (4 years ago)
Author:
oliver.maas
Message:

pt_surface_heating_rate: calculate pt(0) incrementally by using dt_3d instead of calculating it absolutely by using time_since_reference_point, because time_since_reference_point is set to zero for initializing_actions = 'cyclic_fill', add statement for pt_surface_heating_rate in header.f90

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

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

    r4536 r4573  
    2525! -----------------
    2626! $Id$
     27! added statement for pt_surface_heating_rate
     28!
     29! 4536 2020-05-17 17:24:13Z raasch
    2730! output of restart data format added
    2831!
     
    15791582       WRITE ( io, 475 )  pt_surface_initial_change
    15801583    ENDIF
     1584    IF ( pt_surface_heating_rate /= 0.0_wp )  THEN
     1585       WRITE ( io, 476 )  pt_surface_heating_rate
     1586    ENDIF
    15811587    IF ( humidity  .AND.  q_surface_initial_change /= 0.0_wp )  THEN
    1582        WRITE ( io, 476 )  q_surface_initial_change       
     1588       WRITE ( io, 477 )  q_surface_initial_change       
    15831589    ENDIF
    15841590    IF ( passive_scalar  .AND.  q_surface_initial_change /= 0.0_wp )  THEN
    1585        WRITE ( io, 477 )  q_surface_initial_change       
     1591       WRITE ( io, 478 )  q_surface_initial_change       
    15861592    ENDIF
    15871593
     
    19311937            '    the value is negative) by ',F5.2,' K at the beginning of the',&
    19321938                 ' 3D-simulation'/)
    1933 476 FORMAT ('    The surface humidity is increased (or decreased, ',&
     1939476 FORMAT ('    The surface temperature increases (or decreases, ', &
     1940                 'respectively, if'/ &
     1941            '    the value is negative) by ',F8.4,' K/h during the', &
     1942                 ' 3D-simulation'/)
     1943477 FORMAT ('    The surface humidity is increased (or decreased, ',&
    19341944                 'respectively, if the'/ &
    19351945            '    value is negative) by ',E8.1,' kg/kg at the beginning of', &
    19361946                 ' the 3D-simulation'/)
    1937 477 FORMAT ('    The scalar value is increased at the surface (or decreased, ',&
     1947478 FORMAT ('    The scalar value is increased at the surface (or decreased, ',&
    19381948                 'respectively, if the'/ &
    19391949            '    value is negative) by ',E8.1,' kg/m**3 at the beginning of', &
  • palm/trunk/SOURCE/time_integration.f90

    r4565 r4573  
    2525! -----------------
    2626! $Id$
     27! calculate pt(0) incrementally by using dt_3d instead of calculating it absolutely
     28! by using time_since_reference_point, because time_since_reference_point is set
     29! to zero for initializing_actions = 'cyclic_fill'
     30!
     31! 4565 2020-06-15 08:30:38Z oliver.maas
    2732! added new surface temperature forcing method for bc_pt_b = 'dirichlet':
    2833! surface temperature pt(0) can be linearly increased by pt_surface_heating_rate (in K/h)
     
    950955!
    951956!--       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
     957          IF ( pt_surface_heating_rate /= 0.0_wp .AND. intermediate_timestep_count == 1 ) THEN
     958             pt(0,:,:) = pt(0,:,:) + dt_3d * pt_surface_heating_rate / 3600.0_wp
    955959          ENDIF
    956960
Note: See TracChangeset for help on using the changeset viewer.