Changeset 1053 for codes/icosagcm/trunk/src/parallel

Timestamp:

09/24/20 18:17:18 (4 years ago)

Author:

dubos

Message:

trunk : simplify allocate_field -- tested on 4 GPUs (TBC)

File:

• codes/icosagcm/trunk/src/parallel/transfert_mpi.f90 (modified) (1 diff)

codes/icosagcm/trunk/src/parallel/transfert_mpi.f90

@@ -135 +135 @@
     if( request(1)%field_type /= field(1)%field_type ) call dynamico_abort( "init_message : field_type/request mismatch" )
     field_type = request(1)%field_type
-
-    ! Set field%rval4d pointer to always use 4d array
-    do ind = 1, ndomain
-      if( field(ind)%ndim == 2 ) field(ind)%rval4d(1:size(field(ind)%rval2d,1),1:1,1:1) => field(ind)%rval2d
-      ! This is Fortran 2008 : can be avoided by using a subroutine with rval3d as a 1D dummy argument
-      ! (/!\ : using a subroutine might generate a temporary contiguous array)
-      if( field(ind)%ndim == 3 ) field(ind)%rval4d(1:size(field(ind)%rval3d,1), &
-        1:size(field(ind)%rval3d,2), 1:1) => field(ind)%rval3d
-    end do
+    
     dim3 = size(field(1)%rval4d,2)
     dim4 = size(field(1)%rval4d,3)