Por fin he podido hacer algo !!
1.- Repetido el paso de proceso1, pero ahora respetando los NA originales de los datos HadISST. El resultado en el fichero septiembres_proc1.nc
2.- Media de los años 1953 a 1978 en el fichero septiembres_media5378.nc. Y media de los años 1982 a 1999 en el fichero septiembres_media8299.nc.
3.- En el proceso2 tenía una duda: anteriormente solo habíamos utilizado la media para un hemisferio, entre -180 y 0 º de latitud Este, y dejado en el otro hemisferio los valores originales. Pero no tenía claro si preferías utilizar la media para todos los valores, así que he hecho las dos cosas:
En el fichero septiembres_proc2_5378_180.nc, se utiliza la media de los años 1953 a 1978, y se sustituye en aquellos puntos en los que no hubiese datos de AARI y el valor original de HadISST no fuera NA, pero solo entre -180 y 0º. El fichero septiembres_proc2_5378_360.nc hace lo mismo pero para todo el globo. Los correspondientes ficheros septiembres_proc2_8299_180.nc y septiembres_proc2_8299_360.nc son lo mismo pero utilizando la media de 1982 a 1999. Esto se hace con los proceso2b.R y proceso2c.R
4.- Se toman cada uno de estos 4 ficheros nc, y con aunoceros.R se ponen los valores de hielo a 1 si tienen >=0.15, y se ponen a cero los del hemisferio sur. Se generan los ficheros sep_proc3_5378_180.nc, sep_proc3_5378_360.nc, sep_proc3_8299_180.nc, sep_proc3_8299_360.nc.
5.- Los ficheros:
ext_proc1.txt, ext_media8299.txt, ext_media5378.txt, ext_proc2_8299_360.txt, ext_proc2_8299_180.txt, ext_proc2_5378_360.txt, ext_proc2_5378_180.txt, ext_proc3_8299_360.txt, ext_proc3_8299_180.txt, ext_proc3_5378_360.txt, ext_proc3_5378_180.txt
tienen las extensiones correspondientes a cada fichero nc.
Todavía no me ha dado tiempo a revisarlo, así que igual hay algún error gordo, a ver si puedo revisarlo con más tiempo mañana o pasado.
Los ficheros están en la carpeta:
https://drive.google.com/folderview?id=0B0qq1MZw8pX-dEVwNlVwYy1uMmc&usp=sharingPero los ficheros .nc todavía están subiendo, le debe quedar todavía casi media hora para acabar.