Changeset 4783 for palm/trunk/SOURCE
- Timestamp:
- Nov 13, 2020 1:58:45 PM (4 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 5 edited
-
indoor_model_mod.f90 (modified) (1 diff)
-
init_3d_model.f90 (modified) (3 diffs)
-
parin.f90 (modified) (3 diffs)
-
radiation_model_mod.f90 (modified) (1 diff)
-
urban_surface_mod.f90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/indoor_model_mod.f90
r4780 r4783 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! Change parameters for summer_pars and winter_pars (responsible: S. Rissmann)22 ! 23 23 ! 24 24 ! Former revisions: 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Change parameters for summer_pars and winter_pars (responsible: S. Rissmann) 28 ! 29 ! 4780 2020-11-10 11:17:10Z suehring 27 30 ! Enable 3D data output also with 64-bit precision 28 31 ! -
palm/trunk/SOURCE/init_3d_model.f90
r4680 r4783 24 24 ! ----------------- 25 25 ! $Id$ 26 ! bugfix for reading restart data with MPI-I/O (does not work with blockwise I/O) 27 ! 28 ! 4680 2020-09-16 10:20:34Z gronemeier 26 29 ! Add option to fix date or time of the simulation 27 30 ! … … 1086 1089 IF ( TRIM( initializing_actions ) == 'cyclic_fill' ) THEN 1087 1090 1091 ! 1092 !-- Blockwise I/O does not work together with MPI-I/O 1093 IF ( restart_data_format_input(1:3) == 'mpi' ) THEN 1094 CALL rrd_read_parts_of_global 1095 ELSE 1096 DO i = 0, io_blocks-1 1097 IF ( i == io_group ) THEN 1098 CALL rrd_read_parts_of_global 1099 ENDIF 1100 #if defined( __parallel ) 1101 CALL MPI_BARRIER( comm2d, ierr ) 1102 #endif 1103 ENDDO 1104 ENDIF 1105 1106 ENDIF 1107 1108 ! 1109 !-- Read processor specific binary data from restart file. 1110 !-- Blockwise I/O does not work together with MPI-I/O 1111 IF ( restart_data_format_input(1:3) == 'mpi' ) THEN 1112 CALL rrd_local 1113 ELSE 1088 1114 DO i = 0, io_blocks-1 1089 1115 IF ( i == io_group ) THEN 1090 CALL rrd_ read_parts_of_global1116 CALL rrd_local 1091 1117 ENDIF 1092 1118 #if defined( __parallel ) … … 1094 1120 #endif 1095 1121 ENDDO 1096 1097 ENDIF 1098 1099 ! 1100 !-- Read processor specific binary data from restart file 1101 DO i = 0, io_blocks-1 1102 IF ( i == io_group ) THEN 1103 CALL rrd_local 1104 ENDIF 1105 #if defined( __parallel ) 1106 CALL MPI_BARRIER( comm2d, ierr ) 1107 #endif 1108 ENDDO 1122 ENDIF 1109 1123 1110 1124 -
palm/trunk/SOURCE/parin.f90
r4680 r4783 25 25 ! ----------------- 26 26 ! $Id$ 27 ! bugfix for reading restart data with MPI-I/O (does not work with blockwise I/O) 28 ! 29 ! 4680 2020-09-16 10:20:34Z gronemeier 27 30 ! Add option to fix date or time of the simulation 28 31 ! … … 393 396 394 397 CALL location_message( 'reading environment parameters from ENVPAR', 'finished' ) 398 399 CALL location_message( 'reading NAMELIST parameters from PARIN', 'start' ) 400 401 ! 402 !-- Open the NAMELIST-file which is send with this job 403 CALL check_open( 11 ) 404 405 ! 406 !-- Read the control parameters for initialization. 407 !-- The namelist "inipar" must be provided in the NAMELIST-file. 408 READ ( 11, initialization_parameters, ERR=10, END=11 ) 409 GOTO 14 410 411 10 BACKSPACE( 11 ) 412 READ( 11 , '(A)') line 413 CALL parin_fail_message( 'initialization_parameters', line ) 414 415 11 REWIND ( 11 ) 416 READ ( 11, inipar, ERR=12, END=13 ) 417 418 message_string = 'namelist inipar is deprecated and will be ' // & 419 'removed in near future. & Please use namelist ' // & 420 'initialization_parameters instead' 421 CALL message( 'parin', 'PA0017', 0, 1, 0, 6, 0 ) 422 423 GOTO 14 424 425 12 BACKSPACE( 11 ) 426 READ( 11 , '(A)') line 427 CALL parin_fail_message( 'inipar', line ) 428 429 13 message_string = 'no initialization_parameters-namelist found' 430 CALL message( 'parin', 'PA0272', 1, 2, 0, 6, 0 ) 431 432 ! 433 !-- Try to read runtime parameters given by the user for this run 434 !-- (namelist "runtime_parameters"). The namelist "runtime_parmeters" 435 !-- can be omitted. In that case default values are used for the 436 !-- parameters. 437 14 line = ' ' 438 439 REWIND ( 11 ) 440 line = ' ' 441 DO WHILE ( INDEX( line, '&runtime_parameters' ) == 0 ) 442 READ ( 11, '(A)', END=16 ) line 443 ENDDO 444 BACKSPACE ( 11 ) 445 446 ! 447 !-- Read namelist 448 READ ( 11, runtime_parameters, ERR = 15 ) 449 GOTO 18 450 451 15 BACKSPACE( 11 ) 452 READ( 11 , '(A)') line 453 CALL parin_fail_message( 'runtime_parameters', line ) 454 455 16 REWIND ( 11 ) 456 line = ' ' 457 DO WHILE ( INDEX( line, '&d3par' ) == 0 ) 458 READ ( 11, '(A)', END=18 ) line 459 ENDDO 460 BACKSPACE ( 11 ) 461 462 ! 463 !-- Read namelist 464 READ ( 11, d3par, ERR = 17, END = 18 ) 465 466 message_string = 'namelist d3par is deprecated and will be ' // & 467 'removed in near future. &Please use namelist ' // & 468 'runtime_parameters instead' 469 CALL message( 'parin', 'PA0487', 0, 1, 0, 6, 0 ) 470 471 GOTO 18 472 473 17 BACKSPACE( 11 ) 474 READ( 11 , '(A)') line 475 CALL parin_fail_message( 'd3par', line ) 476 477 18 CONTINUE 478 479 ! 480 !-- Check for module namelists and read them 481 CALL module_interface_parin 482 395 483 ! 396 484 !-- Calculate the number of groups into which parallel I/O is split. … … 424 512 io_group = MOD( global_id+1, io_blocks ) 425 513 426 CALL location_message( 'reading NAMELIST parameters from PARIN', 'start' ) 427 ! 428 !-- Data is read in parallel by groups of PEs 429 DO i = 0, io_blocks-1 430 IF ( i == io_group ) THEN 431 432 ! 433 !-- Open the NAMELIST-file which is send with this job 434 CALL check_open( 11 ) 435 436 ! 437 !-- Read the control parameters for initialization. 438 !-- The namelist "inipar" must be provided in the NAMELIST-file. 439 READ ( 11, initialization_parameters, ERR=10, END=11 ) 440 GOTO 14 441 442 10 BACKSPACE( 11 ) 443 READ( 11 , '(A)') line 444 CALL parin_fail_message( 'initialization_parameters', line ) 445 446 11 REWIND ( 11 ) 447 READ ( 11, inipar, ERR=12, END=13 ) 448 449 message_string = 'namelist inipar is deprecated and will be ' // & 450 'removed in near future. & Please use namelist ' // & 451 'initialization_parameters instead' 452 CALL message( 'parin', 'PA0017', 0, 1, 0, 6, 0 ) 453 454 GOTO 14 455 456 12 BACKSPACE( 11 ) 457 READ( 11 , '(A)') line 458 CALL parin_fail_message( 'inipar', line ) 459 460 13 message_string = 'no initialization_parameters-namelist found' 461 CALL message( 'parin', 'PA0272', 1, 2, 0, 6, 0 ) 462 463 ! 464 !-- Try to read runtime parameters given by the user for this run 465 !-- (namelist "runtime_parameters"). The namelist "runtime_parmeters" 466 !-- can be omitted. In that case default values are used for the 467 !-- parameters. 468 14 line = ' ' 469 470 REWIND ( 11 ) 471 line = ' ' 472 DO WHILE ( INDEX( line, '&runtime_parameters' ) == 0 ) 473 READ ( 11, '(A)', END=16 ) line 514 ! 515 !-- If required, read control parameters from restart file (produced by 516 !-- a prior run). All PEs are reading from file created by PE0 (see 517 !-- check_open) 518 IF ( TRIM( initializing_actions ) == 'read_restart_data' ) THEN 519 520 ! 521 !-- If not set by the user in the runtime parameters, the data format for restart 522 !-- input needs to be set now! This is normally done later in check parameters. 523 IF ( TRIM( restart_data_format ) /= 'fortran_binary' .AND. & 524 TRIM( restart_data_format ) /= 'mpi' .AND. & 525 TRIM( restart_data_format ) /= 'mpi_shared_memory' ) THEN 526 message_string = 'illegal restart data format "' // TRIM( restart_data_format ) // '"' 527 CALL message( 'parin', 'PA0724', 1, 2, 0, 6, 0 ) 528 ENDIF 529 IF ( TRIM( restart_data_format_input ) == 'undefined' ) THEN 530 restart_data_format_input = restart_data_format 531 ENDIF 532 533 ! 534 !-- Blockwise I/O does not work together with MPI-I/O 535 IF ( restart_data_format_input(1:3) == 'mpi' .AND. io_blocks /= 1 ) THEN 536 CALL rrd_global 537 ELSE 538 ! 539 !-- Data is read in parallel by groups of PEs 540 DO i = 0, io_blocks-1 541 IF ( i == io_group ) THEN 542 CALL rrd_global 543 ENDIF 544 #if defined( __parallel ) 545 CALL MPI_BARRIER( comm2d, ierr ) 546 #endif 474 547 ENDDO 475 BACKSPACE ( 11 ) 476 477 ! 478 !-- Read namelist 479 READ ( 11, runtime_parameters, ERR = 15 ) 480 GOTO 18 481 482 15 BACKSPACE( 11 ) 483 READ( 11 , '(A)') line 484 CALL parin_fail_message( 'runtime_parameters', line ) 485 486 16 REWIND ( 11 ) 487 line = ' ' 488 DO WHILE ( INDEX( line, '&d3par' ) == 0 ) 489 READ ( 11, '(A)', END=18 ) line 490 ENDDO 491 BACKSPACE ( 11 ) 492 493 ! 494 !-- Read namelist 495 READ ( 11, d3par, ERR = 17, END = 18 ) 496 497 message_string = 'namelist d3par is deprecated and will be ' // & 498 'removed in near future. &Please use namelist ' // & 499 'runtime_parameters instead' 500 CALL message( 'parin', 'PA0487', 0, 1, 0, 6, 0 ) 501 502 GOTO 18 503 504 17 BACKSPACE( 11 ) 505 READ( 11 , '(A)') line 506 CALL parin_fail_message( 'd3par', line ) 507 508 18 CONTINUE 509 510 ! 511 !-- Check for module namelists and read them 512 CALL module_interface_parin 513 514 ! 515 !-- If required, read control parameters from restart file (produced by 516 !-- a prior run). All PEs are reading from file created by PE0 (see 517 !-- check_open) 518 IF ( TRIM( initializing_actions ) == 'read_restart_data' ) THEN 519 520 ! 521 !-- If not set by the user in the runtime parameters, the data format for restart 522 !-- input needs to be set now! This is normally done later in check parameters. 523 IF ( TRIM( restart_data_format ) /= 'fortran_binary' .AND. & 524 TRIM( restart_data_format ) /= 'mpi' .AND. & 525 TRIM( restart_data_format ) /= 'mpi_shared_memory' ) THEN 526 message_string = 'illegal restart data format "' // TRIM( restart_data_format ) // '"' 527 CALL message( 'parin', 'PA0724', 1, 2, 0, 6, 0 ) 528 ENDIF 529 IF ( TRIM( restart_data_format_input ) == 'undefined' ) THEN 530 restart_data_format_input = restart_data_format 531 ENDIF 532 533 CALL rrd_global 534 ! 535 !-- Increment the run count 536 runnr = runnr + 1 537 ! 538 !-- In case of a restart run, the number of user-defined profiles on 539 !-- the restart file (already stored in max_pr_user) has to match the 540 !-- one given for the current run. max_pr_user_tmp is calculated in 541 !-- user_parin and max_pr_user is read in via rrd_global. 542 IF ( max_pr_user /= max_pr_user_tmp ) THEN 543 WRITE( message_string, * ) 'the number of user-defined ', & 544 'profiles given in data_output_pr (', max_pr_user_tmp, & 545 ') does not match the one ', & 546 'found in the restart file (', max_pr_user, ')' 547 CALL message( 'user_parin', 'UI0009', 1, 2, 0, 6, 0 ) 548 ENDIF 549 ELSE 550 max_pr_user = max_pr_user_tmp 548 ENDIF 549 550 551 ! 552 !-- Increment the run count 553 runnr = runnr + 1 554 ! 555 !-- In case of a restart run, the number of user-defined profiles on 556 !-- the restart file (already stored in max_pr_user) has to match the 557 !-- one given for the current run. max_pr_user_tmp is calculated in 558 !-- user_parin and max_pr_user is read in via rrd_global. 559 IF ( max_pr_user /= max_pr_user_tmp ) THEN 560 WRITE( message_string, * ) 'the number of user-defined ', & 561 'profiles given in data_output_pr (', max_pr_user_tmp, & 562 ') does not match the one ', & 563 'found in the restart file (', max_pr_user, ')' 564 CALL message( 'user_parin', 'UI0009', 1, 2, 0, 6, 0 ) 565 ENDIF 566 ELSE 567 max_pr_user = max_pr_user_tmp 568 ENDIF 569 ! 570 !-- Activate spinup 571 IF ( land_surface .OR. urban_surface ) THEN 572 IF ( spinup_time > 0.0_wp ) THEN 573 coupling_start_time = spinup_time 574 time_since_reference_point = simulated_time - coupling_start_time 575 IF ( spinup_pt_mean == 9999999.9_wp ) THEN 576 spinup_pt_mean = pt_surface 551 577 ENDIF 552 ! 553 !-- Activate spinup 554 IF ( land_surface .OR. urban_surface ) THEN 555 IF ( spinup_time > 0.0_wp ) THEN 556 coupling_start_time = spinup_time 557 time_since_reference_point = simulated_time - coupling_start_time 558 IF ( spinup_pt_mean == 9999999.9_wp ) THEN 559 spinup_pt_mean = pt_surface 560 ENDIF 561 end_time = end_time + spinup_time 562 IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN 563 spinup = .TRUE. 564 ELSEIF ( TRIM( initializing_actions ) == 'read_restart_data' & 565 .AND. time_since_reference_point > 0.0_wp ) THEN 566 data_output_during_spinup = .FALSE. !< required for correct ntdim calculation 567 !< in check_parameters for restart run 568 ENDIF 569 ENDIF 578 end_time = end_time + spinup_time 579 IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN 580 spinup = .TRUE. 581 ELSEIF ( TRIM( initializing_actions ) == 'read_restart_data' & 582 .AND. time_since_reference_point > 0.0_wp ) THEN 583 data_output_during_spinup = .FALSE. !< required for correct ntdim calculation 584 !< in check_parameters for restart run 570 585 ENDIF 571 572 ! 573 !-- In case of nested runs, explicitly set nesting boundary conditions. 574 !-- This will overwrite the user settings and basic defaults. 575 !-- bc_lr and bc_ns always need to be cyclic for vertical nesting. 576 IF ( nested_run ) THEN 577 IF ( nesting_mode == 'vertical' ) THEN 578 IF (bc_lr /= 'cyclic' .OR. bc_ns /= 'cyclic' ) THEN 579 WRITE ( message_string, *) 'bc_lr and bc_ns were set to ,', & 580 'cyclic for vertical nesting' 581 CALL message( 'parin', 'PA0428', 0, 0, 0, 6, 0 ) 582 bc_lr = 'cyclic' 583 bc_ns = 'cyclic' 584 ENDIF 585 IF ( child_domain ) THEN 586 bc_uv_t = 'nested' 587 bc_pt_t = 'nested' 588 bc_q_t = 'nested' 589 bc_s_t = 'nested' 590 bc_cs_t = 'nested' 591 bc_p_t = 'neumann' 592 ENDIF 593 ! 594 !-- For other nesting modes only set boundary conditions for 595 !-- nested domains. 596 ELSE 597 IF ( child_domain ) THEN 598 bc_lr = 'nested' 599 bc_ns = 'nested' 600 bc_uv_t = 'nested' 601 bc_pt_t = 'nested' 602 bc_q_t = 'nested' 603 bc_s_t = 'nested' 604 bc_cs_t = 'nested' 605 bc_p_t = 'neumann' 606 ENDIF 607 ENDIF 586 ENDIF 587 ENDIF 588 589 ! 590 !-- In case of nested runs, explicitly set nesting boundary conditions. 591 !-- This will overwrite the user settings and basic defaults. 592 !-- bc_lr and bc_ns always need to be cyclic for vertical nesting. 593 IF ( nested_run ) THEN 594 IF ( nesting_mode == 'vertical' ) THEN 595 IF (bc_lr /= 'cyclic' .OR. bc_ns /= 'cyclic' ) THEN 596 WRITE ( message_string, *) 'bc_lr and bc_ns were set to ,', & 597 'cyclic for vertical nesting' 598 CALL message( 'parin', 'PA0428', 0, 0, 0, 6, 0 ) 599 bc_lr = 'cyclic' 600 bc_ns = 'cyclic' 608 601 ENDIF 609 ! 610 !-- Set boundary conditions also in case the model is offline-nested in 611 !-- larger-scale models. 612 IF ( nesting_offline ) THEN 613 bc_lr = 'nesting_offline' 614 bc_ns = 'nesting_offline' 615 bc_uv_t = 'nesting_offline' 616 bc_pt_t = 'nesting_offline' 617 bc_q_t = 'nesting_offline' 618 ! bc_s_t = 'nesting_offline' ! scalar boundary condition is not clear yet 619 ! bc_cs_t = 'nesting_offline' ! same for chemical species 602 IF ( child_domain ) THEN 603 bc_uv_t = 'nested' 604 bc_pt_t = 'nested' 605 bc_q_t = 'nested' 606 bc_s_t = 'nested' 607 bc_cs_t = 'nested' 620 608 bc_p_t = 'neumann' 621 609 ENDIF 622 623 ! 624 !-- In case of nested runs, make sure that initializing_actions = 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 629 !-- are set properly. An exception is made in case of restart runs and 630 !-- if user decides to do everything by its own. 631 IF ( child_domain .AND. .NOT. ( & 632 TRIM( initializing_actions ) == 'read_restart_data' .OR. & 633 TRIM( initializing_actions ) == 'set_constant_profiles' .OR. & 634 TRIM( initializing_actions ) == 'by_user' ) ) THEN 635 message_string = 'initializing_actions = ' // & 636 TRIM( initializing_actions ) // ' has been ' // & 637 'changed to set_constant_profiles in child ' // & 638 'domain.' 639 CALL message( 'parin', 'PA0492', 0, 0, 0, 6, 0 ) 640 641 initializing_actions = 'set_constant_profiles' 610 ! 611 !-- For other nesting modes only set boundary conditions for 612 !-- nested domains. 613 ELSE 614 IF ( child_domain ) THEN 615 bc_lr = 'nested' 616 bc_ns = 'nested' 617 bc_uv_t = 'nested' 618 bc_pt_t = 'nested' 619 bc_q_t = 'nested' 620 bc_s_t = 'nested' 621 bc_cs_t = 'nested' 622 bc_p_t = 'neumann' 642 623 ENDIF 643 ! 644 !-- Check validity of lateral boundary conditions. This has to be done 645 !-- here because they are already used in init_pegrid and init_grid and 646 !-- therefore cannot be check in check_parameters 647 IF ( bc_lr /= 'cyclic' .AND. bc_lr /= 'dirichlet/radiation' .AND. & 648 bc_lr /= 'radiation/dirichlet' .AND. bc_lr /= 'nested' .AND. & 649 bc_lr /= 'nesting_offline' ) THEN 650 message_string = 'unknown boundary condition: bc_lr = "' // & 651 TRIM( bc_lr ) // '"' 652 CALL message( 'parin', 'PA0049', 1, 2, 0, 6, 0 ) 653 ENDIF 654 IF ( bc_ns /= 'cyclic' .AND. bc_ns /= 'dirichlet/radiation' .AND. & 655 bc_ns /= 'radiation/dirichlet' .AND. bc_ns /= 'nested' .AND. & 656 bc_ns /= 'nesting_offline' ) THEN 657 message_string = 'unknown boundary condition: bc_ns = "' // & 658 TRIM( bc_ns ) // '"' 659 CALL message( 'parin', 'PA0050', 1, 2, 0, 6, 0 ) 660 ENDIF 661 ! 662 !-- Set internal variables used for speed optimization in if clauses 663 IF ( bc_lr /= 'cyclic' ) bc_lr_cyc = .FALSE. 664 IF ( bc_lr == 'dirichlet/radiation' ) bc_lr_dirrad = .TRUE. 665 IF ( bc_lr == 'radiation/dirichlet' ) bc_lr_raddir = .TRUE. 666 IF ( bc_ns /= 'cyclic' ) bc_ns_cyc = .FALSE. 667 IF ( bc_ns == 'dirichlet/radiation' ) bc_ns_dirrad = .TRUE. 668 IF ( bc_ns == 'radiation/dirichlet' ) bc_ns_raddir = .TRUE. 669 ! 670 !-- Radiation-Dirichlet conditions are allowed along one of the horizontal directions only. 671 !-- In general, such conditions along x and y may work, but require a) some code changes 672 !-- (e.g. concerning allocation of c_u, c_v ... arrays), and b) a careful model setup by the 673 !-- user, in order to guarantee that there is a clearly defined outflow at two sides. 674 !-- Otherwise, the radiation condition may produce wrong results. 675 IF ( ( bc_lr_dirrad .OR. bc_lr_raddir ) .AND. ( bc_ns_dirrad .OR. bc_ns_raddir ) ) THEN 676 message_string = 'bc_lr = "' // TRIM( bc_lr ) // '" and bc_ns = "' // & 677 TRIM( bc_ns ) // '" are not allowed to be set at the same time' 678 CALL message( 'parin', 'PA0589', 1, 2, 0, 6, 0 ) 679 ENDIF 680 ! 681 !-- Check in case of initial run, if the grid point numbers are well 682 !-- defined and allocate some arrays which are already needed in 683 !-- init_pegrid or check_parameters. During restart jobs, these arrays 684 !-- will be allocated in rrd_global. All other arrays are allocated 685 !-- in init_3d_model. 686 IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN 687 688 IF ( nx <= 0 ) THEN 689 WRITE( message_string, * ) 'no value or wrong value given', & 690 ' for nx: nx=', nx 691 CALL message( 'parin', 'PA0273', 1, 2, 0, 6, 0 ) 692 ENDIF 693 IF ( ny <= 0 ) THEN 694 WRITE( message_string, * ) 'no value or wrong value given', & 695 ' for ny: ny=', ny 696 CALL message( 'parin', 'PA0274', 1, 2, 0, 6, 0 ) 697 ENDIF 698 IF ( nz <= 0 ) THEN 699 WRITE( message_string, * ) 'no value or wrong value given', & 700 ' for nz: nz=', nz 701 CALL message( 'parin', 'PA0275', 1, 2, 0, 6, 0 ) 702 ENDIF 703 704 ! 705 !-- As a side condition, routine flow_statistics require at least 14 706 !-- vertical grid levels (they are used to store time-series data) 707 !> @todo Remove this restriction 708 IF ( nz < 14 ) THEN 709 WRITE( message_string, * ) 'nz >= 14 is required' 710 CALL message( 'parin', 'PA0362', 1, 2, 0, 6, 0 ) 711 ENDIF 712 713 ! 714 !-- ATTENTION: in case of changes to the following statement please 715 !-- also check the allocate statement in routine rrd_global 716 ALLOCATE( pt_init(0:nz+1), q_init(0:nz+1), s_init(0:nz+1), & 717 ref_state(0:nz+1), sa_init(0:nz+1), ug(0:nz+1), & 718 u_init(0:nz+1), v_init(0:nz+1), vg(0:nz+1), & 719 hom(0:nz+1,2,pr_palm+max_pr_user+max_pr_cs+max_pr_salsa,0:statistic_regions), & 720 hom_sum(0:nz+1,pr_palm+max_pr_user+max_pr_cs+max_pr_salsa,0:statistic_regions) ) 721 722 hom = 0.0_wp 723 724 ENDIF 725 726 ! 727 !-- NAMELIST-file is not needed anymore 728 CALL close_file( 11 ) 729 730 ENDIF 731 #if defined( __parallel ) 732 CALL MPI_BARRIER( MPI_COMM_WORLD, ierr ) 733 #endif 734 ENDDO 624 ENDIF 625 ENDIF 626 ! 627 !-- Set boundary conditions also in case the model is offline-nested in 628 !-- larger-scale models. 629 IF ( nesting_offline ) THEN 630 bc_lr = 'nesting_offline' 631 bc_ns = 'nesting_offline' 632 bc_uv_t = 'nesting_offline' 633 bc_pt_t = 'nesting_offline' 634 bc_q_t = 'nesting_offline' 635 ! bc_s_t = 'nesting_offline' ! scalar boundary condition is not clear yet 636 ! bc_cs_t = 'nesting_offline' ! same for chemical species 637 bc_p_t = 'neumann' 638 ENDIF 639 640 ! 641 !-- In case of nested runs, make sure that initializing_actions = 642 !-- 'set_constant_profiles' even though the constant-profiles 643 !-- initializations for the prognostic variables will be overwritten 644 !-- by pmci_child_initialize and pmci_parent_initialize. This is, 645 !-- however, important e.g. to make sure that diagnostic variables 646 !-- are set properly. An exception is made in case of restart runs and 647 !-- if user decides to do everything by its own. 648 IF ( child_domain .AND. .NOT. ( & 649 TRIM( initializing_actions ) == 'read_restart_data' .OR. & 650 TRIM( initializing_actions ) == 'set_constant_profiles' .OR. & 651 TRIM( initializing_actions ) == 'by_user' ) ) THEN 652 message_string = 'initializing_actions = ' // & 653 TRIM( initializing_actions ) // ' has been ' // & 654 'changed to set_constant_profiles in child ' // & 655 'domain.' 656 CALL message( 'parin', 'PA0492', 0, 0, 0, 6, 0 ) 657 658 initializing_actions = 'set_constant_profiles' 659 ENDIF 660 ! 661 !-- Check validity of lateral boundary conditions. This has to be done 662 !-- here because they are already used in init_pegrid and init_grid and 663 !-- therefore cannot be check in check_parameters 664 IF ( bc_lr /= 'cyclic' .AND. bc_lr /= 'dirichlet/radiation' .AND. & 665 bc_lr /= 'radiation/dirichlet' .AND. bc_lr /= 'nested' .AND. & 666 bc_lr /= 'nesting_offline' ) THEN 667 message_string = 'unknown boundary condition: bc_lr = "' // & 668 TRIM( bc_lr ) // '"' 669 CALL message( 'parin', 'PA0049', 1, 2, 0, 6, 0 ) 670 ENDIF 671 IF ( bc_ns /= 'cyclic' .AND. bc_ns /= 'dirichlet/radiation' .AND. & 672 bc_ns /= 'radiation/dirichlet' .AND. bc_ns /= 'nested' .AND. & 673 bc_ns /= 'nesting_offline' ) THEN 674 message_string = 'unknown boundary condition: bc_ns = "' // & 675 TRIM( bc_ns ) // '"' 676 CALL message( 'parin', 'PA0050', 1, 2, 0, 6, 0 ) 677 ENDIF 678 ! 679 !-- Set internal variables used for speed optimization in if clauses 680 IF ( bc_lr /= 'cyclic' ) bc_lr_cyc = .FALSE. 681 IF ( bc_lr == 'dirichlet/radiation' ) bc_lr_dirrad = .TRUE. 682 IF ( bc_lr == 'radiation/dirichlet' ) bc_lr_raddir = .TRUE. 683 IF ( bc_ns /= 'cyclic' ) bc_ns_cyc = .FALSE. 684 IF ( bc_ns == 'dirichlet/radiation' ) bc_ns_dirrad = .TRUE. 685 IF ( bc_ns == 'radiation/dirichlet' ) bc_ns_raddir = .TRUE. 686 ! 687 !-- Radiation-Dirichlet conditions are allowed along one of the horizontal directions only. 688 !-- In general, such conditions along x and y may work, but require a) some code changes 689 !-- (e.g. concerning allocation of c_u, c_v ... arrays), and b) a careful model setup by the 690 !-- user, in order to guarantee that there is a clearly defined outflow at two sides. 691 !-- Otherwise, the radiation condition may produce wrong results. 692 IF ( ( bc_lr_dirrad .OR. bc_lr_raddir ) .AND. ( bc_ns_dirrad .OR. bc_ns_raddir ) ) THEN 693 message_string = 'bc_lr = "' // TRIM( bc_lr ) // '" and bc_ns = "' // & 694 TRIM( bc_ns ) // '" are not allowed to be set at the same time' 695 CALL message( 'parin', 'PA0589', 1, 2, 0, 6, 0 ) 696 ENDIF 697 ! 698 !-- Check in case of initial run, if the grid point numbers are well 699 !-- defined and allocate some arrays which are already needed in 700 !-- init_pegrid or check_parameters. During restart jobs, these arrays 701 !-- will be allocated in rrd_global. All other arrays are allocated 702 !-- in init_3d_model. 703 IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN 704 705 IF ( nx <= 0 ) THEN 706 WRITE( message_string, * ) 'no value or wrong value given', & 707 ' for nx: nx=', nx 708 CALL message( 'parin', 'PA0273', 1, 2, 0, 6, 0 ) 709 ENDIF 710 IF ( ny <= 0 ) THEN 711 WRITE( message_string, * ) 'no value or wrong value given', & 712 ' for ny: ny=', ny 713 CALL message( 'parin', 'PA0274', 1, 2, 0, 6, 0 ) 714 ENDIF 715 IF ( nz <= 0 ) THEN 716 WRITE( message_string, * ) 'no value or wrong value given', & 717 ' for nz: nz=', nz 718 CALL message( 'parin', 'PA0275', 1, 2, 0, 6, 0 ) 719 ENDIF 720 721 ! 722 !-- As a side condition, routine flow_statistics require at least 14 723 !-- vertical grid levels (they are used to store time-series data) 724 !> @todo Remove this restriction 725 IF ( nz < 14 ) THEN 726 WRITE( message_string, * ) 'nz >= 14 is required' 727 CALL message( 'parin', 'PA0362', 1, 2, 0, 6, 0 ) 728 ENDIF 729 730 ! 731 !-- ATTENTION: in case of changes to the following statement please 732 !-- also check the allocate statement in routine rrd_global 733 ALLOCATE( pt_init(0:nz+1), q_init(0:nz+1), s_init(0:nz+1), & 734 ref_state(0:nz+1), sa_init(0:nz+1), ug(0:nz+1), & 735 u_init(0:nz+1), v_init(0:nz+1), vg(0:nz+1), & 736 hom(0:nz+1,2,pr_palm+max_pr_user+max_pr_cs+max_pr_salsa,0:statistic_regions), & 737 hom_sum(0:nz+1,pr_palm+max_pr_user+max_pr_cs+max_pr_salsa,0:statistic_regions) ) 738 739 hom = 0.0_wp 740 741 ENDIF 742 743 ! 744 !-- NAMELIST-file is not needed anymore 745 CALL close_file( 11 ) 735 746 736 747 CALL location_message( 'reading NAMELIST parameters from PARIN', 'finished' ) -
palm/trunk/SOURCE/radiation_model_mod.f90
r4780 r4783 23 23 ! Current revisions: 24 24 ! ------------------ 25 ! Albedo parameter list extended - building facades added25 ! 26 26 ! 27 27 ! Former revisions: 28 28 ! ----------------- 29 29 ! $Id$ 30 ! Albedo parameter list extended - building facades added 31 ! 32 ! 4780 2020-11-10 11:17:10Z suehring 30 33 ! Enable 3D data output also with 64-bit precision 31 34 ! -
palm/trunk/SOURCE/urban_surface_mod.f90
r4780 r4783 22 22 ! Current revisions: 23 23 ! ----------------- 24 ! 25 ! 26 ! Former revisions: 27 ! ----------------- 28 ! $Id$ 24 29 ! - Default building parameters updated and extended (responsible: S. Rissmann) 25 30 ! - First and second wall layer initialized with individual building properties rather than with 26 31 ! with the same properties (heat capacity and conductivity) 27 32 ! 28 ! Former revisions: 29 ! ----------------- 30 ! $Id$ 33 ! 4780 2020-11-10 11:17:10Z suehring 31 34 ! Enable 3D data output also with 64-bit precision 32 35 !
Note: See TracChangeset
for help on using the changeset viewer.
