Home


About us
PALM
PALM-4U
Installation
Release notes
Documentation
Getting Help
Gallery
Developers area

Changeset 2365 for palm/trunk/SOURCE


Ignore:
Timestamp:
Aug 21, 2017 2:59:59 PM (7 years ago)
Author:
kanani
Message:

Vertical nesting implemented (SadiqHuq?)

Location:
palm/trunk/SOURCE
Files:
1 added
13 edited

Legend:

Unmodified
Added
Removed

  • palm/trunk/SOURCE/Makefile

    r2359 r2365  
    2525# -----------------
    2626# $Id$
     27# Added dependencies for vertical_nesting_mod
     28#
     29# 2339 2017-08-07 13:55:26Z gronemeier
    2730# corrected timestamp in header
    2831#
    … …  
    409412        user_lpm_init.f90 user_lpm_set_attributes.f90 user_module.f90 \
    410413        user_parin.f90 user_read_restart_data.f90 \
    411         user_spectra.f90 user_statistics.f90 virtual_flight_mod.f90 \
     414        user_spectra.f90 user_statistics.f90 vertical_nesting_mod.f90 virtual_flight_mod.f90 \
    412415        wind_turbine_model_mod.f90 write_3d_binary.f90 write_var_list.f90
    413416
    … …  
    425428# The following line is needed for palm_simple_install, don't remove it!
    426429#to_be_replaced_by_include
    427 
    428 #BOUNDS="-Rbc"  # Array bounds checking. Compromises performance seriously.
    429 
    430 F90 = ftn
    431 #COPT = -DMPI_REAL=MPI_DOUBLE_PRECISION -DMPI_2REAL=MPI_2DOUBLE_PRECISION -D__lc -D__parallel -D__fftw \
    432 #             -D__netcdf -D__netcdf4 -D__netcdf4_parallel -D__nopointer
    433 COPT = -DMPI_REAL=MPI_DOUBLE_PRECISION -DMPI_2REAL=MPI_2DOUBLE_PRECISION -D__lc -D__parallel -D__fftw \
    434              -D__netcdf -D__netcdf4 -D__netcdf4_parallel
    435 
    436 #F90FLAGS = -O2 -e Fm -G2 -rm $(BOUNDS) -hnoomp
    437 F90FLAGS = -O2 -e F -G2 -rm $(BOUNDS) -hnoomp
    438 LDFLAGS = $(F90FLAGS) -dynamic
    439 
    440 # End of include
    441430
    442431.SUFFIXES:
    … …  
    469458average_3d_data.o: modules.o cpulog_mod.o mod_kinds.o exchange_horiz_2d.o land_surface_model_mod.o \
    470459        radiation_model_mod.o urban_surface_mod.o
    471 boundary_conds.o: modules.o mod_kinds.o pmc_interface_mod.o surface_mod.o
     460boundary_conds.o: modules.o mod_kinds.o pmc_interface_mod.o surface_mod.o vertical_nesting_mod.o
    472461buoyancy.o: modules.o mod_kinds.o
    473462calc_mean_profile.o: modules.o mod_kinds.o
    … …  
    480469        microphysics_mod.o model_1d_mod.o netcdf_interface_mod.o plant_canopy_model_mod.o pmc_interface_mod.o \
    481470        radiation_model_mod.o spectra_mod.o subsidence_mod.o synthetic_turbulence_generator_mod.o \
    482         wind_turbine_model_mod.o urban_surface_mod.o
     471        wind_turbine_model_mod.o urban_surface_mod.o vertical_nesting_mod.o
    483472close_file.o: modules.o mod_kinds.o netcdf_interface_mod.o
    484473compute_vpt.o: modules.o mod_kinds.o
    … …  
    508497diffusion_v.o: modules.o mod_kinds.o surface_mod.o
    509498diffusion_w.o: modules.o mod_kinds.o surface_mod.o
    510 diffusivities.o: modules.o mod_kinds.o surface_mod.o
     499diffusivities.o: modules.o mod_kinds.o surface_mod.o vertical_nesting_mod.o
    511500disturb_field.o: modules.o cpulog_mod.o mod_kinds.o random_function_mod.o \
    512501                 random_generator_parallel_mod.o
    … …  
    530519init_advec.o: modules.o mod_kinds.o
    531520init_cloud_physics.o: modules.o mod_kinds.o
    532 init_coupling.o: modules.o mod_kinds.o
     521init_coupling.o: modules.o mod_kinds.o vertical_nesting_mod.o
    533522init_dvrp.o: modules.o mod_kinds.o
    534 init_grid.o: modules.o mod_kinds.o advec_ws.o netcdf_interface_mod.o surface_mod.o
     523init_grid.o: modules.o mod_kinds.o advec_ws.o netcdf_interface_mod.o surface_mod.o vertical_nesting_mod.o
    535524init_masks.o: modules.o mod_kinds.o netcdf_interface_mod.o
    536525init_ocean.o: modules.o eqn_state_seawater.o mod_kinds.o
    537 init_pegrid.o: modules.o mod_kinds.o spectra_mod.o synthetic_turbulence_generator_mod.o
     526init_pegrid.o: modules.o mod_kinds.o spectra_mod.o synthetic_turbulence_generator_mod.o \
     527   vertical_nesting_mod.o
    538528init_pt_anomaly.o: modules.o mod_kinds.o
    539529init_rankine.o: modules.o mod_kinds.o
    … …  
    637627surface_layer_fluxes_mod.o: modules.o mod_kinds.o exchange_horiz_2d.o land_surface_model_mod.o \
    638628        urban_surface_mod.o surface_mod.o
    639 surface_mod.o: modules.o mod_kinds.o init_pegrid.o model_1d_mod.o
     629surface_mod.o: modules.o mod_kinds.o model_1d_mod.o
    640630swap_timelevel.o: modules.o cpulog_mod.o mod_kinds.o land_surface_model_mod.o \
    641631   pmc_interface_mod.o urban_surface_mod.o
    … …  
    648638        spectra_mod.o user_actions.o microphysics_mod.o synthetic_turbulence_generator_mod.o \
    649639        surface_layer_fluxes_mod.o surface_mod.o \
    650         urban_surface_mod.o virtual_flight_mod.o wind_turbine_model_mod.o
     640        urban_surface_mod.o vertical_nesting_mod.o virtual_flight_mod.o wind_turbine_model_mod.o
    651641time_integration_spinup.o: cpulog_mod.o modules.o land_surface_model_mod.o mod_kinds.o radiation_model_mod.o surface_layer_fluxes_mod.o urban_surface_mod.o
    652642time_to_string.o: mod_kinds.o
    … …  
    689679user_spectra.o: modules.o mod_kinds.o spectra_mod.o user_module.o
    690680user_statistics.o: modules.o mod_kinds.o netcdf_interface_mod.o user_module.o
     681vertical_nesting_mod.o: modules.o mod_kinds.o surface_mod.o
    691682virtual_flight_mod.o: modules.o cpulog_mod.o mod_kinds.o netcdf_interface_mod.o user_init_flight.o user_flight.o
    692683wind_turbine_model_mod.o: modules.o cpulog_mod.o mod_kinds.o

  • palm/trunk/SOURCE/boundary_conds.f90

    r2320 r2365  
    2525! -----------------
    2626! $Id$
     27! Vertical grid nesting implemented: exclude setting vertical velocity to zero
     28! on fine grid (SadiqHuq)
     29!
     30! 2320 2017-07-21 12:47:43Z suehring
    2731! Remove unused control parameter large_scale_forcing from only-list
    2832!
    … …  
    166170    USE control_parameters,                                                    &
    167171        ONLY:  bc_pt_t_val, bc_q_t_val, bc_s_t_val, constant_diffusion,        &
    168                cloud_physics, dt_3d, humidity,                                 &
     172               cloud_physics, coupling_mode, dt_3d, humidity,                  &
    169173               ibc_pt_b, ibc_pt_t, ibc_q_b, ibc_q_t, ibc_s_b, ibc_s_t,         &
    170174               ibc_sa_t, ibc_uv_b, ibc_uv_t, inflow_l, inflow_n, inflow_r,     &
    … …  
    237241    ENDIF
    238242
    239     IF ( .NOT. nest_domain )  THEN
    240        w_p(nzt:nzt+1,:,:) = 0.0_wp  ! nzt is not a prognostic level (but cf. pres)
     243!
     244!-- Vertical nesting: Vertical velocity not zero at the top of the fine grid
     245    IF (  .NOT.  nest_domain  .AND.                                            &
     246                 TRIM(coupling_mode) /= 'vnested_fine' )  THEN
     247       w_p(nzt:nzt+1,:,:) = 0.0_wp  !< nzt is not a prognostic level (but cf. pres)
    241248    ENDIF
    242249

  • palm/trunk/SOURCE/check_parameters.f90

    r2354 r2365  
    2525! -----------------
    2626! $Id$
     27! Vertical grid nesting implemented: Check coupling mode. Generate file header
     28! (SadiqHuq)
     29!
     30! 2354 2017-08-17 10:49:36Z schwenkel
    2731! Bugfix correlated to lsm_check_data_output_pr.
    2832! If-statement for following checks is essential, otherwise units for lsm output
    … …  
    561565    USE wind_turbine_model_mod,                                                &
    562566        ONLY:  wtm_check_parameters, wind_turbine
     567    USE vertical_nesting_mod,                                                  &
     568        ONLY:  vnested, vnest_check_parameters
    563569
    564570
    … …  
    603609!> @todo Check if any queries for other coupling modes (e.g. precursor_ocean) are missing
    604610    IF ( coupling_mode /= 'uncoupled'            .AND.  &
     611         coupling_mode /= 'vnested_crse'         .AND.  &
     612         coupling_mode /= 'vnested_fine'         .AND.  &
    605613         coupling_mode /= 'atmosphere_to_ocean'  .AND.  &
    606614         coupling_mode /= 'ocean_to_atmosphere' )  THEN
    … …  
    611619!
    612620!-- Check dt_coupling, restart_time, dt_restart, end_time, dx, dy, nx and ny
    613     IF ( coupling_mode /= 'uncoupled')  THEN
     621    IF ( coupling_mode /= 'uncoupled'  .AND.                                   &
     622         coupling_mode(1:8) /= 'vnested_' )  THEN
    614623
    615624       IF ( dt_coupling == 9999999.9_wp )  THEN
    … …  
    809818    IF ( coupling_mode == 'uncoupled' )  THEN
    810819       coupling_string = ''
     820    ELSEIF ( coupling_mode == 'vnested_crse' )  THEN
     821       coupling_string = ' nested (coarse)'
     822    ELSEIF ( coupling_mode == 'vnested_fine' )  THEN
     823       coupling_string = ' nested (fine)'
    811824    ELSEIF ( coupling_mode == 'atmosphere_to_ocean' )  THEN
    812825       coupling_string = ' coupled (atmosphere)'
    … …  
    38583871    ENDIF
    38593872
     3873!-- Vertical nesting: check fine and coarse grid compatibility for data exchange
     3874    IF ( vnested )  CALL vnest_check_parameters
     3875
    38603876    CALL location_message( 'finished', .TRUE. )
    38613877!

  • palm/trunk/SOURCE/diffusivities.f90

    r2233 r2365  
    2525! -----------------
    2626! $Id$
     27! Vertical grid nesting implemented (SadiqHuq)
     28!
     29! 2233 2017-05-30 18:08:54Z suehring
    2730!
    2831! 2232 2017-05-30 17:47:52Z suehring
    … …  
    108111        ONLY :  bc_h
    109112
     113    USE vertical_nesting_mod,                                                  &
     114        ONLY:  vnest_boundary_conds_khkm, vnest_init
     115
    110116    IMPLICIT NONE
    111117
    … …  
    266272       kh(:,nyn+1,:) = kh(:,nyn,:)
    267273    ENDIF
     274!
     275!-- Vertical nesting: set fine grid eddy viscosity top boundary condition
     276    IF ( vnest_init )  CALL vnest_boundary_conds_khkm
    268277
    269278 END SUBROUTINE diffusivities

  • palm/trunk/SOURCE/init_coupling.f90

    r2298 r2365  
    2525! ------------------
    2626! $Id$
     27! Vertical nesting implemented (SadiqHuq)
     28!
     29! 2298 2017-06-29 09:28:18Z raasch
    2730! MPI2 coupling removed
    2831!
    … …  
    6770    USE pegrid
    6871
     72    USE vertical_nesting_mod
     73
    6974    IMPLICIT NONE
    7075
    … …  
    96101       IF ( TRIM( coupling_mode ) == 'coupled_run' )  THEN
    97102          i = 1
     103       ELSEIF ( TRIM( coupling_mode ) == 'vnested_twi' )  THEN
     104          i = 9
    98105       ELSE
    99106          i = 0
    … …  
    111118    i = bc_data(0)
    112119
    113     IF ( i == 0 ) THEN
     120    IF ( i == 0 )  THEN
    114121       coupling_mode = 'uncoupled'
    115122!
    … …  
    121128          CLOSE ( 90 )
    122129       ENDIF
     130    ELSEIF ( i == 9 )  THEN
     131
     132!
     133!--    Set a flag to identify runs with vertical nesting
     134       vnested = .TRUE.
     135       
     136       comm_inter = MPI_COMM_WORLD
     137       
     138!
     139!--    Split the total available PE's into two groups
     140!--    numprocs for Coarse and Fine Grid are specified via mrun argument -N
     141       IF ( myid < bc_data(1) )  THEN
     142          inter_color     = 0
     143          numprocs        = bc_data(1)
     144          coupling_mode   = 'vnested_crse'
     145       ELSE
     146          inter_color     = 1
     147          numprocs        = bc_data(2)
     148          coupling_mode   = 'vnested_fine'
     149       ENDIF
     150       
     151       CALL MPI_COMM_SPLIT( MPI_COMM_WORLD, inter_color, 0, comm_palm, ierr )
     152       comm2d = comm_palm
     153       
     154       OPEN( 90, FILE='VNESTING_PORT_OPENED', FORM='FORMATTED' )
     155       WRITE ( 90, '(''TRUE'')' )
     156       CLOSE ( 90 )
     157     
    123158    ELSE
    124159       comm_inter = MPI_COMM_WORLD
    … …  
    154189    ENDIF
    155190
     191    IF (  TRIM( coupling_mode ) == 'vnested_fine' )  THEN
     192!
     193!-- Set file extension for vertical nesting
     194       coupling_char = '_N'
     195    ENDIF
     196
    156197 END SUBROUTINE init_coupling

  • palm/trunk/SOURCE/init_grid.f90

    r2319 r2365  
    2525! -----------------
    2626! $Id$
     27! Vertical nesting implemented (SadiqHuq)
     28!
     29! 2319 2017-07-20 17:33:17Z suehring
    2730! Remove print statements
    2831!
    … …  
    250253               canyon_height, canyon_wall_left, canyon_wall_south,             &
    251254               canyon_width_x, canyon_width_y, constant_flux_layer,            &
    252                coupling_char, dp_level_ind_b, dz, dz_max, dz_stretch_factor,   &
     255               coupling_char, coupling_mode,                                   &
     256               dp_level_ind_b, dz, dz_max, dz_stretch_factor,                  &
    253257               dz_stretch_level, dz_stretch_level_index, grid_level, ibc_uv_b, &
    254258               io_blocks, io_group, inflow_l, inflow_n, inflow_r, inflow_s,    &
    … …  
    285289    USE surface_mod,                                                           &
    286290        ONLY:  get_topography_top_index, init_bc
     291
     292    USE vertical_nesting_mod,                                                  &
     293        ONLY:  vnested, vnest_init_grid
    287294
    288295    IMPLICIT NONE
    … …  
    20002007    CALL exchange_horiz( l_wall, nbgp )     
    20012008
     2009!
     2010!-- Vertical nesting: communicate vertical grid level arrays between fine and
     2011!-- coarse grid
     2012    IF ( vnested )  CALL vnest_init_grid
    20022013
    20032014 END SUBROUTINE init_grid

  • palm/trunk/SOURCE/init_pegrid.f90

    r2300 r2365  
    2525! -----------------
    2626! $Id$
     27! Vertical nesting implemented (SadiqHuq)
     28!
     29! 2300 2017-06-29 13:31:14Z raasch
    2730! host-specific settings removed
    2831!
    … …  
    228231        ONLY:  nxl_y, nxl_yd, nxl_z, nxr_y, nxr_yd, nxr_z, nyn_x, nyn_z, nys_x,&
    229232               nys_z, nzb_x, nzb_y, nzb_yd, nzt_x, nzt_yd, nzt_y
     233
     234    USE vertical_nesting_mod,                                                  &
     235        ONLY:  vnested, vnest_init_pegrid_domain, vnest_init_pegrid_rank
    230236
    231237    IMPLICIT NONE
    … …  
    335341
    336342!
     343!-- Vertical nesting: store four lists that identify partner ranks to exchange
     344!-- data
     345    IF ( vnested )  CALL vnest_init_pegrid_rank
     346
     347!
    337348!-- Determine sub-topologies for transpositions
    338349!-- Transposition from z to x:
    … …  
    642653    CALL MPI_TYPE_COMMIT( type_xy, ierr )
    643654
    644     IF ( TRIM( coupling_mode ) /= 'uncoupled' )  THEN
     655    IF ( TRIM( coupling_mode ) /= 'uncoupled' .AND. .NOT. vnested )  THEN
    645656   
    646657!
    … …  
    743754    ENDIF
    744755
     756!
     757!-- Store partner grid point co-ordinates as lists.
     758!-- Create custom MPI vector datatypes for contiguous data transfer
     759    IF ( vnested )  CALL vnest_init_pegrid_domain
    745760
    746761#endif

  • palm/trunk/SOURCE/local_stop.f90

    r2101 r2365  
    2525! -----------------
    2626! $Id$
     27! Vertical nesting implemented (SadiqHuq)
     28!
     29! 2101 2017-01-05 16:42:31Z suehring
    2730!
    2831! 2000 2016-08-20 18:09:15Z knoop
    … …  
    7578        ONLY:  nested_run
    7679
     80    USE vertical_nesting_mod,                                                  &
     81        ONLY:  vnested
     82
    7783#if defined( __parallel )
    7884    IF ( coupling_mode == 'uncoupled' )  THEN
    … …  
    9298          ENDIF
    9399       ENDIF
     100    ELSEIF ( coupling_mode(1:8) == 'vnested_' )  THEN
     101
     102       PRINT*, '+++ local_stop:'
     103       PRINT*, '     model "', TRIM( coupling_mode ), '" terminated'
     104!
     105!--    Abort both coarse and fine grid
     106       CALL MPI_ABORT( MPI_COMM_WORLD, 9999, ierr )
    94107    ELSE
    95108

  • palm/trunk/SOURCE/parin.f90

    r2339 r2365  
    2525! -----------------
    2626! $Id$
     27! Vertical grid nesting: add vnest_start_time to d3par (SadiqHuq)
     28!
     29! 2339 2017-08-07 13:55:26Z gronemeier
    2730! corrected timestamp in header
    2831!
    … …  
    354357        ONLY:  wtm_parin
    355358
     359    USE vertical_nesting_mod,                                                  &
     360        ONLY:  vnest_start_time
    356361
    357362    IMPLICIT NONE
    … …  
    445450             skip_time_do2d_xy, skip_time_do2d_xz, skip_time_do2d_yz,          &
    446451             skip_time_do3d, skip_time_domask, synchronous_exchange,           &
    447              termination_time_needed, z_max_do2d
     452             termination_time_needed, vnest_start_time, z_max_do2d
    448453
    449454

  • palm/trunk/SOURCE/read_var_list.f90

    r2339 r2365  
    2525! -----------------
    2626! $Id$
     27! Vertical grid nesting implemented (SadiqHuq)
     28!
     29! 2339 2017-08-07 13:55:26Z gronemeier
    2730! corrected timestamp in header
    2831!
    … …  
    251254    USE synthetic_turbulence_generator_mod,                                    &
    252255        ONLY:  stg_read_restart_data
     256
     257    USE vertical_nesting_mod,                                                  &
     258        ONLY:  vnest_init
    253259
    254260    IMPLICIT NONE
    … …  
    785791             READ ( 13 )  vg_vertical_gradient_level_ind
    786792          CASE ( 'virtual_flight' )
    787              READ ( 13 ) virtual_flight
     793             READ ( 13 )  virtual_flight
     794          CASE ( 'vnest_init' )
     795             READ ( 13 )  vnest_init
    788796          CASE ( 'volume_flow_area' )
    789797             READ ( 13 )  volume_flow_area

  • palm/trunk/SOURCE/time_integration.f90

    r2320 r2365  
    2525! -----------------
    2626! $Id$
     27! Vertical grid nesting implemented (SadiqHuq)
     28!
     29! 2320 2017-07-21 12:47:43Z suehring
    2730! Set bottom boundary conditions after nesting interpolation and anterpolation
    2831!
    … …  
    399402        ONLY:  wind_turbine, wtm_forces
    400403
     404    USE vertical_nesting_mod,                                                  &
     405        ONLY:  vnested, vnest_anterpolate, vnest_anterpolate_e,                &
     406               vnest_boundary_conds, vnest_boundary_conds_khkm,                &
     407               vnest_deallocate, vnest_init, vnest_init_fine,                  &
     408               vnest_start_time
     409
    401410    IMPLICIT NONE
    402411
    … …  
    429438!-- Data exchange between coupled models in case that a call has been omitted
    430439!-- at the end of the previous run of a job chain.
    431     IF ( coupling_mode /= 'uncoupled'  .AND.  run_coupled )  THEN
     440    IF ( coupling_mode /= 'uncoupled'  .AND.  run_coupled .AND. .NOT. vnested)  THEN
    432441!
    433442!--    In case of model termination initiated by the local model the coupler
    … …  
    458467       CALL cpu_log( log_point_s(10), 'timesteps', 'start' )
    459468!
     469!--    Vertical nesting: initialize fine grid
     470       IF ( vnested ) THEN
     471          IF ( .NOT. vnest_init  .AND.  simulated_time >= vnest_start_time )  THEN
     472             CALL cpu_log( log_point(80), 'vnest_init', 'start' )
     473             CALL vnest_init_fine
     474             vnest_init = .TRUE.
     475             CALL cpu_log( log_point(80), 'vnest_init', 'stop' )
     476          ENDIF
     477       ENDIF
     478!
    460479!--    Determine ug, vg and w_subs in dependence on data from external file
    461480!--    LSF_DATA
    … …  
    624643!--       Swap the time levels in preparation for the next time step.
    625644          CALL swap_timelevel
     645
     646!
     647!--       Vertical nesting: Interpolate fine grid data to the coarse grid
     648          IF ( vnest_init ) THEN
     649             CALL cpu_log( log_point(81), 'vnest_anterpolate', 'start' )
     650             CALL vnest_anterpolate
     651             CALL cpu_log( log_point(81), 'vnest_anterpolate', 'stop' )
     652          ENDIF
    626653
    627654          IF ( nested_run )  THEN
    … …  
    736763          IF ( intermediate_timestep_count == 1  .OR. &
    737764                call_psolver_at_all_substeps )  THEN
    738              CALL pres
     765
     766             IF (  vnest_init ) THEN
     767!
     768!--             Compute pressure in the CG, interpolate top boundary conditions
     769!--             to the FG and then compute pressure in the FG
     770                IF ( coupling_mode == 'vnested_crse' )  CALL pres
     771
     772                CALL cpu_log( log_point(82), 'vnest_bc', 'start' )
     773                CALL vnest_boundary_conds
     774                CALL cpu_log( log_point(82), 'vnest_bc', 'stop' )
     775 
     776                IF ( coupling_mode == 'vnested_fine' )  CALL pres
     777
     778!--             Anterpolate TKE, satisfy Germano Identity
     779                CALL cpu_log( log_point(83), 'vnest_anter_e', 'start' )
     780                CALL vnest_anterpolate_e
     781                CALL cpu_log( log_point(83), 'vnest_anter_e', 'stop' )
     782
     783             ELSE
     784
     785                CALL pres
     786
     787             ENDIF
     788
    739789          ENDIF
    740790
    … …  
    910960!
    911961!--    Data exchange between coupled models
    912        IF ( coupling_mode /= 'uncoupled'  .AND.  run_coupled )  THEN
     962       IF ( coupling_mode /= 'uncoupled'  .AND.  run_coupled                   &
     963                                          .AND. .NOT. vnested )  THEN
    913964          time_coupling = time_coupling + dt_3d
    914965
    … …  
    11541205    ENDDO   ! time loop
    11551206
     1207!-- Vertical nesting: Deallocate variables initialized for vertical nesting   
     1208    IF ( vnest_init )  CALL vnest_deallocate
     1209
    11561210    IF ( myid == 0 )  CALL finish_progress_bar
    11571211

  • palm/trunk/SOURCE/timestep.f90

    r2258 r2365  
    2525! -----------------
    2626! $Id$
     27! Vertical grid nesting: Sync fine and coarse grid timestep (SadiqHuq)
     28!
     29! 2258 2017-06-08 07:55:13Z suehring
    2730! Bugfix, add pre-preprocessor directives to enable non-parrallel mode
    2831!
    … …  
    137140        ONLY:  flow_statistics_called, hom, u_max, u_max_ijk, v_max, v_max_ijk,&
    138141               w_max, w_max_ijk
     142
     143    USE vertical_nesting_mod,                                                  &
     144        ONLY:  vnested, vnest_timestep_sync
    139145
    140146    IMPLICIT NONE
    … …  
    376382    ENDIF
    377383
     384!
     385!-- Vertical nesting: coarse and fine grid timestep has to be identical   
     386    IF ( vnested )  CALL vnest_timestep_sync
     387
    378388    CALL cpu_log( log_point(12), 'calculate_timestep', 'stop' )
    379389

  • palm/trunk/SOURCE/write_var_list.f90

    r2339 r2365  
    2525! -----------------
    2626! $Id$
     27! Vertical nesting implemented (SadiqHuq)
     28!
     29! 2339 2017-08-07 13:55:26Z gronemeier
    2730! corrected timestamp in header
    2831!
    … …  
    224227               v_max_ijk, w_max, w_max_ijk
    225228
     229    USE vertical_nesting_mod,                                                  &
     230        ONLY:  vnest_init
    226231
    227232    IMPLICIT NONE
    … …  
    695700    WRITE ( 14 )  'virtual_flight                '
    696701    WRITE ( 14 )  virtual_flight
     702    WRITE ( 14 )  'vnest_init                    '
     703    WRITE ( 14 )  vnest_init
    697704    WRITE ( 14 )  'volume_flow_area              '
    698705    WRITE ( 14 )  volume_flow_area
Note: See TracChangeset for help on using the changeset viewer.
                                                                                                                                                                                                                                                                                                                                                                               
  | Impressum | ©Leibniz Universität Hannover |