1 | program Interpol |
---|
2 | |
---|
3 | !Programa que a partir de un file en ijS de 5km lo tranforma en un file xyS y luego |
---|
4 | ! usando el comando paste podremos tener los resultados en ijxylonlat S |
---|
5 | !Solo usarlo para datos en 5km |
---|
6 | |
---|
7 | implicit none |
---|
8 | |
---|
9 | integer::l,x,y,n |
---|
10 | |
---|
11 | real,dimension(1406160)::I |
---|
12 | real,dimension(1406160)::J |
---|
13 | real,dimension(1406160)::K |
---|
14 | |
---|
15 | !introducir el nombre del file en ij a transformar examples: |
---|
16 | !open(12,file='/pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/prueba_giu1') |
---|
17 | !open(12,file='/pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/grille-surface-1143x1143-ijS') |
---|
18 | !open(12,file='/pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/grille-epaisseur-1143x1143-ijE') |
---|
19 | !open(12,file='/pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/grille-bedelev-1143x1143-ijB') |
---|
20 | !open(12,file='/pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/grille-water-1143x1143-ijW') |
---|
21 | !open(12,file='/pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/grille-bathy-1143x1143-ijBa') |
---|
22 | open(12,file='/pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/grille-ground-1143x1143-ijGr') |
---|
23 | |
---|
24 | |
---|
25 | do l=1,1406160 |
---|
26 | read(12,*)I(l),J(l),K(l) |
---|
27 | end do |
---|
28 | close(12) |
---|
29 | |
---|
30 | n=0 |
---|
31 | do y=-2720000,2855000,5000 |
---|
32 | do x=-2985000,3310000,5000 |
---|
33 | |
---|
34 | n=n+1 |
---|
35 | write(13,*) x, y, K(n) |
---|
36 | end do |
---|
37 | end do |
---|
38 | |
---|
39 | end program Interpol |
---|
40 | |
---|
41 | |
---|
42 | |
---|
43 | ! despues hacer la transformacion de xy en lonlat para ello hacer en la consola: |
---|
44 | ! 1) cp fort.13 grille-surface-1143x1143-xyS |
---|
45 | ! 2) awk '{print $1,$2}' /pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/grille-surface-1143x1143-xyS > xy-5k |
---|
46 | ! 3) awk '{print $1,$2}' /pcgn/home/navas/Grisli/Transfor_datos_giu/Prueba/grille-surface-1143x1143-ijS > ij-5k |
---|
47 | ! 4) mapproject xy-5k -Js0/-90/-71/1:1000 -Fk -C -R-180/+180/-90/0 -I > lonlat-5k |
---|
48 | ! 5) paste ij-5k xy-5k lonlat-5k > ijxylonlat-5km |
---|