Changeset 4680 for palm/trunk/SOURCE/parin.f90
- Timestamp:
- Sep 16, 2020 10:20:34 AM (4 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/parin.f90 (modified) (15 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/parin.f90
r4565 r4680 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 23 ! 22 ! 23 ! 24 24 ! Former revisions: 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Add option to fix date or time of the simulation 28 ! 29 ! 4565 2020-06-15 08:30:38Z oliver.maas 27 30 ! added pt_surface_heating_rate 28 ! 31 ! 29 32 ! 4564 2020-06-12 14:03:36Z raasch 30 33 ! Vertical nesting method of Huq et al. (2019) removed 31 ! 34 ! 32 35 ! 4536 2020-05-17 17:24:13Z raasch 33 36 ! bugfix for restart data format query 34 ! 37 ! 35 38 ! 4505 2020-04-20 15:37:15Z schwenkel 36 39 ! Add flag for saturation check 37 ! 40 ! 38 41 ! 4495 2020-04-13 20:11:20Z raasch 39 42 ! restart data handling with MPI-IO added 40 ! 43 ! 41 44 ! 4360 2020-01-07 11:25:50Z suehring 42 45 ! removed recycling_yshift 43 ! 46 ! 44 47 ! 4227 2019-09-10 18:04:34Z gronemeier 45 48 ! implement new palm_date_time_mod 46 ! 49 ! 47 50 ! 4146 2019-08-07 07:47:36Z gronemeier 48 51 ! added rotation_angle to initialization_parameters 49 ! 52 ! 50 53 ! 4191 2019-08-27 15:45:07Z gronemeier 51 54 ! bugfix: add recycling_method_for_thermodynamic_quantities to inipar namelist 52 ! 55 ! 53 56 ! 4183 2019-08-23 07:33:16Z oliver.maas 54 57 ! replaced recycle_absolute_quantities by recycling_method_for_thermodynamic_quantities 55 ! 58 ! 56 59 ! 4182 2019-08-22 15:20:23Z scharf 57 60 ! Corrected "Former revisions" section 58 ! 61 ! 59 62 ! 4176 2019-08-20 14:10:41Z oliver.maas 60 63 ! added recycle_absolute_quantities to initialization_parameters namelist 61 ! 64 ! 62 65 ! 4173 2019-08-20 12:04:06Z gronemeier 63 66 ! add vdi_internal_controls 64 ! 67 ! 65 68 ! 4131 2019-08-02 11:06:18Z monakurppa 66 69 ! Allocate hom and hom_sum to allow profile output for salsa variables. 67 ! 70 ! 68 71 ! 4079 2019-07-09 18:04:41Z suehring 69 72 ! +monotonic_limiter_z 70 ! 73 ! 71 74 ! 4022 2019-06-12 11:52:39Z suehring 72 75 ! Change default top boundary condition for pressure to Neumann in offline 73 76 ! nesting case 74 ! 77 ! 75 78 ! 4017 2019-06-06 12:16:46Z schwenkel 76 79 ! Introduce alternative switch for debug output during timestepping 77 ! 80 ! 78 81 ! 3885 2019-04-11 11:29:34Z kanani 79 ! Changes related to global restructuring of location messages and introduction 82 ! Changes related to global restructuring of location messages and introduction 80 83 ! of additional debug messages 81 ! 84 ! 82 85 ! 3806 2019-03-21 12:45:50Z raasch 83 86 ! additional check for lateral boundary conditions added 84 ! 87 ! 85 88 ! 3747 2019-02-16 15:15:23Z gronemeier 86 89 ! removed setting of parameter region 87 ! 90 ! 88 91 ! 3746 2019-02-16 12:41:27Z gronemeier 89 92 ! Removed most_method 90 ! 93 ! 91 94 ! 3649 2019-01-02 16:52:21Z suehring 92 95 ! Delete debug-print statements … … 103 106 !------------------------------------------------------------------------------! 104 107 SUBROUTINE parin 105 108 106 109 107 110 USE arrays_3d, & … … 156 159 IMPLICIT NONE 157 160 158 CHARACTER (LEN=80) :: line !< dummy string that contains the current line of the parameter file 161 CHARACTER (LEN=80) :: line !< dummy string that contains the current line of the parameter file 159 162 160 163 INTEGER(iwp) :: global_id !< process id with respect to MPI_COMM_WORLD … … 224 227 turbulent_inflow, turbulent_outflow, & 225 228 use_subsidence_tendencies, ug_surface, ug_vertical_gradient, & 229 use_fixed_date, use_fixed_time, & 226 230 use_free_convection_scaling, & 227 231 ug_vertical_gradient_level, use_surface_fluxes, use_cmax, & … … 298 302 ug_vertical_gradient_level, use_surface_fluxes, use_cmax, & 299 303 use_top_fluxes, use_ug_for_galilei_tr, use_upstream_for_tke, & 304 use_fixed_date, use_fixed_time, & 300 305 use_free_convection_scaling, & 301 306 uv_heights, u_bulk, u_profile, vdi_checks, & … … 304 309 wall_adjustment, wall_heatflux, wall_humidityflux, & 305 310 wall_scalarflux, y_shift, zeta_max, zeta_min, z0h_factor 306 311 307 312 NAMELIST /d3par/ averaging_interval, averaging_interval_pr, & 308 313 cpu_log_barrierwait, create_disturbances, & … … 418 423 io_blocks = global_procs / maximum_parallel_io_streams 419 424 io_group = MOD( global_id+1, io_blocks ) 420 425 421 426 CALL location_message( 'reading NAMELIST parameters from PARIN', 'start' ) 422 427 ! … … 434 439 READ ( 11, initialization_parameters, ERR=10, END=11 ) 435 440 GOTO 14 436 441 437 442 10 BACKSPACE( 11 ) 438 443 READ( 11 , '(A)') line … … 441 446 11 REWIND ( 11 ) 442 447 READ ( 11, inipar, ERR=12, END=13 ) 443 448 444 449 message_string = 'namelist inipar is deprecated and will be ' // & 445 450 'removed in near future. & Please use namelist ' // & 446 451 'initialization_parameters instead' 447 452 CALL message( 'parin', 'PA0017', 0, 1, 0, 6, 0 ) 448 453 449 454 GOTO 14 450 455 451 456 12 BACKSPACE( 11 ) 452 457 READ( 11 , '(A)') line … … 458 463 ! 459 464 !-- Try to read runtime parameters given by the user for this run 460 !-- (namelist "runtime_parameters"). The namelist "runtime_parmeters" 461 !-- can be omitted. In that case default values are used for the 465 !-- (namelist "runtime_parameters"). The namelist "runtime_parmeters" 466 !-- can be omitted. In that case default values are used for the 462 467 !-- parameters. 463 468 14 line = ' ' … … 533 538 !-- In case of a restart run, the number of user-defined profiles on 534 539 !-- the restart file (already stored in max_pr_user) has to match the 535 !-- one given for the current run. max_pr_user_tmp is calculated in 540 !-- one given for the current run. max_pr_user_tmp is calculated in 536 541 !-- user_parin and max_pr_user is read in via rrd_global. 537 542 IF ( max_pr_user /= max_pr_user_tmp ) THEN … … 566 571 567 572 ! 568 !-- In case of nested runs, explicitly set nesting boundary conditions. 573 !-- In case of nested runs, explicitly set nesting boundary conditions. 569 574 !-- This will overwrite the user settings and basic defaults. 570 575 !-- bc_lr and bc_ns always need to be cyclic for vertical nesting. … … 584 589 bc_s_t = 'nested' 585 590 bc_cs_t = 'nested' 586 bc_p_t = 'neumann' 591 bc_p_t = 'neumann' 587 592 ENDIF 588 593 ! 589 !-- For other nesting modes only set boundary conditions for 594 !-- For other nesting modes only set boundary conditions for 590 595 !-- nested domains. 591 ELSE 596 ELSE 592 597 IF ( child_domain ) THEN 593 598 bc_lr = 'nested' … … 616 621 ENDIF 617 622 618 ! 623 ! 619 624 !-- In case of nested runs, make sure that initializing_actions = 620 !-- 'set_constant_profiles' even though the constant-profiles 621 !-- initializations for the prognostic variables will be overwritten 622 !-- by pmci_child_initialize and pmci_parent_initialize. This is, 623 !-- however, important e.g. to make sure that diagnostic variables 625 !-- 'set_constant_profiles' even though the constant-profiles 626 !-- initializations for the prognostic variables will be overwritten 627 !-- by pmci_child_initialize and pmci_parent_initialize. This is, 628 !-- however, important e.g. to make sure that diagnostic variables 624 629 !-- are set properly. An exception is made in case of restart runs and 625 630 !-- if user decides to do everything by its own. … … 631 636 TRIM( initializing_actions ) // ' has been ' // & 632 637 'changed to set_constant_profiles in child ' // & 633 'domain.' 638 'domain.' 634 639 CALL message( 'parin', 'PA0492', 0, 0, 0, 6, 0 ) 635 640 636 641 initializing_actions = 'set_constant_profiles' 637 ENDIF 642 ENDIF 638 643 ! 639 644 !-- Check validity of lateral boundary conditions. This has to be done
Note: See TracChangeset
for help on using the changeset viewer.
