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

Timestamp:

Mar 13, 2018 4:24:40 PM (6 years ago)

Author:

maronga

Message:

bugfix in land surface model and new option in spinup

File:

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

palm/trunk/SOURCE/time_integration_spinup.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Added flag for switching on/off calculation of soil moisture
+! 
+! 2818 2018-02-19 16:42:36Z maronga
 ! Velocity components near walls/ground are now set to the profiles stored in
 ! u_init and v_init. Activated soil moisture calculation during spinup.
@@ -83 +86 @@
 
     USE control_parameters,                                                    &
-        ONLY:  averaging_interval_pr, constant_diffusion, constant_flux_layer, &
+        ONLY:  averaging_interval_pr, calc_soil_moisture_during_spinup,        &
+               constant_diffusion, constant_flux_layer,                        &
                coupling_start_time, current_timestep_number,                   &
                data_output_during_spinup, disturbance_created, dopr_n, do_sum, &
@@ -348 +352 @@
 !--             Call for horizontal upward-facing surfaces
                 CALL lsm_energy_balance( .TRUE., -1 )
-                CALL lsm_soil_model( .TRUE., -1, .TRUE. )
+                CALL lsm_soil_model( .TRUE., -1, calc_soil_moisture_during_spinup )
 !
 !--             Call for northward-facing surfaces
                 CALL lsm_energy_balance( .FALSE., 0 )
-                CALL lsm_soil_model( .FALSE., 0, .TRUE. )
+                CALL lsm_soil_model( .FALSE., 0, calc_soil_moisture_during_spinup )
 !
 !--             Call for southward-facing surfaces
                 CALL lsm_energy_balance( .FALSE., 1 )
-                CALL lsm_soil_model( .FALSE., 1, .TRUE. )
+                CALL lsm_soil_model( .FALSE., 1, calc_soil_moisture_during_spinup )
 !
 !--             Call for eastward-facing surfaces
                 CALL lsm_energy_balance( .FALSE., 2 )
-                CALL lsm_soil_model( .FALSE., 2, .TRUE. )
+                CALL lsm_soil_model( .FALSE., 2, calc_soil_moisture_during_spinup )
 !
 !--             Call for westward-facing surfaces
                 CALL lsm_energy_balance( .FALSE., 3 )
-                CALL lsm_soil_model( .FALSE., 3, .TRUE. )
+                CALL lsm_soil_model( .FALSE., 3, calc_soil_moisture_during_spinup )
 
                 CALL cpu_log( log_point(54), 'land_surface', 'stop' )