program Interpol

!Programa que a partir de un file en ijS de 5km lo tranforma en un file xyS y luego
! usando el comando paste podremos tener los resultados en ijxylonlat S
!Solo usarlo para datos en 5km

implicit none

integer::l,x,y,n

real,dimension(1406160)::I
real,dimension(1406160)::J
real,dimension(1406160)::K

!introducir el nombre del file en ij a transformar examples:
!open(12,file='/pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/prueba_giu1')
!open(12,file='/pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/grille-surface-1143x1143-ijS')
!open(12,file='/pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/grille-epaisseur-1143x1143-ijE')
!open(12,file='/pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/grille-bedelev-1143x1143-ijB')
!open(12,file='/pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/grille-water-1143x1143-ijW')
!open(12,file='/pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/grille-bathy-1143x1143-ijBa')
open(12,file='/pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/grille-ground-1143x1143-ijGr')


do l=1,1406160
   read(12,*)I(l),J(l),K(l)
end do
close(12)

n=0
do y=-2720000,2855000,5000
do x=-2985000,3310000,5000

n=n+1
write(13,*) x, y, K(n)
end do
end do

end program Interpol



! despues hacer la transformacion de xy en lonlat para ello hacer en la consola:
! 1) cp fort.13 grille-surface-1143x1143-xyS
! 2) awk '{print $1,$2}' /pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/grille-surface-1143x1143-xyS > xy-5k
! 3) awk '{print $1,$2}' /pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/grille-surface-1143x1143-ijS > ij-5k
! 4) mapproject xy-5k -Js0/-90/-71/1:1000 -Fk -C -R-180/+180/-90/0 -I > lonlat-5k
! 5) paste ij-5k  xy-5k  lonlat-5k > ijxylonlat-5km