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

Timestamp:

Sep 9, 2020 8:27:58 PM (4 years ago)

Author:

pavelkrc

Message:

Radiative transfer model RTM version 4.1

File:

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

palm/trunk/SOURCE/time_integration.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Implementation of downward facing USM and LSM surfaces
+!
+! 4669 2020-09-09 13:43:47Z pavelkrc
 ! - Fix missing call of radiation after spinup
 ! - Fix calculation of force_radiation_call
@@ -339 +342 @@
 
     USE land_surface_model_mod,                                                                    &
-        ONLY:  lsm_boundary_condition, lsm_energy_balance, lsm_soil_model, skip_time_do_lsm
+        ONLY:  lsm_boundary_condition, lsm_energy_balance, skip_time_do_lsm
 
     USE lagrangian_particle_model_mod,                                                             &
@@ -423 +426 @@
 
     USE urban_surface_mod,                                                                         &
-        ONLY:  usm_boundary_condition, usm_material_heat_model, usm_material_model,                &
-               usm_surface_energy_balance, usm_green_heat_model
+        ONLY:  usm_boundary_condition,                                                             &
+               usm_energy_balance
 
     USE virtual_measurement_mod,                                                                   &
@@ -1108 +1111 @@
 
                 CALL cpu_log( log_point(54), 'land_surface', 'start' )
-!
-!--             Call for horizontal upward-facing surfaces
-                CALL lsm_energy_balance( .TRUE., -1 )
-                CALL lsm_soil_model( .TRUE., -1, .TRUE. )
-!
-!--             Call for northward-facing surfaces
-                CALL lsm_energy_balance( .FALSE., 0 )
-                CALL lsm_soil_model( .FALSE., 0, .TRUE. )
-!
-!--             Call for southward-facing surfaces
-                CALL lsm_energy_balance( .FALSE., 1 )
-                CALL lsm_soil_model( .FALSE., 1, .TRUE. )
-!
-!--             Call for eastward-facing surfaces
-                CALL lsm_energy_balance( .FALSE., 2 )
-                CALL lsm_soil_model( .FALSE., 2, .TRUE. )
-!
-!--             Call for westward-facing surfaces
-                CALL lsm_energy_balance( .FALSE., 3 )
-                CALL lsm_soil_model( .FALSE., 3, .TRUE. )
+                CALL lsm_energy_balance( .FALSE. )
 
 !
@@ -1143 +1127 @@
              IF (urban_surface) THEN
                 CALL cpu_log( log_point(74), 'urban_surface', 'start' )
-
-                CALL usm_surface_energy_balance( .FALSE. )
-                IF ( usm_material_model )  THEN
-                   CALL usm_green_heat_model
-                   CALL usm_material_heat_model ( .FALSE. )
-                ENDIF
+                CALL usm_energy_balance( .FALSE. )
 
 !
@@ -1655 +1634 @@
              CALL cpu_log( log_point(23), 'disturb_heatflux', 'stop' )
           ENDIF
-          IF ( surf_lsm_h%ns    >= 1 )  THEN
+          IF ( surf_lsm_h(0)%ns >= 1 )  THEN
              CALL cpu_log( log_point(23), 'disturb_heatflux', 'start' )
-             CALL disturb_heatflux( surf_lsm_h    )
+             CALL disturb_heatflux( surf_lsm_h(0) )
              CALL cpu_log( log_point(23), 'disturb_heatflux', 'stop' )
           ENDIF
-          IF ( surf_usm_h%ns    >= 1 )  THEN
+          IF ( surf_usm_h(0)%ns >= 1 )  THEN
              CALL cpu_log( log_point(23), 'disturb_heatflux', 'start' )
-             CALL disturb_heatflux( surf_usm_h    )
+             CALL disturb_heatflux( surf_usm_h(0) )
              CALL cpu_log( log_point(23), 'disturb_heatflux', 'stop' )
           ENDIF