PROGRAM find_max_eeta
      USE Rcm_mod_subs
      IMPLICIT NONE
      CHARACTER (LEN=1) :: prefix,answer
      CHARACTER (LEN=10) :: time_string
      CHARACTER (LEN=62) :: info_string
      CHARACTER (LEN=05) :: itrack_string
      CHARACTER (LEN=80) :: dev_name, dev_type
      CHARACTER (LEN=6) :: ut_string
!     
      REAL :: x_bnd(jsize), y_bnd(jsize), eeta_to_plot(isize,jsize), &
              xgsm_max, xgsm_min, ygsm_max, ygsm_min, xgsm_step, ygsm_step, xstep_tmp, ystep_tmp,&
              bright, contra, a1, a2, av, xpoly(4), ypoly(4), x_wedge(2), y_wedge(2), &
              wedge(1000), tr(6), x_plot(2), y_plot(2), x_win(2), y_win(2)
!
      INTEGER :: kc_to_plot,k_to_plot,k_array(ksize),k_count,kc_beg,kc_end , c1, c2, kc
!
      INTEGER :: n, k, i,  pgopen, npts, j, n_rec, pgid, minind, maxind, iv
      LOGICAL :: Do_plot_gridlines = .FALSE.
      INTEGER :: n_contours, color_use
      REAL :: eeta_min, eeta_max(kcsize)
      REAL :: cont_level, cont_step
!
!
      CALL Read_grid ()
      CALL Read_plasma ()
!
      WRITE (*,'(A)',ADVANCE='NO') 'ENTER N-REC: '
      READ (*,*) n_rec
!
      eeta_max = 0.0
      DO n = 1, n_rec
         PRINT*,'N_REC=', n
         CALL Read_array ('rcmeeta', n, label, ARRAY_3D = eeta)
         DO kc = 1, kcsize
            IF (MAXVAL(eeta(:,:,kc)) > eeta_max(kc)) eeta_max(kc)=MAXVAL(eeta(:,:,kc))
         END DO
      END DO
      OPEN (LUN, FILE='eeta_max.txt', STATUS='NEW')
      DO kc = 1, kcsize
         WRITE (LUN,*) kc, eeta_max(kc)
      END DO
      CLOSE (LUN)
      STOP
      END PROGRAM find_max_eeta