Changeset 3469 for palm/trunk/SOURCE

Timestamp:

Oct 30, 2018 8:05:07 PM (6 years ago)

Author:

kanani

Message:

Implement indoor climate and energy demand model

Location:

palm/trunk/SOURCE

Files:

• Makefile (modified) (7 diffs)
• indoor_model_mod.f90 (added)
• init_3d_model.f90 (modified) (3 diffs)
• modules.f90 (modified) (2 diffs)
• parin.f90 (modified) (3 diffs)
• time_integration.f90 (modified) (4 diffs)
• urban_surface_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/Makefile

@@ -25 +25 @@
 # -----------------
 # $Id$
+# Add indoor model (kanani, srissman, tlang),
+# minor formatting
+# 
+# 3467 2018-10-30 19:05:21Z suehring
 # Implementation of a new aerosol module salsa.
 # 
@@ -514 +518 @@
         average_3d_data.f90 \
         basic_constants_and_equations_mod.f90 \
+        biometeorology_mod.f90 \
+        biometeorology_ipt_mod.f90 \
+        biometeorology_pt_mod.f90 \
+        biometeorology_pet_mod.f90 \
+        biometeorology_utci_mod.f90 \
         boundary_conds.f90 \
         buoyancy.f90 \
@@ -553 +562 @@
         gust_mod.f90 \
         header.f90 \
-        biometeorology_mod.f90 \
-        biometeorology_ipt_mod.f90 \
-        biometeorology_pt_mod.f90 \
-        biometeorology_pet_mod.f90 \
-        biometeorology_utci_mod.f90 \
+        indoor_model_mod.f90 \
         inflow_turbulence.f90 \
         init_3d_model.f90 \
@@ -1080 +1085 @@
         surface_mod.o \
         virtual_flight_mod.o
+indoor_model_mod.o: \
+        cpulog_mod.o \
+        date_and_time_mod.o \
+        mod_kinds.o \
+        modules.o \
+        netcdf_data_input_mod.o \
+        surface_mod.o \
+        urban_surface_mod.o
 inflow_turbulence.o: \
         cpulog_mod.o \
@@ -1093 +1106 @@
         disturb_heatflux.o \
         gust_mod.o \
+        indoor_model_mod.o \
         land_surface_model_mod.o \
         large_scale_forcing_nudging_mod.o \
@@ -1415 +1429 @@
         date_and_time_mod.o \
         gust_mod.o \
+        indoor_model_mod.o \
         land_surface_model_mod.o \
         model_1d_mod.o \
@@ -1708 +1723 @@
         disturb_heatflux.o \
         gust_mod.o \
+        indoor_model_mod.o \
         interaction_droplets_ptq.o \
         land_surface_model_mod.o \

palm/trunk/SOURCE/init_3d_model.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Add indoor model (kanani, srissman, tlang)
+! 
+! 3467 2018-10-30 19:05:21Z suehring
 ! Implementation of a new aerosol module salsa.
 ! 
@@ -567 +570 @@
 
     USE indices
-    
+
+    USE indoor_model_mod,                                                      &
+        ONLY:  im_init
+
     USE kinds
 
@@ -2447 +2453 @@
 
 !
+!-- If required, initialize indoor model
+    IF ( indoor_model )  THEN
+       CALL location_message( 'initializing indoor model', .FALSE. )
+       CALL im_init
+       CALL location_message( 'finished', .TRUE. )
+    ENDIF
+
+!
 !-- Initialize the ws-scheme.    
     IF ( ws_scheme_sca .OR. ws_scheme_mom )  CALL ws_init

palm/trunk/SOURCE/modules.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Add indoor model (kanani, srissman, tlang)
+! 
+! 3467 2018-10-30 19:05:21Z suehring
 ! Add biometeorology
 ! 
@@ -1345 +1348 @@
     LOGICAL ::  humidity = .FALSE.                               !< namelist parameter
     LOGICAL ::  humidity_remote = .FALSE.                        !< switch for receiving near-surface humidity flux (atmosphere-ocean coupling)
