Changeset 2512 for palm/trunk/UTIL/combine_plot_fields.f90
- Timestamp:
- Oct 4, 2017 8:26:59 AM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/UTIL/combine_plot_fields.f90
r2365 r2512 25 25 ! ----------------- 26 26 ! $Id$ 27 ! PALM output does not contain ghost layer data any more 28 ! avs- and iso2d-related parts removed, handling of compressed data removed 29 ! 30 ! 2365 2017-08-21 14:59:59Z kanani 27 31 ! Vertical grid nesting implemented (SadiqHuq) 28 32 ! … … 121 125 INTEGER(iwp), DIMENSION(0:1,1000) :: id_var, levels 122 126 123 LOGICAL :: avs_output, compressed, found, iso2d_output, netcdf_output, & 124 netcdf_parallel, netcdf_0, netcdf_1 127 LOGICAL :: found, netcdf_output, netcdf_parallel, netcdf_0, netcdf_1 125 128 LOGICAL :: vnest 126 129 127 130 REAL(wp) :: cpu_start_time, cpu_end_time, dx, simulated_time 128 REAL(wp), DIMENSION(:), ALLOCATABLE :: eta, ho, hu129 131 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: pf, pf_tmp 130 132 REAL(spk), DIMENSION(:,:,:), ALLOCATABLE :: pf3d, pf3d_tmp … … 227 229 228 230 ! 229 !-- Inquire whether an iso2d parameter file exists230 INQUIRE( FILE='PLOT2D_'//modus//'_GLOBAL'//TRIM( model_string ), &231 EXIST=iso2d_output )232 233 !234 231 !-- Inquire whether a NetCDF file exists 235 232 INQUIRE( FILE='DATA_2D_'//modus//'_NETCDF'//TRIM( model_string ), & … … 387 384 !-- File from PE0 contains special information at the beginning, 388 385 !-- concerning the lower and upper indices of the total-domain 389 !-- used in PALM (nxag, nxeg, nyag, nyeg) and the lower and 390 !-- upper indices of the array to be written by this routine 391 !-- (nxa, nxe, nya, nye). Usually in the horizontal directions 392 !-- nxag=-1 and nxa=0 while all other variables have the same 393 !-- value (i.e. nxeg=nxe). 386 !-- used in PALM (nxa, nxe, nya, nye). 394 387 !-- Allocate necessary arrays, open the output file and write 395 388 !-- the coordinate informations needed by ISO2D. 396 389 IF ( id == 0 .AND. fanz(0) == 0 .AND. fanz(1) == 0 ) THEN 397 READ ( id+110 ) nxag, nxeg, nyag, nyeg 390 398 391 READ ( id+110 ) nxa, nxe, nya, nye 399 ALLOCATE ( eta(nya:nye), ho(nxa:nxe), hu(nxa:nxe), & 400 pf(nxag:nxeg,nyag:nyeg) ) 401 READ ( id+110 ) dx, eta, hu, ho 402 403 392 ALLOCATE ( pf(nxa:nxe,nya:nye) ) 404 393 ! 405 394 !-- Set actual domain bounds to total domain … … 407 396 ye_do = nye 408 397 409 IF ( iso2d_output ) THEN410 OPEN ( 2, FILE='PLOT2D_'//modus//TRIM( model_string ),&411 FORM='UNFORMATTED' )412 WRITE ( 2 ) dx, eta, hu, ho413 ENDIF414 398 ENDIF 415 399 ! 416 400 !-- Read output time 417 401 IF ( netcdf_output .AND. id == 0 ) THEN 418 IF ( netcdf_1 ) THEN 419 READ ( id+110, END=998 ) simulated_time, time_step, av 420 ELSE 421 ! 422 !-- For compatibility with earlier PALM versions 423 READ ( id+110, END=998 ) simulated_time, time_step 424 av = 0 425 ENDIF 402 READ ( id+110, END=998 ) simulated_time, time_step, av 426 403 ENDIF 427 404 ! … … 460 437 ENDDO 461 438 ! 462 !-- Write the data of the total domain cross-section 463 IF ( iso2d_output ) WRITE ( 2 ) pf(nxa:nxe,nya:nye) 464 465 ! 466 !-- Write same data in NetCDF format 439 !-- Write data in NetCDF format 467 440 IF ( netcdf_output ) THEN 468 441 #if defined( __netcdf ) … … 544 517 ENDDO 545 518 CLOSE ( 2 ) 546 DEALLOCATE ( eta, ho, hu,pf )519 DEALLOCATE ( pf ) 547 520 548 521 ! … … 598 571 CALL SYSTEM_CLOCK( count, count_rate ) 599 572 cpu_start_time = REAL( count ) / REAL( count_rate ) 600 601 !602 !-- Inquire whether an avs fld file exists603 INQUIRE( FILE='PLOT3D_FLD'//TRIM( model_string ), EXIST=avs_output )604 573 605 574 ! … … 639 608 640 609 ! 641 !-- Combination only works, if data are not compressed. In that case,642 !-- PALM created a flag file (PLOT3D_COMPRESSED)643 INQUIRE ( FILE='PLOT3D_COMPRESSED'//TRIM( model_string ), &644 EXIST=compressed )645 646 !647 610 !-- Find out the number of files and open them 648 DO WHILE ( found .AND. .NOT. compressed)611 DO WHILE ( found ) 649 612 650 613 OPEN ( danz+110, & … … 678 641 PRINT*, ' 3D-data: ', danz, ' file(s) found' 679 642 ELSE 680 IF ( found .AND. compressed ) THEN 681 PRINT*, '+++ no 3D-data processing, since data are compressed' 682 ELSE 683 PRINT*, ' no 3D-data file available' 684 ENDIF 643 PRINT*, ' no 3D-data file available' 685 644 ENDIF 686 645 ENDIF … … 764 723 DO id = 0, danz-1 765 724 ! 766 !-- File from PE0 contains special information at the beginning, 767 !-- concerning the lower and upper indices of the total-domain used 768 !-- in PALM (nxag, nxeg, nyag, nyeg) and the lower and 769 !-- upper indices of the array to be written by this routine (nxa, 770 !-- nxe, nya, nye, nza, nze). Usually nxag=-1 and nxa=0, nyag=-1 771 !-- and nya=0, nzeg=nz and nze=nz_do3d. 772 !-- Allocate necessary array and open the output file. 725 !-- File from PE0 contains at the beginning the index bounds 726 !-- of PALM's total domain. 727 !-- Allocate the array for storing the total domain data 773 728 IF ( id == 0 .AND. fanz(0) == 0 .AND. fanz(1) == 0 ) THEN 774 READ ( id+110 ) nxag, nxeg, nyag, nyeg729 ! READ ( id+110 ) nxag, nxeg, nyag, nyeg 775 730 READ ( id+110 ) nxa, nxe, nya, nye, nza, nze 776 731 ALLOCATE ( pf3d(nxa:nxe,nya:nye,nza:nze) ) 777 IF ( avs_output ) THEN778 OPEN ( 2, FILE='PLOT3D_DATA'//TRIM( model_string ), &779 FORM='UNFORMATTED' )780 ENDIF781 732 ENDIF 782 733 … … 818 769 819 770 ! 820 !-- Write data of the total domain 821 IF ( avs_output ) WRITE ( 2 ) pf3d(nxa:nxe,nya:nye,nza:nze) 822 823 ! 824 !-- Write same data in NetCDF format 771 !-- Write data of the total domain in NetCDF format 825 772 IF ( netcdf_output ) THEN 826 773 #if defined( __netcdf )
Note: See TracChangeset
for help on using the changeset viewer.