viernes, 13 de septiembre de 2019


=============================================
DESARROLLADOR: Ing.Ramirez Quispe, Robert Marlindo.
LENGUAJE               : Fortran 95.                    
PROGRAMA             : Modelo Onda Cinemática .           
                                       Transformar lluvia-escorrentía .   
FECHA                       : 13/09/2019 DD/MM/AA.               
LUGAR                      : U.N.H. Ingeniería Civi-Huancavelica.
CONSULTA               : ramirezquispe1@hotmail.com       
=============================================

EL PROGRAMA EN EL SIGUIENTE ENLACE: DESCARGAR





















domingo, 17 de febrero de 2019



Código fuente Lenguaje HPPPL y Lenguaje C.



Unidad De Post Grado De La Facultad De Ciencias De IngenieríaMaestría En Planificación Estratégica Y Gestión De Proyectos De Ingeniería
ASIGNATURA:
Estadística Aplicada A La Investigación Científica
ESTUDIANTE:
Ramirez Quispe, Robert Marlindo
DOCENTE:
Dr. Alfonso G. Cordero Fernández
I Evaluación Práctica 
DESCARGAR
Huancavelica, 03 de agosto de 2018


sábado, 12 de enero de 2019

jueves, 25 de octubre de 2018



CURSO DE LENGUAJE C++.
Se comparte al público ejemplos realizados
en los cursos realizados de C++
 para los estudiantes de la escuela de
Ing. Civil-Huancavelica.de
la Universidad Nacional de Huancavelica.
                      1._NIVEL 0: DESCARGA
2._NIVEL 1:
3._NIVEL 2:
4._NIVEL 3:
5._NIVEL 4:

viernes, 18 de mayo de 2018

miércoles, 14 de marzo de 2018

DIFERENCIAS FINITAS Y APLICACIÓN DEL MODELO DE ONDA CINEMÁTICA (APROXIMACIÓN DE LA ECUACIÓN DE SAINT-VENANT) RESULTADO OBTENIDOS CON LA AYUDA DEL LENGUAJE C++. 






viernes, 23 de febrero de 2018





Hazen-Williams y Jacobi. Cálculo de caudales en tuberías con la ecuación de Hazen-Williams y aplicación del Método Jacobi. DESCARGAR





jueves, 18 de mayo de 2017

CAUCES NATURALES EN LENGUAJE JULIA


CAUCES NATURALES EN LENGUAJE JULIA
Codificado por: Ramirez Quispe, Robert Marlindo.
Link Descarga: En proceso.
Desarrollado en la Universidad Nacional de Huancavelica-Escuela académico Profesional 
Ingeniera Civil-Hucnavelica






martes, 16 de mayo de 2017

CAUCES NATURALES EN LENGUAJE JULIA


CAUCES NATURALES EN LENGUAJE JULIA
Codificado por: Ramirez Quispe, Robert Marlindo.
Link Descarga: En proceso.
Desarrollado en la Universidad Nacional de Huancavelica-Escuela académico Profesional 
Ingeniera Civil-Hucnavelica

viernes, 20 de enero de 2017

CURVAS IDF EN LENGUAJE JULIA


CURVAS IDF EN LENGUAJE JULIA
Calcula los coeficientes, K,m y n con el Método de Bernard(1932). I =k*T^m/D^n  si quieren ver tabla las tablas, están almacenados en forma de matriz.Solo para poder imprimir.
Link para descargar el programa:

DESCARGAR PROGRAMA

#************* Curvas IDF *************
# Método de Bernard(1932).I =k*T^m/D^n
#Programación       : Ramirez Quispe, Robert Marlindo.
#Estudiante de      : Ing.de Civil- Universidad Nacional de Huancavelica.
#                             ramirezquispe1@hotmail.com
#                             ramirezquispe1@gmail.com
#facebook              : https://www.facebook.com/profile.php?id=100011674563458
#Canal You Tube   : https://www.youtube.com/channel/UCZq0GWDDkHnUKAbUX2HViRw

