Changeset 3857


Ignore:
Timestamp:
Apr 3, 2019 1:00:16 PM (2 years ago)
Author:
knoop
Message:

Introduced 1D arrays for horizontal grid dimensions

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/init_grid.f90

    r3855 r3857  
    369369
    370370    USE arrays_3d,                                                             &
    371         ONLY:  dd2zu, ddzu, ddzu_pres, ddzw, dzu, dzw, zu, zw
     371        ONLY:  dd2zu, ddzu, ddzu_pres, ddzw, dzu, dzw, x, xu, y, yv, zu, zw
    372372       
    373373    USE control_parameters,                                                    &
     
    431431    nysg = nys - nbgp
    432432    nyng = nyn + nbgp
     433
     434!
     435!-- Allocate grid arrays
     436    ALLOCATE( x(0:nx), xu(0:nx) )
     437    DO i = 0, nx
     438       xu(i) = i * dx
     439       x(i)  = i * dx + 0.5_wp * dx
     440    ENDDO
     441
     442    ALLOCATE( y(0:ny), yv(0:ny) )
     443    DO j = 0, ny
     444       yv(j) = j * dy
     445       y(j)  = j * dy + 0.5_wp * dy
     446    ENDDO
    433447
    434448!
  • palm/trunk/SOURCE/modules.f90

    r3747 r3857  
    801801    REAL(wp), DIMENSION(:), ALLOCATABLE ::  v_stokes_zw            !< v-component of Stokes drift velocity at zw levels
    802802    REAL(wp), DIMENSION(:), ALLOCATABLE ::  w_subs                 !< subsidence/ascent velocity
     803    REAL(wp), DIMENSION(:), ALLOCATABLE ::  x                      !< horizontal grid coordinate of v-grid (in m)
     804    REAL(wp), DIMENSION(:), ALLOCATABLE ::  xu                     !< horizontal grid coordinate of u-grid (in m)
     805    REAL(wp), DIMENSION(:), ALLOCATABLE ::  y                      !< horizontal grid coordinate of u-grid (in m)
     806    REAL(wp), DIMENSION(:), ALLOCATABLE ::  yv                     !< horizontal grid coordinate of v-grid (in m)
    803807    REAL(wp), DIMENSION(:), ALLOCATABLE ::  zu                     !< vertical grid coordinate of u-grid (in m)
    804808    REAL(wp), DIMENSION(:), ALLOCATABLE ::  zw                     !< vertical grid coordinate of w-grid (in m)
Note: See TracChangeset for help on using the changeset viewer.