1. Alberto Giménez
  2. GUYRA

Commits

agimenez  committed 21bb59f

- Creación de tablas de diferencias entre proyecciones y la normal 6190.

  • Participants
  • Parent commits 23ff29a
  • Branches default

Comments (0)

Files changed (4)

File resumen_normal_6190.py

View file
  • Ignore whitespace
+#!/usr/bin/env python
+from sys import path
+from numpy import average,std
+import numpy as np
+import locale
+from guyra_lib import *
+
+path.append("/home/clima/cdat-5.2/lib64/python2.6/site-packages/")
+
+import cdms2
+
+LISTA = ["anual",
+        "semestre_EJ",
+        "semestre_JD",
+        "trimestre_EFM",
+        "trimestre_AMJ",
+        "trimestre_JAS",
+        "trimestre_OND"]
+
+def lista_normal():
+    ff = open("tablas/tabla_normales_cli.txt","w")
+    ff.write("\tPRECIPITACION\tTEMPERATURA MEDIA\n")
+    for fDesc in LISTA:
+        ds = cdms2.open("datos/grid_6190_%s.nc" % fDesc)
+        ff.write("\n" + SUBTITULOS[fDesc] + "\n")
+        val1 = ds("prec",squeeze=1,latitude=MI_bnd[1],longitude=MI_bnd[0])
+        val2 = ds("temp",squeeze=1,latitude=MI_bnd[1],longitude=MI_bnd[0])
+        #print val1.shape,val2.shape
+        ff.write(locale.format_string("PROMEDIO\t%.4f\t%.4f\n", (average(val1),average(val2))))
+        ff.write(locale.format_string("DESVIACION STD\t%.4f\t%.4f\n", (std(val1),std(val2))))
+        ds.close()
+    ff.close()
+
+def lista_proyeccion(var):
+    TIT = {"prec":"PRECIPITACION","temp":"TEMPERATURA MEDIA"}
+    ff = open("tablas/tabla_proyeccion_%s.txt" % var,"w")
+    ff.write("%s\n" % TIT[var])
+    ff.write("\tVARIACION A2\t\tVARIACION B2\n")
+    ff.write("\tABSOLUTA (mm/dia)\tPORCENTUAL (%)\tABSOLUTA (mm/dia)\tPORCENTUAL (%)\tDIFERENCIA ABSOLUTA (mm/dia)\n")
+    for fDesc in LISTA:
+        ds = cdms2.open("datos/grid_6190_%s.nc" % fDesc)
+        vds = cdms2.open("datos/grid_%s.nc" % fDesc)
+        ff.write(SUBTITULOS[fDesc] + "\t")
+        val = ds(var,squeeze=1,latitude=MI_bnd[1],longitude=MI_bnd[0])
+        val1 = vds(var+"_a2",squeeze=1,latitude=MI_bnd[1],longitude=MI_bnd[0])
+        val2 = vds(var+"_b2",squeeze=1,latitude=MI_bnd[1],longitude=MI_bnd[0])
+        base = average(val)
+        pval1 = average(val1)
+        pval2 = average(val2)
+        abs1 = pval1 - base
+        por1 = (pval1/base - 1)*100
+        abs2 = pval2 - base
+        por2 = (pval2/base - 1)*100
+        #print pval1,pval2,base
+        ff.write(locale.format_string("%.4f\t%.2f\t%.4f\t%.2f\t%.4f\n",(abs1,por1,abs2,por2,abs1-abs2)))
+        #print val.shape,val1.shape,val2.shape
+        ds.close()
+        vds.close()
+    ff.close()    
+
+if __name__ == '__main__':
+    locale.setlocale(locale.LC_NUMERIC, "es_PY")
+    #lista_normal()
+    lista_proyeccion("prec")
+    #lista_proyeccion("temp")
+

File tablas/tabla_normales_cli.txt

View file
  • Ignore whitespace
+	PRECIPITACION	TEMPERATURA MEDIA
+
+Promedio Anual de Enero a Diciembre
+PROMEDIO	4,6801	21,3932
+DESVIACION STD	0,1611	0,1864
+
+Promedio Semestral de Enero a Junio
+PROMEDIO	4,4668	22,2653
+DESVIACION STD	0,1905	0,1399
+
+Promedio Semestral de Julio a Diciembre
+PROMEDIO	4,7123	19,9567
+DESVIACION STD	0,1963	0,2493
+
+Promedio Trimestral de Enero a Marzo
+PROMEDIO	4,8046	25,5945
+DESVIACION STD	0,1471	0,0994
+
+Promedio Trimestral de Abril a Junio
+PROMEDIO	4,1291	18,9360
+DESVIACION STD	0,2487	0,2728
+
+Promedio Trimestral de Julio a Septiembre
+PROMEDIO	4,0060	17,2982
+DESVIACION STD	0,2458	0,3376
+
+Promedio Trimestral de Octubre a Diciembre
+PROMEDIO	5,7806	23,7442
+DESVIACION STD	0,1854	0,1484

File tablas/tabla_proyeccion_prec.txt

View file
  • Ignore whitespace
+PRECIPITACION
+	VARIACION A2		VARIACION B2
+	ABSOLUTA (mm/dia)	PORCENTUAL (%)	ABSOLUTA (mm/dia)	PORCENTUAL (%)	DIFERENCIA ABSOLUTA (mm/dia)
+Promedio Anual de Enero a Diciembre	0,1802	3,85	0,1682	3,59	0,0120
+Promedio Semestral de Enero a Junio	0,2853	6,39	0,3229	7,23	-0,0376
+Promedio Semestral de Julio a Diciembre	0,1840	3,90	0,1477	3,13	0,0364
+Promedio Trimestral de Enero a Marzo	0,0275	0,57	0,0743	1,55	-0,0468
+Promedio Trimestral de Abril a Junio	0,5432	13,15	0,5715	13,84	-0,0284
+Promedio Trimestral de Julio a Septiembre	0,5871	14,66	0,5826	14,54	0,0045
+Promedio Trimestral de Octubre a Diciembre	-0,4370	-7,56	-0,5557	-9,61	0,1186

File tablas/tabla_proyeccion_temp.txt

View file
  • Ignore whitespace
+TEMPERATURA MEDIA
+	VARIACION A2		VARIACION B2
+	ABSOLUTA (mm/dia)	PORCENTUAL (%)	ABSOLUTA (mm/dia)	PORCENTUAL (%)	DIFERENCIA ABSOLUTA (mm/dia)
+Promedio Anual de Enero a Diciembre	0,9644	4,51	0,8791	4,11	0,0853
+Promedio Semestral de Enero a Junio	0,8305	3,73	0,7083	3,18	0,1223
+Promedio Semestral de Julio a Diciembre	1,0445	5,23	1,0049	5,04	0,0396
+Promedio Trimestral de Enero a Marzo	0,8627	3,37	0,6882	2,69	0,1745
+Promedio Trimestral de Abril a Junio	0,7984	4,22	0,7283	3,85	0,0701
+Promedio Trimestral de Julio a Septiembre	1,2542	7,25	1,1054	6,39	0,1489
+Promedio Trimestral de Octubre a Diciembre	0,9422	3,97	0,9946	4,19	-0,0524