module module_read_ampere !Afegit Santi.
        use hist_module,only: modeltime
        use init_module,only: iyear
        use params_module, only: nmlon, nmlat
        implicit none
        integer :: i, ios, j   !, numlin
!       integer, parameter :: nmlon = 80, nmlat = 97
!       double precision :: lo, la, Jr0, lon(nmlon,nmlat),
!    | lat(nmlon,nmlat)
        real :: Jr_ampere(nmlon+1,nmlat)

      contains
        subroutine read_ampere
          implicit none
          character*100 :: nomarx
          character*3 :: dia
          character*2 :: ut,minut
          character*4 :: any
          write(dia,'(i3.3)') modeltime(1)
          write(ut,'(i2.2)') modeltime(2)
          write(minut,'(i2.2)') modeltime(3)
          write(any,'(i4)') iyear

          j = j + 1
          if (j==1) then
            write(minut,'(i2.2)') modeltime(3) - 2
          endif
          write(nomarx,'(a100)') '/hao/aim2/marsal/tiegcm/
     |AMPERE/Ampere_Jr_data/Jr_'//any//dia//ut//minut//'.dat'
          write(6,*) nomarx
          open(unit=1, file=nomarx, status='old', iostat=ios)
            do i = 1, nmlon+1
              read(1,*,iostat=ios) Jr_ampere(i,1:nmlat)
            enddo  
          close(1)
          write(6,*)'j,Jr_ampere(3,2)', j,Jr_ampere(3,2)

        end subroutine read_ampere

      end module module_read_ampere