Changeset 3014 for palm/trunk/SOURCE/time_integration.f90

Timestamp:

May 9, 2018 8:42:38 AM (6 years ago)

Author:

maronga

Message:

series of bugfixes

File:

• palm/trunk/SOURCE/time_integration.f90 (modified) (4 diffs)

palm/trunk/SOURCE/time_integration.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Fixed bug in IF statement
+! Ensure that the time when calling the radiation to be the time step of the
+! pre-calculated time when first calculate the positions of the sun
+! 
+! 3004 2018-04-27 12:33:25Z Giersch
 ! First call of flow_statistics has been removed. It is already called in 
 ! run_control itself
@@ -491 +496 @@
 
     CHARACTER (LEN=9) ::  time_to_string          !< 
-    INTEGER           ::  n
-    INTEGER           ::  lsp
+    INTEGER(iwp)      ::  it
+    INTEGER(iwp)      ::  lsp
+    INTEGER(iwp)      ::  n
+
 
     REAL(wp) ::  dt_3d_old  !< temporary storage of timestep to be used for
                             !< steering of run control output interval
+    REAL(wp) ::  tsrp_org   !< original value of time_since_reference_point
 !
 !-- At beginning determine the first time step
@@ -1040 +1048 @@
                 ENDIF
 
+!
+!--             Adjust the current_ time to the time step of the radiation model.
+!--             Needed since radiation is pre-calculated and stored only on apparent 
+!--             solar positions
+                it = FLOOR(time_since_reference_point/dt_radiation)
+                tsrp_org = time_since_reference_point
+                time_since_reference_point = REAL(it,wp) * dt_radiation
+
                 CALL radiation_control
 
@@ -1050 +1066 @@
                    CALL cpu_log( log_point(75), 'radiation_interaction', 'stop' )
                 ENDIF
+    
+!
+!--             Return the current time to its original value
+                time_since_reference_point = tsrp_org
 
              ENDIF