Changeset 4230 for palm/trunk/SOURCE/chem_emissions_mod.f90
- Timestamp:
- Sep 11, 2019 1:58:14 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/chem_emissions_mod.f90
r4227 r4230 27 27 ! ----------------- 28 28 ! $Id$ 29 ! Bugfix, consider that time_since_reference_point can be also negative when 30 ! time indices are determined. 31 ! 32 ! 4227 2019-09-10 18:04:34Z gronemeier 29 33 ! implement new palm_date_time_mod 30 34 ! … … 1091 1095 1092 1096 CALL get_date_time( 0.0_wp, second_of_day=time_utc_init ) 1093 CALL get_date_time( time_since_reference_point, hour=hour_of_day ) 1094 1095 days_since_reference_point = INT( FLOOR( ( time_utc_init + time_since_reference_point ) & 1097 CALL get_date_time( MAX( 0.0_wp, time_since_reference_point ), & 1098 hour=hour_of_day ) 1099 1100 days_since_reference_point = INT( FLOOR( ( & 1101 time_utc_init + MAX( 0.0_wp, time_since_reference_point ) ) & 1096 1102 / seconds_per_day ) ) 1097 1103 … … 1124 1130 !-- Update time indices 1125 1131 1126 CALL get_date_time( time_since_reference_point,&1132 CALL get_date_time( MAX( time_since_reference_point, 0.0_wp ), & 1127 1133 day_of_year=day_of_year, hour=hour_of_day ) 1128 1134 index_hh = ( day_of_year - 1_iwp ) * hour_of_day … … 1156 1162 ! 1157 1163 !-- Update time indices 1158 CALL get_date_time( time_since_reference_point, &1159 month=month_of_year, &1160 day=day_of_month, &1161 hour=hour_of_day, &1164 CALL get_date_time( MAX( time_since_reference_point, 0.0_wp ), & 1165 month=month_of_year, & 1166 day=day_of_month, & 1167 hour=hour_of_day, & 1162 1168 day_of_week=day_of_week ) 1163 1169 index_mm = month_of_year … … 1230 1236 ! 1231 1237 !-- Get time-factor for specific hour 1232 CALL get_date_time( time_since_reference_point, hour=hour_of_day ) 1238 CALL get_date_time( MAX( time_since_reference_point, 0.0_wp ), & 1239 hour=hour_of_day ) 1233 1240 1234 1241 index_hh = hour_of_day
Note: See TracChangeset
for help on using the changeset viewer.