Changeset 4735

Timestamp:

Oct 13, 2020 7:00:22 AM (4 years ago)

Author:

raasch

Message:

interface for 3d logical arrays added

File:

• palm/trunk/SOURCE/restart_data_mpi_io_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/restart_data_mpi_io_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! interface for 3d logical arrays added
+!
+! 4694 2020-09-23 15:09:19Z pavelkrc
 ! Fix name of subroutine in message
 !
@@ -334 +337 @@
        MODULE PROCEDURE rrd_mpi_io_int8_3d
        MODULE PROCEDURE rrd_mpi_io_logical
+       MODULE PROCEDURE rrd_mpi_io_logical_3d
        MODULE PROCEDURE rrd_mpi_io_real
        MODULE PROCEDURE rrd_mpi_io_real_2d
@@ -372 +376 @@
        MODULE PROCEDURE wrd_mpi_io_int8_3d
        MODULE PROCEDURE wrd_mpi_io_logical
+       MODULE PROCEDURE wrd_mpi_io_logical_3d
        MODULE PROCEDURE wrd_mpi_io_real
        MODULE PROCEDURE wrd_mpi_io_real_2d
@@ -1939 +1944 @@
 ! Description:
 ! ------------
+!> Read 3D-LOGICAL with MPI-IO
+!--------------------------------------------------------------------------------------------------!
+ SUBROUTINE rrd_mpi_io_logical_3d( name, data )
+
+    IMPLICIT NONE
+
+    CHARACTER(LEN=*), INTENT(IN) ::  name                !<
+
+    INTEGER(iwp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  logical_as_integer  !<
+
+    LOGICAL, INTENT(OUT), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  data  !<
+
+
+    CALL rrd_mpi_io_int4_3d( name, logical_as_integer )
+    data(:,:,:) = ( logical_as_integer(:,:,:) == 1 )
+
+ END SUBROUTINE rrd_mpi_io_logical_3d
+
+
+
+!--------------------------------------------------------------------------------------------------!
+! Description:
+! ------------
 !> Write INTEGER with MPI-IO
 !--------------------------------------------------------------------------------------------------!
@@ -2462 +2490 @@
 
  END SUBROUTINE wrd_mpi_io_logical
+
+
+
+!--------------------------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Write 3D-LOGICAL with MPI-IO
+!--------------------------------------------------------------------------------------------------!
+ SUBROUTINE wrd_mpi_io_logical_3d( name, data )
+
+    IMPLICIT NONE
+
+    CHARACTER(LEN=*), INTENT(IN) ::  name                !<
+
+    INTEGER(iwp) ::  i  !< loop index
+    INTEGER(iwp) ::  j  !< loop index
+    INTEGER(iwp) ::  k  !< loop index
+
+    INTEGER(iwp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  logical_as_integer  !<
+
+    LOGICAL, INTENT(IN), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  data   !<
+
+
+    DO  i = nxlg, nxrg
+       DO  j = nysg, nyng
+          DO  k = nzb, nzt+1
+             IF ( data(k,j,i) )  THEN
+                logical_as_integer(k,j,i) = 1
+             ELSE
+                logical_as_integer(k,j,i) = 0
+             ENDIF
+          ENDDO
+       ENDDO
+    ENDDO
+
+    CALL wrd_mpi_io_int4_3d( name, logical_as_integer )
+
+ END SUBROUTINE wrd_mpi_io_logical_3d