Changeset 3065 for palm/trunk/SOURCE/urban_surface_mod.f90

Timestamp:

Jun 12, 2018 7:03:02 AM (6 years ago)

Author:

Giersch

Message:

New vertical stretching procedure has been introduced

File:

• palm/trunk/SOURCE/urban_surface_mod.f90 (modified) (6 diffs)

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Unused array dxdir was removed, dz was replaced by dzu to consider vertical
+! grid stretching
+! 
+! 3049 2018-05-29 13:52:36Z Giersch
 ! Error messages revised
 ! 
@@ -282 +286 @@
 #if ! defined( __nopointer )
     USE arrays_3d,                                                             &
-        ONLY:  hyp, zu, pt, pt_1, pt_2, p, u, v, w, hyp, tend
+        ONLY:  dzu, hyp, zu, pt, pt_1, pt_2, p, u, v, w, hyp, tend
 #endif
 
@@ -292 +296 @@
     
     USE control_parameters,                                                    &
-        ONLY:  coupling_start_time, dz, topography, dt_3d,                     &
+        ONLY:  coupling_start_time, topography, dt_3d,                         &
                intermediate_timestep_count, initializing_actions,              &
                intermediate_timestep_count_max, simulated_time, end_time,      &
@@ -7207 +7211 @@
         REAL(wp), DIMENSION(nzb:nzt)          :: exn                !< value of the Exner function in layers
         
-        REAL(wp), DIMENSION(0:4)              :: dxdir              !< surface normal direction gridbox length
         REAL(wp)                              :: dtime              !< simulated time of day (in UTC)
         INTEGER(iwp)                          :: dhour              !< simulated hour of day (in UTC)
@@ -7213 +7216 @@
 
 
-        dxdir = (/dz,dy,dy,dx,dx/)
 #if ! defined( __nopointer )
         exn(nzb:nzt) = (hyp(nzb:nzt) / 100000.0_wp )**0.286_wp          !< Exner function
@@ -7634 +7636 @@
                          (dtime/3600.0_wp-REAL(dhour,wp))*aheatprof(k,dhour+1)
                  IF ( aheat(k,j,i) > 0.0_wp )  THEN
-                    pt(k,j,i) = pt(k,j,i) + aheat(k,j,i)*acoef*dt_3d/(exn(k)*rho_cp*dz)
+                    pt(k,j,i) = pt(k,j,i) + aheat(k,j,i)*acoef*dt_3d/(exn(k)*rho_cp*dzu(k))
                  ENDIF
               ENDIF