Changeset 2968 for palm/trunk/SOURCE/land_surface_model_mod.f90

Timestamp:

Apr 13, 2018 11:52:24 AM (6 years ago)

Author:

suehring

Message:

Bugfixes in initalization of land-surface model as well as timeseries data output in case of elevated model surfaces.

File:

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

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix in initialization in case of elevated model surface
+! 
+! 2963 2018-04-12 14:47:44Z suehring
 ! - In initialization of surface types, consider the case that surface_fractions
 !   is not given in static input file.
@@ -2374 +2377 @@
 
        USE indices,                                                            &
-           ONLY:  nx, ny
+           ONLY:  nx, ny, topo_min_level
 
        USE pmc_interface,                                                      &
@@ -2408 +2411 @@
        IF (  .NOT.  cloud_physics )  THEN
           CALL calc_mean_profile( pt, 4 )
-          rho_surface = surface_pressure * 100.0_wp / ( r_d * hom(nzb+1,1,4,0) * exn )
+          rho_surface = surface_pressure * 100.0_wp / ( r_d * hom(topo_min_level+1,1,4,0) * exn )
        ENDIF
 
@@ -2424 +2427 @@
        tt_soil_h_m%var_2d    = 0.0_wp
        tm_soil_h_m%var_2d    = 0.0_wp
-       tm_liq_h_m%var_1d  = 0.0_wp
-       surf_lsm_h%c_liq = 0.0_wp
+       tm_liq_h_m%var_1d     = 0.0_wp
+       surf_lsm_h%c_liq      = 0.0_wp
 
        surf_lsm_h%ghf = 0.0_wp
@@ -2443 +2446 @@
           tt_soil_v_m(l)%var_2d    = 0.0_wp
           tm_soil_v_m(l)%var_2d    = 0.0_wp
-          tm_liq_v_m(l)%var_1d  = 0.0_wp
-          surf_lsm_v(l)%c_liq = 0.0_wp
+          tm_liq_v_m(l)%var_1d     = 0.0_wp
+          surf_lsm_v(l)%c_liq      = 0.0_wp
 
           surf_lsm_v(l)%ghf = 0.0_wp