Changeset 3857

Timestamp:

Apr 3, 2019 1:00:16 PM (5 years ago)

Author:

knoop

Message:

Introduced 1D arrays for horizontal grid dimensions

Location:

palm/trunk/SOURCE

Files:

• init_grid.f90 (modified) (2 diffs)
• modules.f90 (modified) (1 diff)

palm/trunk/SOURCE/init_grid.f90

@@ -369 +369 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  dd2zu, ddzu, ddzu_pres, ddzw, dzu, dzw, zu, zw
+        ONLY:  dd2zu, ddzu, ddzu_pres, ddzw, dzu, dzw, x, xu, y, yv, zu, zw
         
     USE control_parameters,                                                    &
@@ -431 +431 @@
     nysg = nys - nbgp
     nyng = nyn + nbgp
+
+!
+!-- Allocate grid arrays
+    ALLOCATE( x(0:nx), xu(0:nx) )
+    DO i = 0, nx
+       xu(i) = i * dx
+       x(i)  = i * dx + 0.5_wp * dx
+    ENDDO
+
+    ALLOCATE( y(0:ny), yv(0:ny) )
+    DO j = 0, ny
+       yv(j) = j * dy
+       y(j)  = j * dy + 0.5_wp * dy
+    ENDDO
 
 !

palm/trunk/SOURCE/modules.f90

@@ -801 +801 @@
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  v_stokes_zw            !< v-component of Stokes drift velocity at zw levels
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  w_subs                 !< subsidence/ascent velocity
+    REAL(wp), DIMENSION(:), ALLOCATABLE ::  x                      !< horizontal grid coordinate of v-grid (in m)
+    REAL(wp), DIMENSION(:), ALLOCATABLE ::  xu                     !< horizontal grid coordinate of u-grid (in m)
+    REAL(wp), DIMENSION(:), ALLOCATABLE ::  y                      !< horizontal grid coordinate of u-grid (in m)
+    REAL(wp), DIMENSION(:), ALLOCATABLE ::  yv                     !< horizontal grid coordinate of v-grid (in m)
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  zu                     !< vertical grid coordinate of u-grid (in m)
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  zw                     !< vertical grid coordinate of w-grid (in m)