M.Sc. Ing.
Jorge Diaz Zelada
MODELACION DE LA EVAPOTRANSPIRACION POTENCIAL
1) Una vez cambiado el atributo del mapa PyT_Simulados a “Temp”, SE PROCEDE
A CREAR EL MAPA DE PUNTOS DE TEMPERATURAS. Presione botón derecho sobre
el mapa de puntos “PyT_Simulados” y utilice los comandos /Vector Operations/Attribute
Map y con la ayuda de la tabla “TEMP” y las columnas correspondientes de cada mes
(Attribute) crear mapas de puntos (temene, temfeb, temmar, etc).
Utilizar la script T1
2) A partir de los mapas de puntos crear los mapas interpolados de temeperatura
raterizados (temene, temfeb, temmar, etc.): Interpolation/Moving average (utilizar la
georeferencia CuencaCerro)
La macro T2 no corre
3) Cortar los mapas en función al mapa “PunataLimite”: p ej.
temenec:=iff(isundef(Geo_Cerro_Rico),?,temene)
Temfebc:=iff(isundef(Geo_Cerro_Rico),?,temfeb)
temmarc:=iff(isundef(Geo_Cerro_Rico),?,temmar)
temabrc:=iff(isundef(Geo_Cerro_Rico),?,temabr)
temmayc:=iff(isundef(Geo_Cerro_Rico),?,temmay)
temjunc:=iff(isundef(Geo_Cerro_Rico),?,temjun)
temjulc:=iff(isundef(Geo_Cerro_Rico),?,temjul)
temagoc:=iff(isundef(Geo_Cerro_Rico),?,temago)
temsepc:=iff(isundef(Geo_Cerro_Rico),?,temsep)
temoctc:=iff(isundef(Geo_Cerro_Rico),?,temoct)
temnovc:=iff(isundef(Geo_Cerro_Rico),?,temnov)
temdicc:=iff(isundef(Geo_Cerro_Rico),?,temdic)
Hacer correr la macro T3
Sin embargo, estos mapas representan interpolaciones lineales de la temperatura en función
de pocas estaciones y no representan la variación de la temperatura, sobre todo si tomamos
en cuenta la variación topográfica del área de estudio; por tanto, se debe hacer la corrección
de los mismos en función de la relación de 0.46°C de variación de la temperatura por cada
100 metros de variación de altura.
3) A partir de archivo de puntos “PyT_Simulados” y con la ayuda de la columna “COTAS”
de la tabla “TEMP” crear mapas de puntos “temalt” con las alturas de las estaciones.
4) A partir del mapa de puntos “temalt” crear un mapa raster de alturas con los comandos:
Interpolation/Moving average (utilizar la georeferencia CuencaCerro). Este mapa tiene
el mismo nombre.
1
M.Sc. Ing. Jorge Diaz Zelada
5) Cortar este mapa en función al mapa “Geo_Cerro_Rico”: p ej.
temaltc:=iff(isundef(Geo_Cerro_Rico),?,temalt)
6) Obtener un CuencaCerroDEM del área de estudio.
7) Obtener un mapa “difalt”, resultado de la diferencia de alturas entre los mapas “temalt” y
“CuencaCerroDEM”.
Difalt:=temaltc- CuencaCerroDEM
9) Se debe multiplicar el mapa “difalt” con la constante 0.46 y dividirlo por la constante
100, para obtener un mapa “tem_cor” que representa la suma algebraica de la corrección de
temperaturas a los doce mapas de temperatura:
tem_cor:=(0.46*difalt)/100
9) Los mapas de temperatura mensual (temenec, temfebc, etc), deben ser corregidos
sumando (algebraicamente) el mapa “tem_cor”, creando los mapas “temenecor, temfebcor,
etc.)
temenecor:= temenec+tem_cor
temfebcor:= temfebc+tem_cor
temmarcor:= temmarc+tem_cor
temabrcor:= temabrc+tem_cor
temmaycor:= temmayc+tem_cor
temjuncor:= temjunc+tem_cor
temjulcor:= temjulc+tem_cor
temagocor:= temagoc+tem_cor
temsepcor:= temsepc+tem_cor
temoctcor:= temoctc+tem_cor
temnovcor:= temnovc+tem_cor
temdiccor:= temdicc+tem_cor
Hacer correr la macro T4
10) Calcular el índice de calor mensual (“iene”, “ifeb”, etc.)
iene:=(temenecor/5)^1.514
iene:=iff(temenecor<0,0,(temenecor/5)^1.514)
ifeb:=(temfebcor/5)^1.514
ifeb:=iff(temfebcor<0,0,(temfebcor/5)^1.514)
imar:=(temmarcor/5)^1.514
2
M.Sc. Ing. Jorge Diaz Zelada
imar:=iff(temmarcor<0,0,(temmarcor/5)^1.514)
iabr:=(temabrcor/5)^1.514
iabr:=iff(temabrcor<0,0,(temabrcor/5)^1.514)
imay:=(temmaycor/5)^1.514
imay:=iff(temmaycor<0,0,(temmaycor/5)^1.514)
ijun:=(temjuncor/5)^1.514
ijun:=iff(temjuncor<0,0,(temjuncor/5)^1.514)
ijul:=(temjulcor/5)^1.514
ijul:=iff(temjulcor<0,0,(temjulcor/5)^1.514)
iago:=(temagocor/5)^1.514
iago:=iff(temagocor<0,0,(temagocor/5)^1.514)
isep:=(temsepcor/5)^1.514
isep:=iff(temsepcor<0,0,(temsepcor/5)^1.514)
ioct:=(temoctcor/5)^1.514
ioct:=iff(temoctcor<0,0,(temoctcor/5)^1.514)
inov:=(temnovcor/5)^1.514
inov:=iff(temnovcor<0,0,(temnovcor/5)^1.514)
idic:=(temdiccor/5)^1.514
idic:=iff(temdiccor<0,0,(temdiccor/5)^1.514)
Hacer correr la macro T5
11) Calcular el índice de calor anual
I:=iene+ifeb+imar+iabr+imay+ijun+ijul+iago+isep+ioct+inov+idic
12) Crear el mapa del factor de modificación (a)
MapaA:=6.75e-7*I^3-7.71e-5*I^2+1.792e-2*I+0.49239
13) Obtener mapas de horas de luminosidad para cada mes, a partir de la tabla EXCEL
3
M.Sc. Ing. Jorge Diaz Zelada
“TABLAS DE HORAS DE LUMINOSIDAD”. Las constantes mensuales están en función de
la latitud del área de estudio que corresponde a los 17°S. Estos mapas son de constantes en
toda el área de la cuenca para los doce meses del año (dene, dfeb, dmar, etc) de acuerdo a la
tabla de factores de duración posible de luminosidad.
Las horas de luminosidad está en el archivo Excel de “TABLAS US CONSERVATION-
2017”. Con datos 22°
dene:=iff(isundef(Geo_Cerro_Rico),?,1.15)
dfeb:=iff(isundef(Geo_Cerro_Rico),?,1.00)
dmar:=iff(isundef(Geo_Cerro_Rico),?,1.05)
dabr:=iff(isundef(Geo_Cerro_Rico),?,0.97)
dmay:=iff(isundef(Geo_Cerro_Rico),?,0.95)
djun:=iff(isundef(Geo_Cerro_Rico),?,0.90)
djul:=iff(isundef(Geo_Cerro_Rico),?,0.94)
dago:=iff(isundef(Geo_Cerro_Rico),?,0.99)
dsep:=iff(isundef(Geo_Cerro_Rico),?,1.00)
doct:=iff(isundef(Geo_Cerro_Rico),?,1.09)
dnov:=iff(isundef(Geo_Cerro_Rico),?,1.10)
ddic:=iff(isundef(Geo_Cerro_Rico),?,1.16).
Hacer correr la macro T6
14) Se crean las ETP para cada mes aplicando la ecuación de Thornthwaite-Mather
etpene:=16*((10*temenecor)/I)^mapaA*dene
etpene:=iff(temenecor<0,0,16*((10*temenecor)/I)^MapaA*dene)
4
M.Sc. Ing. Jorge Diaz Zelada
etpfeb:=iff(temfebcor<0,0,16*((10*temfebcor)/I)^MapaA*dfeb)
etpmar:=iff(temmarcor<0,0,16*((10*temmarcor)/I)^MapaA*dmar)
etpabr:=iff(temabrcor<0,0,16*((10*temabrcor)/I)^MapaA*dabr)
etpmay:=iff(temmaycor<0,0,16*((10*temmaycor)/I)^MapaA*dmay)
etpjun:=iff(temjuncor<0,0,16*((10*temjuncor)/I)^MapaA*djun)
etpjul:=iff(temjulcor<0,0,16*((10*temjulcor)/I)^MapaA*djul)
etpago:=iff(temagocor<0,0,16*((10*temagocor)/I)^MapaA*dago)
etpsep:=iff(temsepcor<0,0,16*((10*temsepcor)/I)^MapaA*dsep)
etpoct:=iff(temoctcor<0,0,16*((10*temoctcor)/I)^MapaA*doct)
etpnov:=iff(temnovcor<0,0,16*((10*temnovcor)/I)^MapaA*dnov)
etpdic:=iff(temdiccor<0,0,16*((10*temdiccor)/I)^MapaA*ddic)
etpdic:=16*((10*temdiccor)/I)^mapaA*ddic
Hacer correr la macro T7
15) Obtener el mapa de ETP anual sumando todos los mapas de ETP mensual:
)
