Changeset 4654 for palm/trunk/SOURCE/time_integration_spinup.f90

Timestamp:

Aug 28, 2020 1:47:23 PM (4 years ago)

Author:

pavelkrc

Message:

Add possibility to output surface data during spinup

File:

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

palm/trunk/SOURCE/time_integration_spinup.f90

@@ -116 +116 @@
                spinup_pt_mean,                                                                     &
                spinup_time,                                                                        &
+               surface_output,                                                                     &
                timestep_count,                                                                     &
                time_dopr,                                                                          &
@@ -174 +175 @@
     USE statistics,                                                                                &
         ONLY:  flow_statistics_called
+
+    USE surface_data_output_mod,                                                                   &
+        ONLY:  dt_dosurf,                                                                          &
+               skip_time_dosurf,                                                                   &
+               surface_data_output,                                                                &
+               time_dosurf
 
     USE surface_layer_fluxes_mod,                                                                  &
@@ -490 +497 @@
              time_dopr         = time_dopr        + dt_3d
           ENDIF
+          IF ( surface_output )  THEN
+             IF ( simulated_time >= skip_time_dosurf )  THEN
+                time_dosurf    = time_dosurf + dt_3d
+             ENDIF
+          ENDIF
           time_run_control     = time_run_control + dt_3d
 
@@ -545 +557 @@
              time_do3d = MOD( time_do3d, MAX( dt_do3d, dt_3d ) )
           ENDIF
-
+!
+!--       Output of surface data
+          IF ( surface_output )  THEN
+             IF ( time_dosurf >= dt_dosurf )  THEN
+                CALL surface_data_output( 0 )
+                time_dosurf = MOD( time_dosurf, MAX( dt_dosurf, dt_3d ) )
+             ENDIF
+          ENDIF
 
        ENDIF