C $Header: /u/gcmpack/MITgcm/pkg/mdsio/mdsio_rw_slice.F,v 1.1 2009/09/01 19:19:10 jmc Exp $
C $Name:  $

#include "MDSIO_OPTIONS.h"

C--  File mdsio_rw_slice.F: old version of MDS_READ/WRITE_SEC_XZ/YZ S/R with
C    fewer arguments (kept for backward compatibility): call new MDSIO S/R
C    with fixed additional arguments
C--   Contents
C--   o MDSREADFIELDXZ
C--   o MDSREADFIELDYZ
C--   o MDSREADFIELDXZ_LOC
C--   o MDSREADFIELDYZ_LOC
C--   o MDSWRITEFIELDXZ
C--   o MDSWRITEFIELDYZ
C--   o MDSWRITEFIELDXZ_LOC
C--   o MDSWRITEFIELDYZ_LOC

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

      SUBROUTINE MDSREADFIELDXZ(
     I   fName,
     I   filePrec,
     I   arrType,
     I   nNz,
     |   arr,
     I   irecord,
     I   myThid )

C Arguments:
C
C fName      string  :: base name for file to written
C filePrec   integer :: number of bits per word in file (32 or 64)
C arrType    char(2) :: declaration of "arr": either "RS" or "RL"
C nNz        integer :: size of third dimension: normally either 1 or Nr
C arr         RS/RL  :: array to write, arr(:,:,nNz,:,:)
C irecord    integer :: record number to read
C myThid     integer :: thread identifier
C
C Routine now calls MDS_READ_SEC_XZ, just a way to add extra arguments
C to the argument list.
C The 1rst new argument (useCurrentDir=.FALSE.) allows to read files from
C the "mdsioLocalDir" directory (if it is set).
C The 2nd new argument avoid argument array of undefined type (RL/RS).

      IMPLICIT NONE
C Global variables / COMMON blocks
#include "SIZE.h"
c #include "EEPARAMS.h"

C Routine arguments
      CHARACTER*(*) fName
      INTEGER filePrec
      CHARACTER*(2) arrType
      INTEGER nNz
      _RL     arr(*)
      INTEGER irecord
      INTEGER myThid

#ifdef ALLOW_AUTODIFF
C Local variables
      _RL dummyRL(1)
      _RS dummyRS(1)

      IF ( arrType.EQ.'RL' ) THEN
        CALL MDS_READ_SEC_XZ(
     I                fName, filePrec, .FALSE., arrType, nNz,
     O                arr, dummyRS,
     I                irecord, myThid )
      ELSE
        CALL MDS_READ_SEC_XZ(
     I                fName, filePrec, .FALSE., arrType, nNz,
     O                dummyRL, arr,
     I                irecord, myThid )
      ENDIF

#else /* ALLOW_AUTODIFF */
      STOP 'ABNORMAL END: S/R MDSREADFIELDXZ is retired'
#endif /* ALLOW_AUTODIFF */

      RETURN
      END


C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| SUBROUTINE MDSREADFIELDYZ( I fName, I filePrec, I arrType, I nNz, | arr, I irecord, I myThid ) C Arguments: C C fName string :: base name for file to written C filePrec integer :: number of bits per word in file (32 or 64) C arrType char(2) :: declaration of "arr": either "RS" or "RL" C nNz integer :: size of third dimension: normally either 1 or Nr C arr RS/RL :: array to write, arr(:,:,nNz,:,:) C irecord integer :: record number to read C myThid integer :: thread identifier C C Routine now calls MDS_READ_SEC_YZ, just a way to add extra arguments C to the argument list. C The 1rst new argument (useCurrentDir=.FALSE.) allows to read files from C the "mdsioLocalDir" directory (if it is set). C The 2nd new argument avoid argument array of undefined type (RL/RS). IMPLICIT NONE C Global variables / COMMON blocks #include "SIZE.h" c #include "EEPARAMS.h" C Routine arguments CHARACTER*(*) fName INTEGER filePrec CHARACTER*(2) arrType INTEGER nNz _RL arr(*) INTEGER irecord INTEGER myThid #ifdef ALLOW_AUTODIFF C Local variables _RL dummyRL(1) _RS dummyRS(1) IF ( arrType.EQ.'RL' ) THEN CALL MDS_READ_SEC_YZ( I fName, filePrec, .FALSE., arrType, nNz, O arr, dummyRS, I irecord, myThid ) ELSE CALL MDS_READ_SEC_YZ( I fName, filePrec, .FALSE., arrType, nNz, O dummyRL, arr, I irecord, myThid ) ENDIF #else /* ALLOW_AUTODIFF */ STOP 'ABNORMAL END: S/R MDSREADFIELDYZ is retired' #endif /* ALLOW_AUTODIFF */ RETURN END


