cfc_example/code/MDSIO_BUFF_3D.h

C $Header: /u/gcmpack/MITgcm/pkg/mdsio/MDSIO_BUFF_3D.h,v 1.3 2010/12/23 02:42:16 jmc Exp $
C $Name:  $

CBOP
C     !ROUTINE: MDSIO_BUFF_3D.h
C     !INTERFACE:
C     include "MDSIO_BUFF_3D.h"
C
C     !DESCRIPTION:
C     *==========================================================*
C     | MDSIO_BUFF_3D.h
C     | o Shared 3-D Buffers used for I/O
C     *==========================================================*
CEOP

C     size3dBuf  :: buffer 3rd dimension, corresponds to the maximum number
C                   of levels that can be read/written at a time.
C     Note: minimum value = Nr, but in few cases (vertical interpolation,
C           NrPhys from Fizhi, ...)  needs to be larger. Here we pick 2*Nr
C           which should be enough for most applications.
      INTEGER size3dBuf
#ifdef ALLOW_FIZHI
      PARAMETER ( size3dBuf = Nr+NrPhys )
#else
c     PARAMETER ( size3dBuf = 2*Nr )
C     increase buffer size (below) for some pkg/layers output:
      PARAMETER ( size3dBuf = 3*Nr )
#endif

C--   COMMON /MDS_3D_BUFFERS/  3-D Shared Local Buffers
C     Those buffers have be in common block to be shared by all threads;
C     considered to be "owned" by master-thread and any access by other 
C     than master thread needs to be put protected by BARRIER.
C shared3dBuf_rx :: Heap storage buffer to which master thread
C                   read-in/write-from data which all threads copy from
C                   (during read) or copy to (during write).
      COMMON /MDS_3D_BUFFERS/ shared3dBuf_r8, shared3dBuf_r4
      Real*8 shared3dBuf_r8( sNx*sNy*size3dBuf*nSx*nSy )
      Real*4 shared3dBuf_r4( sNx*sNy*size3dBuf*nSx*nSy )

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
1	jmc	1.1	C $Header: /u/gcmpack/MITgcm/pkg/mdsio/MDSIO_BUFF_3D.h,v 1.3 2010/12/23 02:42:16 jmc Exp $
2			C $Name: $
3
4			CBOP
5			C !ROUTINE: MDSIO_BUFF_3D.h
6			C !INTERFACE:
7			C include "MDSIO_BUFF_3D.h"
8			C
9			C !DESCRIPTION:
10			C ==========================================================
11			C \| MDSIO_BUFF_3D.h
12			C \| o Shared 3-D Buffers used for I/O
13			C ==========================================================
14			CEOP
15
16			C size3dBuf :: buffer 3rd dimension, corresponds to the maximum number
17			C of levels that can be read/written at a time.
18			C Note: minimum value = Nr, but in few cases (vertical interpolation,
19			C NrPhys from Fizhi, ...) needs to be larger. Here we pick 2*Nr
20			C which should be enough for most applications.
21			INTEGER size3dBuf
22			#ifdef ALLOW_FIZHI
23			PARAMETER ( size3dBuf = Nr+NrPhys )
24			#else
25			c PARAMETER ( size3dBuf = 2*Nr )
26			C increase buffer size (below) for some pkg/layers output:
27			PARAMETER ( size3dBuf = 3*Nr )
28			#endif
29
30			C-- COMMON /MDS_3D_BUFFERS/ 3-D Shared Local Buffers
31			C Those buffers have be in common block to be shared by all threads;
32			C considered to be "owned" by master-thread and any access by other
33			C than master thread needs to be put protected by BARRIER.
34			C shared3dBuf_rx :: Heap storage buffer to which master thread
35			C read-in/write-from data which all threads copy from
36			C (during read) or copy to (during write).
37			COMMON /MDS_3D_BUFFERS/ shared3dBuf_r8, shared3dBuf_r4
38			Real8 shared3dBuf_r8( sNxsNysize3dBufnSx*nSy )
39			Real4 shared3dBuf_r4( sNxsNysize3dBufnSx*nSy )
40
41			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|