Ignore:
Timestamp:
Oct 4, 2017 8:26:59 AM (4 years ago)
Author:
raasch
Message:

upper bounds of cross section and 3d output changed from nx+1,ny+1 to nx,ny; no output if redundant ghost layer data to NetCDF files

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/check_open.f90

    r2300 r2512  
    2525! -----------------
    2626! $Id$
     27! upper bounds of cross section and 3d output changed from nx+1,ny+1 to nx,ny
     28! no output of ghost layer data
     29! iso2d-related parts removed
     30!
     31! 2300 2017-06-29 13:31:14Z raasch
    2732! -host
    2833!
     
    178183    CHARACTER (LEN=2)   ::  mask_char               !<
    179184    CHARACTER (LEN=2)   ::  suffix                  !<
    180     CHARACTER (LEN=20)  ::  xtext = 'time in s'     !<
    181185    CHARACTER (LEN=30)  ::  filename                !<
    182186    CHARACTER (LEN=80)  ::  rtext                   !<
    183     CHARACTER (LEN=100) ::  batch_scp               !<
    184187    CHARACTER (LEN=100) ::  line                    !<
    185     CHARACTER (LEN=400) ::  command                 !<
    186188
    187189    INTEGER(iwp) ::  av          !<
    188     INTEGER(iwp) ::  numline = 1 !<
    189     INTEGER(iwp) ::  cranz       !<
    190190    INTEGER(iwp) ::  file_id     !<
    191191    INTEGER(iwp) ::  i           !<
    192     INTEGER(iwp) ::  iaddres     !<
    193192    INTEGER(iwp) ::  ioerr       !< IOSTAT flag for IO-commands ( 0 = no error )
    194     INTEGER(iwp) ::  iusern      !<
    195193    INTEGER(iwp) ::  j           !<
    196194    INTEGER(iwp) ::  k           !<
    197     INTEGER(iwp) ::  legpos = 1  !<
    198     INTEGER(iwp) ::  timodex = 1 !<
    199195   
    200     INTEGER(iwp), DIMENSION(10) ::  klist !<
    201 
    202     LOGICAL ::  datleg = .TRUE.               !<
    203     LOGICAL ::  get_filenames                 !<
    204     LOGICAL ::  grid = .TRUE.                 !<
    205     LOGICAL ::  netcdf_extend                 !<
    206     LOGICAL ::  rand = .TRUE.                 !<
    207     LOGICAL ::  swap = .TRUE.                 !<
    208     LOGICAL ::  twoxa = .TRUE.                !<
    209     LOGICAL ::  twoya = .TRUE.                !<
    210 
    211     REAL(wp) ::  ansx = -999.999_wp !<
    212     REAL(wp) ::  ansy = -999.999_wp !<
    213     REAL(wp) ::  gwid = 0.1_wp      !<
    214     REAL(wp) ::  rlegfak = 1.5_wp   !<
    215     REAL(wp) ::  sizex = 250.0_wp   !<
    216     REAL(wp) ::  sizey = 40.0_wp    !<
    217     REAL(wp) ::  texfac = 1.5_wp    !<
    218 
    219     REAL(wp), DIMENSION(:), ALLOCATABLE      ::  eta !<
    220     REAL(wp), DIMENSION(:), ALLOCATABLE      ::  ho  !<
    221     REAL(wp), DIMENSION(:), ALLOCATABLE      ::  hu  !<
    222  
    223 
    224 
    225     NAMELIST /RAHMEN/  numline, cranz, datleg, rtext, swap
    226     NAMELIST /CROSS/   ansx, ansy, grid, gwid, klist, legpos,                  &
    227                        rand, rlegfak, sizex, sizey, texfac,                    &
    228                        timodex, twoxa, twoya, xtext
    229                        
     196    LOGICAL ::  netcdf_extend    !<
    230197
    231198!
     
    422389          IF ( myid == 0  .AND.  .NOT. openfile(file_id)%opened_before )  THEN
    423390!
    424 !--          Output for combine_plot_fields
     391!--          Write index bounds of total domain for combine_plot_fields
    425392             IF ( data_output_2d_on_each_pe  .AND.  myid_char /= '' )  THEN
    426                 WRITE (21)  -nbgp, nx+nbgp, -nbgp, ny+nbgp    ! total array size
    427                 WRITE (21)   0, nx+1,  0, ny+1    ! output part
    428              ENDIF
    429 !
    430 !--          Determine and write ISO2D coordiante header
    431              ALLOCATE( eta(0:ny+1), ho(0:nx+1), hu(0:nx+1) )
    432              hu = 0.0_wp
    433              ho = (ny+1) * dy
    434              DO  i = 1, ny
    435                 eta(i) = REAL( i ) / ( ny + 1.0_wp )
    436              ENDDO
    437              eta(0)    = 0.0_wp
    438              eta(ny+1) = 1.0_wp
    439 
    440              WRITE (21)  dx,eta,hu,ho
    441              DEALLOCATE( eta, ho, hu )
     393                WRITE (21)   0, nx,  0, ny
     394             ENDIF
    442395
    443396          ENDIF
     
    455408          IF ( myid == 0  .AND.  .NOT. openfile(file_id)%opened_before )  THEN
    456409!
    457 !--          Output for combine_plot_fields
     410!--          Write index bounds of total domain for combine_plot_fields
    458411             IF ( data_output_2d_on_each_pe  .AND.  myid_char /= '' )  THEN
    459                 WRITE (22)  -nbgp, nx+nbgp, 0, nz+1    ! total array size
    460                 WRITE (22)   0, nx+1, 0, nz+1    ! output part
    461              ENDIF
    462 !
    463 !--          Determine and write ISO2D coordinate header
    464              ALLOCATE( eta(0:nz+1), ho(0:nx+1), hu(0:nx+1) )
    465              hu = 0.0_wp
    466              ho = zu(nz+1)
    467              DO  i = 1, nz
    468                 eta(i) = REAL( zu(i) ) / zu(nz+1)
    469              ENDDO
    470              eta(0)    = 0.0_wp
    471              eta(nz+1) = 1.0_wp
    472 
    473              WRITE (22)  dx,eta,hu,ho
    474              DEALLOCATE( eta, ho, hu )
     412                WRITE (22)   0, nx, 0, nz+1    ! output part
     413             ENDIF
    475414
    476415          ENDIF
     
    488427          IF ( myid == 0  .AND.  .NOT. openfile(file_id)%opened_before )  THEN
    489428!
    490 !--          Output for combine_plot_fields
     429!--          Write index bounds of total domain for combine_plot_fields
    491430             IF ( data_output_2d_on_each_pe  .AND.  myid_char /= '' )  THEN
    492                 WRITE (23)  -nbgp, ny+nbgp, 0, nz+1    ! total array size
    493                 WRITE (23)   0, ny+1, 0, nz+1    ! output part
    494              ENDIF
    495 !
    496 !--          Determine and write ISO2D coordiante header
    497              ALLOCATE( eta(0:nz+1), ho(0:ny+1), hu(0:ny+1) )
    498              hu = 0.0_wp
    499              ho = zu(nz+1)
    500              DO  i = 1, nz
    501                 eta(i) = REAL( zu(i) ) / zu(nz+1)
    502              ENDDO
    503              eta(0)    = 0.0_wp
    504              eta(nz+1) = 1.0_wp
    505 
    506              WRITE (23)  dx,eta,hu,ho
    507              DEALLOCATE( eta, ho, hu )
     431                WRITE (23)   0, ny, 0, nz+1    ! output part
     432             ENDIF
    508433
    509434          ENDIF
     
    514439                     FORM='UNFORMATTED' )
    515440!
    516 !--       Write coordinate file for AVS
     441!--       Specifications for combine_plot_fields
    517442          IF ( myid == 0 )  THEN
    518443#if defined( __parallel )
    519 !
    520 !--          Specifications for combine_plot_fields
    521              WRITE ( 30 )  -nbgp,nx+nbgp,-nbgp,ny+nbgp
    522              WRITE ( 30 )  0,nx+1,0,ny+1,0,nz_do3d
     444             WRITE ( 30 )  0, nx, 0, ny, 0, nz_do3d
    523445#endif
    524446          ENDIF
Note: See TracChangeset for help on using the changeset viewer.