C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| SUBROUTINE MDSREADFIELDXZ_LOC( I fName, I filePrec, I arrType, I nNz, | arr, I irecord, I myThid ) C Arguments: C C fName string :: base name for file to written C filePrec integer :: number of bits per word in file (32 or 64) C arrType char(2) :: declaration of "arr": either "RS" or "RL" C nNz integer :: size of third dimension: normally either 1 or Nr C arr RS/RL :: array to write, arr(:,:,nNz,:,:) C irecord integer :: record number to read C myThid integer :: thread identifier C C Routine now calls MDS_READ_SEC_XZ, just a way to add extra arguments C to the argument list. C The 1rst new argument (useCurrentDir=.TRUE.) forces to ignore the C "mdsioLocalDir" parameter and to always write to the current directory. C The 2nd new argument avoid argument array of undefined type (RL/RS). IMPLICIT NONE C Global variables / COMMON blocks #include "SIZE.h" c #include "EEPARAMS.h" C Routine arguments CHARACTER*(*) fName INTEGER filePrec CHARACTER*(2) arrType INTEGER nNz _RL arr(*) INTEGER irecord INTEGER myThid #ifdef ALLOW_AUTODIFF C Local variables _RL dummyRL(1) _RS dummyRS(1) IF ( arrType.EQ.'RL' ) THEN CALL MDS_READ_SEC_XZ( I fName, filePrec, .TRUE., arrType, nNz, O arr, dummyRS, I irecord, myThid ) ELSE CALL MDS_READ_SEC_XZ( I fName, filePrec, .TRUE., arrType, nNz, O dummyRL, arr, I irecord, myThid ) ENDIF #else /* ALLOW_AUTODIFF */ STOP 'ABNORMAL END: S/R MDSREADFIELDXZ_LOC is empty' #endif /* ALLOW_AUTODIFF */ RETURN END


C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| SUBROUTINE MDSREADFIELDYZ_LOC( I fName, I filePrec, I arrType, I nNz, | arr, I irecord, I myThid ) C Arguments: C C fName string :: base name for file to written C filePrec integer :: number of bits per word in file (32 or 64) C arrType char(2) :: declaration of "arr": either "RS" or "RL" C nNz integer :: size of third dimension: normally either 1 or Nr C arr RS/RL :: array to write, arr(:,:,nNz,:,:) C irecord integer :: record number to read C myThid integer :: thread identifier C C Routine now calls MDS_READ_SEC_YZ, just a way to add extra arguments C to the argument list. C The 1rst new argument (useCurrentDir=.TRUE.) forces to ignore the C "mdsioLocalDir" parameter and to always write to the current directory. C The 2nd new argument avoid argument array of undefined type (RL/RS). IMPLICIT NONE C Global variables / COMMON blocks #include "SIZE.h" c #include "EEPARAMS.h" C Routine arguments CHARACTER*(*) fName INTEGER filePrec CHARACTER*(2) arrType INTEGER nNz _RL arr(*) INTEGER irecord INTEGER myThid #ifdef ALLOW_AUTODIFF C Local variables _RL dummyRL(1) _RS dummyRS(1) IF ( arrType.EQ.'RL' ) THEN CALL MDS_READ_SEC_YZ( I fName, filePrec, .TRUE., arrType, nNz, O arr, dummyRS, I irecord, myThid ) ELSE CALL MDS_READ_SEC_YZ( I fName, filePrec, .TRUE., arrType, nNz, O dummyRL, arr, I irecord, myThid ) ENDIF #else /* ALLOW_AUTODIFF */ STOP 'ABNORMAL END: S/R MDSREADFIELDYZ_LOC is empty' #endif /* ALLOW_AUTODIFF */ RETURN END


