program RAMS_to_GRADS
!**********************************************************************
! *              .      .    .
! * PROGRAM:    RAMS_to_GRADS
! * PRGMMR:     D. ZUPANSKI     ORG: CIRA/CSU   DATE: 2004-09-23
! *
! * ABSTRACT:   Read an EnsDA history file and convert it to a 
!               format suitable for use in GRADS
! *
! * PROGRAM LOG:
! *
! * 09/09/2003    ..... M. ZUPANSKI: Original 'transform.F'
! * 10/17/2003    ..... M. ZUPANSKI: Innovation normalization
! * 10/22/2003    ..... M. ZUPANSKI: Innovation histogram
! * 10/22/2003    ..... M. ZUPANSKI: plotting
! * 10/24/2003    ..... M. ZUPANSKI: Root-Mean-Squared error calculation
! * 06/25/2004    ..... D. ZUPANSKI: Root-Mean-Squared calculation for RAMS
! * 09/23/2004    ..... D. ZUPANSKI: RAMS_to_GRADS
! *
! **********************************************************************
!-----
     integer,parameter::in_file=21       ! input file #
     integer,parameter::iout_file=51     ! output file #
!-----
     integer :: im,jm,lm
     integer :: histvar_max,kk
     character(len=9)         :: u_name

     real,dimension(:,:,:),allocatable::u
     real,dimension(:,:),allocatable::u_2d
!-----------------------------------------
      histvar_max=6

 100  format(E20.10)
 110  format(3i7)
 300  format(a9)

     rewind in_file
     rewind iout_file

     read(in_file) im,jm,lm

     allocate(u(1:im,1:jm,1:lm))

               do kk=1,histvar_max

     read(in_file) u_name
     read(in_file) u
     write(*,*) " u_name=",u_name," min,max=",minval(u),maxval(u)

     allocate(u_2d(1:im,1:jm))
                  do l=1,lm
     u_2d(:,:)=u(:,:,l)
     write(iout_file) u_2d
                  end do
     deallocate(u_2d)

               end do     !!!! do kk=1,histvar_max

      end program RAMS_to_GRADS