C--   File plot_field.F: Routines for "formatted" I/O in the MITgcm UV 
C--                      implementation.
C--    Contents
C--    o plot_field_xyr8 - Writes a XY Real*8 field
      SUBROUTINE PLOT_FIELD_XYR8( 
     I                            fld, fldNam)

C     /==========================================================\
C     | SUBROUTINE PLOT_FIELD_XYR8                               |
C     | Print out an XY Real 8 field using text map.             |
C     |==========================================================|
C     | This routine references "numerical model" parameters like|
C     | like the integration time. It uses these to create a     |
C     | title for the field before calling a generic execution   |
C     | environment support routine.                             |
C     | This routine can also be edited to cause only some region|
C     | of a field to be printed by default, or every other      |
C     | point etc..                                              |
C     | Other plot formats can also be substituted here.         |
C     \==========================================================/
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "CG2D.h"

C     == Routine arguments ==
      REAL fld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      CHARACTER*(*) fldNam
      INTEGER myThid

C     == Local variables ==
      CHARACTER*(MAX_LEN_MBUF) fldTitle
      INTEGER iStart, iEnd, iStride
      INTEGER jStart, jEnd, jStride
      INTEGER kStart, kEnd, kStride
      INTEGER biStart, biEnd, biStride
      INTEGER bjStart, bjEnd, bjStride

C--    Form name for identifying "plot"
       WRITE(fldTitle,'(A,A,A,I4)') ' Field ', fldNam, ' at iteration ',
     & nIter
C--    Do "plot" using textual contour map "execution environment" routine
C      Substitute other plotting utilities here!
       iStart   =  1-OLx
       iEnd     =  sNx+OLx
       iStride  =  1
       jStart   =  sNy+OLy
       jEnd     =  1-OLy
       jStride  = -1
       kStart   =  1
       kEnd     =  1
       kStride  =  1
       biStart  =  1
       biEnd    =  1
       biStride =  1
       bjStart  =  1
       bjEnd    =  1    
       bjStride = -1
       CALL PRINT_MAPR8(
     I        fld, fldTitle,
     I         1-OLx,sNx+OLx,1-OLy,sNy+OLy,1,1,  1, 1,
     I         iStart,   iEnd,  iStride,
     I         jStart,   jEnd,  jStride,
     I         kStart,   kEnd,  kStride,
     I        biStart,  biEnd, biStride,
     I        bjStart,  bjEnd, bjStride )

      RETURN
      END