C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| SUBROUTINE MDSWRITEFIELDXZ( I fName, I filePrec, I globalFile, I arrType, I nNz, I arr, I irecord, I myIter, I myThid ) C Arguments: C C fName string :: base name for file to write C filePrec integer :: number of bits per word in file (32 or 64) C globalFile logical :: selects between writing a global or tiled file C arrType char(2) :: declaration of "arr": either "RS" or "RL" C nNz integer :: size of third dimension: normally either 1 or Nr C arr RS/RL :: array to write, arr(:,:,nNzdim,:,:) C irecord integer :: record number to write C myIter integer :: time step number C myThid integer :: thread identifier C C Routine now calls MDS_WRITE_REC_XZ, just a way to add extra arguments C to the argument list. C The 1rst new argument (useCurrentDir=.FALSE.) allows to write files to C the "mdsioLocalDir" directory (if it is set). C The 2nd new argument avoid argument array of undefined type (RL/RS). IMPLICIT NONE C Global variables / common blocks #include "SIZE.h" c #include "EEPARAMS.h" C Routine arguments CHARACTER*(*) fName INTEGER filePrec LOGICAL globalFile CHARACTER*(2) arrType INTEGER nNz _RL arr(*) INTEGER irecord INTEGER myIter INTEGER myThid #ifdef ALLOW_AUTODIFF C Local variables _RL dummyRL(1) _RS dummyRS(1) IF ( arrType.EQ.'RL' ) THEN CALL MDS_WRITE_SEC_XZ( I fName, filePrec, globalFile, .FALSE., I arrType, nNz, arr, dummyRS, I irecord, myIter, myThid ) ELSE CALL MDS_WRITE_SEC_XZ( I fName, filePrec, globalFile, .FALSE., I arrType, nNz, dummyRL, arr, I irecord, myIter, myThid ) ENDIF #else /* ALLOW_AUTODIFF */ STOP 'ABNORMAL END: S/R MDSWRITEFIELDXZ is retired' #endif /* ALLOW_AUTODIFF */ RETURN END


C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| SUBROUTINE MDSWRITEFIELDYZ( I fName, I filePrec, I globalFile, I arrType, I nNz, I arr, I irecord, I myIter, I myThid ) C Arguments: C C fName string :: base name for file to write C filePrec integer :: number of bits per word in file (32 or 64) C globalFile logical :: selects between writing a global or tiled file C arrType char(2) :: declaration of "arr": either "RS" or "RL" C nNz integer :: size of third dimension: normally either 1 or Nr C arr RS/RL :: array to write, arr(:,:,nNzdim,:,:) C irecord integer :: record number to write C myIter integer :: time step number C myThid integer :: thread identifier C C Routine now calls MDS_WRITE_REC_YZ, just a way to add extra arguments C to the argument list. C The 1rst new argument (useCurrentDir=.FALSE.) allows to write files to C the "mdsioLocalDir" directory (if it is set). C The 2nd new argument avoid argument array of undefined type (RL/RS). IMPLICIT NONE C Global variables / common blocks #include "SIZE.h" c #include "EEPARAMS.h" C Routine arguments CHARACTER*(*) fName INTEGER filePrec LOGICAL globalFile CHARACTER*(2) arrType INTEGER nNz _RL arr(*) INTEGER irecord INTEGER myIter INTEGER myThid #ifdef ALLOW_AUTODIFF C Local variables _RL dummyRL(1) _RS dummyRS(1) IF ( arrType.EQ.'RL' ) THEN CALL MDS_WRITE_SEC_YZ( I fName, filePrec, globalFile, .FALSE., I arrType, nNz, arr, dummyRS, I irecord, myIter, myThid ) ELSE CALL MDS_WRITE_SEC_YZ( I fName, filePrec, globalFile, .FALSE., I arrType, nNz, dummyRL, arr, I irecord, myIter, myThid ) ENDIF #else /* ALLOW_AUTODIFF */ STOP 'ABNORMAL END: S/R MDSWRITEFIELDYZ is retired' #endif /* ALLOW_AUTODIFF */ RETURN END


