Changeset 2354 for palm

Timestamp:

Aug 17, 2017 10:49:36 AM (7 years ago)

Author:

schwenkel

Message:

bugfix for lsm data output

Location:

palm/trunk/SOURCE

Files:

• check_parameters.f90 (modified) (2 diffs)
• land_surface_model_mod.f90 (modified) (4 diffs)

palm/trunk/SOURCE/check_parameters.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix correlated to lsm_check_data_output_pr.
+! If-statement for following checks is essential, otherwise units for lsm output 
+! are set to 'illegal' and palm will be aborted.
+! 
+! 2348 2017-08-10 10:40:10Z kanani
 ! New: Check for simultaneous use of geostrophic wind and u_profile/v_profile
 ! 
@@ -2903 +2908 @@
                                             dopr_unit(i) )
 
-             CALL lsf_nudging_check_data_output_pr( data_output_pr(i), i,      &
+             IF ( unit == 'illegal' )  THEN
+                CALL lsf_nudging_check_data_output_pr( data_output_pr(i), i,   &
                                                     unit, dopr_unit(i) )
+             ENDIF
 
              IF ( unit == 'illegal' )  THEN

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! minor bugfixes
+! 
+! 2340 2017-08-07 17:11:13Z maronga
 ! Revised root_distribution tabel and implemented a pseudo-generic root fraction
 ! calculation
@@ -1234 +1237 @@
        ENDDO    
     ENDIF
+    nzs = nzt_soil + 1
+
 
     IF ( vegetation_type == 0 )  THEN
@@ -1290 +1295 @@
        
     ddz_soil_center = 1.0_wp / dz_soil_center
-    ddz_soil = 1.0_wp / dz_soil
+    ddz_soil(nzb_soil:nzt_soil) = 1.0_wp / dz_soil(nzb_soil:nzt_soil)
+
 
 
@@ -1431 +1437 @@
 !--    for vegetation type surfaces, the thermal conductivity of the soil is
 !--    needed
+
        IF ( surf%vegetation_surface(m) )  THEN