Changeset 4898


Ignore:
Timestamp:
Mar 4, 2021 3:49:52 PM (4 years ago)
Author:
raasch
Message:

MPI-IO for sky view factors implemented

File:
1 edited

Legend:

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

    r4896 r4898  
    2727! -----------------
    2828! $Id$
     29! MPI-IO for sky view factors implemented
     30!
     31! 4896 2021-03-03 16:10:18Z raasch
    2932! typo in file appendix removed
    3033!
     
    477480               pt_surface,                                                                         &
    478481               read_svf,                                                                           &
     482               restart_data_format_input,                                                          &
    479483               restart_data_format_output,                                                         &
    480484               rho_surface,                                                                        &
     
    506510               nxr,                                                                                &
    507511               nxrg,                                                                               &
     512               nx_on_file,                                                                         &
    508513               ny,                                                                                 &
    509514               nyn,                                                                                &
     
    511516               nys,                                                                                &
    512517               nysg,                                                                               &
     518               ny_on_file,                                                                         &
    513519               nzb,                                                                                &
    514520               nzt,                                                                                &
     
    595601        ONLY:  rd_mpi_io_check_array,                                                              &
    596602               rrd_mpi_io,                                                                         &
    597                wrd_mpi_io
     603               wrd_mpi_io,                                                                         &
     604               rd_mpi_io_open,                                                                     &
     605               rd_mpi_io_surface_filetypes,                                                        &
     606               rrd_mpi_io_surface,                                                                 &
     607               wrd_mpi_io_surface,                                                                 &
     608               rd_mpi_io_close,                                                                    &
     609               tgh
    598610
    599611    USE statistics,                                                                                &
     
    34243436
    34253437!
    3426 !-- Find all discretized apparent solar positions for radiation interaction.
    3427     IF ( radiation_interactions )  CALL radiation_presimulate_solar_pos
    3428 
    3429 !
    3430 !-- If required, read or calculate and write out the SVF
    3431     IF ( radiation_interactions .AND. read_svf )  THEN
    3432 !
    3433 !--    Read sky-view factors and further required data from file
    3434        CALL radiation_read_svf()
    3435 
    3436     ELSEIF ( radiation_interactions .AND. .NOT. read_svf )  THEN
    3437 !
    3438 !--    Calculate SFV and CSF
    3439        CALL radiation_calc_svf()
     3438!-- If required, read or calculate and write out the SVF.
     3439    IF ( radiation_interactions )  THEN
     3440
     3441!
     3442!--    Find all discretized apparent solar positions for radiation interaction.
     3443       CALL radiation_presimulate_solar_pos
     3444
     3445       IF ( read_svf )  THEN
     3446!
     3447!--       Read sky-view factors and further required data from file
     3448          CALL radiation_read_svf()
     3449
     3450       ELSE
     3451!
     3452!--       Calculate svf and csf.
     3453          CALL radiation_calc_svf()
     3454       ENDIF
     3455
     3456       IF ( write_svf )  THEN
     3457!
     3458!--       Write svf, csf svfsurf and csfsurf data to file.
     3459          CALL radiation_write_svf()
     3460       ENDIF
     3461
     3462!
     3463!--    Adjust radiative fluxes.
     3464!--    In case of urban and land surfaces, also call an initial interaction.
     3465       CALL radiation_interaction
     3466
    34403467    ENDIF
    34413468
    3442     IF ( radiation_interactions .AND. write_svf )  THEN
    3443 !
    3444 !--    Write svf, csf svfsurf and csfsurf data to file
    3445        CALL radiation_write_svf()
    3446     ENDIF
    3447 
    3448 !
    3449 !-- Adjust radiative fluxes. In case of urban and land surfaces, also call an initial interaction.
    3450     IF ( radiation_interactions )  THEN
    3451        CALL radiation_interaction
    3452     ENDIF
    3453 
    34543469    IF ( debug_output )  CALL debug_message( 'radiation_init', 'end' )
    3455 
    3456     RETURN ! @Todo: remove, I don't see what we need this for here
    34573470
    34583471 END SUBROUTINE radiation_init
     
    94509463    CHARACTER(rad_version_len) ::  rad_version_field  !<
    94519464
     9465    CHARACTER(LEN=64) ::  rad_version_in  !< rad version on file
     9466
    94529467    INTEGER(iwp) ::  i                      !<
     9468    INTEGER(iwp) ::  ncsfl_tot              !<
    94539469    INTEGER(iwp) ::  ndsidir_from_file = 0  !<
     9470    INTEGER(iwp) ::  ndsidir_tot            !<
    94549471    INTEGER(iwp) ::  npcbl_from_file   = 0  !<
     9472    INTEGER(iwp) ::  npcbl_tot              !<
    94559473    INTEGER(iwp) ::  nsurfl_from_file  = 0  !<
     9474    INTEGER(iwp) ::  nsurfl_tot             !<
     9475    INTEGER(iwp) ::  nsvfl_tot              !< total (sum across all PEs) counter for the different IO variables
    94569476    INTEGER(iwp) ::  nmrtbl_from_file  = 0  !<
    9457 
    9458 
    9459     CALL location_message( 'reading view factors for radiation interaction', 'start' )
    9460 
    9461     DO  i = 0, io_blocks-1
    9462        IF ( i == io_group )  THEN
    9463 
    9464 !
    9465 !--       numprocs_previous_run is only known in case of reading restart data. If a new initial run
    9466 !--       which reads svf data is started the following query will be skipped.
    9467           IF ( initializing_actions == 'read_restart_data' )  THEN
    9468 
    9469              IF ( numprocs_previous_run /= numprocs )  THEN
    9470                 WRITE( message_string, * ) 'A different number of processors between the run ',    &
    9471                                            'that has written the svf data and the one that ',      &
    9472                                            'will read it is not allowed'
    9473                 CALL message( 'check_open', 'PA0491', 1, 2, 0, 6, 0 )
     9477    INTEGER(iwp) ::  nmrtbl_tot             !<
     9478    INTEGER(iwp) ::  nmrtf_tot              !<
     9479
     9480    INTEGER(iwp), DIMENSION(4) ::  global_sum  !<
     9481    INTEGER(iwp), DIMENSION(4) ::  local_sum   !< variables to compute total counter
     9482
     9483    INTEGER(iwp), DIMENSION(nys:nyn,nxl:nxr) ::  global_end    !< global end index
     9484    INTEGER(iwp), DIMENSION(nys:nyn,nxl:nxr) ::  global_start  !< global start index
     9485    INTEGER(iwp), DIMENSION(nys:nyn,nxl:nxr) ::  end_index     !< local end index
     9486    INTEGER(iwp), DIMENSION(nys:nyn,nxl:nxr) ::  start_index   !< local start index
     9487
     9488    LOGICAL ::  data_to_read  !< flag indicating if data is available for current variable
     9489
     9490    REAL(wp), ALLOCATABLE, DIMENSION(:,:) ::  tmp  !<
     9491
     9492
     9493    IF ( TRIM( restart_data_format_input ) == 'fortran_binary' )  THEN
     9494
     9495       CALL location_message( 'reading sky view factors for radiation interaction ' //             &
     9496                              '(Fortran binary format)', 'start' )
     9497
     9498       DO  i = 0, io_blocks-1
     9499          IF ( i == io_group )  THEN
     9500
     9501!
     9502!--          numprocs_previous_run is only known in case of reading restart data. If a new initial
     9503!--          run which reads svf data is started the following query will be skipped.
     9504             IF ( initializing_actions == 'read_restart_data' )  THEN
     9505
     9506                IF ( numprocs_previous_run /= numprocs )  THEN
     9507                   WRITE( message_string, * ) 'A different number of processors between the run ', &
     9508                                              'that has written the svf data and the one that ',   &
     9509                                              'will read it is not allowed'
     9510                   CALL message( 'check_open', 'PA0491', 1, 2, 0, 6, 0 )
     9511                ENDIF
     9512
    94749513             ENDIF
    94759514
     9515!
     9516!--          Open binary file
     9517             CALL check_open( 88 )
     9518
     9519!
     9520!--          Read and check version
     9521             READ ( 88 ) rad_version_field
     9522             IF ( TRIM( rad_version_field ) /= TRIM( rad_version ) )  THEN
     9523                 WRITE( message_string, * ) 'Version of binary SVF file "',                        &
     9524                                            TRIM( rad_version_field ), '" does not match ',        &
     9525                                            'the version of model "', TRIM( rad_version ), '"'
     9526                 CALL message( 'radiation_read_svf', 'PA0482', 1, 2, 0, 6, 0 )
     9527             ENDIF
     9528
     9529!
     9530!--          Read nsvfl, ncsfl, nsurfl, nmrtf
     9531             READ ( 88 ) nsvfl, ncsfl, nsurfl_from_file, npcbl_from_file, ndsidir_from_file,       &
     9532                         nmrtbl_from_file, nmrtf
     9533
     9534             IF ( nsvfl < 0  .OR.  ncsfl < 0 )  THEN
     9535                 WRITE( message_string, * ) 'Wrong number of SVF or CSF'
     9536                 CALL message( 'radiation_read_svf', 'PA0483', 1, 2, 0, 6, 0 )
     9537             ELSE
     9538                 WRITE(debug_string,*) 'Number of SVF, CSF, and nsurfl to read', nsvfl, ncsfl,     &
     9539                                       nsurfl_from_file
     9540                 IF ( debug_output )  CALL debug_message( debug_string, 'info' )
     9541             ENDIF
     9542
     9543             IF ( nsurfl_from_file /= nsurfl )  THEN
     9544                 WRITE( message_string, * ) 'nsurfl from SVF file does not match calculated ',     &
     9545                                            'nsurfl from radiation_interaction_init'
     9546                 CALL message( 'radiation_read_svf', 'PA0490', 1, 2, 0, 6, 0 )
     9547             ENDIF
     9548
     9549             IF ( npcbl_from_file /= npcbl )  THEN
     9550                 WRITE( message_string, * ) 'npcbl from SVF file does not match calculated npcbl', &
     9551                                            ' from radiation_interaction_init'
     9552                 CALL message( 'radiation_read_svf', 'PA0493', 1, 2, 0, 6, 0 )
     9553             ENDIF
     9554
     9555             IF ( ndsidir_from_file /= ndsidir )  THEN
     9556                 WRITE( message_string, * ) 'ndsidir from SVF file does not match calculated ',    &
     9557                                            'ndsidir from radiation_presimulate_solar_pos'
     9558                 CALL message( 'radiation_read_svf', 'PA0494', 1, 2, 0, 6, 0 )
     9559             ENDIF
     9560
     9561             IF ( nmrtbl_from_file /= nmrtbl )  THEN
     9562                 WRITE( message_string, * ) 'nmrtbl from SVF file does not match calculated ',     &
     9563                                            'nmrtbl from radiation_interaction_init'
     9564                 CALL message( 'radiation_read_svf', 'PA0494', 1, 2, 0, 6, 0 )
     9565             ELSE
     9566                 WRITE( debug_string, * ) 'Number of nmrtf to read ', nmrtf
     9567                 IF ( debug_output )  CALL debug_message( debug_string, 'info' )
     9568             ENDIF
     9569
     9570!
     9571!--          Arrays skyvf, skyvft, dsitrans and dsitransc are allready allocated in
     9572!--          radiation_interaction_init and radiation_presimulate_solar_pos
     9573             IF ( nsurfl > 0 )  THEN
     9574                READ( 88 ) skyvf
     9575                READ( 88 ) skyvft
     9576                READ( 88 ) dsitrans
     9577             ENDIF
     9578
     9579             IF ( plant_canopy  .AND.  npcbl > 0 )  THEN
     9580                READ( 88 )  dsitransc
     9581             ENDIF
     9582
     9583!
     9584!--          The allocation of svf, svfsurf, csf, csfsurf, mrtf, mrtft, and mrtfsurf happens in
     9585!--          routine radiation_calc_svf which is not called if the program enters
     9586!--          radiation_read_svf. Therefore these arrays have to be allocated in the following.
     9587             IF ( nsvfl > 0 )  THEN
     9588                ALLOCATE( svf(ndsvf,nsvfl) )
     9589                ALLOCATE( svfsurf(idsvf,nsvfl) )
     9590                READ( 88 ) svf
     9591                READ( 88 ) svfsurf
     9592             ENDIF
     9593
     9594             IF ( plant_canopy  .AND.  ncsfl > 0 )  THEN
     9595                ALLOCATE( csf(ndcsf,ncsfl) )
     9596                ALLOCATE( csfsurf(idcsf,ncsfl) )
     9597                READ( 88 ) csf
     9598                READ( 88 ) csfsurf
     9599             ENDIF
     9600
     9601             IF ( nmrtbl > 0 )  THEN
     9602                READ( 88 ) mrtsky
     9603                READ( 88 ) mrtskyt
     9604                READ( 88 ) mrtdsit
     9605             ENDIF
     9606
     9607             IF ( nmrtf > 0 )  THEN
     9608                ALLOCATE( mrtf(nmrtf) )
     9609                ALLOCATE( mrtft(nmrtf) )
     9610                ALLOCATE( mrtfsurf(2,nmrtf) )
     9611                READ( 88 ) mrtf
     9612                READ( 88 ) mrtft
     9613                READ( 88 ) mrtfsurf
     9614             ENDIF
     9615
     9616!
     9617!--          Close binary file
     9618             CALL close_file( 88 )
     9619
    94769620          ENDIF
    9477 
    9478 !
    9479 !--       Open binary file
    9480           CALL check_open( 88 )
    9481 
    9482 !
    9483 !--       Read and check version
    9484           READ ( 88 ) rad_version_field
    9485           IF ( TRIM( rad_version_field ) /= TRIM( rad_version ) )  THEN
    9486               WRITE( message_string, * ) 'Version of binary SVF file "',                           &
    9487                                          TRIM( rad_version_field ), '" does not match ',           &
    9488                                          'the version of model "', TRIM( rad_version ), '"'
    9489               CALL message( 'radiation_read_svf', 'PA0482', 1, 2, 0, 6, 0 )
     9621#if defined( __parallel )
     9622          CALL MPI_BARRIER( comm2d, ierr )
     9623#endif
     9624       ENDDO
     9625
     9626       CALL location_message( 'reading sky view factors for radiation interaction ' //             &
     9627                              '(Fortran binary format)', 'finished' )
     9628
     9629
     9630    ELSEIF ( restart_data_format_input(1:3) == 'mpi' )  THEN
     9631
     9632!
     9633!--    In case of MPI-IO data is treated like surface data and the respective routines from
     9634!--    restart_data_mpi_io_mod are used for reading. Contrary to restart data, the virtual processor
     9635!--    grid can not be changed between writing and reading svf data.
     9636       CALL location_message( 'reading sky view factors for radiation interaction ' //             &
     9637                              '(MPI-IO format)', 'start' )
     9638!
     9639!--    Open MPI-IO svf file for reading the global data.
     9640       CALL rd_mpi_io_open( 'READ', 'SVFIN' // TRIM( coupling_char ),                              &
     9641                            open_for_global_io_only = .TRUE. )
     9642!
     9643!--    Check general header.
     9644       IF ( tgh%pes_along_x /= npex  .OR.  tgh%pes_along_y /= npey )  THEN
     9645!
     9646!--       Force re-calculation of svfs.
     9647          read_svf = .FALSE.
     9648          WRITE( message_string, '(A,I7,A,I7,A,I7,A,I7,A)' )                                       &
     9649              'virtual PE grid has changed between previous and current run &npex_prev = ',        &
     9650              tgh%pes_along_x, ' npey_prev = ', tgh%pes_along_y, ' npex_new = ', npex,             &
     9651              ' npey_new = ', npey, '&svf wil be re-calculated'
     9652          CALL message( 'radiation_read_svf', 'PA0517', 0, 0, 0, 6, 0 )
     9653          RETURN
     9654       ENDIF
     9655!
     9656!--    Read global variables
     9657       CALL rrd_mpi_io( 'rad_version', rad_version_in )
     9658       CALL rrd_mpi_io( 'nsvfl', nsvfl_tot )
     9659       CALL rrd_mpi_io( 'ncsfl', ncsfl_tot )
     9660       CALL rrd_mpi_io( 'nsurfl', nsurfl_from_file )
     9661       CALL rrd_mpi_io( 'npcbl', npcbl_from_file )
     9662       CALL rrd_mpi_io( 'ndsidir', ndsidir_from_file )
     9663       CALL rrd_mpi_io( 'nmrtbl', nmrtbl_from_file )
     9664       CALL rrd_mpi_io( 'nmrtf', nmrtf_tot )
     9665
     9666       CALL rd_mpi_io_close
     9667!
     9668!--    Compute global values of local counters.
     9669       local_sum(1) = nsurfl
     9670       local_sum(2) = npcbl
     9671       local_sum(3) = ndsidir
     9672       local_sum(4) = nmrtbl
     9673#if defined( __parallel )
     9674       CALL MPI_ALLREDUCE( local_sum, global_sum, SIZE(local_sum), MPI_INTEGER, MPI_SUM, comm2d,   &
     9675                           ierr)
     9676#else
     9677       global_sum = local_sum
     9678#endif
     9679       nsurfl_tot  = global_sum(1)
     9680       npcbl_tot   = global_sum(2)
     9681       ndsidir_tot = global_sum(3)
     9682       nmrtbl_tot  = global_sum(4)
     9683!
     9684!--    Check for errors.
     9685       nx_on_file = tgh%total_nx-1
     9686       ny_on_file = tgh%total_ny-1
     9687
     9688       IF ( nx_on_file /= nx  .OR.  ny_on_file /= ny )  THEN
     9689          WRITE( message_string, '(A,4(A,I7))' )                                                   &
     9690               'total number of grid points along x and y in file SVFIN do not match current run', &
     9691               '&nx_on_file = ', nx_on_file, ' ny_on_file = ', ny_on_file, ' nx = ', nx, ' ny = ', &
     9692               ny
     9693          CALL message( 'radiation_read_svf', 'PA0518', 1, 2, 0, 6, 0 )
     9694       ENDIF
     9695
     9696       IF ( TRIM( rad_version_in ) /= TRIM( rad_version ) )  THEN
     9697          WRITE( message_string, * ) 'Version of binary SVF file "', TRIM( rad_version_field ),    &
     9698                                     '" does not match the version of model "',                    &
     9699                                     TRIM( rad_version ), '"'
     9700          CALL message( 'radiation_read_svf', 'PA0482', 1, 2, 0, 6, 0 )
     9701       ENDIF
     9702
     9703       IF ( nsvfl_tot < 0  .OR.  ncsfl_tot < 0 )  THEN
     9704          WRITE( message_string, * ) 'Wrong number of SVF or CSF'
     9705          CALL message( 'radiation_read_svf', 'PA0483', 1, 2, 0, 6, 0 )
     9706       ELSE
     9707          WRITE(debug_string,*) 'Number of SVF, CSF, and nsurfl to read', nsvfl_tot, ncsfl_tot,  &
     9708                                          nsurfl_tot
     9709          IF ( debug_output )  CALL debug_message( debug_string, 'info' )
     9710       ENDIF
     9711
     9712       IF ( nsurfl_from_file /= nsurfl_tot )  THEN
     9713          WRITE( message_string, * ) 'nsurfl from SVF file does not match calculated ',        &
     9714             'nsurfl from radiation_interaction_init'
     9715          CALL message( 'radiation_read_svf', 'PA0490', 1, 2, 0, 6, 0 )
     9716       ENDIF
     9717
     9718       IF ( npcbl_from_file /= npcbl_tot )  THEN
     9719          WRITE( message_string, * ) 'npcbl from SVF file does not match calculated npcbl ',   &
     9720             'from radiation_interaction_init'
     9721          CALL message( 'radiation_read_svf', 'PA0493', 1, 2, 0, 6, 0 )
     9722       ENDIF
     9723
     9724       IF ( ndsidir_from_file /= ndsidir_tot )  THEN
     9725          WRITE( message_string, * ) 'ndsidir from SVF file does not match calculated ',       &
     9726             'ndsidir from radiation_presimulate_solar_pos'
     9727          CALL message( 'radiation_read_svf', 'PA0494', 1, 2, 0, 6, 0 )
     9728       ENDIF
     9729
     9730       IF ( nmrtbl_from_file /= nmrtbl_tot )  THEN
     9731          WRITE( message_string, * ) 'nmrtbl from SVF file does not match calculated nmrtbl ', &
     9732             'from radiation_interaction_init'
     9733          CALL message( 'radiation_read_svf', 'PA0494', 1, 2, 0, 6, 0 )
     9734       ELSE
     9735          WRITE( debug_string, * ) 'Number of nmrtf to read ', nmrtf_tot
     9736          IF ( debug_output )  CALL debug_message( debug_string, 'info' )
     9737       ENDIF
     9738
     9739!
     9740!--    Open MPI-IO SVF file for read local data.
     9741       CALL rd_mpi_io_open( 'READ', 'SVFIN' // TRIM( coupling_char ) )
     9742
     9743       IF ( nsurfl_tot > 0 )  THEN
     9744!
     9745!--       Read global indices.
     9746          CALL rrd_mpi_io( 'nsurfl_global_start', global_start )
     9747          CALL rrd_mpi_io( 'nsurfl_global_end', global_end )
     9748!
     9749!--       Set file types of variables and compute local indices.
     9750          CALL rd_mpi_io_surface_filetypes( start_index, end_index, data_to_read, global_start,    &
     9751                                            global_end )
     9752          nsurfl = end_index(nyn,nxr)
     9753
     9754          IF ( data_to_read )  THEN
     9755
     9756             CALL rrd_mpi_io_surface( 'skyvf', skyvf )
     9757             CALL rrd_mpi_io_surface( 'skyvft', skyvft )
     9758!
     9759!--          To avoid another overlay of rrd_mpi_io_surface, dsitrans is read as REAL tmp array.
     9760!--          The order of dimensions of dsitrans is different to the order expected by
     9761!--          rrd_mpi_io_surface. Therefor a tranpose of tmp is required.
     9762             ALLOCATE( tmp(SIZE(dsitrans,2),SIZE(dsitrans,1)) )
     9763             CALL rrd_mpi_io_surface( 'dsitrans', tmp )
     9764             dsitrans = TRANSPOSE( tmp )
     9765             DEALLOCATE( tmp )
     9766
    94909767          ENDIF
    94919768
    9492 !
    9493 !--       Read nsvfl, ncsfl, nsurfl, nmrtf
    9494           READ ( 88 ) nsvfl, ncsfl, nsurfl_from_file, npcbl_from_file, ndsidir_from_file,          &
    9495                       nmrtbl_from_file, nmrtf
    9496 
    9497           IF ( nsvfl < 0  .OR.  ncsfl < 0 )  THEN
    9498               WRITE( message_string, * ) 'Wrong number of SVF or CSF'
    9499               CALL message( 'radiation_read_svf', 'PA0483', 1, 2, 0, 6, 0 )
    9500           ELSE
    9501               WRITE(debug_string,*) 'Number of SVF, CSF, and nsurfl to read', nsvfl, ncsfl,        &
    9502                                     nsurfl_from_file
    9503               IF ( debug_output )  CALL debug_message( debug_string, 'info' )
     9769       ENDIF
     9770
     9771       IF ( npcbl_tot > 0 )  THEN
     9772
     9773          CALL rrd_mpi_io( 'npcbl_global_start', global_start )
     9774          CALL rrd_mpi_io( 'npcbl_global_end', global_end )
     9775
     9776          CALL rd_mpi_io_surface_filetypes( start_index, end_index, data_to_read, global_start,    &
     9777                                            global_end )
     9778          npcbl = end_index(nyn,nxr)
     9779
     9780          IF ( data_to_read )  THEN
     9781
     9782             ALLOCATE( tmp(SIZE(dsitransc,2),SIZE(dsitransc,1)) )
     9783             CALL rrd_mpi_io_surface( 'dsitransc', tmp )
     9784             dsitransc = TRANSPOSE( tmp )
     9785             DEALLOCATE( tmp )
     9786
    95049787          ENDIF
    95059788
    9506           IF ( nsurfl_from_file /= nsurfl )  THEN
    9507               WRITE( message_string, * ) 'nsurfl from SVF file does not match calculated ',        &
    9508                                          'nsurfl from radiation_interaction_init'
    9509               CALL message( 'radiation_read_svf', 'PA0490', 1, 2, 0, 6, 0 )
     9789       ENDIF
     9790
     9791       IF ( nsvfl_tot > 0 )  THEN
     9792
     9793          CALL rrd_mpi_io( 'nsvfl_global_start', global_start )
     9794          CALL rrd_mpi_io( 'nsvfl_global_end', global_end )
     9795
     9796          CALL rd_mpi_io_surface_filetypes( start_index, end_index, data_to_read, global_start,    &
     9797                                            global_end )
     9798          nsvfl = end_index(nyn,nxr)
     9799
     9800          IF ( .NOT. ALLOCATED( svf )     )  ALLOCATE( svf(ndsvf,nsvfl) )
     9801          IF ( .NOT. ALLOCATED( svfsurf ) )  ALLOCATE( svfsurf(idsvf,nsvfl) )
     9802
     9803          IF ( data_to_read )  THEN
     9804
     9805             CALL rrd_mpi_io_surface( 'svf', svf )
     9806             ALLOCATE( tmp(SIZE(svfsurf,1),SIZE(svfsurf,2)) )
     9807             CALL rrd_mpi_io_surface( 'svfsurf', tmp )
     9808             svfsurf = tmp
     9809             DEALLOCATE( tmp )
     9810
    95109811          ENDIF
    95119812
    9512           IF ( npcbl_from_file /= npcbl )  THEN
    9513               WRITE( message_string, * ) 'npcbl from SVF file does not match calculated npcbl ',   &
    9514                                          'from radiation_interaction_init'
    9515               CALL message( 'radiation_read_svf', 'PA0493', 1, 2, 0, 6, 0 )
     9813       ENDIF
     9814
     9815       IF ( plant_canopy  )  THEN
     9816
     9817          IF ( ncsfl_tot > 0 )  THEN
     9818
     9819             CALL rrd_mpi_io( 'ncsfl_global_start', global_start )
     9820             CALL rrd_mpi_io( 'ncsfl_global_end', global_end )
     9821
     9822             CALL rd_mpi_io_surface_filetypes( start_index, end_index, data_to_read, global_start, &
     9823                                               global_end )
     9824             ncsfl = end_index(nyn,nxr)
     9825
     9826             IF ( .NOT. ALLOCATED( csf )     )  ALLOCATE( csf(ndcsf,ncsfl) )
     9827             IF ( .NOT. ALLOCATED( csfsurf ) )  ALLOCATE( csfsurf(idcsf,ncsfl) )
     9828
     9829             IF ( data_to_read )  THEN
     9830
     9831                CALL rrd_mpi_io_surface( 'csf', csf )
     9832                ALLOCATE( tmp(SIZE(csfsurf,1),SIZE(csfsurf,2)) )
     9833                CALL rrd_mpi_io_surface( 'csfsurf', tmp )
     9834                csfsurf = tmp
     9835                DEALLOCATE( tmp )
     9836
     9837             ENDIF
     9838
    95169839          ENDIF
    95179840
    9518           IF ( ndsidir_from_file /= ndsidir )  THEN
    9519               WRITE( message_string, * ) 'ndsidir from SVF file does not match calculated ',       &
    9520                                          'ndsidir from radiation_presimulate_solar_pos'
    9521               CALL message( 'radiation_read_svf', 'PA0494', 1, 2, 0, 6, 0 )
     9841       ENDIF
     9842
     9843       IF ( nmrtbl_tot > 0 )  THEN
     9844
     9845          CALL rrd_mpi_io( 'nmrtbl_global_start', global_start )
     9846          CALL rrd_mpi_io( 'nmrtbl_global_end', global_end )
     9847
     9848          CALL rd_mpi_io_surface_filetypes( start_index, end_index, data_to_read, global_start,    &
     9849                                            global_end )
     9850          nmrtbl = end_index(nyn,nxr)
     9851
     9852          IF ( data_to_read )  THEN
     9853
     9854             CALL rrd_mpi_io_surface( 'mrtsky', mrtsky )
     9855             CALL rrd_mpi_io_surface( 'mrtskyt', mrtskyt )
     9856             ALLOCATE( tmp(SIZE(mrtdsit,2),SIZE(mrtdsit,1)) )
     9857             CALL rrd_mpi_io_surface( 'mrtdsit', tmp )
     9858             mrtdsit = TRANSPOSE( tmp )
     9859             DEALLOCATE( tmp )
     9860
    95229861          ENDIF
    9523           IF ( nmrtbl_from_file /= nmrtbl )  THEN
    9524               WRITE( message_string, * ) 'nmrtbl from SVF file does not match calculated nmrtbl ', &
    9525                                          'from radiation_interaction_init'
    9526               CALL message( 'radiation_read_svf', 'PA0494', 1, 2, 0, 6, 0 )
    9527           ELSE
    9528               WRITE( debug_string, * ) 'Number of nmrtf to read ', nmrtf
    9529               IF ( debug_output )  CALL debug_message( debug_string, 'info' )
     9862
     9863       ENDIF
     9864
     9865       IF ( nmrtf_tot > 0 )  THEN
     9866
     9867          CALL rrd_mpi_io( 'nmrtf_global_start', global_start )
     9868          CALL rrd_mpi_io( 'nmrtf_global_end', global_end )
     9869
     9870          CALL rd_mpi_io_surface_filetypes( start_index, end_index, data_to_read, global_start,    &
     9871                                            global_end )
     9872          nmrtf = end_index(nyn,nxr)
     9873
     9874          IF ( .NOT. ALLOCATED( mrtf )     )  ALLOCATE( mrtf(nmrtf) )
     9875          IF ( .NOT. ALLOCATED( mrtft )    )  ALLOCATE( mrtft(nmrtf) )
     9876          IF ( .NOT. ALLOCATED( mrtfsurf ) )  ALLOCATE( mrtfsurf(2,nmrtf) )
     9877
     9878          IF ( data_to_read )  THEN
     9879
     9880             CALL rrd_mpi_io_surface( 'mrtf', mrtf )
     9881             CALL rrd_mpi_io_surface( 'mrtft', mrtft )
     9882             ALLOCATE( tmp(SIZE(mrtfsurf,1),SIZE(mrtfsurf,2)) )
     9883             CALL rrd_mpi_io_surface( 'mrtfsurf', tmp )
     9884             mrtfsurf = tmp
     9885             DEALLOCATE( tmp )
     9886
    95309887          ENDIF
    95319888
    9532 !
    9533 !--       Arrays skyvf, skyvft, dsitrans and dsitransc are allready allocated in
    9534 !--       radiation_interaction_init and radiation_presimulate_solar_pos
    9535           IF ( nsurfl > 0 )  THEN
    9536              READ( 88 ) skyvf
    9537              READ( 88 ) skyvft
    9538              READ( 88 ) dsitrans
    9539           ENDIF
    9540 
    9541           IF ( plant_canopy  .AND.  npcbl > 0 )  THEN
    9542              READ( 88 )  dsitransc
    9543           ENDIF
    9544 
    9545 !
    9546 !--       The allocation of svf, svfsurf, csf, csfsurf, mrtf, mrtft, and mrtfsurf happens in
    9547 !--       routine radiation_calc_svf which is not called if the program enters radiation_read_svf.
    9548 !--       Therefore these arrays have to be allocated in the following.
    9549           IF ( nsvfl > 0 )  THEN
    9550              ALLOCATE( svf(ndsvf,nsvfl) )
    9551              ALLOCATE( svfsurf(idsvf,nsvfl) )
    9552              READ( 88 ) svf
    9553              READ( 88 ) svfsurf
    9554           ENDIF
    9555 
    9556           IF ( plant_canopy  .AND.  ncsfl > 0 )  THEN
    9557              ALLOCATE( csf(ndcsf,ncsfl) )
    9558              ALLOCATE( csfsurf(idcsf,ncsfl) )
    9559              READ( 88 ) csf
    9560              READ( 88 ) csfsurf
    9561           ENDIF
    9562 
    9563           IF ( nmrtbl > 0 )  THEN
    9564              READ( 88 ) mrtsky
    9565              READ( 88 ) mrtskyt
    9566              READ( 88 ) mrtdsit
    9567           ENDIF
    9568 
    9569           IF ( nmrtf > 0 )  THEN
    9570              ALLOCATE( mrtf(nmrtf) )
    9571              ALLOCATE( mrtft(nmrtf) )
    9572              ALLOCATE( mrtfsurf(2,nmrtf) )
    9573              READ( 88 ) mrtf
    9574              READ( 88 ) mrtft
    9575              READ( 88 ) mrtfsurf
    9576           ENDIF
    9577 
    9578 !
    9579 !--       Close binary file
    9580           CALL close_file( 88 )
    9581 
    9582        ENDIF
    9583 #if defined( __parallel )
    9584        CALL MPI_BARRIER( comm2d, ierr )
    9585 #endif
    9586     ENDDO
    9587 
    9588     CALL location_message( 'reading view factors for radiation interaction', 'finished' )
    9589 
     9889       ENDIF
     9890
     9891       CALL rd_mpi_io_close
     9892
     9893       CALL location_message( 'reading sky view factors for radiation interaction ' //             &
     9894                              '(MPI-IO format)', 'finished' )
     9895
     9896    ENDIF
    95909897
    95919898 END SUBROUTINE radiation_read_svf
     
    96049911    IMPLICIT NONE
    96059912
    9606     INTEGER(iwp) ::  i  !<
    9607 
    9608 
    9609     CALL location_message( 'writing view factors for radiation interaction', 'start' )
    9610 
    9611     DO  i = 0, io_blocks-1
    9612        IF ( i == io_group )  THEN
    9613 !
    9614 !--       Open binary file
    9615           CALL check_open( 89 )
    9616 
    9617           WRITE( 89 ) rad_version
    9618           WRITE( 89 ) nsvfl, ncsfl, nsurfl, npcbl, ndsidir, nmrtbl, nmrtf
    9619           IF ( nsurfl > 0 )  THEN
    9620              WRITE( 89 ) skyvf
    9621              WRITE( 89 ) skyvft
    9622              WRITE( 89 ) dsitrans
     9913    INTEGER(iwp), PARAMETER ::  max_i4_value  = 2147483647  !< maximum positive INTEGER(4) value, 2**31 - 1
     9914
     9915    INTEGER(iwp) ::  i            !<
     9916    INTEGER(iwp) ::  ic           !<
     9917    INTEGER(iwp) ::  ierr         !<
     9918    INTEGER(iwp) ::  ind          !<
     9919    INTEGER(iwp) ::  ipcgb        !<
     9920    INTEGER(iwp) ::  isurf        !<
     9921    INTEGER(iwp) ::  j            !<
     9922    INTEGER(iwp) ::  jc           !<
     9923    INTEGER(iwp) ::  ncsfl_tot    !<
     9924    INTEGER(iwp) ::  ndsidir_tot  !<
     9925    INTEGER(iwp) ::  nmrtbl_tot   !<
     9926    INTEGER(iwp) ::  nmrtf_tot    !<
     9927    INTEGER(iwp) ::  npcbl_tot    !<
     9928    INTEGER(iwp) ::  nsurfl_tot   !<
     9929    INTEGER(iwp) ::  nsvfl_tot    !< total (sum over all PEs) counter for the different IO variables
     9930
     9931    INTEGER(idp), DIMENSION(7) ::  local_sum   !< variables to compute total counter
     9932    INTEGER(idp), DIMENSION(7) ::  global_sum  !< INTEGER(idp) to allow check, if total number of values > 2G
     9933
     9934    INTEGER(iwp), DIMENSION(nys:nyn,nxl:nxr) ::  end_index     !< local end index
     9935    INTEGER(iwp), DIMENSION(nys:nyn,nxl:nxr) ::  global_end    !< global end index
     9936    INTEGER(iwp), DIMENSION(nys:nyn,nxl:nxr) ::  global_start  !< global start index
     9937    INTEGER(iwp), DIMENSION(nys:nyn,nxl:nxr) ::  lo_no         !< local number of values
     9938    INTEGER(iwp), DIMENSION(nys:nyn,nxl:nxr) ::  start_index   !< local start index
     9939
     9940    LOGICAL ::  data_to_write  !<  flag indicating if data is available for writing
     9941
     9942    REAL(wp), ALLOCATABLE, DIMENSION(:,:) ::  tmp  !<
     9943
     9944
     9945    IF ( TRIM( restart_data_format_output ) == 'fortran_binary' )  THEN
     9946
     9947       CALL location_message( 'writing sky view factors for radiation interaction ' //             &
     9948                              '(Fortran binary format)', 'start' )
     9949
     9950       DO  i = 0, io_blocks-1
     9951          IF ( i == io_group )  THEN
     9952!
     9953!--          Open binary file
     9954             CALL check_open( 89 )
     9955
     9956             WRITE( 89 ) rad_version
     9957             WRITE( 89 ) nsvfl, ncsfl, nsurfl, npcbl, ndsidir, nmrtbl, nmrtf
     9958             IF ( nsurfl > 0 )  THEN
     9959                WRITE( 89 ) skyvf
     9960                WRITE( 89 ) skyvft
     9961                WRITE( 89 ) dsitrans
     9962             ENDIF
     9963             IF ( npcbl > 0 ) THEN
     9964                WRITE( 89 ) dsitransc
     9965             ENDIF
     9966             IF ( nsvfl > 0 ) THEN
     9967                WRITE( 89 ) svf
     9968                WRITE( 89 ) svfsurf
     9969             ENDIF
     9970             IF ( plant_canopy  .AND.  ncsfl > 0 )  THEN
     9971                 WRITE( 89 ) csf
     9972                 WRITE( 89 ) csfsurf
     9973             ENDIF
     9974             IF ( nmrtbl > 0 )  THEN
     9975                WRITE( 89 ) mrtsky
     9976                WRITE( 89 ) mrtskyt
     9977                WRITE( 89 ) mrtdsit
     9978             ENDIF
     9979             IF ( nmrtf > 0 )  THEN
     9980                 WRITE( 89 ) mrtf
     9981                 WRITE( 89 ) mrtft
     9982                 WRITE( 89 ) mrtfsurf
     9983             ENDIF
     9984!
     9985!--          Close binary file
     9986             CALL close_file( 89 )
     9987
    96239988          ENDIF
    9624           IF ( npcbl > 0 ) THEN
    9625              WRITE( 89 ) dsitransc
     9989#if defined( __parallel )
     9990          CALL MPI_BARRIER( comm2d, ierr )
     9991#endif
     9992       ENDDO
     9993
     9994       CALL location_message( 'writing sky view factors for radiation interaction ' //             &
     9995                              '(Fortran binary format)', 'fnished' )
     9996
     9997    ELSEIF ( restart_data_format_output(1:3) == 'mpi' )  THEN
     9998!
     9999!--    Sky view factor data is treated like surface data and the respective routines from
     10000!--    restart_data_mpi_io_mod are used for writing. In contrary to restart data, the virtual
     10001!--    processor grid can not be changed between writing and reading svf data.
     10002       CALL location_message( 'writing sky view factors for radiation interaction ' //             &
     10003                              '(MPI-IO format)', 'start' )
     10004
     10005!
     10006!--    Open MPI-IO svf file
     10007       CALL rd_mpi_io_open( 'write', 'SVFOUT' // TRIM( coupling_char ) )
     10008!
     10009!--    Write global variables.
     10010       CALL wrd_mpi_io( 'rad_version', rad_version )
     10011!
     10012!--    Sum local number of skyview factor values on all PEs.
     10013       local_sum(1) = nsvfl
     10014       local_sum(2) = ncsfl
     10015       local_sum(3) = nsurfl
     10016       local_sum(4) = npcbl
     10017       local_sum(5) = ndsidir
     10018       local_sum(6) = nmrtbl
     10019       local_sum(7) = nmrtf
     10020#if defined( __parallel )
     10021       CALL MPI_ALLREDUCE( local_sum, global_sum, SIZE( local_sum ), MPI_INTEGER8, MPI_SUM,        &
     10022                           comm2d, ierr)
     10023#else
     10024       global_sum = local_sum
     10025#endif
     10026!
     10027!--    Check, if total number of respective skyview values do not exceed 2**31-1
     10028       DO  i = 1, 7
     10029          IF ( global_sum(i) > max_i4_value )  THEN
     10030             WRITE( message_string, '(A,I1,A)' ) 'number of sky view factor values (', i,          &
     10031                                                 ') > 2**31-1'
     10032             CALL message( 'radiation_write_svf', 'PA0733', 1, 2, 0, 6, 0 )
    962610033          ENDIF
    9627           IF ( nsvfl > 0 ) THEN
    9628              WRITE( 89 ) svf
    9629              WRITE( 89 ) svfsurf
     10034       ENDDO
     10035
     10036       nsvfl_tot   = global_sum(1)
     10037       ncsfl_tot   = global_sum(2)
     10038       nsurfl_tot  = global_sum(3)
     10039       npcbl_tot   = global_sum(4)
     10040       ndsidir_tot = global_sum(5)
     10041       nmrtbl_tot  = global_sum(6)
     10042       nmrtf_tot   = global_sum(7)
     10043!
     10044!--    Write total counters in header section of MPI-IO file.
     10045       CALL wrd_mpi_io( 'nsvfl', nsvfl_tot )
     10046       CALL wrd_mpi_io( 'ncsfl', ncsfl_tot )
     10047       CALL wrd_mpi_io( 'nsurfl', nsurfl_tot )
     10048       CALL wrd_mpi_io( 'npcbl', npcbl_tot )
     10049       CALL wrd_mpi_io( 'ndsidir', ndsidir_tot )
     10050       CALL wrd_mpi_io( 'nmrtbl', nmrtbl_tot )
     10051       CALL wrd_mpi_io( 'nmrtf', nmrtf_tot )
     10052!
     10053!--    Write local data.
     10054!--    All svf values are treated as surface values and use the respective routines from
     10055!--    restart_data_mpi_io_mod.
     10056       IF ( nsurfl > 0 )  THEN
     10057
     10058          lo_no = 0
     10059!
     10060!--       Count surface values on individual grid cells
     10061          DO  isurf = 1, nsurfl
     10062             jc = surfl(iy, isurf)
     10063             ic = surfl(ix, isurf)
     10064             lo_no(jc,ic) = lo_no(jc,ic) + 1
     10065          ENDDO
     10066!
     10067!--       Create local index array similar to surface routines.
     10068          ind = 1
     10069          DO  i = nxl, nxr
     10070             DO  j = nys, nyn
     10071                start_index(j,i) = ind
     10072                end_index(j,i)   = start_index(j,i) + lo_no(j,i) - 1
     10073                ind = ind + lo_no(j,i)
     10074             ENDDO
     10075          ENDDO
     10076
     10077       ELSE
     10078
     10079          start_index = 1
     10080          end_index = 0
     10081
     10082       ENDIF
     10083!
     10084!--    Each PE has to call the next block, therefore nsurfl_tot is used.
     10085!--    This is required to use MPI_FILE_WRITE_All for writing.
     10086       IF ( nsurfl_tot > 0 )  THEN
     10087!
     10088!--       Set file types of variables for this block and compute global indices.
     10089          CALL rd_mpi_io_surface_filetypes( start_index, end_index, data_to_write, global_start,   &
     10090                                            global_end )
     10091
     10092          CALL wrd_mpi_io( 'nsurfl_global_start', global_start )
     10093          CALL wrd_mpi_io( 'nsurfl_global_end', global_end )
     10094
     10095          IF ( data_to_write )  THEN
     10096
     10097             CALL wrd_mpi_io_surface( 'skyvf', skyvf )
     10098             CALL wrd_mpi_io_surface( 'skyvft', skyvft )
     10099!
     10100!--          To avoid another overlay of rrd_mpi_io_surface, dsitrans is written as REAL tmp array.
     10101!--          The order of dimensions of dsitrans is different to the order expected by
     10102!--          rrd_mpi_io_surface. Therefor a tranpose of tmp is required.
     10103             ALLOCATE( tmp(SIZE(dsitrans,2),SIZE(dsitrans,1)) )
     10104             tmp = TRANSPOSE( dsitrans )
     10105             CALL wrd_mpi_io_surface( 'dsitrans', tmp )
     10106             DEALLOCATE( tmp )
     10107
    963010108          ENDIF
    9631           IF ( plant_canopy  .AND.  ncsfl > 0 )  THEN
    9632               WRITE( 89 ) csf
    9633               WRITE( 89 ) csfsurf
     10109
     10110       ENDIF
     10111
     10112       IF ( npcbl > 0 )  THEN
     10113
     10114          lo_no = 0
     10115          DO  isurf = 1, npcbl
     10116             jc = pcbl(iy,isurf)
     10117             ic = pcbl(ix,isurf)
     10118             lo_no(jc,ic) = lo_no(jc,ic) + 1
     10119          ENDDO
     10120
     10121          ind = 1
     10122          DO  i = nxl, nxr
     10123             DO  j = nys, nyn
     10124                start_index(j,i) = ind
     10125                end_index(j,i)   = start_index(j,i) + lo_no(j,i) - 1
     10126                ind = ind+lo_no(j,i)
     10127             ENDDO
     10128          ENDDO
     10129
     10130       ELSE
     10131
     10132          start_index = 1
     10133          end_index = 0
     10134
     10135       ENDIF
     10136
     10137       IF ( npcbl_tot > 0 )  THEN
     10138
     10139          CALL rd_mpi_io_surface_filetypes( start_index, end_index, data_to_write, global_start,   &
     10140                                            global_end )
     10141          CALL wrd_mpi_io( 'npcbl_global_start', global_start )
     10142          CALL wrd_mpi_io( 'npcbl_global_end', global_end )
     10143
     10144          IF ( data_to_write )  THEN
     10145
     10146             ALLOCATE( tmp(SIZE(dsitransc,2),SIZE(dsitransc,1)) )
     10147             tmp = TRANSPOSE( dsitransc )
     10148             CALL wrd_mpi_io_surface( 'dsitransc', tmp )
     10149             DEALLOCATE( tmp )
    963410150          ENDIF
    9635           IF ( nmrtbl > 0 )  THEN
    9636              WRITE( 89 ) mrtsky
    9637              WRITE( 89 ) mrtskyt
    9638              WRITE( 89 ) mrtdsit
     10151       ENDIF
     10152
     10153       IF ( nsvfl > 0 )  THEN
     10154
     10155          lo_no = 0
     10156          DO  j = 1, SIZE( svfsurf, 2 )
     10157             isurf = svfsurf(1,j)
     10158             jc = surfl(iy,isurf)
     10159             ic = surfl(ix,isurf)
     10160             lo_no(jc,ic) = lo_no(jc,ic) + 1
     10161          ENDDO
     10162
     10163          ind = 1
     10164          DO  i = nxl, nxr
     10165             DO  j = nys, nyn
     10166                start_index(j,i) = ind
     10167                end_index(j,i)   = start_index(j,i) + lo_no(j,i) - 1
     10168                ind = ind+lo_no(j,i)
     10169             ENDDO
     10170          ENDDO
     10171
     10172       ELSE
     10173
     10174          start_index = 1
     10175          end_index = 0
     10176       ENDIF
     10177
     10178       IF ( nsvfl_tot > 0 )  THEN
     10179
     10180          CALL rd_mpi_io_surface_filetypes( start_index, end_index, data_to_write, global_start,   &
     10181                                            global_end )
     10182          CALL wrd_mpi_io( 'nsvfl_global_start', global_start )
     10183          CALL wrd_mpi_io( 'nsvfl_global_end', global_end )
     10184
     10185          IF ( data_to_write )  THEN
     10186
     10187             CALL wrd_mpi_io_surface( 'svf', svf )
     10188             ALLOCATE( tmp(SIZE(svfsurf,1),SIZE(svfsurf,2)) )
     10189             tmp(:,:) = svfsurf(:,:)
     10190             CALL wrd_mpi_io_surface( 'svfsurf', tmp )
     10191             DEALLOCATE( tmp )
     10192
    963910193          ENDIF
    9640           IF ( nmrtf > 0 )  THEN
    9641               WRITE( 89 ) mrtf
    9642               WRITE( 89 ) mrtft
    9643               WRITE( 89 ) mrtfsurf
     10194
     10195       ENDIF
     10196
     10197       IF ( plant_canopy  )  THEN
     10198
     10199          lo_no = 0
     10200          IF ( ncsfl > 0 )  THEN
     10201
     10202             DO  j = 1, ncsfl
     10203                ipcgb = csfsurf(1, j)
     10204                jc = pcbl(iy,ipcgb)
     10205                ic = pcbl(ix,ipcgb)
     10206                lo_no(jc,ic) = lo_no(jc,ic) + 1
     10207             ENDDO
     10208
     10209             ind = 1
     10210             DO  i = nxl, nxr
     10211                DO  j = nys, nyn
     10212                   start_index(j,i) = ind
     10213                   end_index(j,i)   = start_index(j,i) + lo_no(j,i) - 1
     10214                   ind = ind+lo_no(j,i)
     10215                ENDDO
     10216             ENDDO
     10217
     10218          ELSE
     10219
     10220             start_index = 1
     10221             end_index = -1
     10222
    964410223          ENDIF
    9645 !
    9646 !--       Close binary file
    9647           CALL close_file( 89 )
    9648 
    9649        ENDIF
    9650 #if defined( __parallel )
    9651        CALL MPI_BARRIER( comm2d, ierr )
    9652 #endif
    9653     ENDDO
    9654 
    9655     CALL location_message( 'writing view factors for radiation interaction', 'finished' )
    9656 
     10224
     10225          IF ( ncsfl_tot > 0 )  THEN
     10226
     10227             CALL rd_mpi_io_surface_filetypes( start_index, end_index, data_to_write, global_start,&
     10228                                               global_end )
     10229             CALL wrd_mpi_io( 'ncsfl_global_start', global_start )
     10230             CALL wrd_mpi_io( 'ncsfl_global_end', global_end )
     10231
     10232             IF ( data_to_write )  THEN
     10233
     10234                IF ( ALLOCATED( csf ) )  THEN
     10235                   CALL wrd_mpi_io_surface( 'csf', csf )
     10236                ELSE
     10237                   ALLOCATE( tmp(ndcsf,0) )
     10238                   CALL wrd_mpi_io_surface( 'csf', tmp )
     10239                   DEALLOCATE( tmp )
     10240                ENDIF
     10241
     10242                IF ( ALLOCATED( csfsurf ) )  THEN
     10243                   ALLOCATE( tmp(SIZE(csfsurf,1),SIZE(csfsurf,2)) )
     10244                   tmp(:,:) = csfsurf(:,:)
     10245                   CALL wrd_mpi_io_surface( 'csfsurf', tmp )
     10246                   DEALLOCATE( tmp )
     10247                ELSE
     10248                   ALLOCATE( tmp(idcsf,ncsfl) )
     10249                   CALL wrd_mpi_io_surface( 'csfsurf', tmp )
     10250                   DEALLOCATE( tmp )
     10251                ENDIF
     10252
     10253             ENDIF
     10254
     10255          ENDIF
     10256
     10257       ENDIF
     10258
     10259       IF ( nmrtbl > 0 )  THEN
     10260
     10261          lo_no = 0
     10262          DO  j = 1, nmrtbl
     10263             jc = mrtbl(iy,j)
     10264             ic = mrtbl(ix,j)
     10265             lo_no(jc,ic) = lo_no(jc,ic) + 1
     10266          ENDDO
     10267
     10268          ind = 1
     10269          DO  i = nxl, nxr
     10270             DO  j = nys, nyn
     10271                start_index(j,i) = ind
     10272                end_index(j,i)   = start_index(j,i) + lo_no(j,i) - 1
     10273                ind = ind+lo_no(j,i)
     10274             ENDDO
     10275          ENDDO
     10276
     10277       ELSE
     10278
     10279          start_index = 1
     10280          end_index = 0
     10281       ENDIF
     10282
     10283       IF ( nmrtbl_tot > 0 )  THEN
     10284
     10285          CALL rd_mpi_io_surface_filetypes( start_index, end_index, data_to_write, global_start,   &
     10286                                            global_end )
     10287          CALL wrd_mpi_io( 'nmrtbl_global_start', global_start )
     10288          CALL wrd_mpi_io( 'nmrtbl_global_end', global_end )
     10289
     10290          IF ( data_to_write )  THEN
     10291
     10292             CALL wrd_mpi_io_surface( 'mrtsky', mrtsky )
     10293             CALL wrd_mpi_io_surface( 'mrtskyt', mrtskyt )
     10294             ALLOCATE( tmp(SIZE(mrtdsit,2),SIZE(mrtdsit,1)) )
     10295             tmp = TRANSPOSE( mrtdsit )
     10296             CALL wrd_mpi_io_surface( 'mrtdsit', tmp )
     10297             DEALLOCATE( tmp )
     10298
     10299          ENDIF
     10300
     10301       ENDIF
     10302
     10303       IF ( nmrtf > 0 )  THEN
     10304
     10305          lo_no = 0
     10306          DO  j = 1, nmrtf
     10307             isurf = mrtfsurf(1,j)
     10308             jc = surfl(iy,isurf)
     10309             ic = surfl(ix,isurf)
     10310             lo_no(jc,ic) = lo_no(jc,ic) + 1
     10311          ENDDO
     10312
     10313          ind = 1
     10314          DO  i = nxl, nxr
     10315             DO  j = nys, nyn
     10316                start_index(j,i) = ind
     10317                end_index(j,i)   = start_index(j,i) + lo_no(j,i) - 1
     10318                ind = ind+lo_no(j,i)
     10319             ENDDO
     10320          ENDDO
     10321
     10322       ELSE
     10323          start_index = 1
     10324          end_index = 0
     10325       ENDIF
     10326
     10327       IF ( nmrtf_tot > 0 )  THEN
     10328
     10329          CALL rd_mpi_io_surface_filetypes( start_index, end_index, data_to_write, global_start,   &
     10330                                            global_end)
     10331          CALL wrd_mpi_io( 'nmrtf_global_start', global_start )
     10332          CALL wrd_mpi_io( 'nmrtf_global_end', global_end )
     10333
     10334          IF ( data_to_write )  THEN
     10335
     10336             CALL wrd_mpi_io_surface ( 'mrtf', mrtf )
     10337             CALL wrd_mpi_io_surface ( 'mrtft', mrtft )
     10338
     10339             ALLOCATE( tmp(SIZE(mrtfsurf,1),SIZE(mrtfsurf,2)) )
     10340             tmp = mrtfsurf
     10341             CALL wrd_mpi_io_surface( 'mrtfsurf', tmp )
     10342             DEALLOCATE( tmp )
     10343
     10344          ENDIF
     10345
     10346       ENDIF
     10347
     10348       CALL rd_mpi_io_close
     10349
     10350       CALL location_message( 'writing sky view factors for radiation interaction ' //             &
     10351                              '(MPI-IO format)', 'finished' )
     10352    ENDIF
    965710353
    965810354 END SUBROUTINE radiation_write_svf
Note: See TracChangeset for help on using the changeset viewer.