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

Timestamp:

Apr 10, 2019 9:31:22 AM (6 years ago)

Author:

suehring

Message:

Make level 3 initialization of urban-surfaces consistent to input data standard; revise flagging of ground-floor level surfaces at buidlings; bugfixes in level 3 initialization of albedo; bugfix in OpenMP directive; check for zero output timestep in surface output; assure that Obukhov lenght does not become zero

File:

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

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix in level 3 initialization of pavement albedo type and pavement 
+! emissivity
+! 
+! 3868 2019-04-08 11:52:36Z suehring
 ! More strict limitation of roughness length when it is in the order of the 
 ! vertical grid spacing 
@@ -4186 +4190 @@
                      pavement_pars_f%fill )                                    &
                    surf_lsm_h%albedo_type(ind_pav_green,m) =                   &
-                                              INT( pavement_pars(ind_p_at,st) )
+                                   INT( pavement_pars_f%pars_xy(ind_p_at,j,i) )
                 IF ( pavement_pars_f%pars_xy(ind_p_emis,j,i) /=                &
                      pavement_pars_f%fill )                                    &
                    surf_lsm_h%emissivity(ind_pav_green,m)  =                   &
-                                              pavement_pars(ind_p_emis,st) 
+                                   pavement_pars_f%pars_xy(ind_p_emis,j,i)
              ENDIF 
 
@@ -4246 +4250 @@
                         pavement_pars_f%fill )                                 &
                       surf_lsm_v(l)%albedo_type(ind_pav_green,m) =             &
-                                            INT( pavement_pars(ind_p_at,st) )
+                                   INT( pavement_pars_f%pars_xy(ind_p_at,j,i) )
 
                    IF ( pavement_pars_f%pars_xy(ind_p_emis,j,i) /=             &
                         pavement_pars_f%fill )                                 &
                       surf_lsm_v(l)%emissivity(ind_pav_green,m)  =             &
-                                            pavement_pars(ind_p_emis,st)  
+                                   pavement_pars_f%pars_xy(ind_p_emis,j,i)  
                 ENDIF 
              ENDDO