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

Timestamp:

Jun 12, 2019 4:50:15 PM (4 years ago)

Author:

suehring

Message:

Improve albedo initialization in land-surface model by setting always some default albedo types, in order to be independent on any user settings; Mask topography at boundary grid cells in nesting mass conservation

File:

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

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Further revision of last commit in order to avoid any side effects when 
+! albedo type is not set in namelist and default albedo type changes. 
+! 
+! 4024 2019-06-12 14:06:46Z suehring
 ! Bugfix in albedo initialization, caused crashes in rrtmg calls
 ! 
@@ -2792 +2796 @@
        ALLOCATE ( surf_lsm_h%albedo_type(0:2,1:surf_lsm_h%ns) )
        ALLOCATE ( surf_lsm_h%emissivity(0:2,1:surf_lsm_h%ns) )
-       surf_lsm_h%albedo_type = albedo_type     
+!
+!--    Initialize albedo type according to its default type, in order to set values
+!--    independent on default albedo_type in radiation model. 
+       surf_lsm_h%albedo_type(ind_veg_wall,:)  =                               &
+                             INT( vegetation_pars(ind_v_at,vegetation_type) )
+       surf_lsm_h%albedo_type(ind_wat_win,:)   =                               &
+                             INT( water_pars(ind_w_at,water_type)           )
+       surf_lsm_h%albedo_type(ind_pav_green,:) =                               &
+                             INT( pavement_pars(ind_p_at,pavement_type)     )
        surf_lsm_h%emissivity  = emissivity
        DO  l = 0, 3
           ALLOCATE ( surf_lsm_v(l)%albedo_type(0:2,1:surf_lsm_v(l)%ns) )
           ALLOCATE ( surf_lsm_v(l)%emissivity(0:2,1:surf_lsm_v(l)%ns)  )
-          surf_lsm_v(l)%albedo_type = albedo_type
+!
+!--       Initialize albedo type according to its default type, in order to 
+!--       set values independent on default albedo_type in radiation model. 
+          surf_lsm_v(l)%albedo_type(ind_veg_wall,:)  =                         &
+                             INT( vegetation_pars(ind_v_at,vegetation_type) )
+          surf_lsm_v(l)%albedo_type(ind_wat_win,:)   =                         &
+                             INT( water_pars(ind_w_at,water_type)           )
+          surf_lsm_v(l)%albedo_type(ind_pav_green,:) =                         &
+                             INT( pavement_pars(ind_p_at,pavement_type)     )
           surf_lsm_v(l)%emissivity  = emissivity
        ENDDO