Changeset 440 for IOIPSL/trunk/src/mathelp.f90

Timestamp:

11/26/08 11:58:38 (15 years ago)

Author:

bellier

Message:

Removing trans_buff

File:

• IOIPSL/trunk/src/mathelp.f90 (modified) (2 diffs)

IOIPSL/trunk/src/mathelp.f90

@@ -10 +10 @@
 !-
   PRIVATE
-  PUBLIC :: mathop,moycum,trans_buff,buildop
+  PUBLIC :: mathop,moycum,buildop
 !-
   INTERFACE mathop
@@ -3113 +3113 @@
 END SUBROUTINE moycum
 !===
-SUBROUTINE trans_buff (ox,sx,oy,sy,oz,sz,xsz,ysz,zsz,v3d,sl,v1d)
-!---------------------------------------------------------------------
-!- This subroutine extracts from the full 3D variable the slab of
-!- data that will be used later. Perhaps there are hardware routines
-!- for this task on some computers. This routine will be obsolete in
-!- a F90 environnement
-!-
-!-sc   V1d = reshape(V3id (ox:ox-1+sx,oy:oy-1+sy,oz:oz-1+sz),
-!-sc                 SHAPE= (/sx*sy*sz/) )
-!-
-!- INPUT
-!- ox  : Origin of slab of data in X
-!- sx  : Size of slab in X
-!- oy  : Origin of slab of data in Y
-!- sy  : Size of slab in Y
-!- oz  : Origin of slab of data in Z
-!- sz  : Size of slab in Z
-!- xsz,ysz,zsz : 3 sizes of full variable v3d
-!- v3d : The full 3D variable
-!- sl  : size of variable v1d
-!- v1d : The 1D variable containing the slab
-!-
-!- VERSION
-!-
-!---------------------------------------------------------------------
-  IMPLICIT NONE
-!-
-  INTEGER :: ox,sx,oy,sy,oz,sz
-  INTEGER :: xsz,ysz,zsz
-  INTEGER :: sl
-  REAL :: v3d(xsz,ysz,zsz)
-  REAL :: v1d(sl)
-!-
-  INTEGER :: ix,iy,iz,ic
-!---------------------------------------------------------------------
-!-
-! We have to consider the case where the zoom starts at (1,1,1)
-! but does not go over the full size.
-!-
-  IF (     (ox .EQ. 1).AND.(oy.EQ. 1).AND.(oz.EQ. 1) &
- &    .AND.(sx .EQ. xsz).AND.(sy .EQ. ysz).AND.(sz .EQ. zsz) ) THEN 
-    DO ic=1,MAX(sx,1)*MAX(sy,1)*MAX(sz,1) 
-      v1d(ic) = v3d(ic,1,1)
-    ENDDO 
-  ELSE IF (     (ox .EQ. 1).AND.(oy .EQ. 1) &
- &         .AND.(sx .EQ. xsz).AND.(sy .EQ. ysz)) THEN 
-    DO iz=oz,(oz-1+sz) 
-      DO ic=1,MAX(sx,1)*MAX(sy,1) 
-        v1d(ic) = v3d(ic,1,iz)
-      ENDDO
-    ENDDO
-  ELSE 
-    ic = 0
-    DO iz=oz,(oz-1+sz)
-      DO iy=oy,(oy-1+sy)
-        DO ix=ox,(ox-1+sx)
-          ic = ic+1
-          v1d(ic) = v3d(ix, iy, iz)
-        ENDDO
-      ENDDO
-    ENDDO
-  ENDIF
-!------------------------
-END SUBROUTINE trans_buff
-!===
 !-----------------
 END MODULE mathelp