Changeset 4226 for palm/trunk/SOURCE/nesting_offl_mod.f90
- Timestamp:
- Sep 10, 2019 5:03:24 PM (5 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/nesting_offl_mod.f90 (modified) (25 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/nesting_offl_mod.f90
r4182 r4226 25 25 ! ----------------- 26 26 ! $Id$ 27 ! - Data input moved into nesting_offl_mod 28 ! - check rephrased 29 ! - time variable is now relative to time_utc_init 30 ! - Define module specific data type for offline nesting in nesting_offl_mod 31 ! 32 ! 4182 2019-08-22 15:20:23Z scharf 27 33 ! Corrected "Former revisions" section 28 34 ! … … 99 105 100 106 USE arrays_3d, & 101 ONLY: dzw, e, diss, pt, pt_init, q, q_init, s, u, u_init, ug, v, & 102 v_init, vg, w, zu, zw 107 ONLY: dzw, & 108 e, & 109 diss, & 110 pt, & 111 pt_init, & 112 q, & 113 q_init, & 114 rdf, & 115 rdf_sc, & 116 s, & 117 u, & 118 u_init, & 119 ug, & 120 v, & 121 v_init, & 122 vg, & 123 w, & 124 zu, & 125 zw 126 127 USE basic_constants_and_equations_mod, & 128 ONLY: g, & 129 pi 103 130 104 131 USE chem_modules, & … … 111 138 bc_dirichlet_r, & 112 139 bc_dirichlet_s, & 140 coupling_char, & 113 141 dt_3d, & 114 142 dz, & … … 132 160 133 161 USE cpulog, & 134 ONLY: cpu_log, log_point 162 ONLY: cpu_log, & 163 log_point, & 164 log_point_s 165 166 USE date_and_time_mod, & 167 ONLY: time_utc_init 135 168 136 169 USE grid_variables … … 145 178 146 179 USE netcdf_data_input_mod, & 147 ONLY: nest_offl 148 180 ONLY: check_existence, & 181 close_input_file, & 182 get_dimension_length, & 183 get_variable, & 184 get_variable_pr, & 185 input_pids_dynamic, & 186 inquire_num_variables, & 187 inquire_variable_names, & 188 input_file_dynamic, & 189 num_var_pids, & 190 open_read_file, & 191 pids_id 192 149 193 USE pegrid 150 194 195 IMPLICIT NONE 196 197 ! 198 !-- Define data type for nesting in larger-scale models like COSMO. 199 !-- Data type comprises u, v, w, pt, and q at lateral and top boundaries. 200 TYPE nest_offl_type 201 202 CHARACTER(LEN=16) :: char_l = 'ls_forcing_left_' !< leading substring for variables at left boundary 203 CHARACTER(LEN=17) :: char_n = 'ls_forcing_north_' !< leading substring for variables at north boundary 204 CHARACTER(LEN=17) :: char_r = 'ls_forcing_right_' !< leading substring for variables at right boundary 205 CHARACTER(LEN=17) :: char_s = 'ls_forcing_south_' !< leading substring for variables at south boundary 206 CHARACTER(LEN=15) :: char_t = 'ls_forcing_top_' !< leading substring for variables at top boundary 207 208 CHARACTER(LEN=100), DIMENSION(:), ALLOCATABLE :: var_names !< list of variable in dynamic input file 209 CHARACTER(LEN=100), DIMENSION(:), ALLOCATABLE :: var_names_chem_l !< names of mesoscale nested chemistry variables at left boundary 210 CHARACTER(LEN=100), DIMENSION(:), ALLOCATABLE :: var_names_chem_n !< names of mesoscale nested chemistry variables at north boundary 211 CHARACTER(LEN=100), DIMENSION(:), ALLOCATABLE :: var_names_chem_r !< names of mesoscale nested chemistry variables at right boundary 212 CHARACTER(LEN=100), DIMENSION(:), ALLOCATABLE :: var_names_chem_s !< names of mesoscale nested chemistry variables at south boundary 213 CHARACTER(LEN=100), DIMENSION(:), ALLOCATABLE :: var_names_chem_t !< names of mesoscale nested chemistry variables at top boundary 214 215 INTEGER(iwp) :: nt !< number of time levels in dynamic input file 216 INTEGER(iwp) :: nzu !< number of vertical levels on scalar grid in dynamic input file 217 INTEGER(iwp) :: nzw !< number of vertical levels on w grid in dynamic input file 218 INTEGER(iwp) :: tind !< time index for reference time in mesoscale-offline nesting 219 INTEGER(iwp) :: tind_p !< time index for following time in mesoscale-offline nesting 220 221 LOGICAL :: init = .FALSE. !< flag indicating that offline nesting is already initialized 222 223 LOGICAL, DIMENSION(:), ALLOCATABLE :: chem_from_file_l !< flags inidicating whether left boundary data for chemistry is in dynamic input file 224 LOGICAL, DIMENSION(:), ALLOCATABLE :: chem_from_file_n !< flags inidicating whether north boundary data for chemistry is in dynamic input file 225 LOGICAL, DIMENSION(:), ALLOCATABLE :: chem_from_file_r !< flags inidicating whether right boundary data for chemistry is in dynamic input file 226 LOGICAL, DIMENSION(:), ALLOCATABLE :: chem_from_file_s !< flags inidicating whether south boundary data for chemistry is in dynamic input file 227 LOGICAL, DIMENSION(:), ALLOCATABLE :: chem_from_file_t !< flags inidicating whether top boundary data for chemistry is in dynamic input file 228 229 REAL(wp), DIMENSION(:), ALLOCATABLE :: surface_pressure !< time dependent surface pressure 230 REAL(wp), DIMENSION(:), ALLOCATABLE :: time !< time levels in dynamic input file 231 REAL(wp), DIMENSION(:), ALLOCATABLE :: zu_atmos !< vertical levels at scalar grid in dynamic input file 232 REAL(wp), DIMENSION(:), ALLOCATABLE :: zw_atmos !< vertical levels at w grid in dynamic input file 233 234 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: ug !< domain-averaged geostrophic component 235 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: vg !< domain-averaged geostrophic component 236 237 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: u_left !< u-component at left boundary 238 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: v_left !< v-component at left boundary 239 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: w_left !< w-component at left boundary 240 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: q_left !< mixing ratio at left boundary 241 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: pt_left !< potentital temperautre at left boundary 242 243 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: u_north !< u-component at north boundary 244 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: v_north !< v-component at north boundary 245 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: w_north !< w-component at north boundary 246 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: q_north !< mixing ratio at north boundary 247 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: pt_north !< potentital temperautre at north boundary 248 249 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: u_right !< u-component at right boundary 250 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: v_right !< v-component at right boundary 251 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: w_right !< w-component at right boundary 252 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: q_right !< mixing ratio at right boundary 253 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: pt_right !< potentital temperautre at right boundary 254 255 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: u_south !< u-component at south boundary 256 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: v_south !< v-component at south boundary 257 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: w_south !< w-component at south boundary 258 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: q_south !< mixing ratio at south boundary 259 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: pt_south !< potentital temperautre at south boundary 260 261 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: u_top !< u-component at top boundary 262 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: v_top !< v-component at top boundary 263 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: w_top !< w-component at top boundary 264 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: q_top !< mixing ratio at top boundary 265 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: pt_top !< potentital temperautre at top boundary 266 267 REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: chem_left !< chemical species at left boundary 268 REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: chem_north !< chemical species at left boundary 269 REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: chem_right !< chemical species at left boundary 270 REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: chem_south !< chemical species at left boundary 271 REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: chem_top !< chemical species at left boundary 272 273 END TYPE nest_offl_type 274 275 REAL(wp) :: fac_dt !< interpolation factor 151 276 REAL(wp) :: zi_ribulk = 0.0_wp !< boundary-layer depth according to bulk Richardson criterion, i.e. the height where Ri_bulk exceeds the critical 152 !< bulk Richardson number of 0.25 277 !< bulk Richardson number of 0.2 278 279 TYPE(nest_offl_type) :: nest_offl !< data structure for data input at lateral and top boundaries (provided by Inifor) 153 280 154 281 SAVE … … 159 286 nesting_offl_calc_zi, & 160 287 nesting_offl_check_parameters, & 288 nesting_offl_geostrophic_wind, & 161 289 nesting_offl_header, & 162 290 nesting_offl_init, & 291 nesting_offl_input, & 292 nesting_offl_interpolation_factor, & 163 293 nesting_offl_mass_conservation, & 164 294 nesting_offl_parin … … 178 308 MODULE PROCEDURE nesting_offl_check_parameters 179 309 END INTERFACE nesting_offl_check_parameters 310 311 INTERFACE nesting_offl_geostrophic_wind 312 MODULE PROCEDURE nesting_offl_geostrophic_wind 313 END INTERFACE nesting_offl_geostrophic_wind 180 314 181 315 INTERFACE nesting_offl_header … … 186 320 MODULE PROCEDURE nesting_offl_init 187 321 END INTERFACE nesting_offl_init 322 323 INTERFACE nesting_offl_input 324 MODULE PROCEDURE nesting_offl_input 325 END INTERFACE nesting_offl_input 326 327 INTERFACE nesting_offl_interpolation_factor 328 MODULE PROCEDURE nesting_offl_interpolation_factor 329 END INTERFACE nesting_offl_interpolation_factor 188 330 189 331 INTERFACE nesting_offl_mass_conservation … … 197 339 CONTAINS 198 340 341 !------------------------------------------------------------------------------! 342 ! Description: 343 ! ------------ 344 !> Reads data at lateral and top boundaries derived from larger-scale model. 345 !------------------------------------------------------------------------------! 346 SUBROUTINE nesting_offl_input 347 348 INTEGER(iwp) :: n !< running index for chemistry variables 349 INTEGER(iwp) :: t !< running index time dimension 350 351 ! 352 !-- Initialize INIFOR forcing in first call. 353 IF ( .NOT. nest_offl%init ) THEN 354 #if defined ( __netcdf ) 355 ! 356 !-- Open file in read-only mode 357 CALL open_read_file( TRIM( input_file_dynamic ) // & 358 TRIM( coupling_char ), pids_id ) 359 ! 360 !-- At first, inquire all variable names. 361 CALL inquire_num_variables( pids_id, num_var_pids ) 362 ! 363 !-- Allocate memory to store variable names. 364 ALLOCATE( nest_offl%var_names(1:num_var_pids) ) 365 CALL inquire_variable_names( pids_id, nest_offl%var_names ) 366 ! 367 !-- Read time dimension, allocate memory and finally read time array 368 CALL get_dimension_length( pids_id, nest_offl%nt, 'time' ) 369 370 IF ( check_existence( nest_offl%var_names, 'time' ) ) THEN 371 ALLOCATE( nest_offl%time(0:nest_offl%nt-1) ) 372 CALL get_variable( pids_id, 'time', nest_offl%time ) 373 ENDIF 374 ! 375 !-- Read vertical dimension of scalar und w grid 376 CALL get_dimension_length( pids_id, nest_offl%nzu, 'z' ) 377 CALL get_dimension_length( pids_id, nest_offl%nzw, 'zw' ) 378 379 IF ( check_existence( nest_offl%var_names, 'z' ) ) THEN 380 ALLOCATE( nest_offl%zu_atmos(1:nest_offl%nzu) ) 381 CALL get_variable( pids_id, 'z', nest_offl%zu_atmos ) 382 ENDIF 383 IF ( check_existence( nest_offl%var_names, 'zw' ) ) THEN 384 ALLOCATE( nest_offl%zw_atmos(1:nest_offl%nzw) ) 385 CALL get_variable( pids_id, 'zw', nest_offl%zw_atmos ) 386 ENDIF 387 ! 388 !-- Read surface pressure 389 IF ( check_existence( nest_offl%var_names, & 390 'surface_forcing_surface_pressure' ) ) THEN 391 ALLOCATE( nest_offl%surface_pressure(0:nest_offl%nt-1) ) 392 CALL get_variable( pids_id, & 393 'surface_forcing_surface_pressure', & 394 nest_offl%surface_pressure ) 395 ENDIF 396 ! 397 !-- Close input file 398 CALL close_input_file( pids_id ) 399 #endif 400 ENDIF 401 ! 402 !-- Check if dynamic driver data input is required. 403 IF ( nest_offl%time(nest_offl%tind_p) <= & 404 MAX( time_since_reference_point, 0.0_wp) + time_utc_init .OR. & 405 .NOT. nest_offl%init ) THEN 406 CONTINUE 407 ! 408 !-- Return otherwise 409 ELSE 410 RETURN 411 ENDIF 412 ! 413 !-- CPU measurement 414 CALL cpu_log( log_point_s(86), 'NetCDF input forcing', 'start' ) 415 416 ! 417 !-- Obtain time index for current point in time. Note, the time coordinate 418 !-- in the input file is relative to time_utc_init. Since time_since_... 419 !-- is negativ when spinup is used, use MAX function to obtain correct 420 !-- time at the beginning. 421 nest_offl%tind = MINLOC( ABS( nest_offl%time - ( & 422 time_utc_init + & 423 MAX( time_since_reference_point, 0.0_wp) )& 424 ), DIM = 1 ) - 1 425 nest_offl%tind_p = nest_offl%tind + 1 426 ! 427 !-- Open file in read-only mode 428 #if defined ( __netcdf ) 429 CALL open_read_file( TRIM( input_file_dynamic ) // & 430 TRIM( coupling_char ), pids_id ) 431 ! 432 !-- Read geostrophic wind components 433 DO t = nest_offl%tind, nest_offl%tind_p 434 CALL get_variable_pr( pids_id, 'ls_forcing_ug', t+1, & 435 nest_offl%ug(t-nest_offl%tind,nzb+1:nzt) ) 436 CALL get_variable_pr( pids_id, 'ls_forcing_vg', t+1, & 437 nest_offl%vg(t-nest_offl%tind,nzb+1:nzt) ) 438 ENDDO 439 ! 440 !-- Read data at lateral and top boundaries. Please note, at left and 441 !-- right domain boundary, yz-layers are read for u, v, w, pt and q. 442 !-- For the v-component, the data starts at nysv, while for the other 443 !-- quantities the data starts at nys. This is equivalent at the north 444 !-- and south domain boundary for the u-component. 445 !-- Note, lateral data is also accessed by parallel IO, which is the reason 446 !-- why different arguments are passed depending on the boundary control 447 !-- flags. Cores that do not belong to the respective boundary just make 448 !-- a dummy read with count = 0, just in order to participate the collective 449 !-- operation. 450 !-- Read data for western boundary 451 CALL get_variable( pids_id, 'ls_forcing_left_u', & 452 nest_offl%u_left, & ! array to be read 453 MERGE( nys+1, 1, bc_dirichlet_l), & ! start index y direction 454 MERGE( nzb+1, 1, bc_dirichlet_l), & ! start index z direction 455 MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), & ! start index time dimension 456 MERGE( nyn-nys+1, 0, bc_dirichlet_l), & ! number of elements along y 457 MERGE( nest_offl%nzu, 0, bc_dirichlet_l), & ! number of elements alogn z 458 MERGE( 2, 0, bc_dirichlet_l), & ! number of time steps (2 or 0) 459 .TRUE. ) ! parallel IO when compiled accordingly 460 461 CALL get_variable( pids_id, 'ls_forcing_left_v', & 462 nest_offl%v_left, & 463 MERGE( nysv, 1, bc_dirichlet_l), & 464 MERGE( nzb+1, 1, bc_dirichlet_l), & 465 MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), & 466 MERGE( nyn-nysv+1, 0, bc_dirichlet_l), & 467 MERGE( nest_offl%nzu, 0, bc_dirichlet_l), & 468 MERGE( 2, 0, bc_dirichlet_l), & 469 .TRUE. ) 470 471 CALL get_variable( pids_id, 'ls_forcing_left_w', & 472 nest_offl%w_left, & 473 MERGE( nys+1, 1, bc_dirichlet_l), & 474 MERGE( nzb+1, 1, bc_dirichlet_l), & 475 MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), & 476 MERGE( nyn-nys+1, 0, bc_dirichlet_l), & 477 MERGE( nest_offl%nzw, 0, bc_dirichlet_l), & 478 MERGE( 2, 0, bc_dirichlet_l), & 479 .TRUE. ) 480 481 IF ( .NOT. neutral ) THEN 482 CALL get_variable( pids_id, 'ls_forcing_left_pt', & 483 nest_offl%pt_left, & 484 MERGE( nys+1, 1, bc_dirichlet_l), & 485 MERGE( nzb+1, 1, bc_dirichlet_l), & 486 MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), & 487 MERGE( nyn-nys+1, 0, bc_dirichlet_l), & 488 MERGE( nest_offl%nzu, 0, bc_dirichlet_l), & 489 MERGE( 2, 0, bc_dirichlet_l), & 490 .TRUE. ) 491 ENDIF 492 493 IF ( humidity ) THEN 494 CALL get_variable( pids_id, 'ls_forcing_left_qv', & 495 nest_offl%q_left, & 496 MERGE( nys+1, 1, bc_dirichlet_l), & 497 MERGE( nzb+1, 1, bc_dirichlet_l), & 498 MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), & 499 MERGE( nyn-nys+1, 0, bc_dirichlet_l), & 500 MERGE( nest_offl%nzu, 0, bc_dirichlet_l), & 501 MERGE( 2, 0, bc_dirichlet_l), & 502 .TRUE. ) 503 ENDIF 504 505 IF ( air_chemistry ) THEN 506 DO n = 1, UBOUND(nest_offl%var_names_chem_l, 1) 507 IF ( check_existence( nest_offl%var_names, & 508 nest_offl%var_names_chem_l(n) ) ) THEN 509 CALL get_variable( pids_id, & 510 TRIM( nest_offl%var_names_chem_l(n) ), & 511 nest_offl%chem_left(:,:,:,n), & 512 MERGE( nys+1, 1, bc_dirichlet_l), & 513 MERGE( nzb+1, 1, bc_dirichlet_l), & 514 MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), & 515 MERGE( nyn-nys+1, 0, bc_dirichlet_l), & 516 MERGE( nest_offl%nzu, 0, bc_dirichlet_l), & 517 MERGE( 2, 0, bc_dirichlet_l), & 518 .TRUE. ) 519 nest_offl%chem_from_file_l(n) = .TRUE. 520 ENDIF 521 ENDDO 522 ENDIF 523 ! 524 !-- Read data for eastern boundary 525 CALL get_variable( pids_id, 'ls_forcing_right_u', & 526 nest_offl%u_right, & 527 MERGE( nys+1, 1, bc_dirichlet_r), & 528 MERGE( nzb+1, 1, bc_dirichlet_r), & 529 MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), & 530 MERGE( nyn-nys+1, 0, bc_dirichlet_r), & 531 MERGE( nest_offl%nzu, 0, bc_dirichlet_r), & 532 MERGE( 2, 0, bc_dirichlet_r), & 533 .TRUE. ) 534 535 CALL get_variable( pids_id, 'ls_forcing_right_v', & 536 nest_offl%v_right, & 537 MERGE( nysv, 1, bc_dirichlet_r), & 538 MERGE( nzb+1, 1, bc_dirichlet_r), & 539 MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), & 540 MERGE( nyn-nysv+1, 0, bc_dirichlet_r), & 541 MERGE( nest_offl%nzu, 0, bc_dirichlet_r), & 542 MERGE( 2, 0, bc_dirichlet_r), & 543 .TRUE. ) 544 545 CALL get_variable( pids_id, 'ls_forcing_right_w', & 546 nest_offl%w_right, & 547 MERGE( nys+1, 1, bc_dirichlet_r), & 548 MERGE( nzb+1, 1, bc_dirichlet_r), & 549 MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), & 550 MERGE( nyn-nys+1, 0, bc_dirichlet_r), & 551 MERGE( nest_offl%nzw, 0, bc_dirichlet_r), & 552 MERGE( 2, 0, bc_dirichlet_r), & 553 .TRUE. ) 554 555 IF ( .NOT. neutral ) THEN 556 CALL get_variable( pids_id, 'ls_forcing_right_pt', & 557 nest_offl%pt_right, & 558 MERGE( nys+1, 1, bc_dirichlet_r), & 559 MERGE( nzb+1, 1, bc_dirichlet_r), & 560 MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), & 561 MERGE( nyn-nys+1, 0, bc_dirichlet_r), & 562 MERGE( nest_offl%nzu, 0, bc_dirichlet_r), & 563 MERGE( 2, 0, bc_dirichlet_r), & 564 .TRUE. ) 565 ENDIF 566 567 IF ( humidity ) THEN 568 CALL get_variable( pids_id, 'ls_forcing_right_qv', & 569 nest_offl%q_right, & 570 MERGE( nys+1, 1, bc_dirichlet_r), & 571 MERGE( nzb+1, 1, bc_dirichlet_r), & 572 MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), & 573 MERGE( nyn-nys+1, 0, bc_dirichlet_r), & 574 MERGE( nest_offl%nzu, 0, bc_dirichlet_r), & 575 MERGE( 2, 0, bc_dirichlet_r), & 576 .TRUE. ) 577 ENDIF 578 579 IF ( air_chemistry ) THEN 580 DO n = 1, UBOUND(nest_offl%var_names_chem_r, 1) 581 IF ( check_existence( nest_offl%var_names, & 582 nest_offl%var_names_chem_r(n) ) ) THEN 583 CALL get_variable( pids_id, & 584 TRIM( nest_offl%var_names_chem_r(n) ), & 585 nest_offl%chem_right(:,:,:,n), & 586 MERGE( nys+1, 1, bc_dirichlet_r), & 587 MERGE( nzb+1, 1, bc_dirichlet_r), & 588 MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), & 589 MERGE( nyn-nys+1, 0, bc_dirichlet_r), & 590 MERGE( nest_offl%nzu, 0, bc_dirichlet_r), & 591 MERGE( 2, 0, bc_dirichlet_r), & 592 .TRUE. ) 593 nest_offl%chem_from_file_r(n) = .TRUE. 594 ENDIF 595 ENDDO 596 ENDIF 597 ! 598 !-- Read data for northern boundary 599 CALL get_variable( pids_id, 'ls_forcing_north_u', & ! array to be read 600 nest_offl%u_north, & ! start index x direction 601 MERGE( nxlu, 1, bc_dirichlet_n ), & ! start index z direction 602 MERGE( nzb+1, 1, bc_dirichlet_n ), & ! start index time dimension 603 MERGE( nest_offl%tind+1, 1, bc_dirichlet_n ), & ! number of elements along x 604 MERGE( nxr-nxlu+1, 0, bc_dirichlet_n ), & ! number of elements alogn z 605 MERGE( nest_offl%nzu, 0, bc_dirichlet_n ), & ! number of time steps (2 or 0) 606 MERGE( 2, 0, bc_dirichlet_n ), & ! parallel IO when compiled accordingly 607 .TRUE. ) 608 609 CALL get_variable( pids_id, 'ls_forcing_north_v', & ! array to be read 610 nest_offl%v_north, & ! start index x direction 611 MERGE( nxl+1, 1, bc_dirichlet_n ), & ! start index z direction 612 MERGE( nzb+1, 1, bc_dirichlet_n ), & ! start index time dimension 613 MERGE( nest_offl%tind+1, 1, bc_dirichlet_n ), & ! number of elements along x 614 MERGE( nxr-nxl+1, 0, bc_dirichlet_n ), & ! number of elements alogn z 615 MERGE( nest_offl%nzu, 0, bc_dirichlet_n ), & ! number of time steps (2 or 0) 616 MERGE( 2, 0, bc_dirichlet_n ), & ! parallel IO when compiled accordingly 617 .TRUE. ) 618 619 CALL get_variable( pids_id, 'ls_forcing_north_w', & ! array to be read 620 nest_offl%w_north, & ! start index x direction 621 MERGE( nxl+1, 1, bc_dirichlet_n ), & ! start index z direction 622 MERGE( nzb+1, 1, bc_dirichlet_n ), & ! start index time dimension 623 MERGE( nest_offl%tind+1, 1, bc_dirichlet_n ), & ! number of elements along x 624 MERGE( nxr-nxl+1, 0, bc_dirichlet_n ), & ! number of elements alogn z 625 MERGE( nest_offl%nzw, 0, bc_dirichlet_n ), & ! number of time steps (2 or 0) 626 MERGE( 2, 0, bc_dirichlet_n ), & ! parallel IO when compiled accordingly 627 .TRUE. ) 628 629 IF ( .NOT. neutral ) THEN 630 CALL get_variable( pids_id, 'ls_forcing_north_pt', & ! array to be read 631 nest_offl%pt_north, & ! start index x direction 632 MERGE( nxl+1, 1, bc_dirichlet_n ), & ! start index z direction 633 MERGE( nzb+1, 1, bc_dirichlet_n ), & ! start index time dimension 634 MERGE( nest_offl%tind+1, 1, bc_dirichlet_n ), & ! number of elements along x 635 MERGE( nxr-nxl+1, 0, bc_dirichlet_n ), & ! number of elements alogn z 636 MERGE( nest_offl%nzu, 0, bc_dirichlet_n ), & ! number of time steps (2 or 0) 637 MERGE( 2, 0, bc_dirichlet_n ), & ! parallel IO when compiled accordingly 638 .TRUE. ) 639 ENDIF 640 IF ( humidity ) THEN 641 CALL get_variable( pids_id, 'ls_forcing_north_qv', & ! array to be read 642 nest_offl%q_north, & ! start index x direction 643 MERGE( nxl+1, 1, bc_dirichlet_n ), & ! start index z direction 644 MERGE( nzb+1, 1, bc_dirichlet_n ), & ! start index time dimension 645 MERGE( nest_offl%tind+1, 1, bc_dirichlet_n ), & ! number of elements along x 646 MERGE( nxr-nxl+1, 0, bc_dirichlet_n ), & ! number of elements alogn z 647 MERGE( nest_offl%nzu, 0, bc_dirichlet_n ), & ! number of time steps (2 or 0) 648 MERGE( 2, 0, bc_dirichlet_n ), & ! parallel IO when compiled accordingly 649 .TRUE. ) 650 ENDIF 651 652 IF ( air_chemistry ) THEN 653 DO n = 1, UBOUND(nest_offl%var_names_chem_n, 1) 654 IF ( check_existence( nest_offl%var_names, & 655 nest_offl%var_names_chem_n(n) ) ) THEN 656 CALL get_variable( pids_id, & 657 TRIM( nest_offl%var_names_chem_n(n) ), & 658 nest_offl%chem_north(:,:,:,n), & 659 MERGE( nxl+1, 1, bc_dirichlet_n ), & 660 MERGE( nzb+1, 1, bc_dirichlet_n ), & 661 MERGE( nest_offl%tind+1, 1, bc_dirichlet_n ), & 662 MERGE( nxr-nxl+1, 0, bc_dirichlet_n ), & 663 MERGE( nest_offl%nzu, 0, bc_dirichlet_n ), & 664 MERGE( 2, 0, bc_dirichlet_n ), & 665 .TRUE. ) 666 nest_offl%chem_from_file_n(n) = .TRUE. 667 ENDIF 668 ENDDO 669 ENDIF 670 ! 671 !-- Read data for southern boundary 672 CALL get_variable( pids_id, 'ls_forcing_south_u', & ! array to be read 673 nest_offl%u_south, & ! start index x direction 674 MERGE( nxlu, 1, bc_dirichlet_s ), & ! start index z direction 675 MERGE( nzb+1, 1, bc_dirichlet_s ), & ! start index time dimension 676 MERGE( nest_offl%tind+1, 1, bc_dirichlet_s ), & ! number of elements along x 677 MERGE( nxr-nxlu+1, 0, bc_dirichlet_s ), & ! number of elements alogn z 678 MERGE( nest_offl%nzu, 0, bc_dirichlet_s ), & ! number of time steps (2 or 0) 679 MERGE( 2, 0, bc_dirichlet_s ), & ! parallel IO when compiled accordingly 680 .TRUE. ) 681 682 CALL get_variable( pids_id, 'ls_forcing_south_v', & ! array to be read 683 nest_offl%v_south, & ! start index x direction 684 MERGE( nxl+1, 1, bc_dirichlet_s ), & ! start index z direction 685 MERGE( nzb+1, 1, bc_dirichlet_s ), & ! start index time dimension 686 MERGE( nest_offl%tind+1, 1, bc_dirichlet_s ), & ! number of elements along x 687 MERGE( nxr-nxl+1, 0, bc_dirichlet_s ), & ! number of elements alogn z 688 MERGE( nest_offl%nzu, 0, bc_dirichlet_s ), & ! number of time steps (2 or 0) 689 MERGE( 2, 0, bc_dirichlet_s ), & ! parallel IO when compiled accordingly 690 .TRUE. ) 691 692 CALL get_variable( pids_id, 'ls_forcing_south_w', & ! array to be read 693 nest_offl%w_south, & ! start index x direction 694 MERGE( nxl+1, 1, bc_dirichlet_s ), & ! start index z direction 695 MERGE( nzb+1, 1, bc_dirichlet_s ), & ! start index time dimension 696 MERGE( nest_offl%tind+1, 1, bc_dirichlet_s ), & ! number of elements along x 697 MERGE( nxr-nxl+1, 0, bc_dirichlet_s ), & ! number of elements alogn z 698 MERGE( nest_offl%nzw, 0, bc_dirichlet_s ), & ! number of time steps (2 or 0) 699 MERGE( 2, 0, bc_dirichlet_s ), & ! parallel IO when compiled accordingly 700 .TRUE. ) 701 702 IF ( .NOT. neutral ) THEN 703 CALL get_variable( pids_id, 'ls_forcing_south_pt', & ! array to be read 704 nest_offl%pt_south, & ! start index x direction 705 MERGE( nxl+1, 1, bc_dirichlet_s ), & ! start index z direction 706 MERGE( nzb+1, 1, bc_dirichlet_s ), & ! start index time dimension 707 MERGE( nest_offl%tind+1, 1, bc_dirichlet_s ), & ! number of elements along x 708 MERGE( nxr-nxl+1, 0, bc_dirichlet_s ), & ! number of elements alogn z 709 MERGE( nest_offl%nzu, 0, bc_dirichlet_s ), & ! number of time steps (2 or 0) 710 MERGE( 2, 0, bc_dirichlet_s ), & ! parallel IO when compiled accordingly 711 .TRUE. ) 712 ENDIF 713 IF ( humidity ) THEN 714 CALL get_variable( pids_id, 'ls_forcing_south_qv', & ! array to be read 715 nest_offl%q_south, & ! start index x direction 716 MERGE( nxl+1, 1, bc_dirichlet_s ), & ! start index z direction 717 MERGE( nzb+1, 1, bc_dirichlet_s ), & ! start index time dimension 718 MERGE( nest_offl%tind+1, 1, bc_dirichlet_s ), & ! number of elements along x 719 MERGE( nxr-nxl+1, 0, bc_dirichlet_s ), & ! number of elements alogn z 720 MERGE( nest_offl%nzu, 0, bc_dirichlet_s ), & ! number of time steps (2 or 0) 721 MERGE( 2, 0, bc_dirichlet_s ), & ! parallel IO when compiled accordingly 722 .TRUE. ) 723 ENDIF 724 725 IF ( air_chemistry ) THEN 726 DO n = 1, UBOUND(nest_offl%var_names_chem_s, 1) 727 IF ( check_existence( nest_offl%var_names, & 728 nest_offl%var_names_chem_s(n) ) ) THEN 729 CALL get_variable( pids_id, & 730 TRIM( nest_offl%var_names_chem_s(n) ), & 731 nest_offl%chem_south(:,:,:,n), & 732 MERGE( nxl+1, 1, bc_dirichlet_s ), & 733 MERGE( nzb+1, 1, bc_dirichlet_s ), & 734 MERGE( nest_offl%tind+1, 1, bc_dirichlet_s ), & 735 MERGE( nxr-nxl+1, 0, bc_dirichlet_s ), & 736 MERGE( nest_offl%nzu, 0, bc_dirichlet_s ), & 737 MERGE( 2, 0, bc_dirichlet_s ), & 738 .TRUE. ) 739 nest_offl%chem_from_file_s(n) = .TRUE. 740 ENDIF 741 ENDDO 742 ENDIF 743 ! 744 !-- Top boundary 745 CALL get_variable( pids_id, 'ls_forcing_top_u', & 746 nest_offl%u_top(0:1,nys:nyn,nxlu:nxr), & 747 nxlu, nys+1, nest_offl%tind+1, & 748 nxr-nxlu+1, nyn-nys+1, 2, .TRUE. ) 749 750 CALL get_variable( pids_id, 'ls_forcing_top_v', & 751 nest_offl%v_top(0:1,nysv:nyn,nxl:nxr), & 752 nxl+1, nysv, nest_offl%tind+1, & 753 nxr-nxl+1, nyn-nysv+1, 2, .TRUE. ) 754 755 CALL get_variable( pids_id, 'ls_forcing_top_w', & 756 nest_offl%w_top(0:1,nys:nyn,nxl:nxr), & 757 nxl+1, nys+1, nest_offl%tind+1, & 758 nxr-nxl+1, nyn-nys+1, 2, .TRUE. ) 759 760 IF ( .NOT. neutral ) THEN 761 CALL get_variable( pids_id, 'ls_forcing_top_pt', & 762 nest_offl%pt_top(0:1,nys:nyn,nxl:nxr), & 763 nxl+1, nys+1, nest_offl%tind+1, & 764 nxr-nxl+1, nyn-nys+1, 2, .TRUE. ) 765 ENDIF 766 IF ( humidity ) THEN 767 CALL get_variable( pids_id, 'ls_forcing_top_qv', & 768 nest_offl%q_top(0:1,nys:nyn,nxl:nxr), & 769 nxl+1, nys+1, nest_offl%tind+1, & 770 nxr-nxl+1, nyn-nys+1, 2, .TRUE. ) 771 ENDIF 772 773 IF ( air_chemistry ) THEN 774 DO n = 1, UBOUND(nest_offl%var_names_chem_t, 1) 775 IF ( check_existence( nest_offl%var_names, & 776 nest_offl%var_names_chem_t(n) ) ) THEN 777 CALL get_variable( pids_id, & 778 TRIM( nest_offl%var_names_chem_t(n) ), & 779 nest_offl%chem_top(0:1,nys:nyn,nxl:nxr,n), & 780 nxl+1, nys+1, nest_offl%tind+1, & 781 nxr-nxl+1, nyn-nys+1, 2, .TRUE. ) 782 nest_offl%chem_from_file_t(n) = .TRUE. 783 ENDIF 784 ENDDO 785 ENDIF 786 787 ! 788 !-- Close input file 789 CALL close_input_file( pids_id ) 790 #endif 791 ! 792 !-- Set control flag to indicate that boundary data has been initially 793 !-- input. 794 nest_offl%init = .TRUE. 795 ! 796 !-- End of CPU measurement 797 CALL cpu_log( log_point_s(86), 'NetCDF input forcing', 'stop' ) 798 799 END SUBROUTINE nesting_offl_input 800 199 801 200 802 !------------------------------------------------------------------------------! … … 210 812 SUBROUTINE nesting_offl_mass_conservation 211 813 212 IMPLICIT NONE213 214 814 INTEGER(iwp) :: i !< grid index in x-direction 215 815 INTEGER(iwp) :: j !< grid index in y-direction … … 281 881 #if defined( __parallel ) 282 882 IF ( collective_wait ) CALL MPI_BARRIER( comm2d, ierr ) 283 CALL MPI_ALLREDUCE( volume_flow_l, volume_flow, 3, MPI_REAL, MPI_SUM, 883 CALL MPI_ALLREDUCE( volume_flow_l, volume_flow, 3, MPI_REAL, MPI_SUM, & 284 884 comm2d, ierr ) 285 885 #else … … 314 914 SUBROUTINE nesting_offl_bc 315 915 316 IMPLICIT NONE317 318 916 INTEGER(iwp) :: i !< running index x-direction 319 917 INTEGER(iwp) :: j !< running index y-direction 320 918 INTEGER(iwp) :: k !< running index z-direction 321 919 INTEGER(iwp) :: n !< running index for chemical species 322 323 REAL(wp) :: fac_dt !< interpolation factor324 920 325 921 REAL(wp), DIMENSION(nzb:nzt+1) :: pt_ref !< reference profile for potential temperature … … 347 943 u_ref_l = 0.0_wp 348 944 v_ref_l = 0.0_wp 349 !350 !-- Determine interpolation factor and limit it to 1. This is because351 !-- t+dt can slightly exceed time(tind_p) before boundary data is updated352 !-- again.353 fac_dt = ( time_since_reference_point - nest_offl%time(nest_offl%tind) &354 + dt_3d ) / &355 ( nest_offl%time(nest_offl%tind_p) - nest_offl%time(nest_offl%tind) )356 fac_dt = MIN( 1.0_wp, fac_dt )357 945 ! 358 946 !-- Set boundary conditions of u-, v-, w-component, as well as q, and pt. … … 863 1451 CALL nesting_offl_calc_zi 864 1452 CALL adjust_sponge_layer 865 1453 1454 CALL cpu_log( log_point(58), 'offline nesting', 'stop' ) 1455 1456 IF ( debug_output_timestep ) CALL debug_message( 'nesting_offl_bc', 'end' ) 1457 1458 1459 END SUBROUTINE nesting_offl_bc 1460 1461 !------------------------------------------------------------------------------! 1462 ! Description: 1463 !------------------------------------------------------------------------------! 1464 !> Update of the geostrophic wind components. 1465 !> @todo: update geostrophic wind also in the child domains (should be done 1466 !> in the nesting. 1467 !------------------------------------------------------------------------------! 1468 SUBROUTINE nesting_offl_geostrophic_wind 1469 1470 INTEGER(iwp) :: k 866 1471 ! 867 1472 !-- Update geostrophic wind components from dynamic input file. … … 874 1479 ug(nzt+1) = ug(nzt) 875 1480 vg(nzt+1) = vg(nzt) 876 877 CALL cpu_log( log_point(58), 'offline nesting', 'stop' ) 878 879 IF ( debug_output_timestep ) CALL debug_message( 'nesting_offl_bc', 'end' ) 880 881 882 END SUBROUTINE nesting_offl_bc 1481 1482 END SUBROUTINE nesting_offl_geostrophic_wind 1483 1484 !------------------------------------------------------------------------------! 1485 ! Description: 1486 !------------------------------------------------------------------------------! 1487 !> Determine the interpolation constant for time interpolation. The 1488 !> calculation is separated from the nesting_offl_bc and 1489 !> nesting_offl_geostrophic_wind in order to be independent on the order 1490 !> of calls. 1491 !------------------------------------------------------------------------------! 1492 SUBROUTINE nesting_offl_interpolation_factor 1493 ! 1494 !-- Determine interpolation factor and limit it to 1. This is because 1495 !-- t+dt can slightly exceed time(tind_p) before boundary data is updated 1496 !-- again. 1497 fac_dt = ( time_utc_init + time_since_reference_point & 1498 - nest_offl%time(nest_offl%tind) + dt_3d ) / & 1499 ( nest_offl%time(nest_offl%tind_p) - nest_offl%time(nest_offl%tind) ) 1500 1501 fac_dt = MIN( 1.0_wp, fac_dt ) 1502 1503 END SUBROUTINE nesting_offl_interpolation_factor 883 1504 884 1505 !------------------------------------------------------------------------------! … … 889 1510 !------------------------------------------------------------------------------! 890 1511 SUBROUTINE nesting_offl_calc_zi 891 892 USE basic_constants_and_equations_mod, &893 ONLY: g894 895 USE kinds896 897 IMPLICIT NONE898 1512 899 1513 INTEGER(iwp) :: i !< loop index in x-direction … … 1070 1684 !------------------------------------------------------------------------------! 1071 1685 SUBROUTINE adjust_sponge_layer 1072 1073 USE arrays_3d, &1074 ONLY: rdf, rdf_sc, zu1075 1076 USE basic_constants_and_equations_mod, &1077 ONLY: pi1078 1079 USE kinds1080 1081 IMPLICIT NONE1082 1686 1083 1687 INTEGER(iwp) :: k !< loop index in z-direction … … 1116 1720 !------------------------------------------------------------------------------! 1117 1721 SUBROUTINE nesting_offl_check_parameters 1118 1119 IMPLICIT NONE1120 1722 ! 1121 1723 !-- Check if offline nesting is applied in nested child domain. … … 1123 1725 message_string = 'Offline nesting is only applicable in root model.' 1124 1726 CALL message( 'offline_nesting_check_parameters', 'PA0622', 1, 2, 0, 6, 0 ) 1125 ENDIF 1727 ENDIF 1126 1728 1127 1729 END SUBROUTINE nesting_offl_check_parameters … … 1133 1735 !------------------------------------------------------------------------------! 1134 1736 SUBROUTINE nesting_offl_parin 1135 1136 IMPLICIT NONE1137 1737 1138 1738 CHARACTER (LEN=80) :: line !< dummy string that contains the current line of the parameter file … … 1174 1774 SUBROUTINE nesting_offl_header ( io ) 1175 1775 1176 IMPLICIT NONE1177 1178 1776 INTEGER(iwp), INTENT(IN) :: io !< Unit of the output file 1179 1777 … … 1199 1797 !------------------------------------------------------------------------------! 1200 1798 SUBROUTINE nesting_offl_init 1201 1202 USE netcdf_data_input_mod, & 1203 ONLY: netcdf_data_input_offline_nesting 1204 1205 IMPLICIT NONE 1206 1799 1207 1800 INTEGER(iwp) :: n !< running index for chemical species 1208 1801 … … 1340 1933 ENDIF 1341 1934 ! 1935 !-- Before initial data input is initiated, check if dynamic input file is 1936 !-- present. 1937 IF ( .NOT. input_pids_dynamic ) THEN 1938 message_string = 'nesting_offline = .TRUE. requires dynamic ' // & 1939 'input file ' // & 1940 TRIM( input_file_dynamic ) // TRIM( coupling_char ) 1941 CALL message( 'nesting_offl_init', 'PA0546', 1, 2, 0, 6, 0 ) 1942 ENDIF 1943 ! 1342 1944 !-- Read COSMO data at lateral and top boundaries 1343 CALL ne tcdf_data_input_offline_nesting1945 CALL nesting_offl_input 1344 1946 ! 1345 1947 !-- Check if sufficient time steps are provided to cover the entire … … 1347 1949 !-- not for the soil/wall spinup. However, as the spinup time is added 1348 1950 !-- to the end_time, this must be considered here. 1349 IF ( end_time - spinup_time > nest_offl%time(nest_offl%nt-1) ) THEN 1350 message_string = 'end_time > provided time in offline nesting.' 1351 CALL message( 'offline_nesting_check_parameters', 'PA0183', & 1352 1, 2, 0, 6, 0 ) 1951 IF ( end_time - spinup_time > & 1952 nest_offl%time(nest_offl%nt-1) - time_utc_init ) THEN 1953 message_string = 'end_time of the simulation exceeds the ' // & 1954 'time dimension in the dynamic input file.' 1955 CALL message( 'nesting_offl_init', 'PA0183', 1, 2, 0, 6, 0 ) 1956 ENDIF 1957 1958 IF ( nest_offl%time(0) /= time_utc_init ) THEN 1959 message_string = 'Offline nesting: time dimension must start at ' // & 1960 ' time_utc_init.' 1961 CALL message( 'nesting_offl_init', 'PA0676', 1, 2, 0, 6, 0 ) 1353 1962 ENDIF 1354 1963 ! … … 1466 2075 !------------------------------------------------------------------------------! 1467 2076 FUNCTION interpolate_in_time( var_t1, var_t2, fac ) 1468 1469 USE kinds1470 1471 IMPLICIT NONE1472 2077 1473 2078 REAL(wp) :: interpolate_in_time !< time-interpolated boundary value
Note: See TracChangeset
for help on using the changeset viewer.
