************************************************************************
* >Fl /4Yb /FPi RADFLCC.FOR FLURAD CORTOR CORTOA CORTOG RADDXF CORSUB SUBIMP
************************************************************************
*       NOMBRE:         RADFLCC
*       ORIGEN:         CATEDRA DE SISTEMAS DE POTENCIA U.N.L.P
*       FECHA:          1997
*       AUTOR:          A.RIFALDI
*       DESCRIPCION:    REDES RADIALES, FLUJO, CORTOCIRCUITO
*       LENGUAJE:       FORTRAN 77
***********************************************************************
* FLUJO DE CARGAS EN REDES RADIALES Y OTROS:
* FLURAD, CORTOR, CORTOG, CORTOA
* ---------------------------------------------------------------------
*
* PRIMERA TARJETA
*     NOMBRE DE IDENTIFICACION DEL TRABAJO 30,10,30,10 LETRAS
*aaaaaaaaaaaaaaaaaaaaaaaaaaaa><bbbbbbbb><cccccccccccccccccccccccccccc><ddddd
*
* UN LOTE DE TARJETAS DE TEXTO A IMAGEN QUE SE IMPRIMEN, LA PRIMERA
*     COLUMNA CON 1, Y EL TEXTO LIBRE, PARA ACABAR 0 EN COLUMNA 1
*     AL MENOS LA ULTIMA TARJETA
* un lote de  tarjetas de explicacion: la primera  columna con 1, y  el texto
* libre 79 caracteres, para acabar 0 en columna 1, al menos la ultima tarjeta
*
* TODAS LAS TARJETAS SIGUIENTES EN FORMATO DE UN DATO CADA 10 COLUMNAS
*     SE REGISTRAN HASTA 7 DATOS NUMERICOS INDICANDO EL PUNTO DECIMAL
*     Y LUEGO HASTA 36 CARACTERES DESDE COLUMNA 71
*
* SEGUNDA TARJETA
*     PB POTENCIA BASE EN MVA
*     VB TENSION BASE EN KV
*     FB FRECUENCIA BASE EN HZ
*     R3 FACTOR DE FASE (1.730508 PARA TRIFASICO)
*     DT DESCRIPCION DE TOPOLOGIA 1 NIVELES, 2 RAMAS
*   PB   ><   VB   ><   FB   ><   R3   ><   DT   ><        ><        ><
*
* TERCERA TARJETA - PRIMERA RAMA - CIRCUITO DE ALIMENTACION
*     PCC POTENCIA DE CORTOCIRCUITO EN MVA O POTENCIA DEL GENERADOR EN MVA
*     ZAL IMPEDANCIA FUNCIONAL DEL GENERADOR, SI PCC ES LA POTENCIA DE
*         CORTOCIRCUITO, PONER 0 EL PROGRAMA ADOPTA 100%, SI EN CAMBIO
*         SE TRATA DE GENERADOR, PONER SU IMPEDANCIA EN %.
*     XSR RELACION ENTRE X Y R DE LA PCC O DEL GENERADOR
*     EE  TENSION EN EL GENERADOR DE THEVENIN EN VALOR RELATIVO (PU)
*     R13 RELACION DE POTENCIA DE CORTOCIRCUITO MONOFASICO A TRIFASICO
*         (PARA CORTOG)
*     XSR1 RELACION ENTRE X Y R DE FALLA MONOFASICA (PARA CORTOG)
*     SIMULT FACTOR DE REDUCCION DE LAS CARGAS, SIMULTANEIDAD, O EL SIGNIFI-
*         CADO QUE QUIERA DARSELE, SIRVE PARA REDUCIR TODAS LAS CARGAS AL
*         MISMO TIEMPO POR EJEMPLO AL 0.5 (SIRVE PARA FLUCAR)
*     NOMBAR NOMBRE DE LA BARRA 36 CARACTERES DESDE COLUMNA 71, A VECES
*         CONVIENE ESCRIBIR LOS 36 CARACTERES, PONENDO EN PATICULAR UN
*         SIMBOLO EN COLUMNA 107 PARA QUE NO LEVANTE BASURA.
*  PCC   ><  ZAL   ><  XSR   ><   EE   ><   R13  ><  XSR1  >< SIMULT ><NOMBAR
*
* LOTE DE TARJETAS -  DATOS DE CADA ELEMENTO DE UNA RAMA
*     FORMADO GRUPOS DE CUATRO TARJETAS AA, BB, CC Y DD LA TARJETA AA TIENE
*     DATOS DE LA IMPEDANCIA DE LA RAMA, BB DATOS DE LA IMPEDANCIA DE
*     SECUENCIA CERO DE LA RAMA, LOS TRANSFORMADORES DE DISTINTOS GRUPOS,
*     CC DATOS DE LA CARGA CONCENTRADA
*     EN NODO FINAL (NUNCA DEBE CONSIDERARSE LA CARGA DE OTRA LINEA)
*     Y DD DATOS DE LA CARGA CONCENTRADA
*     EN NODO FINAL Y QUE REPERESENTA EL APORTE AL CORTO DE ELLA, SI SE
*     TRATA DE MOTORES SU CORRIENTE DE ARRANQUE (NUNCA DEBE CONSIDERARSE
*     LA CARGA DE OTRA LINEA)
*     LOS NODOS SE NUMERAN EN NOTACION POSICIONAL, DE MANERA QUE SU
*     VALOR INDICA LA DISTANCIA AL ORIGEN DEL GRAFO, SE DEBEN RECORRER
*     EN RIGUROSO ORDEN PARA DESCRIBIR EL GRAFO, LAS RAMAS SE
*     IDENTIFICAN CON NUMEROS CORRELATIVOS Y SUS EXTREMOS SE IDENTIFICAN
*     CON EXTREMO INICIAL NUMERO DE LA RAMA ANTERIOR, FINAL NOMERO DE
*     LA RAMA OBJETO.
*     TARJETA AA
*     SON EQUIVALENTES LOS SIGUIENTES REGISTROS
*     R(1)   PT   HL   POTENCIA DEL TRANSFORMADOR EN MVA
*                      O LONGITUD DE LA LINEA EN KM
*     R(2)   ZP   RL   IMPEDANCIA DEL TRANSFORMADOR EN 0/0
*                      O RESISTENCIA DE LA LINEA EN OHM/KM
*     R(3)   PP   XL   PERDIDAS DEL TRANSFORMADOR EN 0/0
*                      O REACTANCIA DE LA LINEA EN OHM/KM
*     R(4)   U1   NP   TENSION PRIMARIA DEL TRANSFORMADOR EN KV
*                      O NUMERO DE CABLES EN PARALELO
*     R(5)   U2   TENSION SECUNDARIA DEL TRANSFORMADOR EN KV
*                      PARA LA LINEA R(5) SE CONSIDERA NULO
*                      Y HAY QUE LLENARLO CON '0.'
*            TAP  PARA TRANSF. DERIVACION +/- en  %; U2 * (1 + TAP / 100)
*                      PARA LA LINEA RELLENAR CON '0.'
*     NOD   NUMERO DE NODO EN DONDE NACE LA RAMA - HAY 2 CRITERIOS
*           DEFINIDOS POR LA VARIABLE DT (DESCRIPCION TOPOLOGIA)
*           DT <= 1 EL PROGRAMA NUMERA LOS NODOS EN BASE A LA DISTANCIA
*                   AL INICIAL ORIGEN, POR LO TANTO HAY VARIAS HOJAS
*                   CON IGUAL NUMERO DE NODO, LOS NODOS ESTAN NUMERADOS
*                   POR NIVELES, CUANDO SE QUIERE VOLVER ATRAS SE
*                   INDICA EL NIVEL DE ARRANQUE, EL PROGRAMA VA AL ULTIMO
*                   NODO CARGADO DE ESE NIVEL
*           DT => 2 LOS NODOS ESTAN NUMERADOS NATURALMENTE (NO SE REPITEN
*                   NUMEROS), EL NODO ORIGEN DE CADA RAMA ES EL ANTERIOR,
*                   CUANDO SE QUIERE VOLVER ATRAS, SE CITA COMO ORIGEN EL
*                   NODO DONDE LA RAMA SE CONECTA (CUIDADO LOS NUMEROS DE
*                   CON PUNTO PUES EL PROGRAMA LEE UN REAL)
*                   NOTA: EL NODO FINAL Y EL NUMERO DE RAMA SON COINCIDENTES
*     NOMBAR NOMBRE DE LA BARRA 36 CARACTERES DESDE COLUMNA 71
*                   EL NOMBRE DE LA BARRA ES EL DEL NODO FINAL DE LA RAMA
*  PT HL ><  ZP RL ><  PP XL ><  UP NP ><   U2   ><   TAP  ><   NOD  ><NOMBAR
*
*     LA SEGUNDA TARJETA DEL GRUPO, BB CONTIENE LOS SIGUIENTES REGISTROS
*     QUE CORRESPONDEN A PARAMETROS DE SECUENCIA CERO (SOLO SON USADOS
*     POR CORTOG)
*     S(1)   CONX      conexion 0 Dy, 1 Yy, 2 Yd, 3 Dd notese que 2 y 3 dan 
*                      sistema aislado y corriente a tierra nula
*     S(2)   Z0   R0   IMPEDANCIA DEL TRANSFORMADOR EN 0/0
*                      O RESISTENCIA DE LA LINEA CONSIDERANDO
*                      LOS PARALELOS DE LAS FASES
*     S(3)   P0   X0   PERDIDAS DEL TRANSFORMADOR EN 0/0
*                      O REACTANCIA DE LA LINEA EN OHM/KM
*     S(4)   ZG   IMPEDANCIA DE PUESTA A TIERRA EN OHM
*     S(5)   XSRG RELACION REACTANCIA RESISTENCIA
*  CON   ><  Z0 R0 ><  P0 X0 ><    ZG  ><  XSRG  >
*
*     TARJETA CC (DATOS USADOS POR FLURAD)
*     PCAR  POTENCIA ACTIVA EN EL NODO
*     QCAR  POTENCIA REACTIVA EN EL NODO (+ INDUCTIVA - CAPACITIVA)
*     GCAR  CONDUCTANCIA DERIVACION EN EL NODO DADA COMO POTENCIA
*           EN EL NODO A TENSION NOMINAL
*     BCAR  SUSCEPTANCIA DERIVACION EN EL NODO (+ INDUCTIVA - CAPACITIVA)
*     UCAR  TENSION DE REFERENCIA DEL NODO
*     U1    TENSION ESPERADA EN VALOR RELATIVO
*   PCAR ><  QCAR  ><  GCAR  ><  BCAR  ><  UCAR  ><   U1   >
*
*     TARJETA DD (DATOS USADOS POR CORTOA, TIE es usado por CORTOR y CORTOG)
*     PCAR  RELLENO
*     QCAR  RELLENO
*     GCAR  POTENCIA ACTIVA QUE REPRESENTA EL ARRANQUE DE LOS MOTORES
*           EN EL NODO A TENSION NOMINAL
*     BCAR  POTENCIA REACTIVA
*     UCAR  TENSION DE REFERENCIA DEL NODO
*     U1    TENSION ESPERADA EN VALOR RELATIVO
*     TIE   TIEMPO DE INTERRUPCION EN SEG (USADO POR CORTOR, CORTOG, CORTOA)
*  00000 >< 000000 ><  GCAR  ><  BCAR  ><  UCAR  ><   U1   ><   TIE  >
*
* LA EJECUCION REINICIA CON UNA TARJETA EN BLANCO
*              FINALIZA POR FALTA DE DATOS
*
* el programa inicia pidiendo se 
* indique el calculo que desea hacer:
*     1 - FLUJO DE CARGA    (FLURAD)
*     2 - CORTO RADIAL      (CORTOR)
*     3 - CORTO CON APORTES (CORTOA)
*     4 - CORTO A TIERRA    (CORTOG)
*     5 - TABLAS DE DATOS   (REDRAD)
*     6 - TERMINAR....   
*
* luego el programa pide el nombre del archivo de datos
* los dos archivos de salida (resultados), y luego          
* Fija el tipo de informe:
*     0 tabla con valores separados por punto y coma (extension CSV)
*     1 texto encabezado solo en la primera hoja
*     2 texto encabezado en todas las hojas
* ---------------------------------------------------------------------