Changeset 147 for palm/trunk/SOURCE/read_3d_binary.f90
- Timestamp:
- Feb 1, 2008 12:41:46 PM (16 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/read_3d_binary.f90 (modified) (49 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/read_3d_binary.f90
r146 r147 5 5 ! ----------------- 6 6 ! Files from which restart data are to be read are determined and subsequently 7 ! opened, 7 ! opened. The total domain on the restart file is allowed to be smaller than 8 ! the current total domain. In this case it will be periodically mapped on the 9 ! current domain (needed for recycling method). 8 10 ! +call of user_read_restart_data, -dopr_time_count, 9 11 ! hom_sum, volume_flow_area, volume_flow_initial moved to read_var_list, … … 62 64 63 65 INTEGER :: files_to_be_opened, i, idum1, j, myid_on_file, & 64 numprocs_on_file, nxlc, nxl_on_file, nxrc, nxr_on_file, nync, & 65 nyn_on_file, nysc, nys_on_file, nzb_on_file, nzt_on_file 66 67 INTEGER, DIMENSION(:), ALLOCATABLE :: file_list 66 numprocs_on_file, nxlc, nxlf, nxlpr, nxl_on_file, nxrc, nxrf, & 67 nxrpr, nxr_on_file, nync, nynf, nynpr, nyn_on_file, nysc, & 68 nysf, nyspr, nys_on_file, nzb_on_file, nzt_on_file, offset_x, & 69 offset_y 70 71 INTEGER, DIMENSION(numprocs_previous_run*4) :: file_list, nxlfa, nxrfa, & 72 nynfa, nysfa, offset_xa, offset_ya 68 73 69 74 REAL, DIMENSION(:,:), ALLOCATABLE :: tmp_2d … … 80 85 !-- of this PE 81 86 files_to_be_opened = 0 82 ALLOCATE( file_list(numprocs_previous_run) )83 87 84 88 DO i = 1, numprocs_previous_run 85 89 90 nxlpr = hor_index_bounds_previous_run(1,i-1) 91 nxrpr = hor_index_bounds_previous_run(2,i-1) 92 nyspr = hor_index_bounds_previous_run(3,i-1) 93 nynpr = hor_index_bounds_previous_run(4,i-1) 94 95 ! 96 !-- Determine the offsets. They may be non-zero in case that the total domain 97 !-- on file is smaller than the current total domain. 98 offset_x = ( nxl / ( nx_on_file + 1 ) ) * ( nx_on_file + 1 ) 99 offset_y = ( nys / ( ny_on_file + 1 ) ) * ( ny_on_file + 1 ) 100 86 101 ! 87 102 !-- Only data which overlap with the current subdomain have to be read 88 IF ( hor_index_bounds_previous_run(1,i-1) <= nxr .AND. & 89 hor_index_bounds_previous_run(2,i-1) >= nxl .AND. & 90 hor_index_bounds_previous_run(3,i-1) <= nyn .AND. & 91 hor_index_bounds_previous_run(4,i-1) >= nys ) THEN 103 IF ( nxlpr <= nxr-offset_x .AND. nxrpr >= nxl-offset_x .AND. & 104 nyspr <= nyn-offset_y .AND. nynpr >= nys-offset_y ) THEN 92 105 93 106 files_to_be_opened = files_to_be_opened + 1 94 107 file_list(files_to_be_opened) = i-1 108 ! 109 !-- Index bounds of overlapping data 110 nxlfa(files_to_be_opened) = MAX( nxl-offset_x, nxlpr ) 111 nxrfa(files_to_be_opened) = MIN( nxr-offset_x, nxrpr ) 112 nysfa(files_to_be_opened) = MAX( nys-offset_y, nyspr ) 113 nynfa(files_to_be_opened) = MIN( nyn-offset_y, nynpr ) 114 115 WRITE (9,*) '*** reading from file: ', i 116 WRITE (9,*) ' index bounds on file:' 117 WRITE (9,*) ' nxlpr=', nxlpr, ' nxrpr=', nxrpr 118 WRITE (9,*) ' nyspr=', nyspr, ' nynpr=', nynpr 119 WRITE (9,*) ' index bounds of current subdmain:' 120 WRITE (9,*) ' nxl =', nxl, ' nxr =', nxr 121 WRITE (9,*) ' nys =', nys, ' nyn =', nyn 122 WRITE (9,*) ' offset used:' 123 WRITE (9,*) ' offset_x=', offset_x, ' offset_y=', offset_y 124 WRITE (9,*) ' bounds of overlapping data:' 125 WRITE (9,*) ' nxlfa=', nxlfa(files_to_be_opened), ' nxrfa=', nxrfa(files_to_be_opened) 126 WRITE (9,*) ' nysfa=', nysfa(files_to_be_opened), ' nynfa=', nynfa(files_to_be_opened) 127 CALL local_flush( 9 ) 128 offset_xa(files_to_be_opened) = offset_x 129 offset_ya(files_to_be_opened) = offset_y 95 130 96 131 ENDIF 97 132 133 ! 134 !-- If the total domain on file is smaller than the current total domain, 135 !-- and if the current subdomain extends beyond the limits of the total 136 !-- domain of file, the respective file may be opened again (three times 137 !-- maximum) to read the still missing parts, which are then added 138 !-- "cyclically". 139 !-- Overlap along x: 140 IF ( ( nxr - offset_x ) > nx_on_file ) THEN 141 142 offset_x = offset_x + ( nx_on_file + 1 ) 143 144 IF ( nxlpr <= nxr-offset_x .AND. nxrpr >= nxl-offset_x .AND. & 145 nyspr <= nyn-offset_y .AND. nynpr >= nys-offset_y ) THEN 146 147 files_to_be_opened = files_to_be_opened + 1 148 file_list(files_to_be_opened) = i-1 149 ! 150 !-- Index bounds of overlapping data 151 nxlfa(files_to_be_opened) = MAX( nxl-offset_x, nxlpr ) 152 nxrfa(files_to_be_opened) = MIN( nxr-offset_x, nxrpr ) 153 nysfa(files_to_be_opened) = MAX( nys-offset_y, nyspr ) 154 nynfa(files_to_be_opened) = MIN( nyn-offset_y, nynpr ) 155 156 WRITE (9,*) '*** reading from file: ', i 157 WRITE (9,*) ' index bounds on file:' 158 WRITE (9,*) ' nxlpr=', nxlpr, ' nxrpr=', nxrpr 159 WRITE (9,*) ' nyspr=', nyspr, ' nynpr=', nynpr 160 WRITE (9,*) ' index bounds of current subdmain:' 161 WRITE (9,*) ' nxl =', nxl, ' nxr =', nxr 162 WRITE (9,*) ' nys =', nys, ' nyn =', nyn 163 WRITE (9,*) ' offset used:' 164 WRITE (9,*) ' offset_x=', offset_x, ' offset_y=', offset_y 165 WRITE (9,*) ' bounds of overlapping data:' 166 WRITE (9,*) ' nxlfa=', nxlfa(files_to_be_opened), ' nxrfa=', nxrfa(files_to_be_opened) 167 WRITE (9,*) ' nysfa=', nysfa(files_to_be_opened), ' nynfa=', nynfa(files_to_be_opened) 168 CALL local_flush( 9 ) 169 offset_xa(files_to_be_opened) = offset_x 170 offset_ya(files_to_be_opened) = offset_y 171 172 ENDIF 173 174 offset_x = offset_x - ( nx_on_file + 1 ) 175 176 ENDIF 177 178 179 ! 180 !-- Overlap along y: 181 IF ( ( nyn - offset_y ) > ny_on_file ) THEN 182 183 offset_y = offset_y + ( ny_on_file + 1 ) 184 185 IF ( nxlpr <= nxr-offset_x .AND. nxrpr >= nxl-offset_x .AND. & 186 nyspr <= nyn-offset_y .AND. nynpr >= nys-offset_y ) THEN 187 188 files_to_be_opened = files_to_be_opened + 1 189 file_list(files_to_be_opened) = i-1 190 ! 191 !-- Index bounds of overlapping data 192 nxlfa(files_to_be_opened) = MAX( nxl-offset_x, nxlpr ) 193 nxrfa(files_to_be_opened) = MIN( nxr-offset_x, nxrpr ) 194 nysfa(files_to_be_opened) = MAX( nys-offset_y, nyspr ) 195 nynfa(files_to_be_opened) = MIN( nyn-offset_y, nynpr ) 196 197 WRITE (9,*) '*** reading from file: ', i 198 WRITE (9,*) ' index bounds on file:' 199 WRITE (9,*) ' nxlpr=', nxlpr, ' nxrpr=', nxrpr 200 WRITE (9,*) ' nyspr=', nyspr, ' nynpr=', nynpr 201 WRITE (9,*) ' index bounds of current subdmain:' 202 WRITE (9,*) ' nxl =', nxl, ' nxr =', nxr 203 WRITE (9,*) ' nys =', nys, ' nyn =', nyn 204 WRITE (9,*) ' offset used:' 205 WRITE (9,*) ' offset_x=', offset_x, ' offset_y=', offset_y 206 WRITE (9,*) ' bounds of overlapping data:' 207 WRITE (9,*) ' nxlfa=', nxlfa(files_to_be_opened), ' nxrfa=', nxrfa(files_to_be_opened) 208 WRITE (9,*) ' nysfa=', nysfa(files_to_be_opened), ' nynfa=', nynfa(files_to_be_opened) 209 CALL local_flush( 9 ) 210 offset_xa(files_to_be_opened) = offset_x 211 offset_ya(files_to_be_opened) = offset_y 212 213 ENDIF 214 215 offset_y = offset_y - ( ny_on_file + 1 ) 216 217 ENDIF 218 219 !-- Overlap along x and y: 220 IF ( ( nxr - offset_x ) > nx_on_file .AND. & 221 ( nyn - offset_y ) > ny_on_file ) THEN 222 223 offset_x = offset_x + ( nx_on_file + 1 ) 224 offset_y = offset_y + ( ny_on_file + 1 ) 225 226 IF ( nxlpr <= nxr-offset_x .AND. nxrpr >= nxl-offset_x .AND. & 227 nyspr <= nyn-offset_y .AND. nynpr >= nys-offset_y ) THEN 228 229 files_to_be_opened = files_to_be_opened + 1 230 file_list(files_to_be_opened) = i-1 231 ! 232 !-- Index bounds of overlapping data 233 nxlfa(files_to_be_opened) = MAX( nxl-offset_x, nxlpr ) 234 nxrfa(files_to_be_opened) = MIN( nxr-offset_x, nxrpr ) 235 nysfa(files_to_be_opened) = MAX( nys-offset_y, nyspr ) 236 nynfa(files_to_be_opened) = MIN( nyn-offset_y, nynpr ) 237 238 WRITE (9,*) '*** reading from file: ', i 239 WRITE (9,*) ' index bounds on file:' 240 WRITE (9,*) ' nxlpr=', nxlpr, ' nxrpr=', nxrpr 241 WRITE (9,*) ' nyspr=', nyspr, ' nynpr=', nynpr 242 WRITE (9,*) ' index bounds of current subdmain:' 243 WRITE (9,*) ' nxl =', nxl, ' nxr =', nxr 244 WRITE (9,*) ' nys =', nys, ' nyn =', nyn 245 WRITE (9,*) ' offset used:' 246 WRITE (9,*) ' offset_x=', offset_x, ' offset_y=', offset_y 247 WRITE (9,*) ' bounds of overlapping data:' 248 WRITE (9,*) ' nxlfa=', nxlfa(files_to_be_opened), ' nxrfa=', nxrfa(files_to_be_opened) 249 WRITE (9,*) ' nysfa=', nysfa(files_to_be_opened), ' nynfa=', nynfa(files_to_be_opened) 250 CALL local_flush( 9 ) 251 offset_xa(files_to_be_opened) = offset_x 252 offset_ya(files_to_be_opened) = offset_y 253 254 ENDIF 255 256 offset_x = offset_x - ( nx_on_file + 1 ) 257 offset_y = offset_y - ( ny_on_file + 1 ) 258 259 ENDIF 260 98 261 ENDDO 99 262 263 ! 264 !-- Save the id-string of the current process, since myid_char may now be used 265 !-- to open files created by PEs with other id. 100 266 myid_char_save = myid_char 101 267 268 ! 269 !-- Test output (remove later) 102 270 DO i = 1, numprocs_previous_run 103 271 WRITE (9,*) 'i=',i-1, ' ibs= ',hor_index_bounds_previous_run(1:4,i-1) … … 128 296 !-- first. 129 297 CALL check_open( 13 ) 298 WRITE (9,*) 'before skipping' 299 CALL local_flush( 9 ) 130 300 IF ( j == 0 ) CALL skip_var_list 301 WRITE (9,*) 'skipping done' 302 CALL local_flush( 9 ) 131 303 132 304 ! … … 228 400 229 401 ! 230 !-- Determine the index range of those gridpoints to be copied from 231 !-- the subdomains on the restart files to the current subdomain 232 nxlc = MAX( nxl, nxl_on_file ) 233 nxrc = MIN( nxr, nxr_on_file ) 234 nysc = MAX( nys, nys_on_file ) 235 nync = MIN( nyn, nyn_on_file ) 402 !-- Get the index range of the subdomain on file which overlap with the 403 !-- current subdomain 404 nxlf = nxlfa(i) 405 nxlc = nxlfa(i) + offset_xa(i) 406 nxrf = nxrfa(i) 407 nxrc = nxrfa(i) + offset_xa(i) 408 nysf = nysfa(i) 409 nysc = nysfa(i) + offset_ya(i) 410 nynf = nynfa(i) 411 nync = nynfa(i) + offset_ya(i) 236 412 237 413 ! … … 244 420 DO WHILE ( TRIM( field_chr ) /= '*** end ***' ) 245 421 422 WRITE (9,*) 'var = ', field_chr 423 CALL local_flush( 9 ) 246 424 SELECT CASE ( TRIM( field_chr ) ) 247 425 … … 249 427 READ ( 13 ) tmp_3d 250 428 e(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 251 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)429 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 252 430 253 431 CASE ( 'e_av' ) … … 257 435 READ ( 13 ) tmp_3d 258 436 e_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 259 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)437 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 260 438 261 439 CASE ( 'e_m' ) 262 440 READ ( 13 ) tmp_3d 263 441 e_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 264 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)442 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 265 443 266 444 CASE ( 'iran' ) ! matching random numbers is still unresolved issue … … 270 448 READ ( 13 ) tmp_3d 271 449 kh(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 272 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)450 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 273 451 274 452 CASE ( 'kh_m' ) 275 453 READ ( 13 ) tmp_3d 276 454 kh_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 277 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)455 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 278 456 279 457 CASE ( 'km' ) 280 458 READ ( 13 ) tmp_3d 281 459 km(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 282 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)460 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 283 461 284 462 CASE ( 'km_m' ) 285 463 READ ( 13 ) tmp_3d 286 464 km_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 287 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)465 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 288 466 289 467 CASE ( 'lwp_av' ) … … 293 471 READ ( 13 ) tmp_2d 294 472 lwp_av(nysc-1:nync+1,nxlc-1:nxrc+1) = & 295 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)473 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 296 474 297 475 CASE ( 'p' ) 298 476 READ ( 13 ) tmp_3d 299 477 p(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 300 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)478 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 301 479 302 480 CASE ( 'p_av' ) … … 306 484 READ ( 13 ) tmp_3d 307 485 p_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 308 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)486 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 309 487 310 488 CASE ( 'pc_av' ) … … 314 492 READ ( 13 ) tmp_3d 315 493 pc_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 316 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)494 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 317 495 318 496 CASE ( 'pr_av' ) … … 322 500 READ ( 13 ) tmp_3d 323 501 pr_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 324 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)502 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 325 503 326 504 CASE ( 'precipitation_amount' ) 327 505 READ ( 13 ) tmp_2d 328 506 precipitation_amount(nysc-1:nync+1,nxlc-1:nxrc+1) = & 329 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)507 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 330 508 331 509 CASE ( 'precipitation_rate_a' ) … … 335 513 READ ( 13 ) tmp_2d 336 514 precipitation_rate_av(nysc-1:nync+1,nxlc-1:nxrc+1) = & 337 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)515 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 338 516 339 517 CASE ( 'pt' ) 340 518 READ ( 13 ) tmp_3d 341 519 pt(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 342 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)520 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 343 521 344 522 CASE ( 'pt_av' ) … … 348 526 READ ( 13 ) tmp_3d 349 527 pt_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 350 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)528 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 351 529 352 530 CASE ( 'pt_m' ) 353 531 READ ( 13 ) tmp_3d 354 532 pt_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 355 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)533 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 356 534 357 535 CASE ( 'q' ) 358 536 READ ( 13 ) tmp_3d 359 537 q(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 360 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)538 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 361 539 362 540 CASE ( 'q_av' ) … … 366 544 READ ( 13 ) tmp_3d 367 545 q_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 368 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)546 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 369 547 370 548 CASE ( 'q_m' ) 371 549 READ ( 13 ) tmp_3d 372 550 q_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 373 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)551 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 374 552 375 553 CASE ( 'ql' ) 376 554 READ ( 13 ) tmp_3d 377 555 ql(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 378 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)556 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 379 557 380 558 CASE ( 'ql_av' ) … … 384 562 READ ( 13 ) tmp_3d 385 563 ql_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 386 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)564 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 387 565 388 566 CASE ( 'ql_c_av' ) … … 392 570 READ ( 13 ) tmp_3d 393 571 ql_c_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 394 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)572 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 395 573 396 574 CASE ( 'ql_v_av' ) … … 400 578 READ ( 13 ) tmp_3d 401 579 ql_v_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 402 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)580 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 403 581 404 582 CASE ( 'ql_vp_av' ) … … 408 586 READ ( 13 ) tmp_3d 409 587 ql_vp_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 410 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)588 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 411 589 412 590 CASE ( 'qs' ) 413 591 READ ( 13 ) tmp_2d 414 592 qs(nysc-1:nync+1,nxlc-1:nxrc+1) = & 415 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)593 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 416 594 417 595 CASE ( 'qsws' ) 418 596 READ ( 13 ) tmp_2d 419 597 qsws(nysc-1:nync+1,nxlc-1:nxrc+1) = & 420 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)598 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 421 599 422 600 CASE ( 'qsws_m' ) 423 601 READ ( 13 ) tmp_2d 424 602 qsws_m(nysc-1:nync+1,nxlc-1:nxrc+1) = & 425 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)603 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 426 604 427 605 CASE ( 'qswst' ) 428 606 READ ( 13 ) tmp_2d 429 607 qswst(nysc-1:nync+1,nxlc-1:nxrc+1) = & 430 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)608 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 431 609 432 610 CASE ( 'qswst_m' ) 433 611 READ ( 13 ) tmp_2d 434 612 qswst_m(nysc-1:nync+1,nxlc-1:nxrc+1) = & 435 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)613 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 436 614 437 615 CASE ( 'qv_av' ) … … 441 619 READ ( 13 ) tmp_3d 442 620 qv_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 443 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)621 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 444 622 445 623 CASE ( 'random_iv' ) ! still unresolved issue … … 453 631 READ ( 13 ) tmp_3d 454 632 rho_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 455 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)633 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 456 634 457 635 CASE ( 'rif' ) 458 636 READ ( 13 ) tmp_2d 459 637 rif(nysc-1:nync+1,nxlc-1:nxrc+1) = & 460 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)638 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 461 639 462 640 CASE ( 'rif_m' ) 463 641 READ ( 13 ) tmp_2d 464 642 rif_m(nysc-1:nync+1,nxlc-1:nxrc+1) = & 465 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)643 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 466 644 467 645 CASE ( 'rif_wall' ) … … 470 648 READ ( 13 ) tmp_4d 471 649 rif_wall(:,nysc-1:nync+1,nxlc-1:nxrc+1,:) = & 472 tmp_4d(:,nys c-1:nync+1,nxlc-1:nxrc+1,:)650 tmp_4d(:,nysf-1:nynf+1,nxlf-1:nxrf+1,:) 473 651 DEALLOCATE( tmp_4d ) 474 652 … … 479 657 READ ( 13 ) tmp_3d 480 658 s_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 481 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)659 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 482 660 483 661 CASE ( 'sa' ) 484 662 READ ( 13 ) tmp_3d 485 663 sa(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 486 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)664 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 487 665 488 666 CASE ( 'sa_av' ) … … 492 670 READ ( 13 ) tmp_3d 493 671 sa_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 494 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)672 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 495 673 496 674 CASE ( 'saswsb' ) 497 675 READ ( 13 ) tmp_2d 498 676 saswsb(nysc-1:nync+1,nxlc-1:nxrc+1) = & 499 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)677 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 500 678 501 679 CASE ( 'saswst' ) 502 680 READ ( 13 ) tmp_2d 503 681 saswst(nysc-1:nync+1,nxlc-1:nxrc+1) = & 504 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)682 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 505 683 506 684 CASE ( 'shf' ) 507 685 READ ( 13 ) tmp_2d 508 686 shf(nysc-1:nync+1,nxlc-1:nxrc+1) = & 509 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)687 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 510 688 511 689 CASE ( 'shf_m' ) 512 690 READ ( 13 ) tmp_2d 513 691 shf_m(nysc-1:nync+1,nxlc-1:nxrc+1) = & 514 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)692 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 515 693 516 694 CASE ( 'spectrum_x' ) … … 523 701 READ ( 13 ) tmp_2d 524 702 ts(nysc-1:nync+1,nxlc-1:nxrc+1) = & 525 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)703 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 526 704 527 705 CASE ( 'ts_av' ) … … 531 709 READ ( 13 ) tmp_2d 532 710 ts_av(nysc-1:nync+1,nxlc-1:nxrc+1) = & 533 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)711 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 534 712 535 713 CASE ( 'tswst' ) 536 714 READ ( 13 ) tmp_2d 537 715 tswst(nysc-1:nync+1,nxlc-1:nxrc+1) = & 538 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)716 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 539 717 540 718 CASE ( 'tswst_m' ) 541 719 READ ( 13 ) tmp_2d 542 720 tswst_m(nysc-1:nync+1,nxlc-1:nxrc+1) = & 543 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)721 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 544 722 545 723 CASE ( 'u' ) 546 724 READ ( 13 ) tmp_3d 547 725 u(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 548 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)726 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 549 727 550 728 CASE ( 'u_av' ) … … 554 732 READ ( 13 ) tmp_3d 555 733 u_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 556 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)734 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 557 735 558 736 CASE ( 'u_m' ) 559 737 READ ( 13 ) tmp_3d 560 738 u_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 561 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)739 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 562 740 563 741 CASE ( 'u_m_l' ) … … 565 743 READ ( 13 ) tmp_3dw 566 744 IF ( outflow_l ) THEN 567 u_m_l(:,nysc-1:nync+1,:) = tmp_3dw(:,nys c-1:nync+1,:)745 u_m_l(:,nysc-1:nync+1,:) = tmp_3dw(:,nysf-1:nynf+1,:) 568 746 ENDIF 569 747 DEALLOCATE( tmp_3dw ) … … 574 752 READ ( 13 ) tmp_3dw 575 753 IF ( outflow_n ) THEN 576 u_m_n(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxl c-1:nxrc+1)754 u_m_n(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxlf-1:nxrf+1) 577 755 ENDIF 578 756 DEALLOCATE( tmp_3dw ) … … 583 761 READ ( 13 ) tmp_3dw 584 762 IF ( outflow_r ) THEN 585 u_m_r(:,nysc-1:nync+1,:) = tmp_3dw(:,nys c-1:nync+1,:)763 u_m_r(:,nysc-1:nync+1,:) = tmp_3dw(:,nysf-1:nynf+1,:) 586 764 ENDIF 587 765 DEALLOCATE( tmp_3dw ) … … 592 770 READ ( 13 ) tmp_3dw 593 771 IF ( outflow_s ) THEN 594 u_m_s(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxl c-1:nxrc+1)772 u_m_s(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxlf-1:nxrf+1) 595 773 ENDIF 596 774 DEALLOCATE( tmp_3dw ) … … 599 777 READ ( 13 ) tmp_2d 600 778 us(nysc-1:nync+1,nxlc-1:nxrc+1) = & 601 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)779 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 602 780 603 781 CASE ( 'usws' ) 604 782 READ ( 13 ) tmp_2d 605 783 usws(nysc-1:nync+1,nxlc-1:nxrc+1) = & 606 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)784 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 607 785 608 786 CASE ( 'uswst' ) 609 787 READ ( 13 ) tmp_2d 610 788 uswst(nysc-1:nync+1,nxlc-1:nxrc+1) = & 611 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)789 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 612 790 613 791 CASE ( 'usws_m' ) 614 792 READ ( 13 ) tmp_2d 615 793 usws_m(nysc-1:nync+1,nxlc-1:nxrc+1) = & 616 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)794 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 617 795 618 796 CASE ( 'uswst_m' ) 619 797 READ ( 13 ) tmp_2d 620 798 uswst_m(nysc-1:nync+1,nxlc-1:nxrc+1) = & 621 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)799 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 622 800 623 801 CASE ( 'us_av' ) … … 627 805 READ ( 13 ) tmp_2d 628 806 us_av(nysc-1:nync+1,nxlc-1:nxrc+1) = & 629 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)807 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 630 808 631 809 CASE ( 'v' ) 632 810 READ ( 13 ) tmp_3d 633 811 v(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 634 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)812 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 635 813 636 814 CASE ( 'v_av' ) … … 640 818 READ ( 13 ) tmp_3d 641 819 v_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 642 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)820 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 643 821 644 822 CASE ( 'v_m' ) 645 823 READ ( 13 ) tmp_3d 646 824 v_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 647 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)825 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 648 826 649 827 CASE ( 'v_m_l' ) … … 651 829 READ ( 13 ) tmp_3dw 652 830 IF ( outflow_l ) THEN 653 v_m_l(:,nysc-1:nync+1,:) = tmp_3dw(:,nys c-1:nync+1,:)831 v_m_l(:,nysc-1:nync+1,:) = tmp_3dw(:,nysf-1:nynf+1,:) 654 832 ENDIF 655 833 DEALLOCATE( tmp_3dw ) … … 660 838 READ ( 13 ) tmp_3dw 661 839 IF ( outflow_n ) THEN 662 v_m_n(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxl c-1:nxrc+1)840 v_m_n(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxlf-1:nxrf+1) 663 841 ENDIF 664 842 DEALLOCATE( tmp_3dw ) … … 669 847 READ ( 13 ) tmp_3dw 670 848 IF ( outflow_r ) THEN 671 v_m_r(:,nysc-1:nync+1,:) = tmp_3dw(:,nys c-1:nync+1,:)849 v_m_r(:,nysc-1:nync+1,:) = tmp_3dw(:,nysf-1:nynf+1,:) 672 850 ENDIF 673 851 DEALLOCATE( tmp_3dw ) … … 678 856 READ ( 13 ) tmp_3dw 679 857 IF ( outflow_s ) THEN 680 v_m_s(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxl c-1:nxrc+1)858 v_m_s(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxlf-1:nxrf+1) 681 859 ENDIF 682 860 DEALLOCATE( tmp_3dw ) … … 685 863 READ ( 13 ) tmp_3d 686 864 vpt(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 687 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)865 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 688 866 689 867 CASE ( 'vpt_av' ) … … 693 871 READ ( 13 ) tmp_3d 694 872 vpt_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 695 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)873 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 696 874 697 875 CASE ( 'vpt_m' ) 698 876 READ ( 13 ) tmp_3d 699 877 vpt_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 700 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)878 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 701 879 702 880 CASE ( 'vsws' ) 703 881 READ ( 13 ) tmp_2d 704 882 vsws(nysc-1:nync+1,nxlc-1:nxrc+1) = & 705 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)883 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 706 884 707 885 CASE ( 'vswst' ) 708 886 READ ( 13 ) tmp_2d 709 887 vswst(nysc-1:nync+1,nxlc-1:nxrc+1) = & 710 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)888 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 711 889 712 890 CASE ( 'vsws_m' ) 713 891 READ ( 13 ) tmp_2d 714 892 vsws_m(nysc-1:nync+1,nxlc-1:nxrc+1) = & 715 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)893 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 716 894 717 895 CASE ( 'vswst_m' ) 718 896 READ ( 13 ) tmp_2d 719 897 vswst_m(nysc-1:nync+1,nxlc-1:nxrc+1) = & 720 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)898 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 721 899 722 900 CASE ( 'w' ) 723 901 READ ( 13 ) tmp_3d 724 902 w(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 725 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)903 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 726 904 727 905 CASE ( 'w_av' ) … … 731 909 READ ( 13 ) tmp_3d 732 910 w_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 733 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)911 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 734 912 735 913 CASE ( 'w_m' ) 736 914 READ ( 13 ) tmp_3d 737 915 w_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & 738 tmp_3d(:,nys c-1:nync+1,nxlc-1:nxrc+1)916 tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) 739 917 740 918 CASE ( 'w_m_l' ) … … 742 920 READ ( 13 ) tmp_3dw 743 921 IF ( outflow_l ) THEN 744 w_m_l(:,nysc-1:nync+1,:) = tmp_3dw(:,nys c-1:nync+1,:)922 w_m_l(:,nysc-1:nync+1,:) = tmp_3dw(:,nysf-1:nynf+1,:) 745 923 ENDIF 746 924 DEALLOCATE( tmp_3dw ) … … 751 929 READ ( 13 ) tmp_3dw 752 930 IF ( outflow_n ) THEN 753 w_m_n(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxl c-1:nxrc+1)931 w_m_n(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxlf-1:nxrf+1) 754 932 ENDIF 755 933 DEALLOCATE( tmp_3dw ) … … 760 938 READ ( 13 ) tmp_3dw 761 939 IF ( outflow_r ) THEN 762 w_m_r(:,nysc-1:nync+1,:) = tmp_3dw(:,nys c-1:nync+1,:)940 w_m_r(:,nysc-1:nync+1,:) = tmp_3dw(:,nysf-1:nynf+1,:) 763 941 ENDIF 764 942 DEALLOCATE( tmp_3dw ) … … 769 947 READ ( 13 ) tmp_3dw 770 948 IF ( outflow_s ) THEN 771 w_m_s(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxl c-1:nxrc+1)949 w_m_s(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxlf-1:nxrf+1) 772 950 ENDIF 773 951 DEALLOCATE( tmp_3dw ) … … 776 954 READ ( 13 ) tmp_2d 777 955 z0(nysc-1:nync+1,nxlc-1:nxrc+1) = & 778 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)956 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 779 957 780 958 CASE ( 'z0_av' ) … … 784 962 READ ( 13 ) tmp_2d 785 963 z0_av(nysc-1:nync+1,nxlc-1:nxrc+1) = & 786 tmp_2d(nys c-1:nync+1,nxlc-1:nxrc+1)964 tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1) 787 965 788 966 CASE DEFAULT … … 801 979 ! 802 980 !-- Read user-defined restart data 803 CALL user_read_restart_data 981 CALL user_read_restart_data( nxlc, nxlf, nxl_on_file, nxrc, nxrf, & 982 nxr_on_file, nync, nynf, nyn_on_file, & 983 nysc, nysf, nys_on_file, tmp_2d, tmp_3d ) 804 984 805 985 !
Note: See TracChangeset
for help on using the changeset viewer.
