Changeset 688 for codes/icosagcm/devel/src/unstructured/transfer_unstructured.F90

Timestamp:

04/09/18 15:22:10 (6 years ago)

Author:

dubos

Message:

devel/unstructured : select double or single precision for physical quantities

File:

• codes/icosagcm/devel/src/unstructured/transfer_unstructured.F90 (modified) (7 diffs)

codes/icosagcm/devel/src/unstructured/transfer_unstructured.F90

@@ -4 +4 @@
   SAVE
   PRIVATE
+
+#include "unstructured.h90"
 
   INCLUDE 'mpif.h'
@@ -47 +49 @@
   SUBROUTINE copy_in(lldata, ndata, llbuf, nbuf, cells, data, buf)
     INTEGER :: ndata, lldata, nbuf, llbuf, cells(nbuf), i
-    REAL(C_DOUBLE) :: data(lldata,ndata), buf(llbuf,nbuf)
+    NUM :: data(lldata,ndata), buf(llbuf,nbuf)
     !$OMP DO
     DO i=1,nbuf
@@ -57 +59 @@
   SUBROUTINE copy_out(lldata, ndata, llbuf, nbuf, cells, data, buf)
     INTEGER :: ndata, lldata, nbuf, llbuf, cells(nbuf), i
-    REAL(C_DOUBLE) :: data(lldata,ndata), buf(llbuf,nbuf)
+    NUM :: data(lldata,ndata), buf(llbuf,nbuf)
     !$OMP DO SCHEDULE(STATIC)
     DO i=1,nbuf
@@ -67 +69 @@
   SUBROUTINE update_halo(index, data)
     INTEGER :: index ! transfer_XXX, cf data_unstructured.F90
-    REAL(C_DOUBLE) :: data(:,:)    
+    NUM :: data(:,:)    
     TYPE(Halo_transfer), POINTER :: halo
-    REAL(C_DOUBLE), POINTER :: buf(:,:)
+    NUM, POINTER :: buf(:,:)
     INTEGER :: req(send_info(index)%ranks+recv_info(index)%ranks), &
          status(MPI_STATUS_SIZE, send_info(index)%ranks+recv_info(index)%ranks)
@@ -85 +87 @@
        num = halo%num(i)
        ireq = ireq+1
-       call MPI_ISEND(buf(1,istart), llbuf*num, MPI_REAL8, halo%rank(i), &
+       call MPI_ISEND(buf(1,istart), llbuf*num, MPI_NUM, halo%rank(i), &
             42, MPI_COMM_WORLD, req(ireq), ierr)
        istart = istart+num
@@ -98 +100 @@
        num = halo%num(i)
        ireq = ireq+1
-       call MPI_IRECV(buf(1,istart), llbuf*num, MPI_REAL8, halo%rank(i), &
+       call MPI_IRECV(buf(1,istart), llbuf*num, MPI_NUM, halo%rank(i), &
             42, MPI_COMM_WORLD, req(ireq), ierr)
        istart = istart+num
@@ -111 +113 @@
   SUBROUTINE update_halo_c(index, lldata, ndata, data) BIND(C, name='dynamico_update_halo')
     INTEGER, VALUE :: index, lldata, ndata
-    REAL(C_DOUBLE) :: data(lldata, ndata)
+    NUM :: data(lldata, ndata)
     CALL update_halo(index, data)
   END SUBROUTINE update_halo_c