#DATOS INGRESAR
#*******************
# Dmm = Duracion en Minutos
Dmn=[  5 ,     10,      15,     20,      25,   30,      60   ];

# Intensidades(mm/hora),para cada duración.
#********************************************
In=[
      32.992   24.692   20.125  17.22   15.182  13.658   8.967 ;
      42.53    31.831   25.944  22.198  19.571  17.607   11.56 ;
      49.745   37.231   30.345  25.964  22.891  20.594   13.521;
      59.284   44.37    36.164  30.942  27.28   24.543   16.113;
      66.499   49.77    40.565  34.708  30.6    27.53    18.074;
      73.714   55.17    44.967  38.474  33.921  30.517   20.036;
      90.468   67.709   55.187  47.219  41.63   37.453   24.589]
#********************************************
#********************************************
n=size(In)
nf=n[1]
nc=n[2]
MoR=zeros(nf,nc)
TrO=ones(nf)    
LisC=ones(nf)
for I=1:nc
    LisC=LisC
  for J=1:nf
    LisC[J]=In[J,I]
  end
    LisC=reverse(sort(LisC))
   for K=1:nf
    MoR[K,I]=LisC[K]
   end
   TrO[I]=(nf+1)/I
end
Cfil=nc*nc
MFIN=zeros(Cfil,8)
LDmn=log10(Dmn)
LTrO=log10(TrO)
CoF=0
sy=0  
sxy1=0
sxy2=0
SX1=0
SX2=0
SSX2=0
SXX12=0
SXX22=0
for I=1:nc
  for J=1:nc
   MFIN[J+CoF,1]=LDmn[I]                    
   MFIN[J+CoF,2]=LTrO[J]                    
   MFIN[J+CoF,3]=log10(MoR[J,I])            
   MFIN[J+CoF,4]=MFIN[J+CoF,2]*MFIN[J+CoF,3]
   MFIN[J+CoF,5]=MFIN[J+CoF,1]*MFIN[J+CoF,3]
   MFIN[J+CoF,6]=MFIN[J+CoF,2]^2            
   MFIN[J+CoF,7]=MFIN[J+CoF,1]^2            
   MFIN[J+CoF,8]=MFIN[J+CoF,1]*MFIN[J+CoF,2]
   sy=sy+MFIN[J+CoF,3]
   sxy1=sxy1+MFIN[J+CoF,2]*MFIN[J+CoF,3]
   sxy2=sxy2+MFIN[J+CoF,1]*MFIN[J+CoF,3]
   SX1=SX1+MFIN[J+CoF,2]
   SX2=SX2+MFIN[J+CoF,1]
   SSX2=SSX2+MFIN[J+CoF,2]^2
   SXX12=SXX12+MFIN[J+CoF,1]*MFIN[J+CoF,2]
   SXX22=SXX22+MFIN[J+CoF,7]  
  end
  CoF=CoF+nc
end

rEsU=inv( [Cfil SX1     SX2;SX1  SSX2   SXX12;SX2  SXX12 SXX22])*[sy;sxy1;sxy2]
println("K = ",round(10^rEsU[1,1],3))
println("m = ",round(rEsU[2,1],3))
println("n = ",round(-rEsU[2,1],3))










domingo, 15 de enero de 2017

PROGRAMA EN PYTHON IILA-SENAMHI-UNI


link de descarga de programa: DESCARGAR PROGRAMA IILA-SENAMHI


"""
                   IILA-SENAMHI-UNI
Programación       : Ramirez Quispe, Robert Marlindo.
Docente            : Ing. Ivan Ayala Bizarro.
Estudiante de      : Ing.de Civil- Universidad Nacional de Huancavelica.
                     ramirezquispe1@hotmail.com
                     ramirezquispe1@gmail.com
facebook           : https://www.facebook.com/profile.php?id=100011674563458
Canal You Tube     : https://www.youtube.com/channel/UCZq0GWDDkHnUKAbUX2HViRw
Agradecimientos al : Ing. Reyder Bendezu, Boza. Por haberme dado, las bases adecuadas sobre programación.

Está totalmente prohibido, cambiar el nombre del autor. 
"""