+    LOGICAL ::  indoor_model = .FALSE.                           !< switch for indoor-climate and energy-demand model
     LOGICAL ::  large_scale_forcing = .FALSE.                    !< namelist parameter
     LOGICAL ::  large_scale_subsidence = .FALSE.                 !< namelist parameter

palm/trunk/SOURCE/parin.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Add indoor model (kanani, srissman, tlang),
+! minor formatting
+! 
+! 3467 2018-10-30 19:05:21Z suehring
 ! Implementation of a new aerosol module salsa.
 ! 
@@ -484 +488 @@
     USE indices,                                                               &
         ONLY:  nx, ny, nz
+
+    USE indoor_model_mod,                                                      &
+        ONLY:  im_parin
 
     USE kinds
@@ -917 +924 @@
           CALL chem_parin
           CALL uvem_parin
-!
-!--       Check if SALSA processes should be carried out and read &salsa_par
-!--       if required
-          CALL salsa_parin          
+          CALL im_parin
+          CALL salsa_parin
 !
 !--       Read user-defined variables

palm/trunk/SOURCE/time_integration.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Add indoor model (kanani, srissman, tlang)
+! 
+! 3467 2018-10-30 19:05:21Z suehring
 ! Implementation of a new aerosol module salsa.
 ! 
@@ -433 +436 @@
                dt_dopr_listing, dt_dots, dt_dvrp, dt_run_control, end_time,    &
                first_call_lpm, first_call_mas, galilei_transformation,         &
-               humidity, intermediate_timestep_count,                          &
+               humidity, indoor_model, intermediate_timestep_count,            &
                intermediate_timestep_count_max,                                &
                land_surface, large_scale_forcing,                              &
@@ -471 +474 @@
     USE indices,                                                               &
         ONLY:  nbgp, nx, nxl, nxlg, nxr, nxrg, nzb, nzt
+
+    USE indoor_model_mod,                                                      &
+        ONLY:  dt_indoor, im_main_heatcool, skip_time_do_indoor, time_indoor
 
     USE interaction_droplets_ptq_mod,                                          &
@@ -1248 +1254 @@
        ENDIF
 !
+!--    If required, calculate indoor temperature, waste heat, heat flux
+!--    through wall, etc.
+!--    dt_indoor steers the frequency of the indoor model calculations
+       IF (       indoor_model                                                    &
+            .AND. intermediate_timestep_count == intermediate_timestep_count_max  &
+            .AND. time_since_reference_point > skip_time_do_indoor )  THEN
+
+          time_indoor = time_indoor + dt_3d
+
+          IF ( time_indoor >= dt_indoor )  THEN
+
+             time_indoor = time_indoor - dt_indoor
+
+             CALL cpu_log( log_point(76), 'indoor_model', 'start' )
+             CALL im_main_heatcool
+             CALL cpu_log( log_point(76), 'indoor_model', 'stop' )
+
+          ENDIF
+       ENDIF
+!
 !--    Increase simulation time and output times
        nr_timesteps_this_run      = nr_timesteps_this_run + 1

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Add missing PUBLIC variables for new indoor model
+! 
+! 3449 2018-10-29 19:36:56Z suehring
 ! Bugfix: Fix average arrays allocations in usm_average_3d_data (J.Resler)
 ! Bugfix: Fix reading wall temperatures (J.Resler)
@@ -1251 +1254 @@
 !-- Public parameters, constants and initial values
     PUBLIC usm_anthropogenic_heat, usm_material_model, usm_wall_mod,           &
-           usm_green_heat_model, usm_temperature_near_surface, building_pars
+           usm_green_heat_model, usm_temperature_near_surface, building_pars,  &
+           nzt_wall, t_surf_10cm_h, t_surf_10cm_v, t_wall_h, t_wall_v,         &
+           t_window_h, t_window_v, building_type