Ignore:
Timestamp:
Jun 28, 2017 7:53:56 AM (4 years ago)
Author:
maronga
Message:

added new spinup mechanism for surface/radiation models

File:
1 edited

Legend:

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

    r2287 r2296  
    2626! -----------------
    2727! $Id$
     28! Bugfix: adjust output of shf. Added support for spinups
     29!
     30! 2287 2017-06-15 16:46:30Z suehring
    2831! Bugfix in determination topography-top index
    2932!
     
    139142   
    140143    USE control_parameters,                                                    &
    141         ONLY:  dz, topography, dt_3d, intermediate_timestep_count,             &
    142                initializing_actions, intermediate_timestep_count_max,          &
    143                simulated_time, end_time, timestep_scheme, tsc,                 &
    144                coupling_char, io_blocks, io_group, message_string,             &
    145                time_since_reference_point, surface_pressure,                   &
    146                g, pt_surface, large_scale_forcing, lsf_surf,                   &
    147                time_do3d, dt_do3d, average_count_3d, varnamelength,            &
    148                urban_surface
     144        ONLY:  coupling_start_time, dz, topography, dt_3d,                     &
     145               intermediate_timestep_count, initializing_actions,              &
     146               intermediate_timestep_count_max, simulated_time, end_time,      &
     147               timestep_scheme, tsc, coupling_char, io_blocks, io_group,       &
     148               message_string, time_since_reference_point, surface_pressure,   &
     149               g, pt_surface, large_scale_forcing, lsf_surf, spinup,           &
     150               spinup_pt_mean, spinup_time, time_do3d, dt_do3d,                &
     151               average_count_3d, varnamelength, urban_surface
    149152
    150153    USE cpulog,                                                                &
     
    34893492!--    Set flag that indicates that the land surface model is switched on
    34903493       urban_surface = .TRUE.
     3494
     3495!
     3496!--    Activate spinup
     3497       IF ( spinup_time > 0.0_wp )  THEN
     3498          coupling_start_time = spinup_time
     3499          end_time = end_time + spinup_time
     3500          IF ( spinup_pt_mean == 9999999.9_wp )  THEN
     3501             spinup_pt_mean = pt_surface
     3502          ENDIF
     3503          spinup = .TRUE.
     3504       ENDIF
    34913505
    34923506 10    CONTINUE
     
    51615175!--        store kinematic surface heat fluxes for utilization in other processes
    51625176!--        diffusion_s, surface_layer_fluxes,...
    5163            surf_usm_h%shf(m) = surf_usm_h%wshf_eb(m) / rho_cp
     5177           surf_usm_h%shf(m) = surf_usm_h%wshf_eb(m) / cp
    51645178
    51655179       ENDDO
     
    52825296!--           store kinematic surface heat fluxes for utilization in other processes
    52835297!--           diffusion_s, surface_layer_fluxes,...
    5284               surf_usm_v(l)%shf(m) = surf_usm_v(l)%wshf_eb(m) / rho_cp
     5298              surf_usm_v(l)%shf(m) = surf_usm_v(l)%wshf_eb(m) / cp
    52855299
    52865300           ENDDO
     
    53165330!--              Instead of this, we can adjust shf in case AH only at surface
    53175331                 surf_usm_h%shf(m) = surf_usm_h%shf(m) +                       &
    5318                                    aheat(j,i) * acoef * ddx * ddy / rho_cp
     5332                                   aheat(j,i) * acoef * ddx * ddy / cp
    53195333              ENDIF
    53205334           ENDDO
Note: See TracChangeset for help on using the changeset viewer.