Changeset 4654 for palm

Timestamp:

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

Author:

pavelkrc

Message:

Add possibility to output surface data during spinup

Location:

palm/trunk/SOURCE

Files:

• radiation_model_mod.f90 (modified) (2 diffs)
• time_integration_spinup.f90 (modified) (4 diffs)

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Add possibility to output surface data during spinup (author: J. Resler)
+!
+! 4653 2020-08-27 08:54:43Z pavelkrc
 ! RTM version 4.0:
 ! - Support full-3D geometry:
@@ -35 +38 @@
 ! - Remove separate direction indices for urban and land faces
 ! - Bugfix: skyvft output
+! Authors: P. Krc, J. Resler (Institute of Computer Science, Prague)
 !
 ! 4644 2020-08-22 15:05:46Z pavelkrc

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