from math import log10
deci=5
while True:
    print("1._Para t<3 horas ")
    print("2._Para 3<=t<=24 horas ")
    print("3._Salir ")
    opc=eval(input("Ingrese opción = "))
    if opc==1:
        a=eval(input("Parámetro Regional a= "))
        b=eval(input("Parámetro Regional b= "))
        K=eval(input("Parámetro Regional K= "))
        n=eval(input("Parámetro Regional n= "))
        T=eval(input("Periodo de retorno T[años]= "))
        Ti=eval(input("Tiempo inicial Ti[m]= "))
        Tf=eval(input("Tiempo final Tf[m]= "))
        In=eval(input("Incremento de tiempo In[m]= "))
        LT=[]
        LP=[]
        LI=[]
        LA=[]        
        PI=[a*(Ti/60)*(1+K*log10(T))*((Ti/60+b)**(n-1))]
        J=0
        while Ti<=Tf:
           P=a*(Ti/60)*(1+K*log10(T))*((Ti/60+b)**(n-1))
           I=a*(1+K*log10(T))*((Ti/60+b)**(n-1))
           LT.append(Ti)
           LP.append(P)
           LI.append(I)
           PI.append(a*((Ti+In)/60)*(1+K*log10(T))*(((Ti+In)/60+b)**(n-1))-P)
           Ti=Ti+In
           J=J+1
        if J%2==0:
            R=2
        else:
            R=1
        n=0
        for K in range(J):
            if K<J/2:
                LA.append(PI[J-2*K-R])
            else:
                LA.append(PI[2*n+1])
                n=n+1
        print("t[min]".center(8),"P[mm]".center(8),"I[mm/hr]".center(8),"PInc".center(8),"PAlt".center(8))
        for I in range(J):
            print(str(round(LT[I],deci)).center(8),str(round(LP[I],deci)).center(8),str(round(LI[I],deci)).center(8),str(round(PI[I],deci)).center(8),str(round(LA[I],deci)).center(8))
    if opc==2:
        a=eval(input("Parámetro Regional a= "))
        K=eval(input("Parámetro Regional K= "))
        n=eval(input("Parámetro Regional n= "))
        T=eval(input("Periodo de retorno T[años]= "))
        Ti=eval(input("Tiempo inicial Ti[m]= "))
        Tf=eval(input("Tiempo final Tf[m]= "))
        In=eval(input("Incremento de tiempo In[m]= "))
        LT=[]
        LP=[]
        LI=[]
        LA=[]        
        PI=[a*(1+K*log10(T))*((Ti)**n)]
        J=0
        while Ti<=Tf:
           P=a*(1+K*log10(T))*((Ti)**n)
           I=a*(1+K*log10(T))*((Ti)**(n-1))
           LT.append(Ti)
           LP.append(P)
           LI.append(I)
           PI.append(a*(1+K*log10(T))*(((Ti+In))**n)-P)
           Ti=Ti+In
           J=J+1
        if J%2==0:
            R=2
        else:
            R=1
        n=0
        for K in range(J):
            if K<J/2:
                LA.append(PI[J-2*K-R])
            else:
                LA.append(PI[2*n+1])
                n=n+1
        print("t[min]".center(8),"P[mm]".center(8),"I[mm/hr]".center(8),"PInc".center(8),"PAlt".center(8))
        for I in range(J):
            print(str(round(LT[I],deci)).center(8),str(round(LP[I],deci)).center(8),str(round(LI[I],deci)).center(8),str(round(PI[I],deci)).center(8),str(round(LA[I],deci)).center(8))
    if opc==3:
        print("Gracias!!!!!")
        break
   




___________________

EJEMPLO PARA t<3 horas:
__________________________







__________________________

EJEMPLO PARA 3<=t< 24 horas.
__________________________