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

Timestamp:

Dec 2, 2014 5:25:50 PM (9 years ago)

Author:

maronga

Message:

added beta version of a land surface model and a simple radiation model for clear sky conditions

File:

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

palm/trunk/SOURCE/time_integration.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Added calls for the land surface model and radiation scheme
 ! 
 ! Former revisions:
@@ -177 +177 @@
                time_dopr_listing, time_dopts, time_dosp, time_dosp_av,         &
                time_dots, time_do_av, time_do_sla, time_disturb, time_dvrp,    &
-               time_run_control, time_since_reference_point, turbulence,       &
+               time_run_control, time_since_reference_point,                   &
+               turbulence,                                                     &
                turbulent_inflow, use_initial_profile_as_reference,             &
                use_single_reference_value, u_gtrans, v_gtrans, ws_scheme_mom,  &
@@ -194 +195 @@
     USE kinds
 
+    USE land_surface_model_mod,                                                &
+        ONLY:  land_surface, lsm_energy_balance, lsm_soil_model
+
     USE ls_forcing_mod,                                                        &
         ONLY:  ls_forcing_surf, ls_forcing_vert
@@ -214 +218 @@
         ONLY:  prognostic_equations_acc, prognostic_equations_cache,           &
                prognostic_equations_vector
+
+    USE radiation_model_mod,                                                   &
+        ONLY: dt_radiation, lsm_radiation, radiation, time_radiation
 
     USE statistics,                                                            &
@@ -577 +584 @@
 !--       velocities at the outflow in case of a non-cyclic lateral wall)
           CALL boundary_conds
+
+!
+!--       When using the land surface model: 
+!--       1) solve energy balance equation to calculate new skin temperature
+!--       2) run soil model 
+          IF ( land_surface )  THEN
+
+             time_radiation = time_radiation + dt_3d
+
+             IF ( time_radiation >= dt_radiation )  THEN
+
+                CALL cpu_log( log_point(50), 'radiation', 'start' )
+
+                time_radiation = time_radiation - dt_radiation
+                CALL lsm_radiation
+
+                CALL cpu_log( log_point(50), 'radiation', 'stop' )
+             ENDIF
+
+             CALL cpu_log( log_point(54), 'land_surface', 'start' )
+
+             CALL lsm_energy_balance
+             CALL lsm_soil_model
+
+             CALL cpu_log( log_point(54), 'land_surface', 'stop' )
+
+          ENDIF
 
 !