Changeset 4242 for palm

Timestamp:

Sep 27, 2019 12:59:10 PM (5 years ago)

Author:

suehring

Message:

Bugfixes in indoor model and in chemical emissions

Location:

palm/trunk/SOURCE

Files:

• chem_emissions_mod.f90 (modified) (5 diffs)
• indoor_model_mod.f90 (modified) (4 diffs)

palm/trunk/SOURCE/chem_emissions_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Adjust index_hh access to new definition accompanied with new 
+! palm_date_time_mod. Note, this is just a preleminary fix. (E Chan)
+! 
+! 4230 2019-09-11 13:58:14Z suehring
 ! Bugfix, consider that time_since_reference_point can be also negative when 
 ! time indices are determined.
@@ -1142 +1146 @@
 !-- Read-in the correspondant time factor
 
-                time_factor(:) = emt_att%hourly_emis_time_factor(:,index_hh)      
+                time_factor(:) = emt_att%hourly_emis_time_factor(:,index_hh+1)      
 
 !
@@ -1191 +1195 @@
                 time_factor(:) = emt_att%mdh_emis_time_factor(:,index_mm) *    &
                                  emt_att%mdh_emis_time_factor(:,index_dd) *    &
-                                 emt_att%mdh_emis_time_factor(:,index_hh)
+                                 emt_att%mdh_emis_time_factor(:,index_hh+1)
 
 !
@@ -1240 +1244 @@
 
           index_hh = hour_of_day
-          time_factor(1) = par_emis_time_factor(index_hh)
+          time_factor(1) = par_emis_time_factor(index_hh+1)
 
        ENDIF  ! emiss_lod
@@ -1470 +1474 @@
              emis_distribution(1,nys:nyn,nxl:nxr,ispec) =                                &
                        emt(match_spec_input(ispec))%                                     &
-                           preproc_emission_data(index_hh,1,nys+1:nyn+1,nxl+1:nxr+1) *   &
+                           preproc_emission_data(index_hh+1,1,nys+1:nyn+1,nxl+1:nxr+1) * &
                        conversion_factor
 

palm/trunk/SOURCE/indoor_model_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Bugfix in array index
+! 
+! 4238 2019-09-25 16:06:01Z suehring
 ! - Bugfix in determination of minimum facade height and in location message
 ! - Bugfix, avoid division by zero
@@ -949 +952 @@
          building_type ==  3  .OR.  building_type == 10  .OR.                  &
          building_type == 11  .OR.  building_type == 12 )  THEN
-       buildings(nb)%ventilation_int_loads = 1
+       buildings(:)%ventilation_int_loads = 1
 !
 !-- Office, building with large windows
@@ -955 +958 @@
              building_type ==  6  .OR.  building_type ==  7  .OR.              &
              building_type ==  8  .OR.  building_type ==  9)  THEN
-       buildings(nb)%ventilation_int_loads = 2
+       buildings(:)%ventilation_int_loads = 2
 !
 !-- Industry, hospitals
@@ -961 +964 @@
              building_type == 15  .OR.  building_type == 16  .OR.              &
              building_type == 17  .OR.  building_type == 18 )  THEN
-       buildings(nb)%ventilation_int_loads = 3
+       buildings(:)%ventilation_int_loads = 3
     ENDIF
 !