Changeset 2550 for palm


Ignore:
Timestamp:
Oct 16, 2017 5:12:01 PM (7 years ago)
Author:
boeske
Message:

enable simulations with complex terrain

Location:
palm/trunk/SOURCE
Files:
6 edited

Legend:

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

    r2513 r2550  
    2525! -----------------
    2626! $Id$
     27! Added checks for complex terrain simulations
     28!
     29! 2513 2017-10-04 09:24:39Z kanani
    2730! Bugfix for some dopr(_initial)_index values and units connected to
    2831! passive-scalar output
     
    39013904    ENDIF
    39023905
     3906!
    39033907!-- Vertical nesting: check fine and coarse grid compatibility for data exchange
    39043908    IF ( vnested )  CALL vnest_check_parameters
     3909
     3910!
     3911!-- Check if topography is read from file in case of complex terrain simulations
     3912    IF ( complex_terrain  .AND.  TRIM( topography ) /= 'read_from_file' )  THEN
     3913       message_string = 'complex_terrain requires topography' //               &
     3914                        ' = ''read_from_file'''
     3915       CALL message( 'check_parameters', 'PA0472', 1, 2, 0, 6, 0 )
     3916    ENDIF
     3917
     3918!
     3919!-- Check if vertical grid stretching is switched off in case of complex
     3920!-- terrain simulations
     3921    IF ( complex_terrain  .AND.  dz_stretch_level < 100000.0_wp )  THEN
     3922       message_string = 'Vertical grid stretching is not allowed for ' //      &
     3923                        'complex_terrain = .T.'
     3924       CALL message( 'check_parameters', 'PA0473', 1, 2, 0, 6, 0 )
     3925    ENDIF
    39053926
    39063927    CALL location_message( 'finished', .TRUE. )
  • palm/trunk/SOURCE/header.f90

    r2544 r2550  
    2525! -----------------
    2626! $Id$
     27! Added output for complex terrain simulations
     28!
     29! 2544 2017-10-13 18:09:32Z maronga
    2730! Moved initial day of year and time to inipar.
    2831!
     
    403406
    404407    USE surface_mod,                                                           &
    405         ONLY:  surf_def_h
     408        ONLY:  surf_def_h, get_topography_top_index
    406409
    407410    USE synthetic_turbulence_generator_mod,                                    &
     
    972975       ELSEIF ( TRIM( topography_grid_convention ) == 'cell_center' )  THEN
    973976          WRITE ( io, 279 )
     977       ENDIF
     978    ENDIF
     979
     980!-- Complex terrain
     981    IF ( complex_terrain )  THEN
     982       WRITE( io, 280 )
     983       IF ( turbulent_inflow )  THEN
     984          WRITE( io, 281 )  zu( get_topography_top_index( 0, 0, 's' ) )
     985       ENDIF
     986       IF ( TRIM( initializing_actions ) == 'cyclic_fill' )  THEN
     987          WRITE( io, 282 )
    974988       ENDIF
    975989    ENDIF
     
    20512065279 FORMAT (' Topography grid definition convention:'/ &
    20522066            ' cell center (scalar grid points)' /)
     2067280 FORMAT (' Complex terrain simulation is activated.')
     2068281 FORMAT ('    --> Mean inflow profiles are adjusted.' / &
     2069            '    --> Elevation of inflow boundary: ', F7.1, ' m' )
     2070282 FORMAT ('    --> Initial data from 3D-precursor run is shifted' / &
     2071            '        vertically depending on local surface height.')
    20532072300 FORMAT (//' Boundary conditions:'/ &
    20542073             ' -------------------'// &
  • palm/trunk/SOURCE/init_3d_model.f90

    r2513 r2550  
    2525! -----------------
    2626! $Id$
     27! Modifications to cyclic fill method and turbulence recycling method in case of
     28! complex terrain simulations
     29!
     30! 2513 2017-10-04 09:24:39Z kanani
    2731! Bugfix in storing initial scalar profile (wrong index)
    2832!
     
    436440    USE surface_mod,                                                           &
    437441        ONLY :  init_surface_arrays, init_surfaces, surf_def_h, surf_lsm_h,    &
    438                 surf_usm_h
     442                surf_usm_h, get_topography_top_index
    439443   
    440444    USE transpose_indices
     
    477481    REAL(wp), DIMENSION(:), ALLOCATABLE ::  ngp_3d_inner_tmp  !<
    478482
     483    INTEGER(iwp) ::  nz_u_shift   !<
     484    INTEGER(iwp) ::  nz_v_shift   !<
     485    INTEGER(iwp) ::  nz_w_shift   !<
     486    INTEGER(iwp) ::  nz_s_shift   !<
     487    INTEGER(iwp) ::  nz_u_shift_l !<
     488    INTEGER(iwp) ::  nz_v_shift_l !<
     489    INTEGER(iwp) ::  nz_w_shift_l !<
     490    INTEGER(iwp) ::  nz_s_shift_l !<
    479491
    480492    CALL location_message( 'allocating arrays', .FALSE. )
     
    14571469
    14581470!
     1471!--    In case of complex terrain and cyclic fill method as initialization,
     1472!--    shift initial data in the vertical direction for each point in the
     1473!--    x-y-plane depending on local surface height
     1474       IF ( complex_terrain  .AND.                                             &
     1475            TRIM( initializing_actions ) == 'cyclic_fill' )  THEN
     1476          DO  i = nxlg, nxrg
     1477             DO  j = nysg, nyng
     1478                nz_u_shift = get_topography_top_index( j, i, 'u' )
     1479                nz_v_shift = get_topography_top_index( j, i, 'v' )
     1480                nz_w_shift = get_topography_top_index( j, i, 'w' )
     1481                nz_s_shift = get_topography_top_index( j, i, 's' )
     1482
     1483                u(nz_u_shift:nzt+1,j,i)  = u(0:nzt+1-nz_u_shift,j,i)               
     1484
     1485                v(nz_v_shift:nzt+1,j,i)  = v(0:nzt+1-nz_v_shift,j,i)
     1486
     1487                w(nz_w_shift:nzt+1,j,i)  = w(0:nzt+1-nz_w_shift,j,i)
     1488
     1489                e(nz_s_shift:nzt+1,j,i)  =  e(0:nzt+1-nz_s_shift,j,i)
     1490                p(nz_s_shift:nzt+1,j,i)  =  p(0:nzt+1-nz_s_shift,j,i)
     1491                pt(nz_s_shift:nzt+1,j,i) = pt(0:nzt+1-nz_s_shift,j,i)
     1492                km(nz_s_shift:nzt+1,j,i) = km(0:nzt+1-nz_s_shift,j,i)
     1493                kh(nz_s_shift:nzt+1,j,i) = kh(0:nzt+1-nz_s_shift,j,i)
     1494             ENDDO
     1495          ENDDO
     1496       ENDIF
     1497
     1498!
    14591499!--    Initialization of the turbulence recycling method
    14601500       IF ( TRIM( initializing_actions ) == 'cyclic_fill'  .AND.               &
     
    14801520          IF ( passive_scalar )                                                &
    14811521             mean_inflow_profiles(:,7) = hom_sum(:,115,0)   ! s
     1522!
     1523!--       In case of complex terrain, determine vertical displacement at inflow
     1524!--       boundary and adjust mean inflow profiles
     1525          IF ( complex_terrain )  THEN
     1526             IF ( nxlg <= 0 .AND. nxrg >= 0 .AND. nysg <= 0 .AND. nyng >= 0 )  THEN
     1527                nz_u_shift_l = get_topography_top_index( 0, 0, 'u' )
     1528                nz_v_shift_l = get_topography_top_index( 0, 0, 'v' )
     1529                nz_w_shift_l = get_topography_top_index( 0, 0, 'w' )
     1530                nz_s_shift_l = get_topography_top_index( 0, 0, 's' )
     1531             ELSE
     1532                nz_u_shift_l = 0
     1533                nz_v_shift_l = 0
     1534                nz_w_shift_l = 0
     1535                nz_s_shift_l = 0
     1536             ENDIF
     1537
     1538#if defined( __parallel )
     1539             CALL MPI_ALLREDUCE(nz_u_shift_l, nz_u_shift, 1, MPI_INTEGER,      &
     1540                                MPI_MAX, comm2d, ierr)
     1541             CALL MPI_ALLREDUCE(nz_v_shift_l, nz_v_shift, 1, MPI_INTEGER,      &
     1542                                MPI_MAX, comm2d, ierr)
     1543             CALL MPI_ALLREDUCE(nz_w_shift_l, nz_w_shift, 1, MPI_INTEGER,      &
     1544                                MPI_MAX, comm2d, ierr)
     1545             CALL MPI_ALLREDUCE(nz_s_shift_l, nz_s_shift, 1, MPI_INTEGER,      &
     1546                                MPI_MAX, comm2d, ierr)
     1547#else
     1548             nz_u_shift = nz_u_shift_l
     1549             nz_v_shift = nz_v_shift_l
     1550             nz_w_shift = nz_w_shift_l
     1551             nz_s_shift = nz_s_shift_l
     1552#endif
     1553
     1554             mean_inflow_profiles(:,1) = 0.0_wp
     1555             mean_inflow_profiles(nz_u_shift:nzt+1,1) = hom_sum(0:nzt+1-nz_u_shift,1,0)  ! u
     1556
     1557             mean_inflow_profiles(:,2) = 0.0_wp
     1558             mean_inflow_profiles(nz_v_shift:nzt+1,2) = hom_sum(0:nzt+1-nz_v_shift,2,0)  ! v
     1559
     1560             mean_inflow_profiles(nz_s_shift:nzt+1,4) = hom_sum(0:nzt+1-nz_s_shift,4,0)  ! pt
     1561             mean_inflow_profiles(nz_s_shift:nzt+1,5) = hom_sum(0:nzt+1-nz_s_shift,8,0)  ! e
     1562
     1563          ENDIF
    14821564
    14831565!
  • palm/trunk/SOURCE/init_grid.f90

    r2478 r2550  
    2525! -----------------
    2626! $Id$
     27! Set lateral boundary conditions for topography on all three ghost layers
     28!
     29! 2478 2017-09-18 13:37:24Z suehring
    2730! Bugfix, correct flag for use_top
    2831!
     
    10411044         
    10421045          IF ( .NOT. bc_ns_cyc )  THEN
    1043              IF ( nys == 0  )  nzb_local(-1,:)   = nzb_local(0,:)
    1044              IF ( nyn == ny )  nzb_local(ny+1,:) = nzb_local(ny,:)
     1046             IF ( nys == 0  )  THEN
     1047                nzb_local(-1,:)   = nzb_local(0,:)
     1048                nzb_local(-2,:)   = nzb_local(0,:)
     1049                nzb_local(-3,:)   = nzb_local(0,:)
     1050             ELSEIF ( nyn == ny )  THEN
     1051                nzb_local(ny+1,:) = nzb_local(ny,:)
     1052                nzb_local(ny+2,:) = nzb_local(ny,:)
     1053                nzb_local(ny+3,:) = nzb_local(ny,:)
     1054             ENDIF
    10451055          ENDIF
    10461056
    10471057          IF ( .NOT. bc_lr_cyc )  THEN
    1048              IF ( nxl == 0  )  nzb_local(:,-1)   = nzb_local(:,0)
    1049              IF ( nxr == nx )  nzb_local(:,nx+1) = nzb_local(:,nx)         
     1058             IF ( nxl == 0  )  THEN
     1059                nzb_local(:,-1)   = nzb_local(:,0)
     1060                nzb_local(:,-2)   = nzb_local(:,0)
     1061                nzb_local(:,-3)   = nzb_local(:,0)
     1062             ELSEIF ( nxr == nx )  THEN
     1063                nzb_local(:,nx+1) = nzb_local(:,nx)
     1064                nzb_local(:,nx+2) = nzb_local(:,nx)         
     1065                nzb_local(:,nx+3) = nzb_local(:,nx)         
     1066             ENDIF         
    10501067          ENDIF
    10511068
  • palm/trunk/SOURCE/modules.f90

    r2508 r2550  
    2525! -----------------
    2626! $Id$
     27! complex_terrain namelist parameter added
     28!
     29! 2508 2017-10-02 08:57:09Z suehring
    2730! Change default value for pt/q/s/sa_vertical_gradient_level
    2831!
     
    11541157    LOGICAL ::  cloud_physics = .FALSE.                      !< namelist parameter
    11551158    LOGICAL ::  cloud_top_radiation = .FALSE.                !< namelist parameter
     1159    LOGICAL ::  complex_terrain = .FALSE.                    !< namelist parameter
    11561160    LOGICAL ::  conserve_volume_flow = .FALSE.               !< namelist parameter
    11571161    LOGICAL ::  constant_diffusion = .FALSE.                 !< diffusion coefficient constant?
     
    19231927
    19241928    CHARACTER (LEN=40) ::  region(0:9)  !< label for statistic region
    1925     
    1926     INTEGER(iwp) ::  pr_palm = 130          !< maximum number of output profiles
     1929 
     1930    INTEGER(iwp) ::  pr_palm = 200          !< maximum number of output profiles
    19271931    INTEGER(iwp) ::  statistic_regions = 0  !< identifier for statistic regions
    1928    
     1932
    19291933    INTEGER(iwp) ::  u_max_ijk(3) = -1  !< index values (i,j,k) of location where u_max occurs
    19301934    INTEGER(iwp) ::  v_max_ijk(3) = -1  !< index values (i,j,k) of location where v_max occurs
  • palm/trunk/SOURCE/parin.f90

    r2544 r2550  
    2525! -----------------
    2626! $Id$
     27! Added complex_terrain
     28!
     29! 2544 2017-10-13 18:09:32Z maronga
    2730! Moved day_of_year_init and time_utc_init to inipar.
    2831!
     
    394397             cloud_physics, cloud_scheme, cloud_top_radiation,                 &
    395398             cloud_water_sedimentation,                                        &
    396              collective_wait, collision_turbulence, conserve_volume_flow,      &
     399             collective_wait, collision_turbulence, complex_terrain,           &
     400             conserve_volume_flow,                                             &
    397401             conserve_volume_flow_mode, constant_flux_layer,                   &
    398402             coupling_start_time, curvature_solution_effects_bulk,             &
Note: See TracChangeset for help on using the changeset viewer.