C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| SUBROUTINE MDSWRITEFIELDXZ_LOC( I fName, I filePrec, I globalFile, I arrType, I nNz, I arr, I irecord, I myIter, I myThid ) C Arguments: C C fName string :: base name for file to write C filePrec integer :: number of bits per word in file (32 or 64) C globalFile logical :: selects between writing a global or tiled file C arrType char(2) :: declaration of "arr": either "RS" or "RL" C nNz integer :: size of third dimension: normally either 1 or Nr C arr RS/RL :: array to write, arr(:,:,nNzdim,:,:) C irecord integer :: record number to write C myIter integer :: time step number C myThid integer :: thread identifier C C Routine now calls MDS_WRITE_REC_XZ, just a way to add extra arguments C to the argument list. C The 1rst new argument (useCurrentDir=.TRUE.) forces to ignore the C "mdsioLocalDir" parameter and to always write to the current directory. C The 2nd new argument avoid argument array of undefined type (RL/RS). IMPLICIT NONE C Global variables / common blocks #include "SIZE.h" c #include "EEPARAMS.h" C Routine arguments CHARACTER*(*) fName INTEGER filePrec LOGICAL globalFile CHARACTER*(2) arrType INTEGER nNz _RL arr(*) INTEGER irecord INTEGER myIter INTEGER myThid #ifdef ALLOW_AUTODIFF C Local variables _RL dummyRL(1) _RS dummyRS(1) IF ( arrType.EQ.'RL' ) THEN CALL MDS_WRITE_SEC_XZ( I fName, filePrec, globalFile, .TRUE., I arrType, nNz, arr, dummyRS, I irecord, myIter, myThid ) ELSE CALL MDS_WRITE_SEC_XZ( I fName, filePrec, globalFile, .TRUE., I arrType, nNz, dummyRL, arr, I irecord, myIter, myThid ) ENDIF #else /* ALLOW_AUTODIFF */ STOP 'ABNORMAL END: S/R MDSWRITEFIELDXZ_LOC is empty' #endif /* ALLOW_AUTODIFF */ RETURN END


C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| SUBROUTINE MDSWRITEFIELDYZ_LOC( I fName, I filePrec, I globalFile, I arrType, I nNz, I arr, I irecord, I myIter, I myThid ) C Arguments: C C fName string :: base name for file to write C filePrec integer :: number of bits per word in file (32 or 64) C globalFile logical :: selects between writing a global or tiled file C arrType char(2) :: declaration of "arr": either "RS" or "RL" C nNz integer :: size of third dimension: normally either 1 or Nr C arr RS/RL :: array to write, arr(:,:,nNzdim,:,:) C irecord integer :: record number to write C myIter integer :: time step number C myThid integer :: thread identifier C C Routine now calls MDS_WRITE_REC_YZ, just a way to add extra arguments C to the argument list. C The 1rst new argument (useCurrentDir=.TRUE.) forces to ignore the C "mdsioLocalDir" parameter and to always write to the current directory. C The 2nd new argument avoid argument array of undefined type (RL/RS). IMPLICIT NONE C Global variables / common blocks #include "SIZE.h" c #include "EEPARAMS.h" C Routine arguments CHARACTER*(*) fName INTEGER filePrec LOGICAL globalFile CHARACTER*(2) arrType INTEGER nNz _RL arr(*) INTEGER irecord INTEGER myIter INTEGER myThid #ifdef ALLOW_AUTODIFF C Local variables _RL dummyRL(1) _RS dummyRS(1) IF ( arrType.EQ.'RL' ) THEN CALL MDS_WRITE_SEC_YZ( I fName, filePrec, globalFile, .TRUE., I arrType, nNz, arr, dummyRS, I irecord, myIter, myThid ) ELSE CALL MDS_WRITE_SEC_YZ( I fName, filePrec, globalFile, .TRUE., I arrType, nNz, dummyRL, arr, I irecord, myIter, myThid ) ENDIF #else /* ALLOW_AUTODIFF */ STOP 'ABNORMAL END: S/R MDSWRITEFIELDYZ_LOC is empty' #endif /* ALLOW_AUTODIFF */ RETURN END