Changeset 2943 for palm

Timestamp:

Apr 3, 2018 4:17:10 PM (7 years ago)

Author:

suehring

Message:

Calculate exner function at all height levels in USM In order to avoid segmentation faults in case radiation-interactions are switched-off, calculate exner function at all height levels in USM .

Location:

palm/trunk/SOURCE

Files:

• radiation_model_mod.f90 (modified) (1 diff)
• urban_surface_mod.f90 (modified) (6 diffs)

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -6278 +6278 @@
                horz_entry = -HUGE(1._wp)
             ELSE
+write(9,*) "igitt", ig, ubound( nzterr ), ig > ubound( nzterr )
+flush(9)
                horz_entry = (nzterr(ig) - origin(1)) / lastdist
             ENDIF

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Calculate exner function at all height levels and remove some un-used
+! variables.
+! 
+! 2932 2018-03-26 09:39:22Z maronga
 ! renamed urban_surface_par to urban_surface_parameters
 ! 
@@ -299 +303 @@
     
     USE radiation_model_mod,                                                   &
-        ONLY:  albedo_type, radiation_interaction, calc_zenith, zenith,                    &
+        ONLY:  albedo_type, radiation_interaction, calc_zenith, zenith,        &
                rad_sw_in, rad_lw_in, rad_sw_out, rad_lw_out,                   &
                sigma_sb, solar_constant, sun_direction, sun_dir_lat,           &
@@ -310 +314 @@
                iup_u, inorth_u, isouth_u, ieast_u, iwest_u, iup_l,             &
                inorth_l, isouth_l, ieast_l, iwest_l, id,                       &
-               iz, iy, ix, idir, jdir, kdir,  nsurf_type, nzub, nzut,          &
-               nzu, pch, nsurf, idsvf, ndsvf,                                  &
+               iz, iy, ix, idir, jdir, kdir,  nsurf_type, nsurf, idsvf, ndsvf, &
                iup_a, idown_a, inorth_a, isouth_a, ieast_a, iwest_a,           &
                idcsf, ndcsf, kdcsf, pct,                                       &
@@ -7156 +7159 @@
         REAL(wp)                              :: lambda_surface_window !< current value of lambda_surface (heat conductivity between air and window)
         REAL(wp)                              :: lambda_surface_green  !< current value of lambda_surface (heat conductivity between air and greeb wall)
-        REAL(wp), DIMENSION(nzub:nzut)        :: exn                !< value of the Exner function in layers
+        REAL(wp), DIMENSION(nzb:nzt)          :: exn                !< value of the Exner function in layers
         
         REAL(wp), DIMENSION(0:4)              :: dxdir              !< surface normal direction gridbox length
@@ -7166 +7169 @@
         dxdir = (/dz,dy,dy,dx,dx/)
 #if ! defined( __nopointer )
-        exn(:) = (hyp(nzub:nzut) / 100000.0_wp )**0.286_wp          !< Exner function
+        exn(nzb:nzt) = (hyp(nzb:nzt) / 100000.0_wp )**0.286_wp          !< Exner function
 #endif
 !        
@@ -7321 +7324 @@
               force_radiation_call_l = .TRUE.
            ENDIF
-!           
-!--        for horizontal surfaces is pt(nzb_s_inner(j,i),j,i) = pt_surf.
-!--        there is no equivalent surface gridpoint for vertical surfaces.
-!--        pt(k,j,i) is calculated for all directions in diffusion_s
-!--        using surface and wall heat fluxes
-#if ! defined( __nopointer )
-!            pt(k-1,j,i) = ( surf_usm_h%frac(0,m) * t_surf_h_p(m)       &
-!                          + surf_usm_h%frac(2,m) * t_surf_window_h_p(m)  &
-!                          + surf_usm_h%frac(1,m) * t_surf_green_h_p(m) ) &
-!                          / exn(k)  ! not for vertical surfaces
-#endif
-
+!
 !--        calculate fluxes
 !--        rad_net_l is never used!