radarsys Commit - r390:5670a70fd8d1 · Jicamarca Repository

Se resolvió plot pattern de ABS

Renato Huallpa -

r390:5670a70fd8d1

parent child

Context file:

r390:5670a70fd8d1

Collapse all files

.env +2 -2

             REDIS_HOST=radarsys-redis
             REDIS_PORT=6300
             POSTGRES_PORT_5432_TCP_ADDR=radarsys-postgres
             POSTGRES_PORT_5432_TCP_PORT=5432
             #POSTGRES_HOST=postgres
             POSTGRES_HOST=localhost
             POSTGRES_USER=docker
             POSTGRES_PASSWORD=docker
             POSTGRES_DB=radarsys
             # DB_NAME=radarsys
             # DB_USER=docker
             # DB_PASSWORD=docker
             PGDATA=/var/lib/postgresql/data
             LC_ALL=C.UTF-8
             TZ=America/Lima
             DOCKER_DATA=/data/dockers/radarsys/
             LOCAL_IP=192.168.1.128
-            MQTT_SERVER=10.10.10.200
+            # MQTT_SERVER=10.10.10.200
-            # MQTT_SERVER = 192.168.100.5
+            MQTT_SERVER = 192.168.100.5
             MQTT_PORT = 1883
             MQTT_KEEPALIVE = 3660
             MQTT_USER_ATRAD=atrad
             MQTT_PASSWORD_ATRAD = atrad
             MQTT_USER = abs
             MQTT_PASSWORD = abs
             MQTT_CLIENT_ID= abs_id
             TOPIC_ABS=abs/beams
             TOPIC_ABS_ACK=abs/beams_ack
             TOPIC_ABS_CHANGE=abs/change_beam

apps/abs/models.py +2 0

             from django.db import models
             from apps.main.models import Configuration , User
             from django.urls import reverse
             from celery.execute import send_task
             from datetime import datetime
             import ast
             import socket
             import json
             import requests
             import struct
             import os, sys, time
             from .mqtt import client as mqtt_client
             from radarsys.socketconfig import sio as sio
             import json
             antenna_default = json.dumps({
                                 "antenna_up": [[0.0,0.0,0.0,0.0,0.5,0.5,0.5,0.5],
                                                 [0.0,0.0,0.0,0.0,0.5,0.5,0.5,0.5],
                                                 [0.0,0.0,0.0,0.0,0.5,0.5,0.5,0.5],
                                                 [0.0,0.0,0.0,0.0,0.5,0.5,0.5,0.5],
                                                 [0.5,0.5,0.5,0.5,1.0,1.0,1.0,1.0],
                                                 [0.5,0.5,0.5,0.5,1.0,1.0,1.0,1.0],
                                                 [0.5,0.5,0.5,0.5,1.0,1.0,1.0,1.0],
                                                 [0.5,0.5,0.5,0.5,1.0,1.0,1.0,1.0]
                                                 ]
                                                ,
                                 "antenna_down": [[0.0,0.0,0.0,0.0,0.5,0.5,0.5,0.5],
                                                 [0.0,0.0,0.0,0.0,0.5,0.5,0.5,0.5],
                                                 [0.0,0.0,0.0,0.0,0.5,0.5,0.5,0.5],
                                                 [0.0,0.0,0.0,0.0,0.5,0.5,0.5,0.5],
                                                 [0.5,0.5,0.5,0.5,3.0,3.0,3.0,3.0],
                                                 [0.5,0.5,0.5,0.5,3.0,3.0,3.0,3.0],
                                                 [0.5,0.5,0.5,0.5,3.0,3.0,3.0,3.0],
                                                 [0.5,0.5,0.5,0.5,3.0,3.0,3.0,3.0]],
                             })
             tx_default = json.dumps({
                             "up": [[1,1,1,1,0,0,0,0],
                                    [1,1,1,1,0,0,0,0],
                                    [1,1,1,1,0,0,0,0],
                                    [1,1,1,1,0,0,0,0],
                                    [0,0,0,0,1,1,1,1],
                                    [0,0,0,0,1,1,1,1],
                                    [0,0,0,0,1,1,1,1],
                                    [0,0,0,0,1,1,1,1]],
                           "down": [[1,1,1,1,0,0,0,0],
                                    [1,1,1,1,0,0,0,0],
                                    [1,1,1,1,0,0,0,0],
                                    [1,1,1,1,0,0,0,0],
                                    [0,0,0,0,1,1,1,1],
                                    [0,0,0,0,1,1,1,1],
                                    [0,0,0,0,1,1,1,1],
                                    [0,0,0,0,1,1,1,1]],
                          })
             rx_default = json.dumps({
                             "up": [[1,1,1,1,0,0,0,0],
                                    [1,1,1,1,0,0,0,0],
                                    [1,1,1,1,0,0,0,0],
                                    [1,1,1,1,0,0,0,0],
                                    [0,0,0,0,1,1,1,1],
                                    [0,0,0,0,1,1,1,1],
                                    [0,0,0,0,1,1,1,1],
                                    [0,0,0,0,1,1,1,1]],
                           "down": [[1,1,1,1,0,0,0,0],
                                    [1,1,1,1,0,0,0,0],
                                    [1,1,1,1,0,0,0,0],
                                    [1,1,1,1,0,0,0,0],
                                    [0,0,0,0,1,1,1,1],
                                    [0,0,0,0,1,1,1,1],
                                    [0,0,0,0,1,1,1,1],
                                    [0,0,0,0,1,1,1,1]],
                          })
             status_default = '0000000000000000000000000000000000000000000000000000000000000000'
             default_messages = {}
             for i in range(1,65):
                 default_messages[str(i)] = "Module "+str(i)
             ues_default = json.dumps({
                             "up": [0.533333,0.00000,1.06667,0.00000],
                             "down": [0.533333,0.00000,1.06667,0.00000]
                         })
             onlyrx_default = json.dumps({
                                 "up": False,
                                 "down": False
                             })
             def up_convertion(cadena):
                 valores = []
                 for c in cadena:
                     if c == 1.0: valores=valores+['000']
                     if c == 2.0: valores=valores+['001']
                     if c == 3.0: valores=valores+['010']
                     if c == 0.0: valores=valores+['011']
                     if c == 0.5: valores=valores+['100']
                     if c == 1.5: valores=valores+['101']
                     if c == 2.5: valores=valores+['110']
                     if c == 3.5: valores=valores+['111']
                 return valores
             def up_conv_bits(value):
                 if value == 1.0: bits="000"
                 if value == 2.0: bits="001"
                 if value == 3.0: bits="010"
                 if value == 0.0: bits="011"
                 if value == 0.5: bits="100"
                 if value == 1.5: bits="101"
                 if value == 2.5: bits="110"
                 if value == 3.5: bits="111"
                 return bits
             def down_convertion(cadena):
                 valores = []
                 for c in cadena:
                     if c == 1.0: valores=valores+['000']
                     if c == 2.0: valores=valores+['001']
                     if c == 3.0: valores=valores+['010']
                     if c == 0.0: valores=valores+['011']
                     if c == 0.5: valores=valores+['100']
                     if c == 1.5: valores=valores+['101']
                     if c == 2.5: valores=valores+['110']
                     if c == 3.5: valores=valores+['111']
                 return valores
             def down_conv_bits(value):
                 if value == 1.0: bits="000"
                 if value == 2.0: bits="001"
                 if value == 3.0: bits="010"
                 if value == 0.0: bits="011"
                 if value == 0.5: bits="100"
                 if value == 1.5: bits="101"
                 if value == 2.5: bits="110"
                 if value == 3.5: bits="111"
                 return bits
             def up_conv_value(bits):
                 if bits == "000": value=1.0
                 if bits == "001": value=2.0
                 if bits == "010": value=3.0
                 if bits == "011": value=0.0
                 if bits == "100": value=0.5
                 if bits == "101": value=1.5
                 if bits == "110": value=2.5
                 if bits == "111": value=3.5
                 return value
             def down_conv_value(bits):
                 if bits == "000": value=1.0
                 if bits == "001": value=2.0
                 if bits == "010": value=3.0
                 if bits == "011": value=0.0
                 if bits == "100": value=0.5
                 if bits == "101": value=1.5
                 if bits == "110": value=2.5
                 if bits == "111": value=3.5
                 return value
             def ip2position(module_number):
                 j=0
                 i=0
                 for x in range(0,module_number-1):
                     j=j+1
                     if j==8:
                         i=i+1
                         j=0
                 pos = [i,j]
                 return pos
             def fromBinary2Char(binary_string):
                 number = int(binary_string, 2)
                 #Plus 33 to avoid more than 1 characters values such as: '\x01'-'\x1f'
                 number = number + 33
                 char = chr(number)
                 return char
             def fromChar2Binary(char):
                 number = ord(char) - 33
                 #Minus 33 to get the real value
                 bits = bin(number)[2:]
                 #To ensure we have a string with 6bits
                 if len(bits) < 6:
                     bits = bits.zfill(6)
                 return bits
             OPERATION_MODES = (
                              (0, 'Manual'),
                              (1, 'Automatic'),
                          )
             class ABSConfiguration(Configuration):
                 active_beam     = models.PositiveSmallIntegerField(verbose_name='Active Beam', default=0)
                 module_status   = models.CharField(verbose_name='Module Status', max_length=10000, default=status_default)
                 operation_mode  = models.PositiveSmallIntegerField(verbose_name='Operation Mode', choices=OPERATION_MODES, default = 0)
                 operation_value = models.FloatField(verbose_name='Periodic (seconds)', default="10", null=True, blank=True)
                 module_messages = models.CharField(verbose_name='Modules Messages', max_length=10000, default=json.dumps(default_messages))
                 class Meta:
                     db_table = 'abs_configurations'
                 def get_absolute_url_plot(self):
                     return reverse('url_plot_abs_patterns', args=[str(self.id)])
                 def parms_to_dict(self):
                     parameters = {}
                     parameters['device_id'] = self.device.id
                     parameters['label']      = self.label
                     parameters['device_type']      = self.device.device_type.name
                     parameters['beams']     = {}
                     beams = ABSBeam.objects.filter(abs_conf=self)
                     b=1
                     for beam in beams:
                         #absbeam = ABSBeam.objects.get(pk=beams[beam])
                         parameters['beams']['beam'+str(b)] = beam.parms_to_dict()#absbeam.parms_to_dict()
                         b+=1
                     return parameters
                 def dict_to_parms(self, parameters):
                     self.label = parameters['label']
                     absbeams  = ABSBeam.objects.filter(abs_conf=self)
                     beams     = parameters['beams']
                     if absbeams:
                         beams_number    = len(beams)
                         absbeams_number = len(absbeams)
                         if beams_number==absbeams_number:
                             i = 1
                             for absbeam in absbeams:
                                 absbeam.dict_to_parms(beams['beam'+str(i)])
                                 i = i+1
                         elif beams_number > absbeams_number:
                             i = 1
                             for absbeam in absbeams:
                                 absbeam.dict_to_parms(beams['beam'+str(i)])
                                 i=i+1
                             for x in range(i,beams_number+1):
                                 new_beam = ABSBeam(
                                            name     =beams['beam'+str(i)]['name'],
                                            antenna  =json.dumps(beams['beam'+str(i)]['antenna']),
                                            abs_conf = self,
                                            tx       =json.dumps(beams['beam'+str(i)]['tx']),
                                            rx       =json.dumps(beams['beam'+str(i)]['rx']),
                                            ues      =json.dumps(beams['beam'+str(i)]['ues']),
                                            only_rx  =json.dumps(beams['beam'+str(i)]['only_rx'])
                                            )
                                 new_beam.save()
                                 i=i+1
                         else: #beams_number < absbeams_number:
                             i = 1
                             for absbeam in absbeams:
                                 if i <= beams_number:
                                     absbeam.dict_to_parms(beams['beam'+str(i)])
                                     i=i+1
                                 else:
                                     absbeam.delete()
                     else:
                         for beam in beams:
                             new_beam = ABSBeam(
                                        name     =beams[beam]['name'],
                                        antenna  =json.dumps(beams[beam]['antenna']),
                                        abs_conf = self,
                                        tx       =json.dumps(beams[beam]['tx']),
                                        rx       =json.dumps(beams[beam]['rx']),
                                        ues      =json.dumps(beams[beam]['ues']),
                                        only_rx  =json.dumps(beams[beam]['only_rx'])
                                        )
                             new_beam.save()
                 def update_from_file(self, parameters):
                     self.dict_to_parms(parameters)
                     self.save()
                 def get_beams(self, **kwargs):
                     '''
                     This function returns ABS Configuration beams
                     '''
                     return ABSBeam.objects.filter(abs_conf=self.pk, **kwargs)
                 def clone(self, **kwargs):
                     beams = self.get_beams()
                     self.pk = None
                     self.id = None
                     for attr, value in kwargs.items():
                         setattr(self, attr, value)
                     self.save()
                     for beam in beams:
                         beam.clone(abs_conf=self)
                     #-----For Active Beam-----
                     new_beams = ABSBeam.objects.filter(abs_conf=self)
                     self.active_beam = new_beams[0].id
                     self.save()
                     #-----For Active Beam-----
                     #-----For Device Status---
                     self.device.status = 3
                     self.device.save()
                     #-----For Device Status---
                     return self
                 def start_device(self):
                     if self.device.status == 3:
                         try:
                             #self.write_device()
                             send_task('task_change_beam', [self.id],)
+                            print("*************************RUNNING ABS**************************",flush=True)
                             self.message = 'ABS running'
                         except Exception as e:
                             self.message = str(e)
                             return False
                         return True
                     else:
                         self.message = 'Please, select Write ABS Device first.'
                         return False
                 def stop_device(self):
                     self.device.status = 2
                     self.device.save()
                     self.message = 'ABS has been stopped.'
+                    print("*************************STOPPED ABS**************************",flush=True)
                     self.save()
                     return True
                 def stop_device_mqtt(self):
                     self.device.status = 2
                     self.device.save()
                     self.message = 'ABS has been stopped.'
                     self.save()
                     mqtt_client.publish(os.environ.get('TOPIC_ABS', 'abs/beams'),"STOP")
                     return True
                 def write_device(self):
                     """
                     This function sends the beams list to every abs module.
                     It needs 'module_conf' function
                     """
                     print("Write 3")
                     beams  = ABSBeam.objects.filter(abs_conf=self)
                     nbeams = len(beams)
                     # Se manda a cero RC para poder realizar cambio de beam
                     if self.experiment is None:
                         confs = []
                     else:
                         confs = Configuration.objects.filter(experiment = self.experiment).filter(type=0)
                     confdds  = ''
                     confjars = ''
                     confrc   = ''
                     #TO STOP DEVICES: DDS-JARS-RC
                     for i in range(0,len(confs)):
                         if i==0:
                             for conf in confs:
                                 if conf.device.device_type.name == 'dds':
                                     confdds = conf
                                     confdds.stop_device()
                                     break
                         if i==1:
                             for conf in confs:
                                 if conf.device.device_type.name == 'jars':
                                     confjars = conf
                                     confjars.stop_device()
                                     break
                         if i==2:
                             for conf in confs:
                                 if conf.device.device_type.name == 'rc':
                                     confrc = conf
                                     confrc.stop_device()
                                     break
                     '''
                     if self.connected_modules() == 0 :
                         print("No encuentra modulos")
                         self.message = "No ABS Module detected."
                         return False
                     '''
                     #-------------Write each abs module-----------
                     if beams:
                         block_id = 0
                         message = 'SNDF{:03d}{:02d}{:02d}'.format(nbeams, nbeams, block_id)
                         for i, status in enumerate(self.module_status):
                             message += ''.join([fromBinary2Char(beam.module_6bits(i)) for beam in beams])
                         status = ['0'] * 64
                         n = 0
                         print("Llega una antes entrar a multicast4")
                         sock = self.send_multicast(message)
                         while True:
                         #for i in range(32):
                             try:
                                 data, address = sock.recvfrom(1024)
                                 print (address, data)
                                 data = data.decode("utf-8")
                                 if data == '1':
                                     status[int(address[0][10:])-1] = '3'
                                     #print (int(address[0][10:])-1)
                                 elif data == '0':
                                     status[int(address[0][10:])-1] = '1'
                             except socket.timeout:
                                 print('Timeout')
                                 break
                             except Exception as e:
                                 print ('Error {}'.format(e))
                                 n += 1
                         sock.close()
                     else:
                         self.message = "ABS Configuration does not have beams"
                         print('No beams')
                         #Start DDS-RC-JARS
                         if confdds:
                             confdds.start_device()
                         if confrc:
                             #print confrc
                             confrc.start_device()
                         if confjars:
                             confjars.start_device()
                         return False
                     if n == 64:
                         self.message = "Could not write ABS Modules"
                         self.device.status = 0
                         self.module_status = ''.join(status)
                         self.save()
                         print('Could not write ABS')
                         #Start DDS-RC-JARS
                         if confdds:
                             confdds.start_device()
                         if confrc:
                             #print confrc
                             confrc.start_device()
                         if confjars:
                             confjars.start_device()
                         return False
                     else:
                         self.message = "ABS Beams List have been sent to ABS Modules"
                         print('ABS beams list sent')
                         self.active_beam = beams[0].pk
                     #Start DDS-RC-JARS
                     if confdds:
                         confdds.start_device()
                     if confrc:
                         #print confrc
                         confrc.start_device()
                     if confjars:
                         confjars.start_device()
                     print('Inicia intento de salvar device.status')
                     self.device.status = 3
                     self.module_status = ''.join(status)
                     #print(status)
                     self.save()
                     print('Estatus salvado')
                     conf_active, __ = ABSActive.objects.get_or_create(pk=1)
                     conf_active.conf = self
                     conf_active.save()
                     return True
                 def write_device_mqtt(self):
                     if self.experiment is None:
                         confs = []
                     else:
                         confs = Configuration.objects.filter(experiment = self.experiment).filter(type=0)
                     confdds  = ''
                     confjars = ''
                     confrc   = ''
                     #TO STOP DEVICES: DDS-JARS-RC
                     for i in range(0,len(confs)):
                         if i==0:
                             for conf in confs:
                                 if conf.device.device_type.name == 'dds':
                                     confdds = conf
                                     confdds.stop_device()
                                     break
                         if i==1:
                             for conf in confs:
                                 if conf.device.device_type.name == 'jars':
                                     confjars = conf
                                     confjars.stop_device()
                                     break
                         if i==2:
                             for conf in confs:
                                 if conf.device.device_type.name == 'rc':
                                     confrc = conf
                                     confrc.stop_device()
                                     break
                     apuntes_up_down=''
                     beams  = ABSBeam.objects.filter(abs_conf=self)
                     inicializacion="{\"beams\":["
                     finalizacion="]}"
                     for beam in beams:
                         beam.antenna=beam.antenna[1:]
                         info="{\"id\":"+str(beam.id)+","+beam.antenna + ","
                         apuntes_up_down=apuntes_up_down+info
                     apuntes_up_down=apuntes_up_down[:len(apuntes_up_down)-1]
                     apuntes_up_down=inicializacion+ apuntes_up_down+finalizacion
                     mqtt_client.publish(os.environ.get('TOPIC_ABS', 'abs/beams'),apuntes_up_down)
                     #Start DDS-RC-JARS
                     if confdds:
                         confdds.start_device()
                     if confrc:
                         #print confrc
                         confrc.start_device()
                     if confjars:
                         confjars.start_device()
                     return True
                 def read_module(self, module):
                     """
                     Read out-bits (up-down) of 1 abs module NOT for Configuration
                     """
                     ip_address  = self.device.ip_address
                     ip_address  = ip_address.split('.')
                     module_seq  = (ip_address[0],ip_address[1],ip_address[2])
                     dot         = '.'
                     module_ip   = dot.join(module_seq)+'.'+str(module)
                     module_port = self.device.port_address
                     read_route = 'http://'+module_ip+':'+str(module_port)+'/read'
                     module_status = json.loads(self.module_status)
                     print(read_route)
                     module_bits = ''
                     try:
                         r_read      = requests.get(read_route, timeout=0.5)
                         answer      = r_read.json()
                         module_bits = answer['allbits']
                     except:
                         return {}
                     return module_bits
                 def read_device(self):
                     parms = {}
                     # Reads active modules.
                     module_status = json.loads(self.module_status)
                     total = 0
                     for status in module_status:
                         if module_status[status] != 0:
                             module_bits = self.read_module(int(status))
                             bits={}
                             if module_bits:
                                 bits = (str(module_bits['um2']) + str(module_bits['um1']) + str(module_bits['um0']) +
                                         str(module_bits['dm2']) + str(module_bits['dm1']) + str(module_bits['dm0']) )
                                 parms[str(status)] = bits
                             total +=1
                     if total==0:
                         self.message = "No ABS Module detected. Please select 'Status'."
                         return False
                     self.message = "ABS Modules have been read"
                     #monitoreo_tx = JROABSClnt_01CeCnMod000000MNTR10
                     return parms
                 def connected_modules(self):
                     """
                     This function returns the number of connected abs-modules without updating.
                     """
                     num = 0
                     print(self.module_status)
                     for i, status in  enumerate(self.module_status):
                         if status != '0':
                             num += 1
                             #print('status {}:{}'.format(i+1, status))
                     return num
                 def send_multicast(self, message):
                     #print("Send multicast")
                     multicast_group = ('224.3.29.71', 10000)
                     # Create the datagram socket
                     sock = socket.socket(socket.AF_INET, socket.SOCK_DGRAM)
                     sock.settimeout(1)
                     local_ip = os.environ.get('LOCAL_IP', '0.0.0.0')
                     local_ip = '0.0.0.0'
                     print("He llegado a IP local")
                     sock.setsockopt(socket.IPPROTO_IP, socket.IP_MULTICAST_IF, socket.inet_aton(local_ip))
                     sock.sendto(message.encode(), multicast_group)
                     print('Sending ' + message)
                     return sock
                 def status_device(self):
                     """
                     This function returns the status of all abs-modules as one.
                     If at least one module is connected, its answer is "1"
                     """
                     print ('Status device')
                     print (self.active_beam)
                     beams = ABSBeam.objects.filter(abs_conf=self)
                     #print beams[self.active_beam-1].module_6bits(0)
                     active = ABSActive.objects.get(pk=1)
                     if active.conf != self:
                         self.message = 'La configuracion actual es la del siguiente enlace %s.' % active.conf.get_absolute_url()
                         self.message +=  "\n"
                         self.message += 'Se debe realizar un write en esta configuracion para luego obtener un status valido.'
                         return False
                     sock = self.send_multicast('MNTR')
                     n = 0
                     status = ['0'] * 64
                     while True:
                     #for i in range(32):
                         #if True:
                         try:
                             print("Recibiendo")
                             address = None
                             data, address = sock.recvfrom(2)
                             print (address, data)
                             print("!!!!")
                             data = data.decode()
                             aux_mon = "1"
                             aux_expected = aux_mon
                             if(len(data)==2):
                                 print ("data[1]: ")
                                 print (data[1])
                                 aux_mon = fromChar2Binary(data[1])
                                 print (aux_mon)
                                 aux_i = (str(address[0]).split('.'))[3]
                                 print (aux_i)
                                 print ('Active beam')
                                 beam_active = ABSBeam.objects.get(pk=self.active_beam)
                                 print (beam_active)
                                 aux_expected = beam_active.module_6bits(int(aux_i)-1)
                                 print (aux_expected)
                             print ("data[0]: ")
                             print (data[0])
                             if data[0] == '1':
                                 status[int(address[0][10:])-1] = '3'
                                 if aux_mon == aux_expected:
                                     print ('Es igual')
                                 else:
                                     print ('Es diferente')
                                     status[int(address[0][10:])-1] = '2'
                             elif data[0] == '0':
                                 status[int(address[0][10:])-1] = '1'
                             n += 1
                             print('Module: {} connected'.format(address))
                         except socket.timeout:
                             print('Timeout')
                             break
                         except:
                             print('Module: {} error'.format(address))
                             pass
                     sock.close()
                     if n > 0:
                         self.message = 'ABS modules Status have been updated.'
                         self.device.status = 1
                     else:
                         self.device.status = 0
                         self.message = 'No ABS module is connected.'
                     self.module_status = ''.join(status)
                     self.save()
                     return self.device.status
                 def send_beam(self, beam_pos):
                     """
                     This function connects to a multicast group and sends the beam number
                     to all abs modules.
                     """
                     print ('Send beam')
                     print (self.active_beam)
                     beams = ABSBeam.objects.filter(abs_conf=self)
                     #print beams[self.active_beam-1].module_6bits(0)
                     active = ABSActive.objects.get(pk=1)
                     if active.conf != self:
                         self.message = 'La configuracion actual es la del siguiente enlace %s.' % active.conf.get_absolute_url()
                         self.message +=  "\n"
                         self.message += 'Se debe realizar un write en esta configuracion para luego obtener un status valido.'
                         return False
                     # Se manda a cero RC para poder realizar cambio de beam
                     if self.experiment is None:
                         confs = []
                     else:
                         confs = Configuration.objects.filter(experiment = self.experiment).filter(type=0)
                     confdds  = ''
                     confjars = ''
                     confrc   = ''
                     #TO STOP DEVICES: DDS-JARS-RC
                     for i in range(0,len(confs)):
                         if i==0:
                             for conf in confs:
                                 if conf.device.device_type.name == 'dds':
                                     confdds = conf
                                     confdds.stop_device()
                                     break
                         if i==1:
                             for conf in confs:
                                 if conf.device.device_type.name == 'jars':
                                     confjars = conf
                                     confjars.stop_device()
                                     break
                         if i==2:
                             for conf in confs:
                                 if conf.device.device_type.name == 'rc':
                                     confrc = conf
                                     confrc.stop_device()
                                     break
                     if beam_pos > 0:
                         beam_pos = beam_pos - 1
                     else:
                         beam_pos = 0
                     #El indice del apunte debe ser menor que el numero total de apuntes
                     #El servidor tcp en el embebido comienza a contar desde 0
                     status = ['0'] * 64
                     message = 'CHGB{}'.format(beam_pos)
                     sock = self.send_multicast(message)
                     while True:
                     #for i in range(32):
                         try:
                             data, address = sock.recvfrom(1024)
                             print (address, data)
                             data = data.decode()
                             if data == '1':
                                 status[int(address[0][10:])-1] = '3'
                             elif data == '0':
                                 status[int(address[0][10:])-1] = '1'
                         except socket.timeout:
                             print('Timeout')
                             break
                         except  Exception as e:
                             print ('Error {}'.format(e))
                             pass
                     sock.close()
                     #Start DDS-RC-JARS
                     if confdds:
                         confdds.start_device()
                     if confrc:
                         #print confrc
                         confrc.start_device()
                     if confjars:
                         confjars.start_device()
                     self.message = "ABS Beam has been changed"
                     self.module_status = ''.join(status)
                     self.save()
                     return True
                 def get_absolute_url_import(self):
                     return reverse('url_import_abs_conf', args=[str(self.id)])
             class ABSActive(models.Model):
                 conf = models.ForeignKey(ABSConfiguration, null=True, verbose_name='ABS Configuration', on_delete=models.CASCADE)
                 class Meta:
                     db_table = 'abs_absactive'
             class ABSBeam(models.Model):
                 name         = models.CharField(max_length=60, default='Beam')
                 antenna      = models.CharField(verbose_name='Antenna', max_length=1000, default=antenna_default)
                 abs_conf     = models.ForeignKey('ABSConfiguration', null=True,
                 verbose_name='ABS Configuration', on_delete=models.CASCADE)
                 tx           = models.CharField(verbose_name='Tx', max_length=1000, default=tx_default)
                 rx           = models.CharField(verbose_name='Rx', max_length=1000, default=rx_default)
                 s_time       = models.TimeField(verbose_name='Star Time', default='00:00:00')
                 e_time       = models.TimeField(verbose_name='End Time', default='23:59:59')
                 ues          = models.CharField(verbose_name='Ues', max_length=100, default=ues_default)
                 only_rx      = models.CharField(verbose_name='Only RX', max_length=40, default=onlyrx_default)
                 class Meta:
                     db_table = 'abs_beams'
                 def __unicode__(self):
                     return u'%s' % (self.name)
                 def parms_to_dict(self):
                     parameters = {}
                     parameters['name']          = self.name
                     parameters['antenna']       = ast.literal_eval(self.antenna)
                     parameters['abs_conf']      = self.abs_conf.name
                     parameters['tx']            = ast.literal_eval(self.tx)
                     parameters['rx']            = ast.literal_eval(self.rx)
                     parameters['s_time']        = self.s_time.strftime("%H:%M:%S")
                     parameters['e_time']        = self.e_time.strftime("%H:%M:%S")
                     parameters['ues']           = ast.literal_eval(self.ues)
                     parameters['only_rx']       = json.loads(self.only_rx)
                     return parameters
                 def dict_to_parms(self, parameters):
                     self.name     = parameters['name']
                     self.antenna  = json.dumps(parameters['antenna'])
                     #self.abs_conf = parameters['abs_conf']
                     self.tx       = json.dumps(parameters['tx'])
                     self.rx       = json.dumps(parameters['rx'])
                     #self.s_time = parameters['s_time']
                     #self.e_time = parameters['e_time']
                     self.ues      = json.dumps(parameters['ues'])
                     self.only_rx     = json.dumps(parameters['only_rx'])
                     self.save()
                 def clone(self, **kwargs):
                     self.pk = None
                     self.id = None
                     for attr, value in kwargs.items():
                         setattr(self, attr, value)
                     self.save()
                     return self
                 def module_6bits(self, module):
                     """
                     This function reads antenna pattern and choose 6bits (upbits-downbits) for one abs module
                     """
                     module += 1
                     if module > 64:
                         beam_bits = ""
                         return beam_bits
                     data      = ast.literal_eval(self.antenna)
                     up_data   = data['antenna_up']
                     down_data = data['antenna_down']
                     pos        = ip2position(module)
                     up_value   = up_data[pos[0]][pos[1]]
                     down_value = down_data[pos[0]][pos[1]]
                     up_bits   = up_conv_bits(up_value)
                     down_bits = down_conv_bits(down_value)
                     beam_bits = up_bits+down_bits
                     return beam_bits
                 @property
                 def get_upvalues(self):
                     """
                     This function reads antenna pattern and show the up-value of one abs module
                     """
                     data      = ast.literal_eval(self.antenna)
                     up_data   = data['antenna_up']
                     up_values = []
                     for data in up_data:
                         for i in range(0,8):
                             up_values.append(data[i])
                     return up_values
                 @property
                 def antenna_upvalues(self):
                     """
                     This function reads antenna pattern and show the up - values of one abs beam
                     in a particular order
                     """
                     data      = ast.literal_eval(self.antenna)
                     up_data   = data['antenna_up']
                     return up_data
                 @property
                 def antenna_downvalues(self):
                     """
                     This function reads antenna pattern and show the down - values of one abs beam
                     in a particular order
                     """
                     data      = ast.literal_eval(self.antenna)
                     down_data   = data['antenna_down']
                     return down_data
                 @property
                 def get_downvalues(self):
                     """
                     This function reads antenna pattern and show the down-value of one abs module
                     """
                     data      = ast.literal_eval(self.antenna)
                     down_data = data['antenna_down']
                     down_values = []
                     for data in down_data:
                         for i in range(0,8):
                             down_values.append(data[i])
                     return down_values
                 @property
                 def get_up_ues(self):
                     """
                     This function shows the up-ues-value of one beam
                     """
                     data      = ast.literal_eval(self.ues)
                     up_ues = data['up']
                     return up_ues
                 @property
                 def get_down_ues(self):
                     """
                     This function shows the down-ues-value of one beam
                     """
                     data      = ast.literal_eval(self.ues)
                     down_ues = data['down']
                     return down_ues
                 @property
                 def get_up_onlyrx(self):
                     """
                     This function shows the up-onlyrx-value of one beam
                     """
                     data      = json.loads(self.only_rx)
                     up_onlyrx = data['up']
                     return up_onlyrx
                 @property
                 def get_down_onlyrx(self):
                     """
                     This function shows the down-onlyrx-value of one beam
                     """
                     data      = json.loads(self.only_rx)
                     down_onlyrx = data['down']
                     return down_onlyrx
                 @property
                 def get_tx(self):
                     """
                     This function shows the tx-values of one beam
                     """
                     data = json.loads(self.tx)
                     return data
                 @property
                 def get_uptx(self):
                     """
                     This function shows the up-tx-values of one beam
                     """
                     data = json.loads(self.tx)
                     up_data = data['up']
                     up_values = []
                     for data in up_data:
                         for i in range(0,8):
                             up_values.append(data[i])
                     return up_values
                 @property
                 def get_downtx(self):
                     """
                     This function shows the down-tx-values of one beam
                     """
                     data = json.loads(self.tx)
                     down_data = data['down']
                     down_values = []
                     for data in down_data:
                         for i in range(0,8):
                             down_values.append(data[i])
                     return down_values
                 @property
                 def get_rx(self):
                     """
                     This function shows the rx-values of one beam
                     """
                     data = json.loads(self.rx)
                     return data
                 @property
                 def get_uprx(self):
                     """
                     This function shows the up-rx-values of one beam
                     """
                     data = json.loads(self.rx)
                     up_data = data['up']
                     up_values = []
                     for data in up_data:
                         for i in range(0,8):
                             up_values.append(data[i])
                     return up_values
                 @property
                 def get_downrx(self):
                     """
                     This function shows the down-rx-values of one beam
                     """
                     data = json.loads(self.rx)
                     down_data = data['down']
                     down_values = []
                     for data in down_data:
                         for i in range(0,8):
                             down_values.append(data[i])
                     return down_values

apps/abs/templates/abs_downpattern_img.html +1 -5

             {% load static %}
             {% load django_bootstrap5 %}
             {% load main_tags %}
-            {% block content %}
+            <img id="imgMainDown" src="{% url 'url_plot_beam' beam.abs_conf.id beam.id 'down' %}" alt="Error ploting..." style="width: 250px;height:250px; background-image: url({% static 'images/loader.gif' %});background-repeat: no-repeat;background-position: 50% 50%;">
-            <div id="PictureOverJRODown" style="float: right">
-              <img id="imgMainDown" src="{% url 'url_plot_beam' beam.abs_conf.id beam.id 'down' %}" alt="Error ploting..." style="width:340px;height:340px; background-image: url({% static 'images/loader.gif' %});background-repeat: no-repeat;background-position: 50% 50%;">
-            </div>
-            {% endblock %}

apps/abs/templates/abs_pattern.html +25 -7

             {% load static %}
             {% load django_bootstrap5 %}
             {% load main_tags %}
             {% block content %}
             <style>
             .abs {
                 border: 2px solid #00334d;
                 vertical-align: center;
                 display: inline-block;
                 font-size: 95%;
             }
             .abs tr{
                 border-bottom: 0px solid #00334d;
             }
             .abs td{
                 border-right: 0px solid #00334d;
                 text-align: center;
                 padding: 4px;
             }
             .pattern {
               border: 2px solid #00334d;
               vertical-align: center;
               font-weight: bold;
               text-align: center;
             }
             .pattern tr{
                 border: 1px solid #ffffff;
                 background-color: #ecf0f1;
             }
             .pattern td{
                 border: 2px solid #e2e2e7;
                 text-align: center;
             }
             .north_quarter{
               border: 2px solid #00334d;
               vertical-align: center;
               font-weight: bold;
             }
             .north_quarter tr{
                 border: 1px solid #ffffff;
                 background-color: #ecf0f1;
             }
             .north_quarter td{
                 border: 2px solid #e2e2e7;
                 text-align: center;
             }
             .east_quarter{
               border: 2px solid #00334d;
               vertical-align: center;
               font-weight: bold;
             }
             .east_quarter tr{
                 border: 1px solid #ffffff;
                 background-color: #ecf0f1;
             }
             .east_quarter td{
                 border: 2px solid #e2e2e7;
                 text-align: center;
             }
             .west_quarter{
               border: 2px solid #00334d;
               vertical-align: center;
               font-weight: bold;
             }
             .west_quarter tr{
                 border: 1px solid #ffffff;
                 background-color: #ecf0f1;
             }
             .west_quarter td{
                 border: 2px solid #e2e2e7;
                 text-align: center;
             }
             .south_quarter{
               border: 2px solid #00334d;
               vertical-align: center;
               font-weight: bold;
             }
             .south_quarter tr{
                 border: 1px solid #ffffff;
                 background-color: #ecf0f1;
             }
             .south_quarter td{
                 border: 2px solid #e2e2e7;
                 text-align: center;
             }
             .abs_tx {
                 border: 2px solid #00334d;
                 vertical-align: center;
                 display: inline-block;
                 vertical-align: top;
                 margin-left: 10px;
                 vertical-align: top;
                 font-size: 95%;
             }
             .abs_tx tr:nth-last-child(1){
                 border-bottom: 0px solid #00334d;
             }
             .abs_tx td {
                 text-align: center;
                 padding: 4px;
             }
             .abs_rx {
                 border: 2px solid #00334d;
                 vertical-align: center;
                 display: inline-block;
                 margin-left: 10px;
                 text-align: center;
                 font-size: 95%;
             }
+            .image_plot {
+              border: 2px solid #00334d;
+              display: inline-block;
+              vertical-align: 50px;
+              margin-left: 10px;
+            }
             .abs_rx tr:nth-last-child(1){
                 border-bottom: 0px solid #00334d;
             }
             .abs_rx td {
                 text-align: center;
                 padding: 4px;
             }
             .tx {
               border: 2px solid #00334d;
               vertical-align: center;
               font-weight: bold;
             }
             .tx tr{
                 border: 1px solid #ffffff;
                 background-color: #ecf0f1;
             }
             .tx td{
                 border: 2px solid #e2e2e7;
                 text-align: center;
             }
             .rx {
               border: 2px solid #00334d;
               vertical-align: center;
               font-weight: bold;
             }
             .rx tr{
                 border: 1px solid #ffffff;
                 background-color: #ecf0f1;
             }
             .rx td{
                 border: 2px solid #e2e2e7;
                 text-align: center;
             }
             .pattern_image{
               display: inline-block;
             }
             </style>
             <div id="UP" class="panel-group">
               <div class="panel panel-default">
                 <div class="panel-heading">UP</div>
                 <div class="panel-body">
                   <table class="abs">
                     <tr>
                       <td> <b>Antenna</b>
                         <table id="antenna_up" class="pattern">
                           <tr>
                             <td> North Quarter
                               <table class="north_quarter">
                                 <tr> <td>{{beam.get_upvalues.0}}</td> <td>{{beam.get_upvalues.1}}</td> <td>{{beam.get_upvalues.2}}</td> <td>{{beam.get_upvalues.3}}</td> </tr>
                                 <tr> <td>{{beam.get_upvalues.8}}</td> <td>{{beam.get_upvalues.9}}</td> <td>{{beam.get_upvalues.10}}</td> <td>{{beam.get_upvalues.11}}</td> </tr>
                                 <tr> <td>{{beam.get_upvalues.16}}</td> <td>{{beam.get_upvalues.17}}</td> <td>{{beam.get_upvalues.18}}</td> <td>{{beam.get_upvalues.19}}</td> </tr>
                                 <tr> <td>{{beam.get_upvalues.24}}</td> <td>{{beam.get_upvalues.25}}</td> <td>{{beam.get_upvalues.26}}</td> <td>{{beam.get_upvalues.27}}</td> </tr>
                               </table>
                             </td>
                             <td> East Quarter
                               <table class="east_quarter">
                                 <tr> <td>{{beam.get_upvalues.4}}</td> <td>{{beam.get_upvalues.5}}</td> <td>{{beam.get_upvalues.6}}</td> <td>{{beam.get_upvalues.7}}</td> </tr>
                                 <tr> <td>{{beam.get_upvalues.12}}</td> <td>{{beam.get_upvalues.13}}</td> <td>{{beam.get_upvalues.14}}</td> <td>{{beam.get_upvalues.15}}</td> </tr>
                                 <tr> <td>{{beam.get_upvalues.20}}</td> <td>{{beam.get_upvalues.21}}</td> <td>{{beam.get_upvalues.22}}</td> <td>{{beam.get_upvalues.23}}</td> </tr>
                                 <tr> <td>{{beam.get_upvalues.28}}</td> <td>{{beam.get_upvalues.29}}</td> <td>{{beam.get_upvalues.30}}</td> <td>{{beam.get_upvalues.31}}</td> </tr>
                               </table>
                             </td>
                           </tr>
                           <tr>
                             <td> West Quarter
                               <table class="west_quarter">
                                 <tr> <td>{{beam.get_upvalues.32}}</td> <td>{{beam.get_upvalues.33}}</td> <td>{{beam.get_upvalues.34}}</td> <td>{{beam.get_upvalues.35}}</td> </tr>
                                 <tr> <td>{{beam.get_upvalues.40}}</td> <td>{{beam.get_upvalues.41}}</td> <td>{{beam.get_upvalues.42}}</td> <td>{{beam.get_upvalues.43}}</td> </tr>
                                 <tr> <td>{{beam.get_upvalues.48}}</td> <td>{{beam.get_upvalues.49}}</td> <td>{{beam.get_upvalues.50}}</td> <td>{{beam.get_upvalues.51}}</td> </tr>
                                 <tr> <td>{{beam.get_upvalues.56}}</td> <td>{{beam.get_upvalues.57}}</td> <td>{{beam.get_upvalues.58}}</td> <td>{{beam.get_upvalues.59}}</td> </tr>
                               </table>
                             </td>
                             <td> South Quarter
                               <table class="south_quarter">
                                 <tr> <td>{{beam.get_upvalues.36}}</td> <td>{{beam.get_upvalues.37}}</td> <td>{{beam.get_upvalues.38}}</td> <td>{{beam.get_upvalues.39}}</td> </tr>
                                 <tr> <td>{{beam.get_upvalues.44}}</td> <td>{{beam.get_upvalues.45}}</td> <td>{{beam.get_upvalues.46}}</td> <td>{{beam.get_upvalues.47}}</td> </tr>
                                 <tr> <td>{{beam.get_upvalues.52}}</td> <td>{{beam.get_upvalues.53}}</td> <td>{{beam.get_upvalues.54}}</td> <td>{{beam.get_upvalues.55}}</td> </tr>
                                 <tr> <td>{{beam.get_upvalues.60}}</td> <td>{{beam.get_upvalues.61}}</td> <td>{{beam.get_upvalues.62}}</td> <td>{{beam.get_upvalues.63}}</td> </tr>
                               </table>
                             </td>
                           </tr>
                         </table>
                       </td>
                     </tr>
                   </table>
                   <table class="abs_tx">
                     <tr>
                       <td> <b>TX</b>
                         <table id="tx_up" class="tx">
                           <tr>
                             <td> North Quarter
                               <table align="center" class="north_quarter">
                                 <tr>
                                   <td>{{beam.get_tx.up.0.0}}</td> <td>{{beam.get_tx.up.0.1}}</td> <td>{{beam.get_tx.up.0.2}}</td> <td>{{beam.get_tx.up.0.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.up.1.0}}</td> <td>{{beam.get_tx.up.1.1}}</td> <td>{{beam.get_tx.up.1.2}}</td> <td>{{beam.get_tx.up.1.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.up.2.0}}</td> <td>{{beam.get_tx.up.2.1}}</td> <td>{{beam.get_tx.up.2.2}}</td> <td>{{beam.get_tx.up.2.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.up.3.0}}</td> <td>{{beam.get_tx.up.3.1}}</td> <td>{{beam.get_tx.up.3.2}}</td> <td>{{beam.get_tx.up.3.3}}</td>
                                 </tr>
                               </table>
                             </td>
                             <td> East Quarter
                               <table align="center" class="east_quarter">
                                 <tr>
                                   <td>{{beam.get_tx.up.0.4}}</td> <td>{{beam.get_tx.up.0.5}}</td> <td>{{beam.get_tx.up.0.6}}</td> <td>{{beam.get_tx.up.0.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.up.1.4}}</td> <td>{{beam.get_tx.up.1.5}}</td> <td>{{beam.get_tx.up.1.6}}</td> <td>{{beam.get_tx.up.1.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.up.2.4}}</td> <td>{{beam.get_tx.up.2.5}}</td> <td>{{beam.get_tx.up.2.6}}</td> <td>{{beam.get_tx.up.2.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.up.3.4}}</td> <td>{{beam.get_tx.up.3.5}}</td> <td>{{beam.get_tx.up.3.6}}</td> <td>{{beam.get_tx.up.3.7}}</td>
                                 </tr>
                               </table>
                             </td>
                           </tr>
                           <tr>
                             <td> West Quarter
                               <table align="center" class="west_quarter">
                                 <tr>
                                   <td>{{beam.get_tx.up.4.0}}</td> <td>{{beam.get_tx.up.4.1}}</td> <td>{{beam.get_tx.up.4.2}}</td> <td>{{beam.get_tx.up.4.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.up.5.0}}</td> <td>{{beam.get_tx.up.5.1}}</td> <td>{{beam.get_tx.up.5.2}}</td> <td>{{beam.get_tx.up.5.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.up.6.0}}</td> <td>{{beam.get_tx.up.6.1}}</td> <td>{{beam.get_tx.up.6.2}}</td> <td>{{beam.get_tx.up.6.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.up.7.0}}</td> <td>{{beam.get_tx.up.7.1}}</td> <td>{{beam.get_tx.up.7.2}}</td> <td>{{beam.get_tx.up.7.3}}</td>
                                 </tr>
                               </table>
                             </td>
                             <td> South Quarter
                               <table align="center" class="south_quarter">
                                 <tr>
                                   <td>{{beam.get_tx.up.4.4}}</td> <td>{{beam.get_tx.up.4.5}}</td> <td>{{beam.get_tx.up.4.6}}</td> <td>{{beam.get_tx.up.4.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.up.5.4}}</td> <td>{{beam.get_tx.up.5.5}}</td> <td>{{beam.get_tx.up.5.6}}</td> <td>{{beam.get_tx.up.5.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.up.6.4}}</td> <td>{{beam.get_tx.up.6.5}}</td> <td>{{beam.get_tx.up.6.6}}</td> <td>{{beam.get_tx.up.6.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.up.7.4}}</td> <td>{{beam.get_tx.up.7.5}}</td> <td>{{beam.get_tx.up.7.6}}</td> <td>{{beam.get_tx.up.7.7}}</td>
                                 </tr>
                               </table>
                             </td>
                           </tr>
                         </table>
                       </td>
                     </tr>
                   </table>
                   <table class="abs_rx">
                     <tr>
                       <td> <b>RX</b>
                         <table id="rx_up" class="rx">
                           <tr>
                             <td> North Quarter
                               <table align="center" class="north_quarter">
                                 <tr>
                                   <td>{{beam.get_rx.up.0.0}}</td> <td>{{beam.get_rx.up.0.1}}</td> <td>{{beam.get_rx.up.0.2}}</td> <td>{{beam.get_rx.up.0.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.up.1.0}}</td> <td>{{beam.get_rx.up.1.1}}</td> <td>{{beam.get_rx.up.1.2}}</td> <td>{{beam.get_rx.up.1.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.up.2.0}}</td> <td>{{beam.get_rx.up.2.1}}</td> <td>{{beam.get_rx.up.2.2}}</td> <td>{{beam.get_rx.up.2.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.up.3.0}}</td> <td>{{beam.get_rx.up.3.1}}</td> <td>{{beam.get_rx.up.3.2}}</td> <td>{{beam.get_rx.up.3.3}}</td>
                                 </tr>
                               </table>
                             </td>
                             <td> East Quarter
                               <table align="center" class="east_quarter">
                                 <tr>
                                   <td>{{beam.get_rx.up.0.4}}</td> <td>{{beam.get_rx.up.0.5}}</td> <td>{{beam.get_rx.up.0.6}}</td> <td>{{beam.get_rx.up.0.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.up.1.4}}</td> <td>{{beam.get_rx.up.1.5}}</td> <td>{{beam.get_rx.up.1.6}}</td> <td>{{beam.get_rx.up.1.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.up.2.4}}</td> <td>{{beam.get_rx.up.2.5}}</td> <td>{{beam.get_rx.up.2.6}}</td> <td>{{beam.get_rx.up.2.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.up.3.4}}</td> <td>{{beam.get_rx.up.3.5}}</td> <td>{{beam.get_rx.up.3.6}}</td> <td>{{beam.get_rx.up.3.7}}</td>
                                 </tr>
                               </table>
                             </td>
                           </tr>
                           <tr>
                             <td> West Quarter
                               <table align="center" class="west_quarter">
                                 <tr>
                                   <td>{{beam.get_rx.up.4.0}}</td> <td>{{beam.get_rx.up.4.1}}</td> <td>{{beam.get_rx.up.4.2}}</td> <td>{{beam.get_rx.up.4.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.up.5.0}}</td> <td>{{beam.get_rx.up.5.1}}</td> <td>{{beam.get_rx.up.5.2}}</td> <td>{{beam.get_rx.up.5.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.up.6.0}}</td> <td>{{beam.get_rx.up.6.1}}</td> <td>{{beam.get_rx.up.6.2}}</td> <td>{{beam.get_rx.up.6.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.up.7.0}}</td> <td>{{beam.get_rx.up.7.1}}</td> <td>{{beam.get_rx.up.7.2}}</td> <td>{{beam.get_rx.up.7.3}}</td>
                                 </tr>
                               </table>
                             </td>
                             <td> South Quarter
                               <table align="center" class="south_quarter">
                                 <tr>
                                   <td>{{beam.get_rx.up.4.4}}</td> <td>{{beam.get_rx.up.4.5}}</td> <td>{{beam.get_rx.up.4.6}}</td> <td>{{beam.get_rx.up.4.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.up.5.4}}</td> <td>{{beam.get_rx.up.5.5}}</td> <td>{{beam.get_rx.up.5.6}}</td> <td>{{beam.get_rx.up.5.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.up.6.4}}</td> <td>{{beam.get_rx.up.6.5}}</td> <td>{{beam.get_rx.up.6.6}}</td> <td>{{beam.get_rx.up.6.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.up.7.4}}</td> <td>{{beam.get_rx.up.7.5}}</td> <td>{{beam.get_rx.up.7.6}}</td> <td>{{beam.get_rx.up.7.7}}</td>
                                 </tr>
                               </table>
                             </td>
                           </tr>
                         </table>
                       </td>
                     </tr>
                   </table>
-                  {% if not edit %}
+                  <table class="image_plot">
-                    {% include "abs_uppattern_img.html" %}
+                    <tr>
-                  {% endif %}
+                      <td>
+                        {% if not edit %}
+                          {% include "abs_uppattern_img.html" %}
+                        {% endif %}
+                        <!-- <img id="imgMains" src="{% url 'url_plot_beam' beam.abs_conf.id beam.id 'up' %}" alt="Error ploting..." style="width: 250px;height:250px; background-image: url({% static 'images/loader.gif' %});background-repeat: no-repeat;background-position: 50% 50%;">   -->
+                      </td>
+                    </tr>
+                  </table>
                   <br>
                   <div id="up_ues" style="display: inline-block">
                     Ues: {{beam.get_up_ues}}
                   </div>
                   <div style="margin-left: 70px; display: inline-block">
                     <input type="checkbox" id="up_onlyrx" {% if beam.get_up_onlyrx == True %} checked="True" {% endif %} disabled>
                       Only RX
                     </input>
                   </div>
                 </div>
               </div>
             </div>
             <div id="DOWN" class="panel-group">
               <div class="panel panel-default">
                 <div class="panel-heading">DOWN</div>
                 <div class="panel-body">
                   <table class="abs">
                     <tr>
                       <td> <b>Antenna</b>
                         <table id="antenna_down" class="pattern">
                           <tr>
                             <td> <b>North Quarter</b>
                               <table class="north_quarter">
                                 <tr> <td>{{beam.get_downvalues.0}}</td> <td>{{beam.get_downvalues.1}}</td> <td>{{beam.get_downvalues.2}}</td> <td>{{beam.get_downvalues.3}}</td> </tr>
                                 <tr> <td>{{beam.get_downvalues.8}}</td> <td>{{beam.get_downvalues.9}}</td> <td>{{beam.get_downvalues.10}}</td> <td>{{beam.get_downvalues.11}}</td> </tr>
                                 <tr> <td>{{beam.get_downvalues.16}}</td> <td>{{beam.get_downvalues.17}}</td> <td>{{beam.get_downvalues.18}}</td> <td>{{beam.get_downvalues.19}}</td> </tr>
                                 <tr> <td>{{beam.get_downvalues.24}}</td> <td>{{beam.get_downvalues.25}}</td> <td>{{beam.get_downvalues.26}}</td> <td>{{beam.get_downvalues.27}}</td> </tr>
                               </table>
                             </td>
                             <td> <b>East Quarter</b>
                               <table class="east_quarter">
                                 <tr> <td>{{beam.get_downvalues.4}}</td> <td>{{beam.get_downvalues.5}}</td> <td>{{beam.get_downvalues.6}}</td> <td>{{beam.get_downvalues.7}}</td> </tr>
                                 <tr> <td>{{beam.get_downvalues.12}}</td> <td>{{beam.get_downvalues.13}}</td> <td>{{beam.get_downvalues.14}}</td> <td>{{beam.get_downvalues.15}}</td> </tr>
                                 <tr> <td>{{beam.get_downvalues.20}}</td> <td>{{beam.get_downvalues.21}}</td> <td>{{beam.get_downvalues.22}}</td> <td>{{beam.get_downvalues.23}}</td> </tr>
                                 <tr> <td>{{beam.get_downvalues.28}}</td> <td>{{beam.get_downvalues.29}}</td> <td>{{beam.get_downvalues.30}}</td> <td>{{beam.get_downvalues.31}}</td> </tr>
                               </table>
                             </td>
                           </tr>
                           <tr>
                             <td> <b>West Quarter</b>
                               <table class="west_quarter">
                                 <tr> <td>{{beam.get_downvalues.32}}</td> <td>{{beam.get_downvalues.33}}</td> <td>{{beam.get_downvalues.34}}</td> <td>{{beam.get_downvalues.35}}</td> </tr>
                                 <tr> <td>{{beam.get_downvalues.40}}</td> <td>{{beam.get_downvalues.41}}</td> <td>{{beam.get_downvalues.42}}</td> <td>{{beam.get_downvalues.43}}</td> </tr>
                                 <tr> <td>{{beam.get_downvalues.48}}</td> <td>{{beam.get_downvalues.49}}</td> <td>{{beam.get_downvalues.50}}</td> <td>{{beam.get_downvalues.51}}</td> </tr>
                                 <tr> <td>{{beam.get_downvalues.56}}</td> <td>{{beam.get_downvalues.57}}</td> <td>{{beam.get_downvalues.58}}</td> <td>{{beam.get_downvalues.59}}</td> </tr>
                               </table>
                             </td>
                             <td> <b>South Quarter</b>
                               <table class="south_quarter">
                                 <tr> <td>{{beam.get_downvalues.36}}</td> <td>{{beam.get_downvalues.37}}</td> <td>{{beam.get_downvalues.38}}</td> <td>{{beam.get_downvalues.39}}</td> </tr>
                                 <tr> <td>{{beam.get_downvalues.44}}</td> <td>{{beam.get_downvalues.45}}</td> <td>{{beam.get_downvalues.46}}</td> <td>{{beam.get_downvalues.47}}</td> </tr>
                                 <tr> <td>{{beam.get_downvalues.52}}</td> <td>{{beam.get_downvalues.53}}</td> <td>{{beam.get_downvalues.54}}</td> <td>{{beam.get_downvalues.55}}</td> </tr>
                                 <tr> <td>{{beam.get_downvalues.60}}</td> <td>{{beam.get_downvalues.61}}</td> <td>{{beam.get_downvalues.62}}</td> <td>{{beam.get_downvalues.63}}</td> </tr>
                               </table>
                             </td>
                           </tr>
                         </table>
                       </td>
                     </tr>
                   </table>
                   <table class="abs_tx">
                     <tr>
                       <td> <b>TX</b>
                         <table id="tx_down" class="tx">
                           <tr>
                             <td> <b>North Quarter</b>
                               <table align="center" class="north_quarter">
                                 <tr>
                                   <td>{{beam.get_tx.down.0.0}}</td> <td>{{beam.get_tx.down.0.1}}</td> <td>{{beam.get_tx.down.0.2}}</td> <td>{{beam.get_tx.down.0.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.down.1.0}}</td> <td>{{beam.get_tx.down.1.1}}</td> <td>{{beam.get_tx.down.1.2}}</td> <td>{{beam.get_tx.down.1.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.down.2.0}}</td> <td>{{beam.get_tx.down.2.1}}</td> <td>{{beam.get_tx.down.2.2}}</td> <td>{{beam.get_tx.down.2.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.down.3.0}}</td> <td>{{beam.get_tx.down.3.1}}</td> <td>{{beam.get_tx.down.3.2}}</td> <td>{{beam.get_tx.down.3.3}}</td>
                                 </tr>
                               </table>
                             </td>
                             <td> <b>East Quarter</b>
                               <table align="center" class="east_quarter">
                                 <tr>
                                   <td>{{beam.get_tx.down.0.4}}</td> <td>{{beam.get_tx.down.0.5}}</td> <td>{{beam.get_tx.down.0.6}}</td> <td>{{beam.get_tx.down.0.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.down.1.4}}</td> <td>{{beam.get_tx.down.1.5}}</td> <td>{{beam.get_tx.down.1.6}}</td> <td>{{beam.get_tx.down.1.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.down.2.4}}</td> <td>{{beam.get_tx.down.2.5}}</td> <td>{{beam.get_tx.down.2.6}}</td> <td>{{beam.get_tx.down.2.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.down.3.4}}</td> <td>{{beam.get_tx.down.3.5}}</td> <td>{{beam.get_tx.down.3.6}}</td> <td>{{beam.get_tx.down.3.7}}</td>
                                 </tr>
                               </table>
                             </td>
                           </tr>
                           <tr>
                             <td> <b>West Quarter</b>
                               <table align="center" class="west_quarter">
                                 <tr>
                                   <td>{{beam.get_tx.down.4.0}}</td> <td>{{beam.get_tx.down.4.1}}</td> <td>{{beam.get_tx.down.4.2}}</td> <td>{{beam.get_tx.down.4.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.down.5.0}}</td> <td>{{beam.get_tx.down.5.1}}</td> <td>{{beam.get_tx.down.5.2}}</td> <td>{{beam.get_tx.down.5.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.down.6.0}}</td> <td>{{beam.get_tx.down.6.1}}</td> <td>{{beam.get_tx.down.6.2}}</td> <td>{{beam.get_tx.down.6.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.down.7.0}}</td> <td>{{beam.get_tx.down.7.1}}</td> <td>{{beam.get_tx.down.7.2}}</td> <td>{{beam.get_tx.down.7.3}}</td>
                                 </tr>
                               </table>
                             </td>
                             <td> <b>South Quarter</b>
                               <table align="center" class="south_quarter">
                                 <tr>
                                   <td>{{beam.get_tx.down.4.4}}</td> <td>{{beam.get_tx.down.4.5}}</td> <td>{{beam.get_tx.down.4.6}}</td> <td>{{beam.get_tx.down.4.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.down.5.4}}</td> <td>{{beam.get_tx.down.5.5}}</td> <td>{{beam.get_tx.down.5.6}}</td> <td>{{beam.get_tx.down.5.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.down.6.4}}</td> <td>{{beam.get_tx.down.6.5}}</td> <td>{{beam.get_tx.down.6.6}}</td> <td>{{beam.get_tx.down.6.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_tx.down.7.4}}</td> <td>{{beam.get_tx.down.7.5}}</td> <td>{{beam.get_tx.down.7.6}}</td> <td>{{beam.get_tx.down.7.7}}</td>
                                 </tr>
                               </table>
                             </td>
                           </tr>
                         </table>
                       </td>
                     </tr>
                   </table>
                   <table class="abs_rx">
                     <tr>
                       <td> <b>RX</b>
                         <table id="rx_down" class="rx">
                           <tr>
                             <td> <b>North Quarter</b>
                               <table align="center" class="north_quarter">
                                 <tr>
                                   <td>{{beam.get_rx.down.0.0}}</td> <td>{{beam.get_rx.down.0.1}}</td> <td>{{beam.get_rx.down.0.2}}</td> <td>{{beam.get_rx.down.0.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.down.1.0}}</td> <td>{{beam.get_rx.down.1.1}}</td> <td>{{beam.get_rx.down.1.2}}</td> <td>{{beam.get_rx.down.1.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.down.2.0}}</td> <td>{{beam.get_rx.down.2.1}}</td> <td>{{beam.get_rx.down.2.2}}</td> <td>{{beam.get_rx.down.2.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.down.3.0}}</td> <td>{{beam.get_rx.down.3.1}}</td> <td>{{beam.get_rx.down.3.2}}</td> <td>{{beam.get_rx.down.3.3}}</td>
                                 </tr>
                               </table>
                             </td>
                             <td> <b>East Quarter</b>
                               <table align="center" class="east_quarter">
                                 <tr>
                                   <td>{{beam.get_rx.down.0.4}}</td> <td>{{beam.get_rx.down.0.5}}</td> <td>{{beam.get_rx.down.0.6}}</td> <td>{{beam.get_rx.down.0.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.down.1.4}}</td> <td>{{beam.get_rx.down.1.5}}</td> <td>{{beam.get_rx.down.1.6}}</td> <td>{{beam.get_rx.down.1.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.down.2.4}}</td> <td>{{beam.get_rx.down.2.5}}</td> <td>{{beam.get_rx.down.2.6}}</td> <td>{{beam.get_rx.down.2.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.down.3.4}}</td> <td>{{beam.get_rx.down.3.5}}</td> <td>{{beam.get_rx.down.3.6}}</td> <td>{{beam.get_rx.down.3.7}}</td>
                                 </tr>
                               </table>
                             </td>
                           </tr>
                           <tr>
                             <td> <b>West Quarter</b>
                               <table align="center" class="west_quarter">
                                 <tr>
                                   <td>{{beam.get_rx.down.4.0}}</td> <td>{{beam.get_rx.down.4.1}}</td> <td>{{beam.get_rx.down.4.2}}</td> <td>{{beam.get_rx.down.4.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.down.5.0}}</td> <td>{{beam.get_rx.down.5.1}}</td> <td>{{beam.get_rx.down.5.2}}</td> <td>{{beam.get_rx.down.5.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.down.6.0}}</td> <td>{{beam.get_rx.down.6.1}}</td> <td>{{beam.get_rx.down.6.2}}</td> <td>{{beam.get_rx.down.6.3}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.down.7.0}}</td> <td>{{beam.get_rx.down.7.1}}</td> <td>{{beam.get_rx.down.7.2}}</td> <td>{{beam.get_rx.down.7.3}}</td>
                                 </tr>
                               </table>
                             </td>
                             <td> <b>South Quarter</b>
                               <table class="south_quarter">
                                 <tr>
                                   <td>{{beam.get_rx.down.4.4}}</td> <td>{{beam.get_rx.down.4.5}}</td> <td>{{beam.get_rx.down.4.6}}</td> <td>{{beam.get_rx.down.4.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.down.5.4}}</td> <td>{{beam.get_rx.down.5.5}}</td> <td>{{beam.get_rx.down.5.6}}</td> <td>{{beam.get_rx.down.5.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.down.6.4}}</td> <td>{{beam.get_rx.down.6.5}}</td> <td>{{beam.get_rx.down.6.6}}</td> <td>{{beam.get_rx.down.6.7}}</td>
                                 </tr>
                                 <tr>
                                   <td>{{beam.get_rx.down.7.4}}</td> <td>{{beam.get_rx.down.7.5}}</td> <td>{{beam.get_rx.down.7.6}}</td> <td>{{beam.get_rx.down.7.7}}</td>
                                 </tr>
                               </table>
                             </td>
                           </tr>
                         </table>
                       </td>
                     </tr>
                   </table>
-                  {% if not edit %}
+                  <table class="image_plot">
-                    {% include "abs_downpattern_img.html" %}
+                    <tr>
-                  {% endif %}
+                      <td>
+                        {% if not edit %}
+                          {% include "abs_downpattern_img.html" %}
+                        {% endif %}
+                      </td>
+                    </tr>
+                  </table>
                   <br>
                   <div id="down_ues" style="display: inline-block">
                     Ues: {{beam.get_down_ues}}
                   </div>
                   <div style="margin-left: 70px; display: inline-block">
                     <input type="checkbox" id="up_onlyrx" {% if beam.get_down_onlyrx == True %} checked="True" {% endif %} disabled>
                       Only RX
                     </input>
                   </div>
                 </div>
               </div>
             </div>
             {% endblock %}
             <script>
             </script>

apps/abs/templates/abs_patterns.html +1 -1

             {% extends "dev_conf.html" %}
             {% load static %}
             {% load django_bootstrap5 %}
             {% load main_tags %}
             {% block content %}
             {% if abs_beams %}
               <div>
                 <h4>Beams:</h4>
                 <div class="container">
                   <div class="btn-group">
                     {% for abs_beam in abs_beams %}
                       <button id="bt_beam{{ forloop.counter }}" type="button"  class="btn btn-default">{{ forloop.counter }}</button>
                     {% endfor %}
                   </div>
                 </div>
               </div>
               <br>
               {% if beam %}
                 {% include "abs_pattern.html" %}
               {% endif %}
             {% else %}
               <div>
                 <h4>Beams:</h4>
                 <p style="color:#b4bcc2; margin-left: 5%;"><i>No Beams...</i></p>
               </div>
             {% endif %}
             <script>
             {% for abs_beam in abs_beams %}
               $("#bt_beam{{ forloop.counter }}").click(function() {
-                document.location = "{% url 'url_plot_abs_patterns' abs_beam.abs_conf.id abs_beam.id %}";
+                document.location = "{% url 'url_plot_abs_patterns2' abs_beam.abs_conf.id abs_beam.id %}";
               });
             {% endfor %}
             //function ChangeColor() {
             //    document.getElementById("button_1").style.backgroundColor = "#2c3e50";
             //    document.getElementById("button_1").style.color = "#ecf0f1";
             //}
             </script>
             {% endblock %}

apps/abs/templates/abs_uppattern_img.html +1 -7

@@ -1,11 +1,5
1	{% load static %}	1	{% load static %}
2	{% load django_bootstrap5 %}	2	{% load django_bootstrap5 %}
3	{% load main_tags %}	3	{% load main_tags %}
4		4
5	{% block content %}	5	<img id="imgMainUp" src="{% url 'url_plot_beam' beam.abs_conf.id beam.id 'up' %}" alt="Error ploting..." style="width: 250px;height:250px; background-image: url({% static 'images/loader.gif' %});background-repeat: no-repeat;background-position: 50% 50%;">
6
7	<div id="PictureOverJROUp" style="float: right">
8	<img id="imgMain" src="{% url 'url_plot_beam' beam.abs_conf.id beam.id 'up' %}" alt="Error ploting..." style="width:340px;height:340px; background-image: url({% static 'images/loader.gif' %});background-repeat: no-repeat;background-position: 50% 50%;">
9	</div>
10
11	{% endblock %}

apps/abs/urls.py +2 -2

             from django.urls import path
             from apps.abs import views
             urlpatterns = (
                 path('<int:id_conf>/', views.abs_conf, name='url_abs_conf'),
                 path('<int:id_conf>/edit/', views.abs_conf_edit, name='url_edit_abs_conf'),
                 path('alert/', views.abs_conf_alert, name='url_alert_abs_conf'),
                 path('<int:id_conf>/import/', views.import_file, name='url_import_abs_conf'),
                 #url(r'^(?P<id_conf>-?\d+)/status/', views.abs_conf, {'status_request':True},name='url_status_abs_conf'),
                 path('<int:id_conf>/change_beam/<int:id_beam>/', views.send_beam, name='url_send_beam'),
                 path('<int:id_conf>/plot/', views.plot_patterns, name='url_plot_abs_patterns'),
-                path('<int:id_conf>/plot/<int:id_beam>/', views.plot_patterns, name='url_plot_abs_patterns'),
+                path('<int:id_conf>/plot/<int:id_beam>/', views.plot_patterns, name='url_plot_abs_patterns2'),
-                path('<int:id_conf>/plot/<int:id_beam>/<int:antenna>/pattern.png', views.plot_pattern, name='url_plot_beam'),
+                path('<int:id_conf>/plot/<int:id_beam>/<slug:antenna>/pattern.png/', views.plot_pattern, name='url_plot_beam'),
                 path('<int:id_conf>/add_beam/', views.add_beam, name='url_add_abs_beam'),
                 path('<int:id_conf>/beam/<int:id_beam>/delete/', views.remove_beam, name='url_remove_abs_beam'),
                 path('<int:id_conf>/beam/<int:id_beam>/edit/', views.edit_beam, name='url_edit_abs_beam'),
             )

apps/abs/utils/Graphics_OverJro.py +12 -12

             """
             The module GRAPHICS_OVERJRO.py gathers classes or/and functions to create graphics from OVER-JRO
             project (e.g. antenna patterns, skynoise, ...).
             MODULES CALLED:
             TIME, NUMPY, MATPLOTLIB, TIMETOOLS
             MODIFICATION HISTORY:
             Created by Ing. Freddy Galindo (frederickgalindo@gmail.com). ROJ Oct 18, 2009.
             """
             import time
             import numpy
             import sys
             import os
             # set HOME environment variable to a directory the httpd server can write to
             #os.environ[ 'HOME' ] = '/usr/local/www/htdocs/overJro/tempReports'
             #os.environ[ 'HOME' ] = '/home/dsuarez/Pictures'
             #os.environ[ 'HOME' ] = '/tmp/'
             import matplotlib
             #if ide==1:
             #    matplotlib.use('Qt4Agg')
             #elif ide==2:
             #    matplotlib.use("Agg")
             #else:
             #    matplotlib.use('TKAgg')
             matplotlib.use("Agg")
             #matplotlib.interactive(1)
             import matplotlib.pyplot
             #import Numeric
             #import scipy
             import scipy.interpolate
             from .Astro_Coords import Equatorial , CelestialBodies
             from .TimeTools import Time , Julian
             from .Graphics_Miscens import ColorTable
             from .Misc_Routines import CoFactors,Vector
             class AntPatternPlot:
                 def __init__(self):
                     """
                     AntPatternPlot creates an object to call methods to plot the antenna pattern.
                     Modification History
                     --------------------
                     Created by Freddy Galindo, ROJ, 06 October 2009.
                     """
                     self.fig = matplotlib.pyplot.figure(figsize=(8,8), facecolor='white')
                     self.ax = self.fig.add_subplot(111)
                 def contPattern(self,iplot=0,gpath='',filename='',mesg='',amp=None ,x=None ,y=None ,getCut=None,title='', save=True):
                     """
                     contPattern plots a contour map of the antenna pattern.
                     Parameters
                     ----------
                     iplot = A integer to specify if the plot is  the first, second, ...  The default va-
                       lue is 0.
                     Examples
                     --------
                     >> Over_Jro.JroPattern(pattern=2).contPattern()
                     Modification history
                     --------------------
                     Converted to Python by Freddy R. Galindo, ROJ, 06 October 2009.
                     """
                     if getCut == 1:
                         return
                     xmax = numpy.max(x)
                     xmin = numpy.min(x)
                     ymax = numpy.max(y)
                     ymin = numpy.min(y)
                     levels = numpy.array([1e-3,1e-2,1e-1,0.5,1.0])
                     tmp = numpy.round(10*numpy.log10(levels),decimals=1)
                     labels = []
-                    for i in numpy.arange(5):labels.append(str(numpy.int(tmp[i])))
+                    for i in numpy.arange(5):labels.append(str(int(tmp[i])))
                     colors = ((0,0,1.),(0,170/255.,0),(127/255.,1.,0),(1.,109/255.,0),(128/255.,0,0))
                     CS = self.ax.contour(x,y,amp.transpose(),levels,colors=colors)
                     fmt = {}
                     for l,s in zip(CS.levels,labels):
                         fmt[l] = s
                     self.ax.annotate('Ng',xy=(-0.05,1.04),xytext=(0.01,0.962),xycoords='axes fraction',arrowprops=dict(facecolor='black', width=1.,shrink=0.2),fontsize=15.)
                     self.ax.annotate(mesg,xy=(0,0),xytext=(0.01,0.01),xycoords='figure fraction')
                     self.ax.clabel(CS,CS.levels,inline=True,fmt=fmt,fontsize=10)
                     self.ax.set_xlim(xmin,xmax)
                     self.ax.set_ylim(ymin,ymax)
                     self.ax.set_title("Total Pattern: " + title)
                     self.ax.set_xlabel("West to South")
                     self.ax.set_ylabel("West to North")
                     self.ax.grid(True)
                     if save:
                         save_fig = os.path.join(gpath,filename)
                         self.fig.savefig(save_fig,format='png')
                 def close(self):
                     matplotlib.pyplot.close(self.fig)
                 def plotRaDec(self,gpath=None,filename=None,jd=2452640.5,ra_obs=None,xg=None,yg=None,x=None,y=None, save=True):
                     """
                     plotRaDec draws right ascension and declination lines on a JRO plane. This function
                     must call after conPattern.
                     Parameters
                     ----------
                     jd = A scalar giving the Julian date.
                     ra_obs = Scalar giving the right ascension of the observatory.
                     xg = A 3-element array to specify ..
                     yg = A 3-element array to specify ..
                     Examples
                     --------
                     >> Over_Jro.JroPattern(pattern=2).contPattern()
                     >> Over_Jro.JroPattern(pattern=2).plotRaDec(jd=jd,ra_obs=ra_obs,xg=xg,yg=yg)
                     Modification history
                     --------------------
                     Converted to Python by Freddy R. Galindo, ROJ, 06 October 2009.
                     """
                     # Finding RA of observatory for a specific date
                     if ra_obs is None:ra_obs = numpy.array([23.37060849])
                     if xg is None:xg = numpy.array([0.62918474,-0.77725579,0.])
                     if yg is None:yg = numpy.array([0.77700346,0.62898048,0.02547905])
                     # Getting HA and DEC axes
                     mindec = -28; maxdec = 4; incdec = 2.
-                    ndec = numpy.int((maxdec - mindec)/incdec) + 1
+                    ndec = int((maxdec - mindec)/incdec) + 1
                     minha = -20; maxha = 20; incha = 2.
-                    nha = numpy.int((maxha - minha)/incha) + 1
+                    nha = int((maxha - minha)/incha) + 1
                     #mcosx = numpy.zeros((nha,ndec))
                     #mcosy = numpy.zeros((nha,ndec))
                     ha_axes = numpy.reshape(numpy.arange(nha)*incha + minha,(nha,1))
                     ones_dec = numpy.reshape(numpy.zeros(ndec) + 1,(ndec,1))
                     ha_axes = numpy.dot(ha_axes,ones_dec.transpose())
                     ha_axes2 = numpy.array(ra_obs - ha_axes)
                     dec_axes = numpy.reshape(numpy.arange(ndec)*incdec + mindec,(ndec,1))
                     ones_ra = numpy.reshape(numpy.zeros(nha) + 1,(nha,1))
                     dec_axes = numpy.dot(ones_ra,dec_axes.transpose())
                     dec_axes2 = numpy.array(dec_axes)
                     ObjHor = Equatorial(ha_axes2,dec_axes2,jd)
                     [alt,az,ha] = ObjHor.change2AltAz()
                     z = numpy.transpose(alt)*CoFactors.d2r  ; z = z.flatten()
                     az = numpy.transpose(az)*CoFactors.d2r  ; az = az.flatten()
                     vect = numpy.array([numpy.cos(z)*numpy.sin(az),numpy.cos(z)*numpy.cos(az),numpy.sin(z)])
                     xg = numpy.atleast_2d(xg)
                     dcosx = numpy.array(numpy.dot(xg,vect))
                     yg = numpy.atleast_2d(yg)
                     dcosy = numpy.array(numpy.dot(yg,vect))
                     mcosx = dcosx.reshape(ndec,nha)
                     mcosy = dcosy.reshape(ndec,nha)
                     # Defining NAN for points outof limits.
                     xmax = numpy.max(x)
                     xmin = numpy.min(x)
                     ymax = numpy.max(y)
                     ymin = numpy.min(y)
                     factor = 1.3
                     noval = numpy.where((mcosx>(xmax*factor)) | (mcosx<(xmin*factor)))
                     if noval[0].size>0:mcosx[noval] = numpy.nan
                     noval = numpy.where((mcosy>(ymax*factor)) | (mcosy<(ymin*factor)))
                     if noval[0].size>0:mcosy[noval] = numpy.nan
                     # Plotting HA and declination grid.
-                    iha0 = numpy.int((0 - minha)/incha)
+                    iha0 = int((0 - minha)/incha)
-                    idec0 = numpy.int((-14 - mindec)/incdec)
+                    idec0 = int((-14 - mindec)/incdec)
                     colorgrid = (1.,109/255.,0)
                     self.ax.plot(mcosx.transpose(),mcosy.transpose(),color=colorgrid,linestyle='--')
                     for idec in numpy.arange(ndec):
                         if idec != idec0:
                             valx = (mcosx[idec,iha0]<=xmax) & (mcosx[idec,iha0]>=xmin)
                             valy = (mcosy[idec,iha0]<=ymax) & (mcosy[idec,iha0]>=ymin)
                             if valx & valy:
-                                text = str(numpy.int(mindec + incdec*idec))+'$^o$'
+                                text = str(int(mindec + incdec*idec))+'$^o$'
                                 self.ax.text(mcosx[idec,iha0],mcosy[idec,iha0],text)
                     matplotlib.pyplot.plot(mcosx,mcosy,color=colorgrid,linestyle='--')
                     for iha in numpy.arange(nha):
                         if iha != iha0:
                             valx = (mcosx[idec0,iha]<=xmax) & (mcosx[idec0,iha]>=xmin)
                             valy = (mcosy[idec0,iha]<=ymax) & (mcosy[idec0,iha]>=ymin)
                             if valx & valy:
-                                text = str(4*numpy.int(minha + incha*iha))+"'"
+                                text = str(4*int(minha + incha*iha))+"'"
                                 self.ax.text(mcosx[idec0,iha],mcosy[idec0,iha],text)
                     if save:
                         save_fig = os.path.join(gpath,filename)
                         matplotlib.pyplot.savefig(save_fig,format='png')
                 def plotBField(self,gpath,filename,dcos,alpha, nlon, nlat, dcosxrange, dcosyrange, heights, alpha_i, save=True):
                     """
                     plotBField draws the magnetic field in a directional cosines plot.
                     Parameters
                     ----------
                     dcos = An 4-dimensional array giving the directional cosines of the magnetic field
                       over the desired place.
                     alpha = An 3-dimensional array giving the angle of the magnetic field over the desi-
                       red place.
                     nlon = An integer to specify the number of elements per longitude.
                     nlat = An integer to specify the number of elements per latitude.
                     dcosxrange = A 2-element array giving the range of the  directional cosines  in the
                       "x" axis.
                     dcosyrange = A 2-element array giving the range of the  directional cosines  in the
                       "y" axis.
                     heights = An array giving the heights (km) where the magnetic field will be modeled               By default the magnetic field will be computed at 100, 500 and 1000km.
                     alpha_i = Angle to interpolate the magnetic field.
                     Modification History
                     --------------------
                     Converted to Python by Freddy R. Galindo, ROJ, 07 October 2009.
                     """
                     handles = []
                     objects = []
                     colors = ['k','m','c','b','g','r','y']
                     marker = ['-+','-*','-D','-x','-s','->','-o','-^']
                     alpha_location = numpy.zeros((nlon,2,heights.size))
                     for ih in numpy.arange(heights.size):
                         alpha_location[:,0,ih] = dcos[:,0,ih,0]
                         for ilon in numpy.arange(nlon):
                             myx = (alpha[ilon,:,ih])[::-1]
                             myy = (dcos[ilon,:,ih,0])[::-1]
                             tck = scipy.interpolate.splrep(myx,myy,s=0)
                             mydcosx = scipy.interpolate.splev(alpha_i,tck,der=0)
                             myx = (alpha[ilon,:,ih])[::-1]
                             myy = (dcos[ilon,:,ih,1])[::-1]
                             tck = scipy.interpolate.splrep(myx,myy,s=0)
                             mydcosy = scipy.interpolate.splev(alpha_i,tck,der=0)
                             alpha_location[ilon,:,ih] = numpy.array([mydcosx, mydcosy])
                         ObjFig, = self.ax.plot(alpha_location[:,0,ih],alpha_location[:,1,ih],
-                            marker[ih % 8],color=colors[numpy.int(ih/8)],ms=4.5,lw=0.5)
+                            marker[ih % 8],color=colors[int(ih/8)],ms=4.5,lw=0.5)
                         handles.append(ObjFig)
-                        objects.append(numpy.str(heights[ih]) + ' km')
+                        objects.append(str(heights[ih]) + ' km')
                     self.ax.legend(handles,objects,loc="lower right", numpoints=1, handlelength=0.3,
                                    handletextpad=0.02, borderpad=0.3, labelspacing=0.1)
                     if save:
                         save_fig = os.path.join(gpath,filename)
                         matplotlib.pyplot.savefig(save_fig,format='png')
             class BFieldPlot:
                 def __init__(self):
                     """
                     BFieldPlot creates an object for drawing magnetic Field lines over Jicamarca.
                     Modification History
                     --------------------
                     Created by Freddy Galindo, ROJ, 07 October 2009.
                     """
                     self.alpha_location = 1
             #        pass
                 def plotBField(self,gpath,filename,dcos,alpha, nlon, nlat, dcosxrange, dcosyrange, heights, alpha_i):
                     """
                     plotBField draws the magnetic field in a directional cosines plot.
                     Parameters
                     ----------
                     dcos = An 4-dimensional array giving the directional cosines of the magnetic field
                       over the desired place.
                     alpha = An 3-dimensional array giving the angle of the magnetic field over the desi-
                       red place.
                     nlon = An integer to specify the number of elements per longitude.
                     nlat = An integer to specify the number of elements per latitude.
                     dcosxrange = A 2-element array giving the range of the  directional cosines  in the
                       "x" axis.
                     dcosyrange = A 2-element array giving the range of the  directional cosines  in the
                       "y" axis.
                     heights = An array giving the heights (km) where the magnetic field will be modeled               By default the magnetic field will be computed at 100, 500 and 1000km.
                     alpha_i = Angle to interpolate the magnetic field.
                     Modification History
                     --------------------
                     Converted to Python by Freddy R. Galindo, ROJ, 07 October 2009.
                     """
                     handles = []
                     objects = []
                     colors = ['k','m','c','b','g','r','y']
                     marker = ['-+','-*','-D','-x','-s','->','-o','-^']
                     alpha_location = numpy.zeros((nlon,2,heights.size))
                     for ih in numpy.arange(heights.size):
                         alpha_location[:,0,ih] = dcos[:,0,ih,0]
                         for ilon in numpy.arange(nlon):
                             myx = (alpha[ilon,:,ih])[::-1]
                             myy = (dcos[ilon,:,ih,0])[::-1]
                             tck = scipy.interpolate.splrep(myx,myy,s=0)
                             mydcosx = scipy.interpolate.splev(alpha_i,tck,der=0)
                             myx = (alpha[ilon,:,ih])[::-1]
                             myy = (dcos[ilon,:,ih,1])[::-1]
                             tck = scipy.interpolate.splrep(myx,myy,s=0)
                             mydcosy = scipy.interpolate.splev(alpha_i,tck,der=0)
                             alpha_location[ilon,:,ih] = numpy.array([mydcosx, mydcosy])
                         ObjFig, = matplotlib.pyplot.plot(alpha_location[:,0,ih],alpha_location[:,1,ih], \
-                            marker[ih % 8],color=colors[numpy.int(ih/8)],ms=4.5,lw=0.5)
+                            marker[ih % 8],color=colors[int(ih/8)],ms=4.5,lw=0.5)
                         handles.append(ObjFig)
-                        objects.append(numpy.str(heights[ih]) + ' km')
+                        objects.append(str(heights[ih]) + ' km')
                     matplotlib.pyplot.xlim(dcosxrange[0],dcosxrange[1])
                     matplotlib.pyplot.ylim(dcosyrange[0],dcosyrange[1])
                     try:
                         ObjlegB = matplotlib.pyplot.legend(handles,objects,loc="lower right", numpoints=1, handlelength=0.3, \
                                             handletextpad=0.02, borderpad=0.3, labelspacing=0.1)
                     except:
                         ObjlegB = matplotlib.pyplot.legend(handles,objects,loc=[0.01,0.75], numpoints=1, handlelength=0, \
                         pad=0.015, handletextsep=0.02,labelsep=0.01)
                     matplotlib.pyplot.setp(ObjlegB.get_texts(),fontsize='small')
                     matplotlib.pyplot.gca().add_artist(ObjlegB)
                     save_fig = os.path.join(gpath,filename)
                     matplotlib.pyplot.savefig(save_fig,format='png')
                     self.alpha_location = alpha_location
             class CelestialObjectsPlot:
                 def __init__(self,jd,dec,tod,maxha_min,show_object=None):
                     self.jd = jd
                     self.dec = dec
                     self.tod = tod
                     self.maxha_min = maxha_min
                     if show_object==None:show_object=numpy.zeros(4)+2
                     self.show_object = show_object
                     self.dcosx_sun = 1
                     self.dcosy_sun = 1
                     self.ha_sun = 1
                     self.time_sun = 1
                     self.dcosx_moon = 1
                     self.dcosy_moon = 1
                     self.ha_moon = 1
                     self.time_moon = 1
                     self.dcosx_hydra = 1
                     self.dcosy_hydra = 1
                     self.ha_hydra = 1
                     self.time_hydra = 1
                     self.dcosx_galaxy = 1
                     self.dcosy_galaxy = 1
                     self.ha_galaxy = 1
                     self.time_galaxy = 1
                 def drawObject(self,glat,glon,xg,yg,dcosxrange,dcosyrange,gpath='',filename=''):
                     jd = self.jd
                     main_dec = self.dec
                     tod = self.tod
                     maxha_min = self.maxha_min
                     mesg = "Drawing celestial objects over Observatory"
             #        print mesg
             #        if textid!=None:textid.append(mesg)
                     maxlev = 24; minlev = 0; maxcol = 39; mincol = 10
                     handles = []
                     objects = ['$Sun$','$Moon$','$Hydra$','$Galaxy$']
                     marker = ['--^','--s','--*','--o']
                     # Getting RGB table to plot celestial object over Jicamarca
                     colortable = ColorTable(table=1).readTable()
                     for io in (numpy.arange(4)+1):
                         if self.show_object[io]!=0:
                             ObjBodies = CelestialBodies()
                             if io==1:
                                 [ra,dec,sunlon,sunobliq] = ObjBodies.sunpos(jd)
                             elif io==2:
                                 [ra,dec,dist,moonlon,moonlat] = ObjBodies.moonpos(jd)
                             elif io==3:
                                 [ra,dec] = ObjBodies.hydrapos()
                             elif io==4:
                                 [maxra,ra] = ObjBodies.skynoise_jro(dec_cut=main_dec)
                                 ra = maxra*15.
                                 dec = main_dec
                             ObjEq = Equatorial(ra,dec,jd,lat=glat,lon=glon)
                             [alt, az, ha] = ObjEq.change2AltAz()
                             vect = numpy.array([az,alt]).transpose()
                             vect = Vector(vect,direction=0).Polar2Rect()
                             dcosx = numpy.array(numpy.dot(vect,xg))
                             dcosy = numpy.array(numpy.dot(vect,yg))
                             wrap = numpy.where(ha>=180.)
                             if wrap[0].size>0:ha[wrap] = ha[wrap] - 360.
                             val = numpy.where((numpy.abs(ha))<=(maxha_min*0.25))
                             if val[0].size>2:
                                 tod_1 = tod*1.
                                 shift_1 = numpy.where(tod>12.)
                                 tod_1[shift_1] = tod_1[shift_1] - 24.
                                 tod_2 = tod*1.
                                 shift_2 = numpy.where(tod<12.)
                                 tod_2[shift_2] = tod_2[shift_2] + 24.
                                 diff0 = numpy.nanmax(tod[val])  -  numpy.nanmin(tod[val])
                                 diff1 = numpy.nanmax(tod_1[val]) - numpy.nanmin(tod_1[val])
                                 diff2 = numpy.nanmax(tod_2[val]) - numpy.nanmin(tod_2[val])
                                 if ((diff0<=diff1) & (diff0<=diff2)):
                                     tod_0 = tod
                                 elif ((diff1<diff0) & (diff1<diff2)):
                                     tod_0 = tod_1
                                 else:
                                     tod_0 = tod_2
                                 if io==1:
                                     self.dcosx_sun = dcosx[val]
                                     self.dcosy_sun = dcosy[val]
                                     self.ha_sun = ha[val]
                                     self.time_sun = numpy.median(tod_0[val])
                                 elif io==2:
                                     self.dcosx_moon = dcosx[val]
                                     self.dcosy_moon = dcosy[val]
                                     self.ha_moon = ha[val]
                                     self.time_moon = numpy.median(tod_0[val])
                                 elif io==3:
                                     self.dcosx_hydra = dcosx[val]
                                     self.dcosy_hydra = dcosy[val]
                                     self.ha_hydra = ha[val]
                                     self.time_hydra = numpy.mean(tod_0[val])
                                 elif io==4:
                                     self.dcosx_galaxy = dcosx[val]
                                     self.dcosy_galaxy = dcosy[val]
                                     self.ha_galaxy = ha[val]
                                     self.time_galaxy = numpy.mean(tod_0[val])
                                 index = numpy.mean(tod_0[val]) - minlev
                                 index = (index*(maxcol - mincol)/(maxlev - minlev)) + mincol
-                                index = numpy.int(index)
+                                index = int(index)
                                 figobjects, = matplotlib.pyplot.plot(dcosx[val],dcosy[val],marker[io-1],\
                                   lw=1,ms=7,mew=0,color=tuple(colortable[:,index]))
                                 handles.append(figobjects)
                     xmax = numpy.max(dcosxrange[1])
                     xmin = numpy.min(dcosxrange[0])
                     ymax = numpy.max(dcosyrange[1])
                     ymin = numpy.min(dcosyrange[0])
                     matplotlib.pyplot.xlim(xmin,xmax)
                     matplotlib.pyplot.ylim(ymin,ymax)
                     val = numpy.where(self.show_object[1:]>0)
                     objects = numpy.array(objects)
                     objects = list(objects[val])
                     try:
                         ObjlegC = matplotlib.pyplot.legend(handles,objects,loc="lower left", numpoints=1, handlelength=0.3, \
                         borderpad=0.3, handletextpad=0.02,labelspacing=0.1)
                     except:
                         ObjlegC = matplotlib.pyplot.legend(handles,objects,loc=[0.01,0.75], numpoints=1, handlelength=0, \
                         pad=0.015, handletextsep=0.02,labelsep=0.01)
                     matplotlib.pyplot.setp(ObjlegC.get_texts(),fontsize='small')
                     ObjlegC.isaxes = False
                     save_fig = os.path.join(gpath,filename)
                     matplotlib.pyplot.savefig(save_fig,format='png')
             class PatternCutPlot:
                 def __init__(self,nsubplots):
                     self.nsubplots = nsubplots
                     self.fig = None
                     self.__plot_width  = 8
                     if self.nsubplots == 5:
                         self.__plot_height = 11
                     if self.nsubplots == 4:
                         self.__plot_height = 9
                     if self.nsubplots == 3:
                         self.__plot_height = 7
                     if self.nsubplots == 2:
                         self.__plot_height = 5
                     if self.nsubplots == 1:
                         self.__plot_height = 3
                     self.fig = matplotlib.pyplot.figure(num = 4,figsize = (self.__plot_width, self.__plot_height))
                     if self.nsubplots < 5:
                         self.__height_inch = 1.1 #altura de los subplots (pulgadas)
                         top_inch = 1.5/2.7 #espacio entre el primer subplot y el limite superior del plot
                         self.__vspace_plot_inch = 1.0#1.5/2 # espacio vertical entre subplots
                         self.__left = 0.1
                     else:
                         self.__height_inch = 1.1 #altura de los subplots (pulgadas)
                         top_inch = 1.5/2.7 #espacio entre el primer subplot y el limite superior del plot
                         self.__vspace_plot_inch = 1.0 # espacio vertical entre subplots
                         self.__left = 0.1
                     self.__bottom_inch = self.__plot_height - (self.__height_inch + top_inch)
                     self.__height = self.__height_inch/self.__plot_height
                     self.__width = 0.8
                 def drawCut(self,io,patterns,npatterns,ha,otitle,subtitle,ptitle):
                     t_cuts = ['B','Sun','Moon','Hydra','Galaxy']
                     self.__bottom = self.__bottom_inch/self.__plot_height
                     subp = self.fig.add_axes([self.__left,self.__bottom,self.__width,self.__height])
                     on_axis_angle = -4.65562
                     for icut in numpy.arange(npatterns):
                         # Getting Antenna cut.
                         pattern = patterns[icut]
                         power = numpy.abs(pattern/numpy.nanmax(pattern))
                         max_power_db = numpy.round(10.*numpy.log10(numpy.nanmax(pattern)),2)
                         bval = numpy.where(power[:,0]==numpy.nanmax(power))
                         beta = -0.25*(ha[bval[0]] + on_axis_angle)
             #            print 'Angle (deg): '+"%f"%(beta)
                         subp.plot(ha,power)
                     xmax = numpy.max(numpy.nanmin(ha))
                     xmin = numpy.min(numpy.nanmax(ha))
                     ymax = numpy.max(1)
                     ymin = numpy.min(0)
                     subp.set_xlim(xmin, xmax)
                     subp.set_ylim(ymin, ymax)
                     subp.set_title(otitle + ' ' + ptitle,size="medium")
                     subp.text(0.5, 1.26,subtitle[0],
                      horizontalalignment='center',
                      verticalalignment='center',
                      transform = subp.transAxes)
                     xlabels = subp.get_xticks()
                     subp.set_xticklabels(xlabels,size="small")
                     ylabels = subp.get_yticks()
                     subp.set_yticklabels(ylabels,size="small")
                     subp.set_xlabel('Hour angle (min) (+ve to West)',size="small")
                     subp.set_ylabel("Power [Max: " + str(max_power_db) + ' dB]',size="small")
                     subp.grid()
                     self.__bottom_inch = self.__bottom_inch - (self.__height_inch + self.__vspace_plot_inch)
             class SkyNoisePlot:
                 def __init__(self,date,powr,time,time_lst):
                     """
                     SkyNoisePlot class creates an object which represents the SkyNoise Object to genera-
                     te a SkyNoise map.
                     Parameters
                     ----------
                     date = A List of 3 elements to define the desired date ([year, month, day]).
                     powr = An array giving the SkyNoise power for the desired time.
                     time = An array giving the number of seconds since 1970 to the desired time.
                     time_lst = Set this input to an array to define the Local Sidereal Time of the desi-
                       red time.
                     Modification History
                     --------------------
                     Created by Freddy Galindo, ROJ, 18 October 2009.
                     """
                     self.date = date
                     self.powr = powr
                     self.time = time
                     self.time_lst = time_lst

apps/abs/utils/overJroShow.py +13 -10

             #!/usr/bin/python
             import sys, os, os.path
             import traceback
             import cgi
             from http import cookies
             import time, datetime
             import types
             import numpy
             import numpy.fft
             import scipy.linalg
             import scipy.special
             from io import StringIO
             #import Misc_Routines
             from .Misc_Routines import CoFactors
             #import TimeTools
             from .TimeTools import Time , Julian ,Doy2Date
             #import JroAntSetup
             from .JroAntSetup import ReturnSetup
             #import Graphics_OverJro
             from .Graphics_OverJro import AntPatternPlot ,BFieldPlot,CelestialObjectsPlot,PatternCutPlot,SkyNoisePlot
             #import Astro_Coords
             from .Astro_Coords import Geodetic ,AltAz ,CelestialBodies
             class JroPattern():
                 def __init__(self,pattern=0,path=None,filename=None,nptsx=101,nptsy=101,maxphi=5,fftopt=0, \
                     getcut=0,dcosx=None,dcosy=None,eomwl=6,airwl=4, **kwargs):
                     """
                     JroPattern class creates an object to represent the useful parameters for beam mode-
                     lling of the Jicamarca VHF radar.
                     Parameters
                     ----------
                     pattern = An integer (See JroAntSetup to know the available values) to load a prede-
                       fined configuration. The default value  is 0. To use a  user-defined configuration
                       pattern  must be None.
                     path = A string giving the directory that contains the user-configuration file. PATH
                       will work if pattern is None.
                     filename = A string giving the  name of the  user-configuration  file. FILENAME will
                       work if pattern is None.
                     nptsx = A scalar to specify the number of points  used to define the angular resolu-
                       tion in the "x" axis. The default value is 101.
                     nptsy = A scalar to specify the number of points  used to define the angular resolu-
                       tion in the "x" axis. The default value is 101.
                     maxphi = A scalar giving the maximum (absolute) angle (in degree) to model the ante-
                       nna pattern. The default value is 5 degrees.
                     fftopt = Set this input  to 1 to model  the beam using FFT.  To model  using antenna
                       theory set to 0 (default value).
                     getcut = Set to 1 to show an antenna cut instead of a contour plot of itself (set to
 ). The defautl value is 0.
                     dcosx = An array giving the directional cosines for the  x-axis. DCOSX  will work if
                       getcut is actived.
                     dcosy = An array giving the directional cosines for the  y-axis. DCOSY  will work if
                       getcut is actived.
                     eomwl = A scalar giving the radar wavelength. The default value is 6m (50 MHZ).
                     airwl = Set this input to float (or intger) to specify the wavelength (in meters) of
                       the transmitted EOM wave in the air. The default value is 4m.
                     Modification History
                     --------------------
                     Converted to Object-oriented Programming by Freddy Galindo, ROJ, 20 September 2009.
                     """
                     # Getting antenna configuration.
                     if filename:
                         setup = ReturnSetup(path=path,filename=filename,pattern=pattern)
                         ues = setup["ues"]
                         phase = setup["phase"]
                         gaintx = setup["gaintx"]
                         gainrx = setup["gainrx"]
                         justrx = setup["justrx"]
                         self.title = setup["title"]
                     else:
                         ues = kwargs["ues"]
                         phase = kwargs["phases"]
                         gaintx = kwargs["gain_tx"]
                         gainrx = kwargs["gain_rx"]
                         justrx = kwargs["just_rx"]
                         self.title = kwargs.get("title", "JRO Pattern")
                     # Defining attributes for JroPattern class.
                     # Antenna configuration
                     self.uestx = ues
                     self.phasetx = phase
                     self.gaintx = gaintx
                     self.uesrx = ues
                     self.phaserx = phase
                     self.gainrx = gainrx
                     self.justrx = justrx
                     # Pattern resolution & method to model
                     self.maxphi = maxphi
                     self.nptsx = nptsx
                     self.nptsy = nptsy
                     self.fftopt = fftopt
                     # To get a cut of the pattern.
                     self.getcut = getcut
                     maxdcos = numpy.sin(maxphi*CoFactors.d2r)
                     if dcosx==None:dcosx = ((numpy.arange(nptsx,dtype=float)/(nptsx-1))-0.5)*2*maxdcos
                     if dcosy==None:dcosy = ((numpy.arange(nptsy,dtype=float)/(nptsy-1))-0.5)*2*maxdcos
                     self.dcosx = dcosx
                     self.dcosy = dcosy
                     self.nx = dcosx.size
                     self.ny = dcosy.size*(getcut==0) + (getcut==1)
                     self.eomwl = eomwl
                     self.airwl = airwl
                     self.kk = 2.*numpy.pi/eomwl
                     self.pattern = None
                     self.meanpos = None
                     self.norpattern = None
                     self.maxpattern = None
                     self.getPattern()
                 def getPattern(self):
                     """
                     getpattern method returns the modeled total antenna pattern and its mean position.
                     Return
                     ------
                     pattern = An array giving the Modelled antenna pattern.
                     mean_pos = A 2-elements array giving the mean position of the main beam.
                     Examples
                     --------
                     >> [pattern, mean_pos] = JroPattern(pattern=2).getPattern()
                     >> print meanpos
                     [  8.08728085e-14  -4.78193873e-14]
                     Modification history
                     --------------------
                     Developed by Jorge L. Chau.
                     Converted to Python by Freddy R. Galindo, ROJ, 20 September 2009.
                     """
                     if (self.fftopt>0) and (self.getcut>0):
                         #print  "Conflict bewteen fftopt and getcut"
                         #print  "To get a cut of the antenna pattern uses ffopt=0"
                         return None, None
                     if (self.fftopt==0):
                         # Getting antenna pattern using the array method
                         self.pattern = self.__usingArray(rx=1)
                         if (self.justrx==0):self.pattern = self.pattern*self.__usingArray(rx=0)
                     elif (self.fftopt>0):
                         # Getting antenna pattern using FFT method
                         self.pattern = self.__usingFFT(rx=1)
                         if (self.justrx==0):self.pattern = self.pattern*self.__usingFFT(rx=0)
                     self.maxpattern = numpy.nanmax(self.pattern)
                     self.norpattern = self.pattern/self.maxpattern
                     if self.getcut==0:self.__getBeamPars()
                 def __usingArray(self,rx):
                     """
                     __usingArray method returns the Jicamarca antenna pattern computed using array model
                     pattern = dipolepattern x modulepattern
                     Parameters
                     ----------
                     rx = Set to 1 to use the Rx information. Otherwise set to 0 for Tx.
                     Return
                     ------
                     pattern = An array giving the modelled antenna pattern using the array model.
                     Modification history
                     --------------------
                     Developed by Jorge L. Chau.
                     Converted to Python by Freddy R. Galindo, ROJ, 20 September 2009.
                     """
                     if rx==1:
                         ues = self.uesrx
                         phase = self.phaserx
                         gain = self.gainrx
                     elif rx==0:
                         ues = self.uestx
                         phase = self.phasetx
                         gain = self.gaintx
                     ues = ues*360./self.airwl
                     phase = phase*360./self.airwl
                     for ii in range(4):
                         if ii==0:dim = numpy.array([4,0,8,4])   # WEST
                         elif ii==1:dim = numpy.array([0,0,4,4]) # NORTH
                         elif ii==2:dim = numpy.array([0,4,4,8]) # EAST
                         elif ii==3:dim = numpy.array([4,4,8,8]) # SOUTH
                         xi = dim[0]; xf = dim[2]; yi = dim[1]; yf = dim[3]
                         phase[xi:xf,yi:yf] = phase[xi:xf,yi:yf] + ues[ii]
                     phase = -phase
                     ar = self.eomwl*numpy.array([[0.5,6., 24.5],[0.5,6.,24.5]])
                     nr = numpy.array([[12.,4.,2.],[12.,4.,2.]])
                     lr = 0.25*self.eomwl*numpy.array([[0,0.,0],[0.,0,0]])
                     # Computing module and dipole patterns.
                     pattern = (numpy.abs(self.__dipPattern(ar,nr,lr)*self.__modPattern(phase,gain)))**2
                     return pattern
                 def __usingFFT(self,rx):
                     """
                     __usingFFT method returns the Jicamarca antenna pattern computed using The Fast Fou-
                     rier Transform.
                     pattern = iFFT(FFT(gain*EXP(j*phase)))
                     Parameters
                     ----------
                     rx = Set to 1 to use the Rx information. Otherwise set to 0 for Tx.
                     Return
                     ------
                     pattern = An array giving the modelled antenna pattern using the array model.
                     Modification history
                     --------------------
                     Developed by Jorge L. Chau.
                     Converted to Python by Freddy R. Galindo, ROJ, 20 September 2009.
+                    Updated to Django 4.1.5 and Python 3.9.16 by Renato Huallpa M. March 2023.
                     """
                     if rx==1:
                         ues = self.uesrx
                         phase = self.phaserx
                         gain = self.gainrx
                     elif rx==0:
                         ues = self.uestx
                         phase = self.phasetx
                         gain = self.gaintx
                     ues = ues*360./self.airwl
                     phase = phase*360./self.airwl
                     for ii in range(4):
                         if ii==0:dim = numpy.array([4,0,8,4])   # WEST
                         elif ii==1:dim = numpy.array([0,0,4,4]) # NORTH
                         elif ii==2:dim = numpy.array([0,4,4,8]) # EAST
                         elif ii==3:dim = numpy.array([4,4,8,8]) # SOUTH
                         xi = dim[0]; xf = dim[2]; yi = dim[1]; yf = dim[3]
                         phase[xi:xf,yi:yf] = phase[xi:xf,yi:yf] + ues[ii]
                     phase = -phase
                     delta_x = self.eomwl/2.
                     delta_y = self.eomwl/2.
                     nxfft = 2048
                     nyfft = 2048
                     dcosx = (numpy.arange(nxfft) - (0.5*nxfft))/(nxfft*delta_x)*self.eomwl
                     dcosy = (numpy.arange(nyfft) - (0.5*nyfft))/(nyfft*delta_y)*self.eomwl
                     fft_gain = numpy.zeros((nxfft,nyfft))
                     fft_phase = numpy.zeros((nxfft,nyfft))
                     nx = 8
                     ny = 8
                     ndx =12
                     ndy =12
                     for iy in numpy.arange(ny):
                         for ix in numpy.arange(nx):
                             ix1 = nxfft/2-self.nx/2*ndx+ix*ndx
                             if ix<(nx/2):ix1 = ix1 - 1
                             if ix>=(nx/2):ix1 = ix1 + 1
                             iy1 = nyfft/2-ny/2*ndx+iy*ndy
                             if iy<(ny/2):iy1 = iy1 - 1
                             if iy>=(ny/2):iy1 = iy1 + 1
                             fft_gain[ix1:ix1+ndx-1,iy1:iy1+ndy-1] = gain[ix,ny-1-iy]
                             fft_phase[ix1:ix1+ndx-1,iy1:iy1+ndy-1] = phase[ix,ny-1-iy]
                     fft_phase = fft_phase*CoFactors.d2r
-                    pattern = numpy.abs(numpy.fft.fft2(fft_gain*numpy.exp(numpy.complex(0,1)*fft_phase)))**2
+                    pattern = numpy.abs(numpy.fft.fft2(fft_gain*numpy.exp(complex(0,1)*fft_phase)))**2
                     pattern = numpy.fft.fftshift(pattern)
                     xvals = numpy.where((dcosx>=(numpy.min(self.dcosx))) & (dcosx<=(numpy.max(self.dcosx))))
                     yvals = numpy.where((dcosy>=(numpy.min(self.dcosy))) & (dcosy<=(numpy.max(self.dcosy))))
                     pattern = pattern[xvals[0][0]:xvals[0][-1],yvals[0][0]:yvals[0][-1]]
                     return pattern
                 def __readAttenuation(self):
                     """
                     _readAttenuation reads the attenuations' file and returns an array giving  these va-
                     lues (dB). The ext file must be in the directory "resource".
                     Return
                     ------
                     attenuation = An array giving attenuation values read from the text file.
                     Modification history
                     --------------------
                     Developed by Jorge L. Chau.
                     Converted to Python by Freddy R. Galindo, ROJ, 20 September 2009.
                     """
             #         attenuation = None
             # #        foldr = sys.path[-1] + os.sep + "resource" + os.sep
             #         base_path = os.path.dirname(os.path.abspath(__file__))
             #         #foldr = './resource'
             #         #filen = "attenuation.txt"
             #         attenuationFile = os.path.join(base_path,"resource","attenuation.txt")
             #         #ff = open(os.path.join(foldr,filen),'r')
             #         ff = open(attenuationFile,'r')
             #         print(ff.read())
             #         exec(ff.read())
             #         ff.close()
                     attenuation = numpy.array([[[-21.25,-15.25,-9.25,-3.25,03.25,09.25,15.25,21.25],
             		         [-21.25,-15.25,-9.25,-3.25,03.25,09.25,15.25,21.25],
                		         [-21.25,-15.25,-9.25,-3.25,03.25,09.25,15.25,21.25],
             		         [-21.25,-15.25,-9.25,-3.25,03.25,09.25,15.25,21.25],
             		         [-21.25,-15.25,-9.25,-3.25,03.25,09.25,15.25,21.25],
             		         [-21.25,-15.25,-9.25,-3.25,03.25,09.25,15.25,21.25],
             		         [-21.25,-15.25,-9.25,-3.25,03.25,09.25,15.25,21.25],
             		         [-21.25,-15.25,-9.25,-3.25,03.25,09.25,15.25,21.25]],
                    		         [[21.25,21.25,21.25,21.25,21.25,21.25,21.25,21.25],
             		         [15.25,15.25,15.25,15.25,15.25,15.25,15.25,15.25],
             		         [09.25,09.25,09.25,09.25,09.25,09.25,09.25,09.25],
             		         [03.25,03.25,03.25,03.25,03.25,03.25,03.25,03.25],
             		         [-03.25,-03.25,-03.25,-03.25,-03.25,-03.25,-03.25,-03.25],
             		         [-09.25,-09.25,-09.25,-09.25,-09.25,-09.25,-09.25,-09.25],
             		         [-15.25,-15.25,-15.25,-15.25,-15.25,-15.25,-15.25,-15.25],
                 		         [-21.25,-21.25,-21.25,-21.25,-21.25,-21.25,-21.25,-21.25]]])
                     return attenuation
                 def __dipPattern(self,ar,nr,lr):
                     """
                     _dipPattern function computes the dipole's pattern to  the Jicamarca radar. The next
                     equation defines the pattern as a function of the mainlobe direction:
                     sincx = SIN(k/2*n0x*(a0x*SIN(phi)*COS(alpha)))/SIN(k/2*(a0x*SIN(phi)*COS(alpha)))
                     sincy = SIN(k/2*n0y*(a0y*SIN(phi)*SIN(alpha)))/SIN(k/2*(a0y*SIN(phi)*SIN(alpha)))
                     A0(phi,alpha) = sincx*sincy
                     Parameters
                     ----------
                     ar = ?
                     nr = ?
                     lr = ?
                     Return
                     ------
                     dipole = An array giving antenna pattern from the dipole point of view..
                     Modification history
                     --------------------
                     Developed by Jorge L. Chau.
                     Converted to Python by Freddy R. Galindo, ROJ, 20 September 2009.
                     """
                     dipole = numpy.zeros((self.nx,self.ny),dtype=complex)
                     for iy in range(self.ny):
                         for ix in range(self.nx):
                             yindex = iy*(self.getcut==0) + ix*(self.getcut==1)
                             argx = ar[0,0]*self.dcosx[ix] - lr[0,0]
                             if argx == 0.0:
                                 junkx = nr[0,0]
                             else:
                                 junkx = numpy.sin(0.5*self.kk*nr[0,0]*argx)/numpy.sin(0.5*self.kk*argx)
                             argy = ar[1,0]*self.dcosy[yindex] - lr[1,0]
                             if argy == 0.0:
                                 junky = nr[1,0]
                             else:
                                 junky = numpy.sin(0.5*self.kk*nr[1,0]*argy)/numpy.sin(0.5*self.kk*argy)
                             dipole[ix,iy] = junkx*junky
                     return dipole
                 def __modPattern(self,phase,gain):
                     """
                     ModPattern computes the module's pattern to the Jicamarca radar.  The next equation
                     defines    the pattern as a function mainlobe direction:
                     phasex = pos(x)*SIN(phi)*COS(alpha)
                     phasey = pos(y)*SIN(phi)*SIN(alpha)
                     A1(phi,alpha) = TOTAL(gain*EXP(COMPLEX(0,k*(phasex+phasey)+phase)))
                     Parameters
                     ----------
                     phase = Bidimensional array (8x8) giving the phase (in meters) of each module.
                     gain  = Bidimensional array (8x8) giving to  define modules  will be active  (ones)
                       and which will not (zeros).
                     Return
                     ------
                     module = An array giving antenna pattern from the module point of view..
                     Modification history
                     --------------------
                     Developed by Jorge L. Chau.
                     Converted to Python by Freddy R. Galindo, ROJ, 20 September 2009.
+                    Updated to Django 4.1.5 and Python 3.9.16 by Renato Huallpa M. March 2023.
                     """
                     pos  = self.eomwl*self.__readAttenuation()
                     posx = pos[0,:,:]
                     posy = pos[1,:,:]
                     phase = phase*CoFactors.d2r
                     module = numpy.zeros((self.nx,self.ny),dtype=complex)
                     for iy in range(self.ny):
                         for ix in range(self.nx):
                             yindex = iy*(self.getcut==0) + ix*(self.getcut==1)
                             phasex = posx*self.dcosx[ix]
                             phasey = posy*self.dcosy[yindex]
-                            tmp = gain*numpy.exp(numpy.complex(0,1.)*(self.kk*(phasex+phasey)+phase))
+                            tmp = gain*numpy.exp(complex(0,1.)*(self.kk*(phasex+phasey)+phase))
                             module[ix,iy] = tmp.sum()
                     return module
                 def __getBeamPars(self):
                     """
                     _getBeamPars computes the main-beam parameters of the antenna.
                     Modification history
                     --------------------
                     Developed by Jorge L. Chau.
                     Converted to Python by Freddy R. Galindo, ROJ, 20 September 2009.
                     """
                     dx = self.dcosx[1] - self.dcosx[0]
                     dy = self.dcosy[1] - self.dcosy[0]
                     amp = self.norpattern
                     xx =  numpy.resize(self.dcosx,(self.nx,self.nx)).transpose()
                     yy =  numpy.resize(self.dcosy,(self.ny,self.ny))
                     mm0 = amp[numpy.where(amp > 0.5)]
                     xx0 = xx[numpy.where(amp > 0.5)]
                     yy0 = yy[numpy.where(amp > 0.5)]
                     xc = numpy.sum(mm0*xx0)/numpy.sum(mm0)
                     yc = numpy.sum(mm0*yy0)/numpy.sum(mm0)
                     rc = numpy.sqrt(mm0.size*dx*dy/numpy.pi)
                     nnx = numpy.where(numpy.abs(self.dcosx - xc) < rc)
                     nny = numpy.where(numpy.abs(self.dcosy - yc) < rc)
                     mm1 = amp[numpy.min(nnx):numpy.max(nnx)+1,numpy.min(nny):numpy.max(nny)+1]
                     xx1 = self.dcosx[numpy.min(nnx):numpy.max(nnx)+1]
                     yy1 = self.dcosy[numpy.min(nny):numpy.max(nny)+1]
                     # fitting data into the main beam.
-                    import gaussfit
+                    #import apps.abs.utils.gaussfit as gaussfit
-                    params = gaussfit.fitgaussian(mm1)
+                    from .gaussfit import fitgaussian
+                    params = fitgaussian(mm1)
                     # Tranforming from indexes to axis' values
                     xcenter = xx1[0] + (((xx1[xx1.size-1] - xx1[0])/(xx1.size -1))*(params[1]))
                     ycenter = yy1[0] + (((yy1[yy1.size-1] - yy1[0])/(yy1.size -1))*(params[2]))
                     xwidth  = ((xx1[xx1.size-1] - xx1[0])/(xx1.size-1))*(params[3])*(1/CoFactors.d2r)
                     ywidth  = ((yy1[yy1.size-1] - yy1[0])/(yy1.size-1))*(params[4])*(1/CoFactors.d2r)
                     meanwx = (xwidth*ywidth)
                     meanpos = numpy.array([xcenter,ycenter])
                     #print  'Position: %f %f' %(xcenter,ycenter)
                     #print  'Widths:   %f %f' %(xwidth, ywidth)
                     #print  'BWHP:     %f' %(2*numpy.sqrt(2*meanwx)*numpy.sqrt(-numpy.log(0.5)))
                     self.meanpos = meanpos
             class BField():
                 def __init__(self,year=None,doy=None,site=1,heights=None,alpha_i=90):
                     """
                     BField class creates an object to get  the Magnetic  field for a  specific date and
                     height(s).
                     Parameters
                     ----------
                     year = A scalar giving the desired year.  If the value is None (default value) then
                       the current year will be used.
                     doy = A scalar giving the desired day of the year. If the value is None (default va-
                       lue) then the current doy will be used.
                     site = An integer to choose the geographic coordinates of the place where the magne-
                       tic field will be computed. The default value is over Jicamarca (site=1)
                     heights = An array giving the heights (km) where the magnetic field will be modeled By default the magnetic field will be computed at 100, 500 and 1000km.
                     alpha_i = Angle to interpolate the magnetic field.
                     Modification History
                     --------------------
                     Converted to Object-oriented Programming by Freddy Galindo, ROJ, 07 October 2009.
                     """
                     tmp = time.localtime()
                     if year==None: year = tmp[0]
                     if doy==None:  doy = tmp[7]
                     self.year = year
                     self.doy = doy
                     self.site = site
                     if heights==None:heights = numpy.array([100,500,1000])
                     self.heights = heights
                     self.alpha_i = alpha_i
                 def getBField(self,maglimits=numpy.array([-7,-7,7,7])):
                     """
                     getBField models the magnetic field for a different heights in a specific date.
                     Parameters
                     ----------
                     maglimits = An 4-elements array giving ..... The default value is [-7,-7,7,7].
                     Return
                     ------
                     dcos = An 4-dimensional array giving the directional cosines of the magnetic field
                       over the desired place.
                     alpha = An 3-dimensional array giving the angle of the magnetic field over the desi-
                       red place.
                     Modification History
                     --------------------
                     Converted to Python by Freddy R. Galindo, ROJ, 07 October 2009.
                     """
                     x_ant = numpy.array([1,0,0])
                     y_ant = numpy.array([0,1,0])
                     z_ant = numpy.array([0,0,1])
                     if self.site==0:
                         title_site = "Magnetic equator"
                         coord_site = numpy.array([-76+52./60.,-11+57/60.,0.5])
                     elif self.site==1:
                         title_site = 'Jicamarca'
                         coord_site = [-76-52./60.,-11-57/60.,0.5]
                         theta = (45+5.35)*numpy.pi/180.      # (50.35 and 1.46 from Fleish Thesis)
                         delta = -1.46*numpy.pi/180
                         x_ant1 = numpy.roll(self.rotvector(self.rotvector(x_ant,1,delta),3,theta),1)
                         y_ant1 = numpy.roll(self.rotvector(self.rotvector(y_ant,1,delta),3,theta),1)
                         z_ant1 = numpy.roll(self.rotvector(self.rotvector(z_ant,1,delta),3,theta),1)
                         ang0 = -1*coord_site[0]*numpy.pi/180.
                         ang1 = coord_site[1]*numpy.pi/180.
                         x_ant = self.rotvector(self.rotvector(x_ant1,2,ang1),3,ang0)
                         y_ant = self.rotvector(self.rotvector(y_ant1,2,ang1),3,ang0)
                         z_ant = self.rotvector(self.rotvector(z_ant1,2,ang1),3,ang0)
                     else:
             #            print "No defined Site. Skip..."
                         return None
                     nhei = self.heights.size
                     pt_intercep = numpy.zeros((nhei,2))
                     nfields = 1
                     grid_res = 0.5
-                    nlon = int(numpy.int(maglimits[2] - maglimits[0])/grid_res + 1)
+                    nlon = int(int(maglimits[2] - maglimits[0])/grid_res + 1)
-                    nlat = int(numpy.int(maglimits[3] - maglimits[1])/grid_res + 1)
+                    nlat = int(int(maglimits[3] - maglimits[1])/grid_res + 1)
                     location = numpy.zeros((nlon,nlat,2))
                     mlon = numpy.atleast_2d(numpy.arange(nlon)*grid_res + maglimits[0])
                     mrep = numpy.atleast_2d(numpy.zeros(nlat) + 1)
                     location0 = numpy.dot(mlon.transpose(),mrep)
                     mlat = numpy.atleast_2d(numpy.arange(nlat)*grid_res + maglimits[1])
                     mrep = numpy.atleast_2d(numpy.zeros(nlon) + 1)
                     location1 = numpy.dot(mrep.transpose(),mlat)
                     location[:,:,0] = location0
                     location[:,:,1] = location1
                     alpha = numpy.zeros((nlon,nlat,nhei))
                     rr = numpy.zeros((nlon,nlat,nhei,3))
                     dcos = numpy.zeros((nlon,nlat,nhei,2))
                     global first_time
                     first_time = None
                     for ilon in numpy.arange(nlon):
                         for ilat in numpy.arange(nlat):
                             outs = self.__bdotk(self.heights,
                                                 self.year + self.doy/366.,
                                                 coord_site[1],
                                                 coord_site[0],
                                                 coord_site[2],
                                                 coord_site[1]+location[ilon,ilat,1],
                                                 location[ilon,ilat,0]*720./180.)
                             alpha[ilon, ilat,:] = outs[1]
                             rr[ilon, ilat,:,:] = outs[3]
                             mrep = numpy.atleast_2d((numpy.zeros(nhei)+1)).transpose()
                             tmp = outs[3]*numpy.dot(mrep,numpy.atleast_2d(x_ant))
                             tmp = tmp.sum(axis=1)
                             dcos[ilon,ilat,:,0] = tmp/numpy.sqrt((outs[3]**2).sum(axis=1))
                             mrep = numpy.atleast_2d((numpy.zeros(nhei)+1)).transpose()
                             tmp = outs[3]*numpy.dot(mrep,numpy.atleast_2d(y_ant))
                             tmp = tmp.sum(axis=1)
                             dcos[ilon,ilat,:,1] = tmp/numpy.sqrt((outs[3]**2).sum(axis=1))
                     return dcos, alpha, nlon, nlat
                 def __bdotk(self,heights,tm,gdlat=-11.95,gdlon=-76.8667,gdalt=0.0,decd=-12.88, ham=-4.61666667):
                     global first_time
                     # Mean Earth radius in Km WGS 84
                     a_igrf = 6371.2
                     bk = numpy.zeros(heights.size)
                     alpha = numpy.zeros(heights.size)
                     bfm = numpy.zeros(heights.size)
                     rr = numpy.zeros((heights.size,3))
                     rgc = numpy.zeros((heights.size,3))
-                    ObjGeodetic = Astro_Coords.Geodetic(gdlat,gdalt)
+                    ObjGeodetic = Geodetic(gdlat,gdalt)
                     [gclat,gcalt] = ObjGeodetic.change2geocentric()
                     gclat = gclat*numpy.pi/180.
                     gclon = gdlon*numpy.pi/180.
                     # Antenna position from center of Earth
                     ca_vector = [numpy.cos(gclat)*numpy.cos(gclon),numpy.cos(gclat)*numpy.sin(gclon),numpy.sin(gclat)]
                     ca_vector = gcalt*numpy.array(ca_vector)
                     dec = decd*numpy.pi/180.
                     # K  vector respect to the center of earth.
                     klon = gclon + ham*numpy.pi/720.
                     k_vector = [numpy.cos(dec)*numpy.cos(klon),numpy.cos(dec)*numpy.sin(klon),numpy.sin(dec)]
                     k_vector = numpy.array(k_vector)
                     for ih in numpy.arange(heights.size):
                         # Vector from Earth's center to volume of interest
                         rr[ih,:] = k_vector*heights[ih]
                         cv_vector = numpy.squeeze(ca_vector) + rr[ih,:]
                         cv_gcalt = numpy.sqrt(numpy.sum(cv_vector**2.))
                         cvxy = numpy.sqrt(numpy.sum(cv_vector[0:2]**2.))
                         radial = cv_vector/cv_gcalt
                         east = numpy.array([-1*cv_vector[1],cv_vector[0],0])/cvxy
                         comp1 = east[1]*radial[2] - radial[1]*east[2]
                         comp2 = east[2]*radial[0] - radial[2]*east[0]
                         comp3 = east[0]*radial[1] - radial[0]*east[1]
                         north = -1*numpy.array([comp1, comp2, comp3])
                         rr_k = cv_vector - numpy.squeeze(ca_vector)
                         u_rr = rr_k/numpy.sqrt(numpy.sum(rr_k**2.))
                         cv_gclat = numpy.arctan2(cv_vector[2],cvxy)
                         cv_gclon = numpy.arctan2(cv_vector[1],cv_vector[0])
                         bhei = cv_gcalt-a_igrf
                         blat = cv_gclat*180./numpy.pi
                         blon = cv_gclon*180./numpy.pi
                         bfield = self.__igrfkudeki(bhei,tm,blat,blon)
                         B = (bfield[0]*north + bfield[1]*east - bfield[2]*radial)*1.0e-5
                         bfm[ih] = numpy.sqrt(numpy.sum(B**2.)) #module
                         bk[ih] = numpy.sum(u_rr*B)
                         alpha[ih] = numpy.arccos(bk[ih]/bfm[ih])*180/numpy.pi
                         rgc[ih,:] = numpy.array([cv_gclon, cv_gclat, cv_gcalt])
                     return bk, alpha, bfm, rr, rgc
                 def __igrfkudeki(self,heights,time,latitude,longitude,ae=6371.2):
                     """
                     __igrfkudeki calculates the International Geomagnetic Reference Field for given in-
                     put conditions based on IGRF2005 coefficients.
                     Parameters
                     ----------
                     heights = Scalar or vector giving the height above the Earth  of the point in ques-
                       tion in kilometers.
                     time = Scalar or vector giving the decimal year of time in question (e.g. 1991.2).
                     latitude = Latitude of point in question in decimal degrees. Scalar or vector.
                     longitude = Longitude of point in question in decimal degrees. Scalar or vector.
                     ae =
                     first_time =
                     Return
                     ------
                     bn =
                     be =
                     bd =
                     bmod =
                     balpha =
                     first_time =
                     Modification History
                     --------------------
                     Converted to Python by Freddy R. Galindo, ROJ, 03 October 2009.
                     """
                     global first_time
                     global gs, hs, nvec, mvec, maxcoef
                     heights = numpy.atleast_1d(heights)
                     time = numpy.atleast_1d(time)
                     latitude = numpy.atleast_1d(latitude)
                     longitude = numpy.atleast_1d(longitude)
                     if numpy.max(latitude)==90:
             #            print "Field calculations are not supported at geographic poles"
                         pass
                     # output arrays
                     bn = numpy.zeros(heights.size)
                     be = numpy.zeros(heights.size)
                     bd = numpy.zeros(heights.size)
                     if first_time==None:first_time=0
                     time0 = time[0]
                     if time!=first_time:
                         #print "Getting coefficients for", time0
                         [periods,g,h ] = self.__readIGRFcoeff()
                         top_year = numpy.max(periods)
                         nperiod = (top_year - 1900)/5 + 1
                         maxcoef = 10
                         if time0>=2000:maxcoef = 12
                         # Normalization array for Schmidt fucntions
                         multer = numpy.zeros((2+maxcoef,1+maxcoef)) + 1
                         for cn in (numpy.arange(maxcoef)+1):
                             for rm in (numpy.arange(cn)+1):
                                 tmp = numpy.arange(2*rm) + cn - rm + 1.
                                 multer[rm+1,cn] = ((-1.)**rm)*numpy.sqrt(2./tmp.prod())
                         schmidt = multer[1:,1:].transpose()
                         # n and m arrays
                         nvec = numpy.atleast_2d(numpy.arange(maxcoef)+2)
                         mvec = numpy.atleast_2d(numpy.arange(maxcoef+1)).transpose()
                         # Time adjusted igrf g and h with Schmidt normalization
                         # IGRF coefficient arrays: g0(n,m), n=1, maxcoeff,m=0, maxcoeff, ...
                         if time0<top_year:
                             dtime = (time0 - 1900) % 5
                             ntime = (time0 - 1900 - dtime)/5
                         else:
                             # Estimating coefficients for times > top_year
                             dtime = (time0 - top_year) + 5
                             ntime = g[:,0,0].size - 2
                         g0 = g[ntime,1:maxcoef+1,:maxcoef+1]
                         h0 = h[ntime,1:maxcoef+1,:maxcoef+1]
                         gdot = g[ntime+1,1:maxcoef+1,:maxcoef+1]-g[ntime,1:maxcoef+1,:maxcoef+1]
                         hdot = h[ntime+1,1:maxcoef+1,:maxcoef+1]-h[ntime,1:maxcoef+1,:maxcoef+1]
                         gs = (g0 + dtime*(gdot/5.))*schmidt[:maxcoef,0:maxcoef+1]
                         hs = (h0 + dtime*(hdot/5.))*schmidt[:maxcoef,0:maxcoef+1]
                         first_time = time0
                     for ii in numpy.arange(heights.size):
                         # Height dependence array rad = (ae/(ae+height))**(n+3)
                         rad = numpy.atleast_2d((ae/(ae + heights[ii]))**(nvec+1))
                         # Sin and Cos of m times longitude phi arrays
                         mphi = mvec*longitude[ii]*numpy.pi/180.
                         cosmphi = numpy.atleast_2d(numpy.cos(mphi))
                         sinmphi = numpy.atleast_2d(numpy.sin(mphi))
                         # Cos of colatitude theta
                         c = numpy.cos((90 - latitude[ii])*numpy.pi/180.)
                         # Legendre functions p(n,m|c)
                         [p,dp]= scipy.special.lpmn(maxcoef+1,maxcoef+1,c)
                         p = p[:,:-1].transpose()
                         s = numpy.sqrt((1. - c)*(1 + c))
                         # Generate derivative array dpdtheta = -s*dpdc
                         dpdtheta = c*p/s
                         for m in numpy.arange(maxcoef+2):    dpdtheta[:,m] = m*dpdtheta[:,m]
                         dpdtheta = dpdtheta + numpy.roll(p,-1,axis=1)
                         # Extracting arrays required for field calculations
                         p = p[1:maxcoef+1,:maxcoef+1]
                         dpdtheta = dpdtheta[1:maxcoef+1,:maxcoef+1]
                         # Weigh p and dpdtheta with gs and hs coefficients.
                         gp = gs*p
                         hp = hs*p
                         gdpdtheta = gs*dpdtheta
                         hdpdtheta = hs*dpdtheta
                         # Calcultate field components
                         matrix0 = numpy.dot(gdpdtheta,cosmphi)
                         matrix1 = numpy.dot(hdpdtheta,sinmphi)
                         bn[ii] = numpy.dot(rad,(matrix0 + matrix1))
                         matrix0 = numpy.dot(hp,(mvec*cosmphi))
                         matrix1 = numpy.dot(gp,(mvec*sinmphi))
                         be[ii] = numpy.dot((-1*rad),((matrix0 - matrix1)/s))
                         matrix0 = numpy.dot(gp,cosmphi)
                         matrix1 = numpy.dot(hp,sinmphi)
                         bd[ii] = numpy.dot((-1*nvec*rad),(matrix0 + matrix1))
                     bmod = numpy.sqrt(bn**2. + be**2. + bd**2.)
                     btheta = numpy.arctan(bd/numpy.sqrt(be**2. + bn**2.))*180/numpy.pi
                     balpha = numpy.arctan(be/bn)*180./numpy.pi
                     #bn : north
                     #be : east
                     #bn : radial
                     #bmod : module
                     return bn, be, bd, bmod, btheta, balpha
                 def str2num(self, datum):
                     try:
                         return int(datum)
                     except:
                         try:
                             return float(datum)
                         except:
                             return datum
                 def __readIGRFfile(self, filename):
                     list_years=[]
                     for i in range(1,24):
                         list_years.append(1895.0 + i*5)
                     epochs=list_years
                     epochs.append(epochs[-1]+5)
                     nepochs = numpy.shape(epochs)
                     gg = numpy.zeros((13,14,nepochs[0]),dtype=float)
                     hh = numpy.zeros((13,14,nepochs[0]),dtype=float)
                     coeffs_file=open(filename)
                     lines=coeffs_file.readlines()
                     coeffs_file.close()
                     for line in lines:
                         items = line.split()
                         g_h = items[0]
                         n = self.str2num(items[1])
                         m = self.str2num(items[2])
                         coeffs = items[3:]
                         for i in range(len(coeffs)-1):
                             coeffs[i] = self.str2num(coeffs[i])
                         #coeffs = numpy.array(coeffs)
                         ncoeffs = numpy.shape(coeffs)[0]
                         if g_h == 'g':
                 #            print n," g ",m
                             gg[n-1,m,:]=coeffs
                         elif g_h=='h':
                 #            print n," h ",m
                             hh[n-1,m,:]=coeffs
                 #        else :
                 #            continue
                 #    Ultimo Reordenamiento para  almacenar .
                     gg[:,:,nepochs[0]-1] = gg[:,:,nepochs[0]-2] + 5*gg[:,:,nepochs[0]-1]
                     hh[:,:,nepochs[0]-1] = hh[:,:,nepochs[0]-2] + 5*hh[:,:,nepochs[0]-1]
             #        return numpy.array([gg,hh])
                     periods = numpy.array(epochs)
                     g = gg
                     h = hh
                     return periods, g, h
                 def __readIGRFcoeff(self,filename="igrf10coeffs.dat"):
                     """
                     __readIGRFcoeff reads the coefficients from  a binary file which is located  in the
                     folder "resource."
                     Parameter
                     ---------
                     filename = A string to specify the name of the file which contains thec coeffs. The
                     default value is "igrf10coeffs.dat"
                     Return
                     ------
                     periods = A lineal array giving...
                     g1 =
                     h1 =
                     Modification History
                     --------------------
                     Converted to Python by Freddy R. Galindo, ROJ, 03 October 2009.
                     """
             # #        igrfile = sys.path[-1] + os.sep + "resource" + os.sep + filename
             #         igrfile = os.path.join('./resource',filename)
             #         f = open(igrfile,'rb')
             #         #f = open(os.getcwd() + os.sep + "resource" + os.sep + filename,'rb')
             #
             #         # Reading SkyNoise Power (lineal scale)
             #         periods = numpy.fromfile(f,numpy.dtype([('var','<f4')]),23)
             #         periods = periods['var']
             #
             #         g = numpy.fromfile(f,numpy.dtype([('var','<f8')]),23*14*14)
             #         g = g['var'].reshape((14,14,23)).transpose()
             #
             #         h = numpy.fromfile(f,numpy.dtype([('var','<f8')]),23*14*14)
             #         h = h['var'].reshape((14,14,23)).transpose()
             #
             #         f.close()
                     base_path = os.path.dirname(os.path.abspath(__file__))
                     filename = os.path.join(base_path,"resource","igrf11coeffs.txt")
                     period_v, g_v, h_v = self.__readIGRFfile(filename)
                     g2 = numpy.zeros((14,14,24))
                     h2 = numpy.zeros((14,14,24))
                     g2[1:14,:,:] = g_v
                     h2[1:14,:,:] = h_v
                     g = numpy.transpose(g2, (2,0,1))
                     h = numpy.transpose(h2, (2,0,1))
                     periods = period_v.copy()
                     return periods, g, h
                 def rotvector(self,vector,axis=1,ang=0):
                     """
                     rotvector function returns the new vector generated rotating the rectagular coords.
                     Parameters
                     ----------
                     vector = A lineal 3-elements array (x,y,z).
                     axis = A integer to specify  the axis used to rotate the coord systems. The default
                       value is 1.
                         axis = 1 -> Around "x"
                         axis = 2 -> Around "y"
                         axis = 3 -> Around "z"
                     ang = Angle of rotation (in radians). The default value is zero.
                     Return
                     ------
                     rotvector = A lineal array of 3 elements giving the new coordinates.
                     Modification History
                     --------------------
                     Converted to Python by Freddy R. Galindo, ROJ, 01 October 2009.
                     """
                     if axis==1:
                         t = [[1,0,0],[0,numpy.cos(ang),numpy.sin(ang)],[0,-numpy.sin(ang),numpy.cos(ang)]]
                     elif axis==2:
                         t = [[numpy.cos(ang),0,-numpy.sin(ang)],[0,1,0],[numpy.sin(ang),0,numpy.cos(ang)]]
                     elif axis==3:
                         t = [[numpy.cos(ang),numpy.sin(ang),0],[-numpy.sin(ang),numpy.cos(ang),0],[0,0,1]]
                     rotvector = numpy.array(numpy.dot(numpy.array(t),numpy.array(vector)))
                     return rotvector
             class overJroShow:
             #    __serverdocspath = '/usr/local/www/htdocs'
             #    __tmpDir = 'overJro/tempReports'
             #     __serverdocspath = '/Users/dsuarez/Pictures'
             #     __tmpDir = 'overjro'
                 __serverdocspath = ''
                 __tmpDir = ''
                 def __init__(self, title=''):
                     self.year = None
                     self.month = None
                     self.dom = None
                     self.pattern = None
                     self.maxphi = None
                     self.heights = None
                     self.filename = None
                     self.showType = None
                     self.path = None
                     self.objects = None
                     self.nptsx = 101
                     self.nptsy = 101
                     self.fftopt = 0
                     self.site = 1
                     self.dcosx = 1
                     self.dcosy = 1
                     self.dcosxrange = None
                     self.dcosyrange = None
                     self.maxha_min= 0.
                     self.show_object = None
                     self.dcosx_mag = None
                     self.dcosy_mag = None
                     self.ha_mag = None
                     self.time_mag = None
                     self.main_dec = None
                     self.ObjC = None
                     self.ptitle = title
                     self.path4plotname = None
                     self.plotname0 = None
                     self.plotname1 = None
                     self.plotname2 = None
                     self.scriptHeaders = 0
                     self.glat = -11.95
                     self.glon = -76.8667
                     self.UT = 5 #timezone
                     self.glat = -11.951481
                     self.glon = -76.874383
             #       self.outputHead('Show Plot')
             #        self.printBody()
                 def setScriptState(self):
                     self.madForm = cgi.FieldStorage()
                     if self.madForm.has_key('serverdocspath'):
                         self.__serverdocspath = self.madForm.getvalue('serverdocspath')#'/usr/local/www/htdocs'
                     if self.madForm.has_key('tmpdir'):
                         self.__tmpDir = self.madForm.getvalue('tmpdir')#'overJro/tempReports'
                     if self.madForm.has_key('showType'):
                         self.showType = int(self.madForm.getvalue('showType'))
                     if self.showType == 0 or self.showType == 1:
             #             if self.madForm.has_key('year') and \
             #                 self.madForm.has_key('month') and \
             #                 self.madForm.has_key('dom') and \
             #                 self.madForm.has_key('pattern') and \
             #                 self.madForm.has_key('maxphi') and \
             #                 self.madForm.has_key('objects') and \
             #                 self.madForm.has_key('heights'):
                         if self.madForm.has_key('year') and \
                             self.madForm.has_key('month') and \
                             self.madForm.has_key('dom') and \
                             self.madForm.has_key('maxphi') and \
                             self.madForm.has_key('objects') and \
                             self.madForm.has_key('heights'):
                                 self.year = int(self.madForm.getvalue('year'))
                                 self.month = int(self.madForm.getvalue('month'))
                                 self.dom = int(self.madForm.getvalue('dom'))
                                 self.maxphi = float(self.madForm.getvalue('maxphi'))
                                 if self.madForm.has_key('pattern'):
                                     tmp_pattern = self.madForm.getvalue('pattern') #pattern es predifinido en listado o definido por el usuario
                                     self.pattern=[]
                                     if tmp_pattern[0] == '[':
                                         tmp_pattern=tmp_pattern[1:]
                                     if tmp_pattern[-1] == ']':
                                         tmp_pattern=tmp_pattern[0:len(tmp_pattern)-1]
                                     for s in tmp_pattern.split(','):
                                         self.pattern.append(float(s))
                                 elif self.madForm.has_key('filename'):
                                     if self.madForm.has_key('filename'):
                                         self.filename = self.madForm.getvalue('filename') # nombre de archivo: patron de radiacion definido por el usuario
                                     if self.madForm.has_key('path'):
                                         self.path = self.madForm.getvalue('path') #path donde se encuentra el archivo: patron de radiacion del usuario
                                 else:
                                     print ("Content-Type: text/html\n")
                                     print ('<h3> This cgi plot script was called without the proper arguments.</h3>')
                                     print ('<p> This is a script used to plot Antenna Cuts over Jicamarca Antenna</p>')
                                     print ('<p> Required arguments:</p>')
                                     print ('<p>    pattern - chekbox indicating objects over jicamarca antenna</p>')
                                     print ('<p> or')
                                     print ('<p>    filename - The pattern defined by users is a file text')
                                     print ('<p>    path - folder with pattern files')
                                     sys.exit(0)
                                 tmp_heights = self.madForm.getvalue('heights')
                                 self.heights=[]
                                 if tmp_heights[0] == '[':
                                     tmp_heights=tmp_heights[1:]
                                 if tmp_heights[-1] == ']':
                                     tmp_heights=tmp_heights[0:len(tmp_heights)-1]
                                 for s in tmp_heights.split(','):
                                     self.heights.append(float(s))
                                 self.heights = numpy.array(self.heights)
                                 tmp_objects = self.madForm.getvalue('objects') #lista con los objetos a graficar en el patron de radiacion
                                 self.objects=[]
                                 if tmp_objects[0] == '[':
                                     tmp_objects=tmp_objects[1:]
                                 if tmp_objects[-1] == ']':
                                     tmp_objects=tmp_objects[0:len(tmp_objects)-1]
                                 for s in tmp_objects.split(','):
                                     self.objects.append(int(s))
                                 if self.showType == 1:
                                     if numpy.sum(self.objects) == 0:
                                         if  self.scriptHeaders == 0:
                                             print ("Content-Type: text/html\n")
                                         print ('<h3> This cgi plot script was called without the proper arguments.</h3>')
                                         print ('<p> This is a script used to plot Antenna Cuts over Jicamarca Antenna</p>')
                                         print ('<p> Required arguments:</p>')
                                         print ('<p>    objects - chekbox indicating objects over jicamarca antenna</p>')
                                         print ('<p> Please, options in "Select Object" must be checked')
                                         sys.exit(0)
                                 #considerar para futura implementacion
                                 if self.madForm.has_key('filename'):
                                     self.filename = self.madForm.getvalue('filename') # nombre de archivo: patron de radiacion definido por el usuario
                                 if self.madForm.has_key('path'):
                                     self.path = self.madForm.getvalue('path') #path donde se encuentra el archivo: patron de radiacion del usuario
                         else:
                             if  self.scriptHeaders == 0:
                                 print ("Content-Type: text/html\n")
                             print ('<h3> This cgi plot script was called without the proper arguments.</h3>')
                             print ('<p> This is a script used to plot Pattern Field and Celestial Objects over Jicamarca Antenna</p>')
                             print ('<p> Required arguments:</p>')
                             print ('<p>    year - year of event</p>')
                             print ('<p>    month - month of event</p>')
                             print ('<p>    dom - day of month</p>')
                             print ('<p>    pattern - pattern is defined by "Select an Experiment" list box</p>')
                             print ('<p>    maxphi - maxphi is defined by "Max Angle" text box</p>')
                             print ('<p>    objects - objects is a list defined by checkbox in "Select Object"</p>')
                             print ('<p>    heights - heights is defined by "Heights" text box, for default heights=[100,500,1000]</p>')
                             print ('<p>    showType - showType is a hidden element for show plot of Pattern&Object or Antenna Cuts or Sky Noise</p>')
                             sys.exit(0)
                     if self.showType == 2:
                         if self.madForm.has_key('year') and \
                             self.madForm.has_key('month') and \
                             self.madForm.has_key('dom'):
                             self.year = int(self.madForm.getvalue('year'))
                             self.month = int(self.madForm.getvalue('month'))
                             self.dom = int(self.madForm.getvalue('dom'))
                         else:
                             if  self.scriptHeaders == 0:
                                 print ("Content-Type: text/html\n")
                             print ('<h3> This cgi plot script was called without the proper arguments.</h3>')
                             print ('<p> This is a script used to plot Sky Noise over Jicamarca Antenna</p>')
                             print ('<p> Required arguments:</p>')
                             print ('<p>    year - year of event</p>')
                             print ('<p>    month - month of event</p>')
                             print ('<p>    dom - day of month</p>')
                             sys.exit(0)
                 def initParameters1(self):
                     gui=1
                     if self.pattern==None:
                         if gui==1: self.filename = self.filename.split(',')
                     pattern = numpy.atleast_1d(self.pattern)
                     filename = numpy.atleast_1d(self.filename)
                     npatterns = numpy.max(numpy.array([pattern.size,filename.size]))
                     self.pattern = numpy.resize(pattern,npatterns)
                     self.filename = numpy.resize(filename,npatterns)
                     self.doy = datetime.datetime(self.year,self.month,self.dom).timetuple().tm_yday
                     if self.objects==None:
                         self.objects=numpy.zeros(5)
                     else:
                         tmp = numpy.atleast_1d(self.objects)
                         self.objects = numpy.zeros(5)
                         self.objects[0:tmp.size] = tmp
                     self.show_object = self.objects
                     self.maxha_min = 4*self.maxphi*numpy.sqrt(2)*1.25
                     if self.heights==None:
                         self.heights = numpy.array([100.,500.,1000.])
                     #ROJ geographic coordinates and time zone
                     self.glat = -11.95
                     self.glon = -76.8667
                     self.UT = 5 #timezone
                     self.glat = -11.951481
                     self.glon = -76.874383
                     self.junkjd = Time(self.year,self.month,self.dom).change2julday()
                     self.junklst = Julian(self.junkjd).change2lst(longitude=self.glon)
                     # Finding RA of observatory for a specific date
                     self.ra_obs = self.junklst*CoFactors.h2d
                 def initParameters(self):
                     # Defining plot filenames
                     self.path4plotname = os.path.join(self.__serverdocspath,self.__tmpDir)
                     self.plotname0 = 'over_jro_0_%i.png'% (time.time()) #plot pattern & objects
                     self.plotname1 = 'over_jro_1_%i.png'% (time.time()) #plot antenna cuts
                     self.plotname2 = 'over_jro_2_%i.png'% (time.time()) #plot sky noise
                     # Defining antenna axes respect to geographic coordinates (See Ochs report).
             #        alfa = 1.46*Misc_Routines.CoFactors.d2r
             #        theta = 51.01*Misc_Routines.CoFactors.d2r
                     alfa = 1.488312*CoFactors.d2r
                     th = 6.166710 + 45.0
                     theta = th*CoFactors.d2r
                     sina = numpy.sin(alfa)
                     cosa = numpy.cos(alfa)
                     MT1 = numpy.array([[1,0,0],[0,cosa,-sina],[0,sina,cosa]])
                     sinb = numpy.sin(theta)
                     cosb = numpy.cos(theta)
                     MT2 = numpy.array([[cosb,sinb,0],[-sinb,cosb,0],[0,0,1]])
                     self.MT3 = numpy.array(numpy.dot(MT2, MT1)).transpose()
                     self.xg = numpy.dot(self.MT3.transpose(),numpy.array([1,0,0]))
                     self.yg = numpy.dot(self.MT3.transpose(),numpy.array([0,1,0]))
                     self.zg = numpy.dot(self.MT3.transpose(),numpy.array([0,0,1]))
                 def plotPattern2(self, date, phases, gain_tx, gain_rx, ues, just_rx):
                     # Plotting Antenna patterns.
                     self.initParameters()
                     self.doy = datetime.datetime(date.year,date.month,date.day).timetuple().tm_yday
                     self.junkjd = Time(self.year,self.month,self.dom).change2julday()
                     self.junklst = Julian(self.junkjd).change2lst(longitude=self.glon)
                     self.ra_obs = self.junklst*CoFactors.h2d
                     date = Time(date.year,date.month,date.day).change2strdate(mode=2)
                     mesg = 'Over Jicamarca: ' + date[0]
                     ObjAnt = JroPattern(pattern=0,
                                         filename=None,
                                         path=None,
                                         nptsx=self.nptsx,
                                         nptsy=self.nptsy,
                                         #maxphi=self.maxphi,
                                         fftopt=self.fftopt,
                                         phases=numpy.array(phases),
                                         gain_tx=numpy.array(gain_tx),
                                         gain_rx=numpy.array(gain_rx),
                                         ues=numpy.array(ues),
                                         just_rx=just_rx
                                         )
                     dum = AntPatternPlot()
                     dum.contPattern(iplot=0,
                                     gpath=self.path4plotname,
                                     filename=self.plotname0,
                                     mesg=mesg,
                                     amp=ObjAnt.norpattern,
                                     x=ObjAnt.dcosx,
                                     y=ObjAnt.dcosy,
                                     getCut=ObjAnt.getcut,
                                     title=self.ptitle,
                                     save=False)
                     dum.plotRaDec(gpath=self.path4plotname,
                                   filename=self.plotname0,
                                   jd=self.junkjd,
                                   ra_obs=self.ra_obs,
                                   xg=self.xg,
                                   yg=self.yg,
                                   x=ObjAnt.dcosx,
                                   y=ObjAnt.dcosy,
                                   save=False)
                     ObjB = BField(self.year,self.doy,1,self.heights)
                     [dcos, alpha, nlon, nlat] = ObjB.getBField()
                     dum.plotBField('', '',dcos,alpha,nlon,nlat,
                                    self.dcosxrange,
                                    self.dcosyrange,
                                    ObjB.heights,
                                    ObjB.alpha_i,
                                    save=False)
                     return dum.fig
                 def plotPattern(self):
                     # Plotting Antenna patterns.
                     npatterns = numpy.size(self.pattern)
                     if npatterns==1:
                         if self.pattern[0] == None: npatterns = self.filename.__len__()
                     date = TimeTools.Time(self.year,self.month,self.dom).change2strdate(mode=2)
                     mesg = 'Over Jicamarca: ' + date[0]
                     title = ''
                     for ii in numpy.arange(npatterns):
                         ObjAnt = JroPattern(pattern=self.pattern[ii],
                                             filename=self.filename[ii],
                                             path=self.path,
                                             nptsx=self.nptsx,
                                             nptsy=self.nptsy,
                                             maxphi=self.maxphi,
                                             fftopt=self.fftopt)
                         title += ObjAnt.title
                         # Plotting Contour Map
                         self.path4plotname = '/home/jespinoza/workspace/radarsys/trunk/webapp/apps/abs/static/images'
                         dum = AntPatternPlot()
                         dum.contPattern(iplot=ii,
                                                                     gpath=self.path4plotname,
                                                                     filename=self.plotname0,
                                                                     mesg=mesg,
                                                                     amp=ObjAnt.norpattern,
                                                                     x=ObjAnt.dcosx,
                                                                     y=ObjAnt.dcosy,
                                                                     getCut=ObjAnt.getcut,
                                                                     title=title)
             #                                                         title=ObjAnt.title)
             #             self.ptitle = ObjAnt.title
                         if ii != (npatterns-1):
                             title += '+'
                         vect_ant = numpy.array([ObjAnt.meanpos[0],ObjAnt.meanpos[1],numpy.sqrt(1-numpy.sum(ObjAnt.meanpos**2.))])
                         vect_geo = numpy.dot(scipy.linalg.inv(self.MT3),vect_ant)
                         vect_polar = Misc_Routines.Vector(numpy.array(vect_geo),direction=1).Polar2Rect()
-                        [ra,dec,ha] = Astro_Coords.AltAz(vect_polar[1],vect_polar[0],self.junkjd).change2equatorial()
+                        [ra,dec,ha] = AltAz(vect_polar[1],vect_polar[0],self.junkjd).change2equatorial()
                         print('Main beam position (HA(min), DEC(degrees)): %f %f')%(ha*4.,dec)
                         self.main_dec = dec
                     self.ptitle = title
                     AntPatternPlot().plotRaDec(gpath=self.path4plotname,
                                                                 filename=self.plotname0,
                                                                 jd=self.junkjd,
                                                                 ra_obs=self.ra_obs,
                                                                 xg=self.xg,
                                                                 yg=self.yg,
                                                                 x=ObjAnt.dcosx,
                                                                 y=ObjAnt.dcosy)
                     self.dcosx = ObjAnt.dcosx
                     self.dcosy = ObjAnt.dcosy
                     self.dcosxrange = [numpy.min(self.dcosx),numpy.max(self.dcosx)]
                     self.dcosyrange = [numpy.min(self.dcosy),numpy.max(self.dcosy)]
                 def plotBfield(self):
                     if self.show_object[0]>0:
                         # Getting B field
                         ObjB = BField(self.year,self.doy,self.site,self.heights)
                         [dcos, alpha, nlon, nlat] = ObjB.getBField()
                         # Plotting B field.
             #            print "Drawing magnetic field over Observatory"
                         Obj = BFieldPlot()
                         Obj.plotBField(self.path4plotname,self.plotname0,dcos,alpha,nlon,nlat,self.dcosxrange,self.dcosyrange,ObjB.heights,ObjB.alpha_i)
                         if self.show_object[0]>1:
                             Bhei = 0
                             dcosx = Obj.alpha_location[:,0,Bhei]
                             dcosy = Obj.alpha_location[:,1,Bhei]
                             vect_ant = [dcosx,dcosy,numpy.sqrt(1.-(dcosx**2. + dcosy**2.))]
                             vect_ant = numpy.array(vect_ant)
                             vect_geo = numpy.dot(scipy.linalg.inv(self.MT3),vect_ant)
                             vect_geo = numpy.array(vect_geo).transpose()
                             vect_polar = Misc_Routines.Vector(vect_geo,direction=1).Polar2Rect()
-                            [ra,dec,ha] = Astro_Coords.AltAz(vect_polar[1,:],vect_polar[0,:],self.junkjd).change2equatorial()
+                            [ra,dec,ha] = AltAz(vect_polar[1,:],vect_polar[0,:],self.junkjd).change2equatorial()
                             val = numpy.where(ha>=180)
                             if val[0].size>0:ha[val] = ha[val] -360.
                             val = numpy.where(numpy.abs(ha)<=self.maxphi)
                             if val[0].size>2:
                                 self.dcosx_mag = dcosx[val]
                                 self.dcosy_mag = dcosy[val]
                                 self.ha_mag = ha[val]
                                 self.time_mag = 0
                 def plotCelestial(self):
                     ntod = 24.*16.
                     tod = numpy.arange(ntod)/ntod*24.
                     [month,dom] = Doy2Date(self.year,self.doy).change2date()
                     jd = Time(self.year,month,dom,tod+self.UT).change2julday()
                     if numpy.sum(self.show_object[1:]>0)!=0:
                         self.ObjC = CelestialObjectsPlot(jd,self.main_dec,tod,self.maxha_min,self.show_object)
                         self.ObjC.drawObject(self.glat,
                                         self.glon,
                                         self.xg,
                                         self.yg,
                                         self.dcosxrange,
                                         self.dcosyrange,
                                         self.path4plotname,
                                         self.plotname0)
                 def plotAntennaCuts(self):
             #        print "Drawing antenna cuts"
                     incha = 0.05 # min
                     nha = numpy.int32(2*self.maxha_min/incha) + 1.
                     newha = numpy.arange(nha)/nha*2.*self.maxha_min - self.maxha_min
                     nha_star = numpy.int32(200./incha)
                     star_ha = (numpy.arange(nha_star) - (nha_star/2))*nha_star
                     #Init ObjCut for PatternCutPlot()
                     view_objects = numpy.where(self.show_object>0)
                     subplots = len(view_objects[0])
                     ObjCut = PatternCutPlot(subplots)
                     for io in (numpy.arange(5)):
                         if self.show_object[io]==2:
                             if io==0:
                                 if self.dcosx_mag.size!=0:
                                     dcosx = self.dcosx_mag
                                     dcosy = self.dcosy_mag
                                     dcosz = 1 - numpy.sqrt(dcosx**2. + dcosy**2.)
                                     # Finding rotation of B respec to antenna coords.
                                     [mm,bb] = scipy.polyfit(dcosx,dcosy,1)
                                     alfa = 0.0
                                     theta = -1.*numpy.arctan(mm)
                                     sina = numpy.sin(alfa); cosa = numpy.cos(alfa)
                                     MT1 = [[1,0,0],[0,cosa,-sina],[0,sina,cosa]]
                                     MT1 = numpy.array(MT1)
                                     sinb = numpy.sin(theta); cosb = numpy.cos(theta)
                                     MT2 = [[cosb,sinb,0],[-sinb,cosb,0],[0,0,1]]
                                     MT2 = numpy.array(MT2)
                                     MT3_mag = numpy.dot(MT2, MT1)
                                     MT3_mag = numpy.array(MT3_mag).transpose()
                                     # Getting dcos respec to B coords
                                     vector = numpy.array([dcosx,dcosy,dcosz])
                                     nvector = numpy.dot(MT3_mag,vector)
                                     nvector = numpy.array(nvector).transpose()
             ##                        print 'Rotation (deg) %f'%(theta/Misc_Routines.CoFactors.d2r)
                                     yoffset = numpy.sum(nvector[:,1])/nvector[:,1].size
             #                        print 'Dcosyoffset %f'%(yoffset)
                                     ha = self.ha_mag*4.
                                     time = self.time_mag
                                     width_star = 0.1 # half width in minutes
                                     otitle = 'B Perp. cut'
             #                    else:
             #                        print "No B perp. over Observatory"
             #
             #
                             elif io==1:
                                 if self.ObjC.dcosx_sun.size!=0:
                                     dcosx = self.ObjC.dcosx_sun
                                     dcosy = self.ObjC.dcosy_sun
                                     ha = self.ObjC.ha_sun*4.0
                                     time = self.ObjC.time_sun
                                     width_star = 2. # half width in minutes
                                     otitle = 'Sun cut'
             #                    else:
             #                        print "Sun is not passing over Observatory"
                             elif io==2:
                                 if self.ObjC.dcosx_moon.size!=0:
                                     dcosx = self.ObjC.dcosx_moon
                                     dcosy = self.ObjC.dcosy_moon
                                     ha = self.ObjC.ha_moon*4
                                     time = self.ObjC.time_moon
                                     m_distance = 404114.6  # distance to the Earth in km
                                     m_diameter = 1734.4    # diameter in km.
                                     width_star = numpy.arctan(m_distance/m_diameter)
                                     width_star = width_star/2./CoFactors.d2r*4.
                                     otitle = 'Moon cut'
             #                    else:
             #                        print "Moon is not passing over Observatory"
                             elif io==3:
                                 if self.ObjC.dcosx_hydra.size!=0:
                                     dcosx = self.ObjC.dcosx_hydra
                                     dcosy = self.ObjC.dcosy_hydra
                                     ha = self.ObjC.ha_hydra*4.
                                     time = self.ObjC.time_hydra
                                     width_star = 0.25 # half width in minutes
                                     otitle = 'Hydra cut'
             #                    else:
             #                        print "Hydra is not passing over Observatory"
                             elif io==4:
                                 if self.ObjC.dcosx_galaxy.size!=0:
                                     dcosx = self.ObjC.dcosx_galaxy
                                     dcosy = self.ObjC.dcosy_galaxy
                                     ha = self.ObjC.ha_galaxy*4.
                                     time = self.ObjC.time_galaxy
                                     width_star = 25. # half width in minutes
                                     otitle = 'Galaxy cut'
             #                    else:
             #                        print "Galaxy center is not passing over Jicamarca"
             #
             #
                             hour = numpy.int32(time)
                             mins = numpy.int32((time - hour)*60.)
                             secs = numpy.int32(((time - hour)*60. - mins)*60.)
                             ObjT = Time(self.year,self.month,self.dom,hour,mins,secs)
                             subtitle = ObjT.change2strdate()
                             star_cut = numpy.exp(-(star_ha/width_star)**2./2.)
                             pol = scipy.polyfit(ha,dcosx,3.)
                             polx = numpy.poly1d(pol); newdcosx = polx(newha)
                             pol = scipy.polyfit(ha,dcosy,3.)
                             poly = numpy.poly1d(pol);newdcosy = poly(newha)
                             patterns = []
                             for icut in numpy.arange(self.pattern.size):
                                 # Getting Antenna cut.
                                 Obj = JroPattern(dcosx=newdcosx,
                                                     dcosy=newdcosy,
                                                     getcut=1,
                                                     pattern=self.pattern[icut],
                                                     path=self.path,
                                                     filename=self.filename[icut])
                                 Obj.getPattern()
                                 patterns.append(Obj.pattern)
                             ObjCut.drawCut(io,
                                            patterns,
                                            self.pattern.size,
                                            newha,
                                            otitle,
                                            subtitle,
                                            self.ptitle)
                     ObjCut.saveFig(self.path4plotname,self.plotname1)
                 def plotSkyNoise(self):
             #       print 'Creating SkyNoise map over Jicamarca'
                     dom = self.dom
                     month = self.month
                     year = self.year
                     julian = Time(year,month,dom).change2julday()
-                    [powr,time, lst] = Astro_Coords.CelestialBodies().skyNoise(julian)
+                    [powr,time, lst] = CelestialBodies().skyNoise(julian)
                     SkyNoisePlot([year,month,dom],powr,time,lst).getPlot(self.path4plotname,self.plotname2)
                 def outputHead(self,title):
                     print ("Content-Type: text/html")
                     print (self).scriptHeaders = 1
                     print ('<html>')
                     print ('<head>')
                     print ('\t<title>' + title + '</title>')
                     print ('<style type="text/css">')
                     print ('body')
                     print ('{')
                     print ('background-color:#ffffff;')
                     print ('}')
                     print ('h1')
                     print ('{')
                     print ('color:black;')
                     print ('font-size:18px;')
                     print ('text-align:center;')
                     print ('}')
                     print ('p')
                     print ('{')
                     print ('font-family:"Arial";')
                     print ('font-size:16px;')
                     print ('color:black;')
                     print ('}')
                     print ('</style>')
             #        self.printJavaScript()
                     print ('</head>')
                 def printJavaScript(self):
                     print
                 def printBody(self):
                     print ('<body>')
             #        print '<h1>Test Input Parms</h1>'
             #        for key in self.madForm.keys():
             #                #print '<p> name=' + str(key)
             #                if type(self.madForm.getvalue(key)) == types.ListType:
             #                    for value in self.madForm.getvalue(key):
             #                        print '<p> name=' + str(key) + \
             #                              ' value=' + value + ''
             #                else:
             #                    print '<p> name=' + str(key) + \
             #                              ' value=' + str(cgi.escape(self.madForm.getvalue(key))) + ''
                     print ('<form name="form1" method="post" target="showFrame">')
                     print ('    <div align="center">')
                     print ('        <table width=98% border="1" cellpadding="1">')
                     print ('            <tr>')
                     print ('                <td colspan="2" align="center">')
                     if self.showType == 0:
                         print ('                    <IMG SRC="%s" BORDER="0" >')%(os.path.join(os.sep + self.__tmpDir,self.plotname0))
                     if self.showType == 1:
                         print ('                    <IMG SRC="%s" BORDER="0" >')%(os.path.join(os.sep + self.__tmpDir,self.plotname1))
                     if self.showType == 2:
                         print ('                    <IMG SRC="%s" BORDER="0" >')%(os.path.join(os.sep + self.__tmpDir,self.plotname2))
                     print ('                </td>')
                     print ('            </tr>')
                     print ('        </table>')
                     print ('    </div>')
                     print ('</form>')
                     print ('</body>')
                     print ('</html>')
                 #def execute(self, serverdocspath, tmpdir, currentdate, finalpath, showType=0, maxphi=5.0, objects="[1,1]", heights="[150,500,1000]"):
                 def setInputParameters(self, serverpath, currentdate, finalpath, showType=0, maxphi=5.0, objects="[1,1]", heights="[150,500,1000]"):
                     self.objects=[]
                     self.heights=[]
                     #self.__serverdocspath = serverdocspath
                     self.__serverdocspath = os.path.split(serverpath)[0]
                     #self.__tmpDir = tmpdir
                     self.__tmpDir = os.path.split(serverpath)[1]
                     self.showType = int(showType)
                     self.year = int(currentdate.strftime("%Y")) # Get year of currentdate
                     self.month = int(currentdate.strftime("%m")) # Get month of currentdate
                     self.dom = int(currentdate.strftime("%d")) # Get day of currentdate
                     self.filename = os.path.split(finalpath)[1]
                     self.path = os.path.split(finalpath)[0]
                     self.maxphi = float(maxphi)
                     tmp_objects = (objects.replace("[","")).replace("]","")
                     for s in tmp_objects.split(','):
                         self.objects.append(int(s))
                     tmp_heights = (heights.replace("[","")).replace("]","")
                     for s in tmp_heights.split(','):
                         self.heights.append(float(s))
                     self.heights = numpy.array(self.heights)
                 def setupParameters(self):
                     self.initParameters()
                 def initParametersCGI(self):
                     self.setScriptState()
                     self.initParameters()
                 def execute(self):
                     if self.showType == 0 or self.showType == 1:
                         self.initParameters1()
                         self.plotPattern()
                         if numpy.sum(self.show_object>0) != 0:
                             self.plotBfield()
                             self.plotCelestial()
                         if numpy.sum(self.show_object>1) != 0:
                             self.plotAntennaCuts()
                     if self.showType == 2:
                         self.plotSkyNoise()
                 def getPlot(self):
                     return os.path.join(self.__serverdocspath,self.__tmpDir,self.plotname0)
             if __name__ == '__main__':
                 # Script overJroShow.py
                 # This script only calls the init function of the class overJroShow()
                 # All work is done by the init function
                 phases = numpy.array([[2.0,0.0,1.5,1.5,1.0,1.0,1.0,0.5],
                         [2.0,2.5,2.5,3.5,0.5,1.0,1.0,1.0],
                         [2.5,2.5,1.0,1.0,0.5,0.5,0.5,0.5],
                         [1.0,1.0,1.0,1.0,0.5,0.5,0.5,1.0],
                         [0.5,0.5,0.5,0.5,0.5,0.0,0.0,0.0],
                         [0.5,0.5,1.0,0.5,0.0,0.0,0.0,0.0],
                         [0.5,0.5,0.5,1.0,0.0,0.0,0.0,0.0],
                         [0.5,0.5,0.5,0.5,0.0,0.0,0.0,0.0]])
                 gain_tx = numpy.array([[0,0,0,0,0,0,0,0],
                         [0,0,0,0,0,0,0,0],
                         [0,0,0,0,0,0,0,0],
                         [0,0,0,0,0,0,0,0],
                         [0,0,0,0,1,1,1,1],
                         [0,0,0,0,0,0,0,0],
                         [0,0,0,0,0,0,0,0],
                         [0,0,0,0,0,0,0,0]])
                 gain_rx = numpy.array([[0,0,0,0,0,0,0,0],
                         [0,0,1,0,0,0,0,0],
                         [0,0,1,0,0,0,0,0],
                         [0,0,0,0,0,0,0,0],
                         [0,0,0,0,0,0,0,0],
                         [0,0,0,0,0,0,0,0],
                         [0,0,0,0,0,0,0,0],
                         [0,0,0,0,0,0,0,0]])
                 jro = overJroShow()
                 fig = jro.plotPattern2(datetime.datetime.today(),
                                        phases=phases,
                                        gain_tx=gain_tx,
                                        gain_rx=gain_rx,
                                        ues=numpy.array([0.0,0.0,0.0,0.0]),
                                        just_rx=0)
                 fig.savefig('./pat.png')

apps/jars/models.py +10 -3

             import json
             import requests
             from django.db import models
             from django.core.validators import MinValueValidator, MaxValueValidator
             from django.urls import reverse
             from apps.main.models import Configuration
             from apps.main.utils import Params
             from .utils import create_jarsfiles
             # Create your models here.
             EXPERIMENT_TYPE = (
                 (0, 'RAW_DATA'),
                 (1, 'PDATA'),
             )
             DATA_TYPE = (
                 (0, 'SHORT'),
                 (1, 'FLOAT'),
             )
             DECODE_TYPE = (
                 (0, 'None'),
                 (1, 'TimeDomain'),
                 (2, 'FreqDomain'),
                 (3, 'InvFreqDomain'),
             )
             FILTER = '{"id":1, "clock": 60, "multiplier": 5, "frequency": 49.92, "f_decimal":	721554506, "fir": 2, "cic_2": 12, "cic_5": 25}'
             class JARSFilter(models.Model):
                 JARS_NBITS  = 32
                 name        = models.CharField(verbose_name='Name', max_length=60, unique=True, default='')
                 clock       = models.FloatField(verbose_name='Clock In (MHz)', validators=[
                                           MinValueValidator(5), MaxValueValidator(75)], null=True, default=60)
                 multiplier  = models.PositiveIntegerField(verbose_name='Multiplier', validators=[
                                                          MinValueValidator(1), MaxValueValidator(20)], default=5)
                 frequency   = models.FloatField(verbose_name='Frequency (MHz)', validators=[
                                               MaxValueValidator(150)], null=True, default=49.9200)
                 f_decimal   = models.BigIntegerField(verbose_name='Frequency (Decimal)', validators=[
                                                    MinValueValidator(-9223372036854775808), MaxValueValidator(2**JARS_NBITS-1)], null=True, default=721554505)
                 cic_2       = models.PositiveIntegerField(verbose_name='CIC2', validators=[
                                                     MinValueValidator(2), MaxValueValidator(100)], default=10)
                 scale_cic_2 = models.PositiveIntegerField(verbose_name='Scale CIC2', validators=[
                                                           MinValueValidator(0), MaxValueValidator(6)], default=1)
                 cic_5       = models.PositiveIntegerField(verbose_name='CIC5', validators=[
                                                     MinValueValidator(1), MaxValueValidator(100)], default=1)
                 scale_cic_5 = models.PositiveIntegerField(verbose_name='Scale CIC5', validators=[
                                                           MinValueValidator(0), MaxValueValidator(20)], default=5)
                 fir         = models.PositiveIntegerField(verbose_name='FIR', validators=[
                                                   MinValueValidator(1), MaxValueValidator(100)], default=6)
                 scale_fir   = models.PositiveIntegerField(verbose_name='Scale FIR', validators=[
                                                         MinValueValidator(0), MaxValueValidator(7)], default=3)
                 number_taps = models.PositiveIntegerField(verbose_name='Number of taps', validators=[
                                                           MinValueValidator(1), MaxValueValidator(256)], default=4)
                 taps        = models.CharField(verbose_name='Taps', max_length=1600, default='0')
                 class Meta:
                     db_table = 'jars_filters'
                 def __unicode__(self):
                     return u'%s' % (self.name)
                 def jsonify(self):
                     data    = {}
                     ignored = ()
                     for field in self._meta.fields:
                         if field.name in ignored:
                             continue
                         data[field.name] = field.value_from_object(self)
                     return data
                 def parms_to_dict(self):
                     parameters = {}
                     parameters['name']       = self.name
                     parameters['clock']      = float(self.clock)
                     parameters['multiplier'] = int(self.multiplier)
                     parameters['frequency']  = float(self.frequency)
                     parameters['f_decimal']  = int(self.frequency)
                     parameters['fir']        = int(self.fir)
                     parameters['cic_2']      = int(self.cic_2)
                     parameters['cic_5']      = int(self.cic_5)
                     return parameters
                 def dict_to_parms(self, parameters):
                     self.name       = parameters['name']
                     self.clock      = parameters['clock']
                     self.multiplier = parameters['multiplier']
                     self.frequency  = parameters['frequency']
                     self.f_decimal  = parameters['f_decimal']
                     self.fir        = parameters['fir']
                     self.cic_2      = parameters['cic_2']
                     self.cic_5      = parameters['cic_5']
                 def dict_to_parms_new(self, parameters):
                     self.name        = parameters['name']
                     self.clock       = parameters['clock']
                     self.multiplier  = parameters['multiplier']
                     self.frequency   = parameters['frequency']
                     self.f_decimal   = parameters['f_decimal']
                     self.fir         = parameters['fir']
                     self.cic_2       = parameters['cic_2']
                     self.cic_5       = parameters['cic_5']
                     self.scale_fir   = parameters['scale_fir']
                     self.scale_cic_2 = parameters['scale_cic_2']
                     self.scale_cic_5 = parameters['scale_cic_5']
                     self.number_taps = parameters['number_taps']
                     self.taps        = parameters['taps']
             class JARSConfiguration(Configuration):
                 ADC_RESOLUTION   = 8
                 PCI_DIO_BUSWIDTH = 32
                 HEADER_VERSION   = 1103
                 BEGIN_ON_START   = True
                 REFRESH_RATE     = 1
                 exp_type         = models.PositiveIntegerField(
                     verbose_name='Experiment Type', choices=EXPERIMENT_TYPE, default=0)
                 cards_number     = models.PositiveIntegerField(verbose_name='Number of Cards', validators=[
                                                            MinValueValidator(1), MaxValueValidator(4)], default=1)
                 channels_number  = models.PositiveIntegerField(verbose_name='Number of Channels', validators=[
                                                               MinValueValidator(1), MaxValueValidator(8)], default=5)
                 channels         = models.CharField(
                     verbose_name='Channels', max_length=15, default='1,2,3,4,5')
                 data_type        = models.PositiveIntegerField(
                     verbose_name='Data Type', choices=DATA_TYPE, default=0)
                 raw_data_blocks  = models.PositiveIntegerField(
                     verbose_name='Raw Data Blocks', validators=[MaxValueValidator(5000)], default=60)
                 profiles_block   = models.PositiveIntegerField(
                     verbose_name='Profiles Per Block', default=400)
                 acq_profiles     = models.PositiveIntegerField(
                     verbose_name='Acquired Profiles', default=400)
                 ftp_interval     = models.PositiveIntegerField(
                     verbose_name='FTP Interval', default=60)
                 fftpoints        = models.PositiveIntegerField(
                     verbose_name='FFT Points', default=16)
                 cohe_integr_str  = models.PositiveIntegerField(
                     verbose_name='Coh. Int. Stride', validators=[MinValueValidator(1)], default=30)
                 cohe_integr      = models.PositiveIntegerField(
                     verbose_name='Coherent Integrations', validators=[MinValueValidator(1)], default=30)
                 incohe_integr    = models.PositiveIntegerField(
                     verbose_name='Incoherent Integrations', validators=[MinValueValidator(1)], default=30)
                 decode_data      = models.PositiveIntegerField(
                     verbose_name='Decode Data', choices=DECODE_TYPE, default=0)
                 post_coh_int     = models.BooleanField(
                     verbose_name='Post Coherent Integration', default=False)
                 spectral_number  = models.PositiveIntegerField(
                     verbose_name='# Spectral Combinations', validators=[MinValueValidator(1)], default=1)
                 spectral         = models.CharField(
                     verbose_name='Combinations', max_length=5000, default='[0, 0],')
                 create_directory = models.BooleanField(
                     verbose_name='Create Directory Per Day', default=True)
                 include_expname  = models.BooleanField(
                     verbose_name='Experiment Name in Directory', default=False)
                 #view_raw_data    = models.BooleanField(verbose_name='View Raw Data', default=True)
                 save_ch_dc       = models.BooleanField(
                     verbose_name='Save Channels DC', default=True)
                 save_data        = models.BooleanField(verbose_name='Save Data', default=True)
                 filter_parms     = models.CharField(
                     max_length=10000, default=FILTER)
                 filter           = models.ForeignKey(
                     'JARSFilter', verbose_name='Filter', null=True, blank=True, on_delete=models.CASCADE)
                 class Meta:
                     db_table = 'jars_configurations'
                 def filter_resolution(self):
                     filter_parms = json.loads(self.filter_parms)
                     clock        = float(filter_parms['clock'])
                     cic_2        = filter_parms['cic_2']
                     cic_5        = filter_parms['cic_5']
                     fir          = filter_parms['fir']
                     resolution   = round((clock/(cic_2*cic_5*fir)), 2)
                     return resolution
                 def dict_to_parms(self, params, id=None):
                     if id is not None:
                         data = Params(params).get_conf(id_conf=id)
                     else:
                         data = Params(params).get_conf(dtype='jars')
                         data['filter_parms'] = params['filter_parms']
                     # self.name = data['name']
                     self.exp_type = data['exp_type']
                     #----PDATA----
                     if self.exp_type == 1:
                         self.incohe_integr   = data['incohe_integr']
                         self.spectral_number = data['spectral_number']
                         self.spectral        = data['spectral']
                         self.fftpoints       = data['fftpoints']
                         self.save_ch_dc      = data['save_ch_dc']
                     else:
                         self.raw_data_blocks = data['raw_data_blocks']
                     #----PDATA----
                     self.cards_number    = data['cards_number']
                     self.channels_number = data['channels_number']
                     self.channels        = data['channels']
                     self.data_type       = data['data_type']
                     self.profiles_block  = data['profiles_block']
                     self.acq_profiles    = data['acq_profiles']
                     self.ftp_interval    = data['ftp_interval']
                     self.cohe_integr_str = data['cohe_integr_str']
                     self.cohe_integr     = data['cohe_integr']
                     #----DECODE----
                     self.decode_data  = data['decode_data']
                     self.post_coh_int = data['post_coh_int']
                     #----DECODE----
                     self.create_directory = data['create_directory']
                     self.include_expname  = data['include_expname']
                     self.save_data        = data['save_data']
                     self.filter_parms     = json.dumps(data['filter_parms'])
                     self.save()
                 def parms_to_text(self, file_format='jars'):
                     data = self.experiment.parms_to_dict()
                     for key in data['configurations']['allIds']:
                         if data['configurations']['byId'][key]['device_type'] in ('dds', 'cgs'):
                             data['configurations']['allIds'].remove(key)
                             data['configurations']['byId'].pop(key)
                         elif data['configurations']['byId'][key]['device_type'] == 'jars':
                             data['configurations']['byId'][key] = self.parms_to_dict(
                             )['configurations']['byId'][str(self.pk)]
                         elif data['configurations']['byId'][key]['device_type'] == 'rc':
                                 data['configurations']['byId'][key]['pulses'] = ''
                                 data['configurations']['byId'][key]['delays'] = ''
                     rc_ids = [pk for pk in data['configurations']['allIds']
                               if data['configurations']['byId'][pk]['device_type'] == 'rc']
                     mix_ids = [pk for pk in rc_ids if data['configurations']
                                ['byId'][pk]['mix']]
                     if mix_ids:
                         params = data['configurations']['byId'][mix_ids[0]]['parameters']
                         rc = data['configurations']['byId'][params.split(
                             '-')[0].split('|')[0]]
                         rc['mix'] = True
                         data['configurations']['byId'][rc['id']] = rc
                     elif len(rc_ids) == 0:
                         self.message = 'File needs RC configuration'
                         return ''
                     json_data = json.dumps(data)
                     racp_file, filter_file = create_jarsfiles(json_data)
                     if file_format == 'racp':
                         return racp_file
                     return filter_file
                 def request(self, cmd, method='get', **kwargs):
                     req = getattr(requests, method)(self.device.url(cmd), **kwargs)
                     payload = req.json()
                     return payload
                 def status_device(self):
                     try:
                         payload = self.request('status',
                                                params={'name': self.experiment.name})
                         self.device.status = payload['status']
                         self.device.save()
                         self.message = payload['message']
                     except Exception as e:
                         self.device.status = 0
                         self.message = str(e)
                         self.device.save()
                         return False
                     return True
                 def stop_device(self):
                     try:
                         payload = self.request('stop', 'post')
                         self.device.status = payload['status']
                         self.device.save()
                         self.message = payload['message']
+                        print("------STOP JARS------",flush=True)
                     except Exception as e:
                         self.device.status = 0
                         self.message = str(e)
                         self.device.save()
                         return False
                     return True
                 def read_device(self):
                     try:
                         payload = self.request(
                             'read', params={'name': self.experiment.name})
                         self.message = 'Configuration loaded'
                     except:
                         self.device.status = 0
                         self.device.save()
                         self.message = 'Could not read JARS configuration.'
                         return False
                     return payload
                 def write_device(self):
                     if self.device.status == 3:
                         self.message = 'Could not configure device. Software Acquisition is running'
                         return False
+                    print("------ JARS WRITING ------",flush=True)
+                    self.device.status = 5  #Busy
+                    print("STATUS:", self.device.status,flush=True)
                     data = self.experiment.parms_to_dict()
                     #print(data)
                     for key in data['configurations']['allIds']:
                         if data['configurations']['byId'][key]['device_type'] in ('dds', 'cgs'):
                             data['configurations']['allIds'].remove(key)
                             data['configurations']['byId'].pop(key)
                         elif data['configurations']['byId'][key]['device_type'] == 'rc':
                                 data['configurations']['byId'][key]['pulses'] = ''
                                 data['configurations']['byId'][key]['delays'] = ''
                     rc_ids = [pk for pk in data['configurations']['allIds']
                               if data['configurations']['byId'][pk]['device_type'] == 'rc']
                     if len(rc_ids) == 0:
                         self.message = 'Missing RC configuration'
                         return False
                     json_data = json.dumps(data)
+                    print("STATUS:", self.device.status,flush=True)
                     try:
                         payload = self.request('write', 'post', json=json_data)
                         self.device.status = payload['status']
                         self.message = payload['message']
                         self.device.save()
                         if self.device.status == 1:
                             return False
+                        print("------ JARS WRITED ------",flush=True)
                     except Exception as e:
                         self.device.status = 0
                         self.message = str(e)
                         self.device.save()
                         return False
+                    print("FINAL STATUS:", self.device.status,flush=True)
                     return True
                 def start_device(self):
                     try:
                         payload = self.request('start', 'post',
                                                json={'name': self.experiment.name})
                         self.device.status = payload['status']
                         self.message = payload['message']
                         self.device.save()
                         if self.device.status == 1:
                             return False
                     except Exception as e:
                         self.device.status = 0
                         self.message = str(e)
                         self.device.save()
                         return False
+                    print("------START JARS------",flush=True)
                     return True
                 def get_log(self):
                     payload = None
                     try:
                         payload = requests.get(self.device.url('get_log'), params={
                                                'name': self.experiment.name})
                     except:
                         self.device.status = 0
                         self.device.save()
                         self.message = 'Jars API is not running.'
                         return False
                     self.message = 'Jars API is running'
                     return payload
                 def update_from_file(self, filename):
                     f = JARSFile(filename)
                     self.dict_to_parms(f.data)
                     self.save()
                 def get_absolute_url_import(self):
                     return reverse('url_import_jars_conf', args=[str(self.id)])
                 def get_absolute_url_read(self):
                     return reverse('url_read_jars_conf', args=[str(self.id)])
                 def get_absolute_url_log(self):
                     return reverse('url_get_jars_log', args=[str(self.id)])

apps/main/management/commands/hello.py +3 -4

             from django.core.management.base import BaseCommand
             from apps.main.models import Campaign, Location
             from datetime import datetime,timedelta
             from apps.main.views import radar_start
             from django.shortcuts import render, redirect,get_object_or_404, HttpResponse
             from django.urls import reverse
             from django.utils.timezone import is_aware
             from django.contrib import messages
             from django.http import HttpResponseRedirect
             from apps.main.views import experiment_start
             from apps.main.models import Experiment, Configuration
             class Command(BaseCommand):
                 """
                 Restart experiment every night at 05:00 am.
                 Example:
                     manage.py restart_experiment
                 """
                 def handle(self, *args, **options):
                     print("\n\n")
                     all_campaigns=Campaign.objects.all()
                     campaigns = Campaign.objects.filter(start_date__lte=datetime.now(),
                                                     end_date__gte=datetime.now()).order_by('-start_date')
                     for campaign in all_campaigns:
                         if campaign.start_date<datetime.now() and campaign.end_date > datetime.now():
                             radar=campaign.get_experiments_by_radar(radar=None)
                             for rad in radar:
                                 # print("RADR", rad)
                                 radar_id=rad["id"]
                                 # print("RADR_",radar_id)
-                                radar_start_scheduler(campaign.id,radar_id)
+                                radar_write_start_scheduler(campaign.id,radar_id)
                             print(campaign.name, "\t\t Campaign already running")
                         else:
                             radar=campaign.get_experiments_by_radar(radar=None)
                             radar_id=radar[0]["id"]
                             if campaign.experiments.all()[0].status !=1:
                                 print(campaign.name, "\t\t Stopping Campaign...")
                                 a=radar_stop_scheduler(campaign.id,radar_id,campaign.experiments.all()[0].id)
                                 print("New Status: ", a)
                             else:
                                 print(campaign.name,"\t\t\t Campaign already stooped")
-            def radar_start_scheduler(id_camp,id_radar):
+            def radar_write_start_scheduler(id_camp,id_radar):
                 campaign    = get_object_or_404(Campaign, pk=id_camp)
                 experiments = campaign.get_experiments_by_radar(id_radar)[0]['experiments']
-                now         = datetime.now()
                 # print(campaign)
                 # print(experiments)
                 for exp in experiments:
                     exp = get_object_or_404(Experiment, pk=exp.id)
                     # print("---------DEBUGG-------------")
                     # print(exp)
                     if exp.status == 2:
                         print('\t\t\t {} \t\t  Experiment already runnnig'.format(exp))
                     else:
                         exp.status = exp.start()
                         if exp.status == 0:
                             print('\t\t\t {} \t\tExperiment not start'.format(exp))
                         if exp.status == 2:
                             print('\t\t\t {} \t\tExperiment started'.format(exp))
                         if exp.status == 4:
-                            print('\t\t\t {} \t\tExperiment with state uknown, please reset'.format(exp))
+                            print('\t\t\t {} \t\tExperiment with state uknown, please reset (Stop and start manually)'.format(exp))
                         exp.save()
             def radar_stop_scheduler(id_camp,id_radar,id_experiment):
                 '''
                 Stop experiments's devices
                 DDS-JARS-RC-CGS-ABS
                 '''
                 exp=get_object_or_404(Experiment,pk=id_experiment)
                 if exp.status == 2:
                     confs = Configuration.objects.filter(experiment=id_experiment,type = 0).order_by('device__device_type__sequence')
                     confs = confs.exclude(device__device_type__name='cgs')
                     try:
                         for conf in confs:
                             # print(conf)
                             conf.stop_device()
                             exp.status= 1
                     except:
                         exp.status= 0
                     exp.save()
                 return exp.status

apps/main/models.py +4 0

             import os
             import json
             import requests
             import time
             from datetime import datetime
             try:
                 from polymorphic.models import PolymorphicModel
             except:
                 from polymorphic import PolymorphicModel
             from django.template.base import kwarg_re
             from django.db import models
             from django.urls import reverse
             from django.core.validators import MinValueValidator, MaxValueValidator
             from django.shortcuts import get_object_or_404
             from django.contrib.auth.models import User
             from django.db.models.signals import post_save
             from django.dispatch import receiver
             from apps.main.utils import Params
             from apps.rc.utils import RCFile
             from apps.jars.utils import RacpFile
             from devices.dds import api as dds_api
             from devices.dds import data as dds_data
             DEV_PORTS = {
                             'rc'    : 2000,
                             'dds'   : 2000,
                             'jars'  : 2000,
                             'usrp'  : 2000,
                             'cgs'   : 8080,
                             'abs'   : 8080,
                             'dds_rest': 80
                         }
             RADAR_STATES = (
                              (0, 'No connected'),
                              (1, 'Connected'),
                              (2, 'Configured'),
                              (3, 'Running'),
                              (4, 'Scheduled'),
                          )
             EXPERIMENT_TYPE = (
                                (0, 'RAW_DATA'),
                                (1, 'PDATA'),
                                )
             DECODE_TYPE = (
                          (0, 'None'),
                          (1, 'TimeDomain'),
                          (2, 'FreqDomain'),
                          (3, 'InvFreqDomain'),
                          )
             DEV_STATES = (
                              (0, 'No connected'),
                              (1, 'Connected'),
                              (2, 'Configured'),
                              (3, 'Running'),
                              (4, 'Unknown'),
+                             (5, 'Busy')
                          )
             DEV_TYPES = (
                             ('', 'Select a device type'),
                             ('rc', 'Radar Controller'),
                             ('dds', 'Direct Digital Synthesizer'),
                             ('jars', 'Jicamarca Radar Acquisition System'),
                             ('usrp', 'Universal Software Radio Peripheral'),
                             ('cgs', 'Clock Generator System'),
                             ('abs', 'Automatic Beam Switching'),
                             ('dds_rest', 'Direct Digital Synthesizer_REST'),
                             ('atrad', 'Transmitter ATRAD'),
                         )
             EXP_STATES = (
                              (0,'Error'),                 #RED
                              (1,'Cancelled'),             #YELLOW
                              (2,'Running'),               #GREEN
                              (3,'Scheduled'),             #BLUE
                              (4,'Unknown'),               #WHITE
                          )
             CONF_TYPES = (
                              (0, 'Active'),
                              (1, 'Historical'),
                          )
             class Profile(models.Model):
                 user = models.OneToOneField(User, on_delete=models.CASCADE)
                 theme = models.CharField(max_length=30, default='spacelab')
             @receiver(post_save, sender=User)
             def create_user_profile(sender, instance, created, **kwargs):
                 if created:
                     Profile.objects.create(user=instance)
             @receiver(post_save, sender=User)
             def save_user_profile(sender, instance, **kwargs):
                 instance.profile.save()
             class Location(models.Model):
                 name = models.CharField(max_length = 30)
                 description = models.TextField(blank=True, null=True)
                 class Meta:
                     db_table = 'db_location'
                 def __str__(self):
                     return u'%s' % self.name
                 def get_absolute_url(self):
                     return reverse('url_location', args=[str(self.id)])
             class DeviceType(models.Model):
                 name = models.CharField(max_length = 10, choices = DEV_TYPES, default = 'dds_rest')
                 sequence = models.PositiveSmallIntegerField(default=55)
                 description = models.TextField(blank=True, null=True)
                 class Meta:
                     db_table = 'db_device_types'
                 def __str__(self):
                     return u'%s' % self.get_name_display()
             class Device(models.Model):
                 device_type  = models.ForeignKey('DeviceType', on_delete=models.CASCADE)
                 location     = models.ForeignKey('Location', on_delete=models.CASCADE)
                 ip_address   = models.GenericIPAddressField(protocol='IPv4', default='0.0.0.0')
                 port_address = models.PositiveSmallIntegerField(default=2000)
                 description  = models.TextField(blank=True, null=True)
                 status       = models.PositiveSmallIntegerField(default=4, choices=DEV_STATES)
                 conf_active  = models.PositiveIntegerField(default=0, verbose_name='Current configuration')
                 class Meta:
                     db_table = 'db_devices'
                 def __str__(self):
                     ret = u'{} [{}]'.format(self.device_type.name.upper(), self.location.name)
                     return ret
                 @property
                 def name(self):
                     return str(self)
                 def get_status(self):
                     return self.status
                 @property
                 def status_color(self):
                     color = 'muted'
                     if self.status == 0:
                         color = "danger"
                     elif self.status == 1:
                         color = "warning"
                     elif self.status == 2:
                         color = "info"
                     elif self.status == 3:
                         color = "success"
                     return color
                 def url(self, path=None):
                     if path:
                         return 'http://{}:{}/{}/'.format(self.ip_address, self.port_address, path)
                     else:
                         return 'http://{}:{}/'.format(self.ip_address, self.port_address)
                 def get_absolute_url(self):
                     return reverse('url_device', args=[str(self.id)])
                 def get_absolute_url_edit(self):
                     return reverse('url_edit_device', args=[str(self.id)])
                 def get_absolute_url_delete(self):
                     return reverse('url_delete_device', args=[str(self.id)])
                 def change_ip(self, ip_address, mask, gateway, dns, **kwargs):
                     if self.device_type.name=='dds':
                         try:
                             answer = dds_api.change_ip(ip = self.ip_address,
                                                        port = self.port_address,
                                                        new_ip = ip_address,
                                                        mask = mask,
                                                        gateway = gateway)
                             if answer[0]=='1':
                                 self.message = '25|DDS - {}'.format(answer)
                                 self.ip_address = ip_address
                                 self.save()
                             else:
                                 self.message = '30|DDS - {}'.format(answer)
                                 return False
                         except Exception as e:
                             self.message = '40|{}'.format(str(e))
                             return False
                     elif self.device_type.name=='rc':
                         headers = {'content-type': "application/json",
                                    'cache-control': "no-cache"}
                         ip      = [int(x) for x in ip_address.split('.')]
                         dns     = [int(x) for x in dns.split('.')]
                         gateway = [int(x) for x in gateway.split('.')]
                         subnet  = [int(x) for x in mask.split('.')]
                         payload = {
                             "ip": ip,
                             "dns": dns,
                             "gateway": gateway,
                             "subnet": subnet
                             }
                         req = requests.post(self.url('changeip'), data=json.dumps(payload), headers=headers)
                         try:
                             answer = req.json()
                             if answer['changeip']=='ok':
                                 self.message    = '25|IP succesfully changed'
                                 self.ip_address = ip_address
                                 self.save()
                             else:
                                 self.message = '30|An error ocuur when changing IP'
                         except Exception as e:
                             self.message = '40|{}'.format(str(e))
                     else:
                         self.message = 'Not implemented'
                         return False
                     return True
             class Campaign(models.Model):
                 template    = models.BooleanField(default=False)
                 name        = models.CharField(max_length=60, unique=True)
                 start_date  = models.DateTimeField(blank=True, null=True)
                 end_date    = models.DateTimeField(blank=True, null=True)
                 tags        = models.CharField(max_length=40, blank=True, null=True)
                 description = models.TextField(blank=True, null=True)
                 experiments = models.ManyToManyField('Experiment', blank=True)
                 author      = models.ForeignKey(User, null=True, blank=True,on_delete=models.CASCADE)
                 class Meta:
                     db_table = 'db_campaigns'
                     ordering = ('name',)
                 def __str__(self):
                     if self.template:
                         return u'{} (template)'.format(self.name)
                     else:
                         return u'{}'.format(self.name)
                 def jsonify(self):
                     data = {}
                     ignored = ('template')
                     for field in self._meta.fields:
                         if field.name in ignored:
                             continue
                         data[field.name] = field.value_from_object(self)
                     data['start_date'] = data['start_date'].strftime('%Y-%m-%d')
                     data['end_date']   = data['end_date'].strftime('%Y-%m-%d')
                     return data
                 def parms_to_dict(self):
                     params = Params({})
                     params.add(self.jsonify(), 'campaigns')
                     for exp in Experiment.objects.filter(campaign = self):
                         params.add(exp.jsonify(), 'experiments')
                         configurations = Configuration.objects.filter(experiment=exp, type=0)
                         for conf in configurations:
                             params.add(conf.jsonify(), 'configurations')
                             if conf.device.device_type.name=='rc':
                                 for line in conf.get_lines():
                                     params.add(line.jsonify(), 'lines')
                     return params.data
                 def dict_to_parms(self, parms, CONF_MODELS):
                     experiments = Experiment.objects.filter(campaign = self)
                     if experiments:
                         for experiment in experiments:
                             experiment.delete()
                     for id_exp in parms['experiments']['allIds']:
                         exp_parms = parms['experiments']['byId'][id_exp]
                         dum = (datetime.now() - datetime(1970, 1, 1)).total_seconds()
                         exp = Experiment(name='{}'.format(dum))
                         exp.save()
                         exp.dict_to_parms(parms, CONF_MODELS, id_exp=id_exp)
                         self.experiments.add(exp)
                     camp_parms = parms['campaigns']['byId'][parms['campaigns']['allIds'][0]]
                     self.name       = '{}-{}'.format(camp_parms['name'], datetime.now().strftime('%y%m%d'))
                     self.start_date = camp_parms['start_date']
                     self.end_date   = camp_parms['end_date']
                     self.tags       = camp_parms['tags']
                     self.save()
                     return self
                 def get_experiments_by_radar(self, radar=None):
                     ret = []
                     if radar:
                         locations = Location.objects.filter(pk=radar)
                     else:
                         locations = set([e.location for e in self.experiments.all()])
                     for loc in locations:
                         dum = {}
                         dum['name']        = loc.name
                         dum['id']          = loc.pk
                         dum['experiments'] = [e for e in self.experiments.all() if e.location==loc]
                         ret.append(dum)
                     return ret
                 def get_absolute_url(self):
                     return reverse('url_campaign', args=[str(self.id)])
                 def get_absolute_url_edit(self):
                     return reverse('url_edit_campaign', args=[str(self.id)])
                 def get_absolute_url_delete(self):
                     return reverse('url_delete_campaign', args=[str(self.id)])
                 def get_absolute_url_export(self):
                     return reverse('url_export_campaign', args=[str(self.id)])
                 def get_absolute_url_import(self):
                     return reverse('url_import_campaign', args=[str(self.id)])
             class RunningExperiment(models.Model):
                 radar              = models.OneToOneField('Location', on_delete=models.CASCADE)
                 running_experiment = models.ManyToManyField('Experiment', blank = True)
                 status             = models.PositiveSmallIntegerField(default=0, choices=RADAR_STATES)
             class Experiment(models.Model):
                 template   = models.BooleanField(default=False)
                 name       = models.CharField(max_length=40, default='', unique=True)
                 location   = models.ForeignKey('Location', null=True, blank=True, on_delete=models.CASCADE)
                 freq       = models.FloatField(verbose_name='Operating Freq. (MHz)', validators=[MinValueValidator(1), MaxValueValidator(10000)], default=49.9200)
                 start_time = models.TimeField(default='00:00:00')
                 end_time   = models.TimeField(default='23:59:59')
                 task       = models.CharField(max_length=36, default='', blank=True, null=True)
                 status     = models.PositiveSmallIntegerField(default=4, choices=EXP_STATES)
                 author     = models.ForeignKey(User, null=True, blank=True,on_delete=models.CASCADE)
                 hash       = models.CharField(default='', max_length=64, null=True, blank=True)
                 class Meta:
                     db_table = 'db_experiments'
                     ordering = ('template', 'name')
                 def __str__(self):
                     if self.template:
                         return u'%s (template)' % (self.name)
                     else:
                         return u'%s' % (self.name)
                 def jsonify(self):
                     data = {}
                     ignored = ('template')
                     for field in self._meta.fields:
                         if field.name in ignored:
                             continue
                         data[field.name] = field.value_from_object(self)
                     data['start_time'] = data['start_time'].strftime('%H:%M:%S')
                     data['end_time']   = data['end_time'].strftime('%H:%M:%S')
                     data['location']   = self.location.name
                     data['configurations'] = ['{}'.format(conf.pk) for
                         conf in Configuration.objects.filter(experiment=self, type=0)]
                     return data
                 @property
                 def radar_system(self):
                     return self.location
                 def clone(self, **kwargs):
                     confs = Configuration.objects.filter(experiment=self, type=0)
                     self.pk = None
                     self.name = '{}_{:%y%m%d%H%M%S}'.format(self.name, datetime.now())
                     for attr, value in kwargs.items():
                         setattr(self, attr, value)
                     self.save()
                     for conf in confs:
                         conf.clone(experiment=self, template=False)
                     return self
                 def start(self):
                     '''
                     Configure and start experiments's devices
                     ABS-CGS-DDS-RC-JARS
                     '''
                     confs = []
                     allconfs = Configuration.objects.filter(experiment=self, type = 0).order_by('-device__device_type__sequence')
                     rc_mix   = [conf for conf in allconfs if conf.device.device_type.name=='rc' and conf.mix]
                     if rc_mix:
                         for conf in allconfs:
                             if conf.device.device_type.name == 'rc' and  not conf.mix:
                                 continue
                             confs.append(conf)
                     else:
                         confs = allconfs
                     print("confs: ",confs)
                     #try:
                     for conf in confs:
                         print("conf->",conf)
                         conf.stop_device()
                         conf.write_device()
                         conf.device.conf_active = conf.pk
                         conf.device.save()
                         conf.start_device()
                         time.sleep(1)
                     #except:
                         #return 0
                     return 2
                 def stop(self):
                     '''
                     Stop experiments's devices
                     DDS-JARS-RC-CGS-ABS
                     '''
                     confs = Configuration.objects.filter(experiment=self, type = 0).order_by('device__device_type__sequence')
                     confs = confs.exclude(device__device_type__name='cgs')
                     try:
                         for conf in confs:
                             conf.stop_device()
                     except:
                         return 0
                     return 1
                 def get_status(self):
                     if self.status == 3:
                         return
                     confs =  Configuration.objects.filter(experiment=self, type=0)
                     for conf in confs:
                         conf.status_device()
                     total = confs.aggregate(models.Sum('device__status'))['device__status__sum']
                     if total==2*confs.count():
                         status = 1
                     elif total == 3*confs.count():
                         status = 2
                     else:
                         status = 0
                     self.status = status
                     self.save()
                 def status_color(self):
                     color = 'muted'
                     if self.status == 0:
                         color = "danger"
                     elif self.status == 1:
                         color = "warning"
                     elif self.status == 2:
                         color = "success"
                     elif self.status == 3:
                         color = "info"
                     return color
                 def parms_to_dict(self):
                     params = Params({})
                     params.add(self.jsonify(), 'experiments')
                     configurations = Configuration.objects.filter(experiment=self, type=0)
                     for conf in configurations:
                         params.add(conf.jsonify(), 'configurations')
                         if conf.device.device_type.name=='rc':
                             for line in conf.get_lines():
                                 params.add(line.jsonify(), 'lines')
                     return params.data
                 def dict_to_parms(self, parms, CONF_MODELS, id_exp=None):
                     configurations = Configuration.objects.filter(experiment=self)
                     if id_exp is not None:
                         exp_parms = parms['experiments']['byId'][id_exp]
                     else:
                         exp_parms = parms['experiments']['byId'][parms['experiments']['allIds'][0]]
                     if configurations:
                         for configuration in configurations:
                             configuration.delete()
                     for id_conf in exp_parms['configurations']:
                         conf_parms = parms['configurations']['byId'][id_conf]
                         device     = Device.objects.filter(device_type__name=conf_parms['device_type'])[0]
                         model      = CONF_MODELS[conf_parms['device_type']]
                         conf = model(
                             experiment = self,
                             device = device,
                             )
                         conf.dict_to_parms(parms, id=id_conf)
                     location, created = Location.objects.get_or_create(name=exp_parms['location'])
                     self.name = '{}-{}'.format(exp_parms['name'], datetime.now().strftime('%y%m%d'))
                     self.location   = location
                     self.start_time = exp_parms['start_time']
                     self.end_time   = exp_parms['end_time']
                     self.save()
                     return self
                 def get_absolute_url(self):
                     return reverse('url_experiment', args=[str(self.id)])
                 def get_absolute_url_edit(self):
                     return reverse('url_edit_experiment', args=[str(self.id)])
                 def get_absolute_url_delete(self):
                     return reverse('url_delete_experiment', args=[str(self.id)])
                 def get_absolute_url_import(self):
                     return reverse('url_import_experiment', args=[str(self.id)])
                 def get_absolute_url_export(self):
                     return reverse('url_export_experiment', args=[str(self.id)])
                 def get_absolute_url_start(self):
                     return reverse('url_start_experiment', args=[str(self.id)])
                 def get_absolute_url_stop(self):
                     return reverse('url_stop_experiment', args=[str(self.id)])
             class Configuration(PolymorphicModel):
                 template        = models.BooleanField(default=False)
                 # name          = models.CharField(verbose_name="Configuration Name", max_length=40, default='')
                 device          = models.ForeignKey('Device', verbose_name='Device', null=True, on_delete=models.CASCADE)
                 label           = models.CharField(verbose_name="Label", max_length=40, default='', blank=True, null=True)
                 experiment      = models.ForeignKey('Experiment', verbose_name='Experiment', null=True, blank=True, on_delete=models.CASCADE)
                 type            = models.PositiveSmallIntegerField(default=0, choices=CONF_TYPES)
                 created_date    = models.DateTimeField(auto_now_add=True)
                 programmed_date = models.DateTimeField(auto_now=True)
                 parameters      = models.TextField(default='{}')
                 author          = models.ForeignKey(User, null=True, blank=True,on_delete=models.CASCADE)
                 hash            = models.CharField(default='', max_length=64, null=True, blank=True)
                 message = ""
                 class Meta:
                     db_table = 'db_configurations'
                     ordering = ('device__device_type__name',)
                 def __str__(self):
                     ret = u'{} '.format(self.device.device_type.name.upper())
                     if 'mix' in [f.name for f in self._meta.get_fields()]:
                         if self.mix:
                             ret = '{} MIX '.format(self.device.device_type.name.upper())
                     if 'label' in [f.name for f in self._meta.get_fields()]:
                         ret += '{}'.format(self.label)
                     if self.template:
                         ret += ' (template)'
                     return ret
                 @property
                 def name(self):
                     return str(self)
                 def jsonify(self):
                     data = {}
                     ignored = ('type', 'polymorphic_ctype', 'configuration_ptr',
                                'created_date', 'programmed_date', 'template', 'device',
                                'experiment')
                     for field in self._meta.fields:
                         if field.name in ignored:
                             continue
                         data[field.name] = field.value_from_object(self)
                     data['device_type'] = self.device.device_type.name
                     if self.device.device_type.name == 'rc':
                         data['lines'] = ['{}'.format(line.pk) for line in self.get_lines()]
                         data['delays'] = self.get_delays()
                         data['pulses'] = self.get_pulses()
                     elif self.device.device_type.name == 'jars':
                         data['decode_type'] = DECODE_TYPE[self.decode_data][1]
                     elif self.device.device_type.name == 'dds':
                         data['frequencyA_Mhz'] = float(data['frequencyA_Mhz'])
                         data['frequencyB_Mhz'] = float(data['frequencyB_Mhz'])
                         data['phaseA'] = dds_data.phase_to_binary(data['phaseA_degrees'])
                         data['phaseB'] = dds_data.phase_to_binary(data['phaseB_degrees'])
                     elif self.device.device_type.name == 'dds_rest':
                         data['frequencyA_Mhz'] = float(data['frequencyA_Mhz'])
                         data['frequencyB_Mhz'] = float(data['frequencyB_Mhz'])
                         data['phaseA'] = dds_data.phase_to_binary(data['phaseA_degrees'])
                         data['phaseB'] = dds_data.phase_to_binary(data['phaseB_degrees'])
                         data['delta_frequency_Mhz'] = float(data['delta_frequency_Mhz'] or 0.00)
                         data['update_clock_Mhz']    = float(data['update_clock_Mhz'] or 0.00)
                         data['ramp_rate_clock_Mhz'] = float(data['ramp_rate_clock_Mhz'] or 0.0)
                     return data
                 def clone(self, **kwargs):
                     self.pk = None
                     self.id = None
                     for attr, value in kwargs.items():
                         setattr(self, attr, value)
                     self.save()
                     return self
                 def parms_to_dict(self):
                     params = Params({})
                     params.add(self.jsonify(), 'configurations')
                     if self.device.device_type.name=='rc':
                         for line in self.get_lines():
                             params.add(line.jsonify(), 'lines')
                     return params.data
                 def parms_to_text(self):
                     raise NotImplementedError("This method should be implemented in %s Configuration model" %str(self.device.device_type.name).upper())
                 def parms_to_binary(self):
                     raise NotImplementedError("This method should be implemented in %s Configuration model" %str(self.device.device_type.name).upper())
                 def dict_to_parms(self, parameters, id=None):
                     params = Params(parameters)
                     if id:
                         data = params.get_conf(id_conf=id)
                     else:
                         data = params.get_conf(dtype=self.device.device_type.name)
                     if data['device_type']=='rc':
                         self.clean_lines()
                         lines = data.pop('lines', None)
                         for line_id in lines:
                             pass
                     for key, value in data.items():
                         if key not in ('id', 'device_type'):
                             setattr(self, key, value)
                     self.save()
                 def export_to_file(self, format="json"):
                     content_type = ''
                     if format == 'racp':
                         content_type = 'text/plain'
                         filename = '%s_%s.%s' %(self.device.device_type.name, self.name, 'racp')
                         content = self.parms_to_text(file_format = 'racp')
                     if format == 'text':
                         content_type = 'text/plain'
                         filename = '%s_%s.%s' %(self.device.device_type.name, self.name, self.device.device_type.name)
                         content = self.parms_to_text()
                     if format == 'binary':
                         content_type = 'application/octet-stream'
                         filename = '%s_%s.bin' %(self.device.device_type.name, self.name)
                         content = self.parms_to_binary()
                     if not content_type:
                         content_type = 'application/json'
                         filename = '%s_%s.json' %(self.device.device_type.name, self.name)
                         content = json.dumps(self.parms_to_dict(), indent=2)
                     fields = {'content_type':content_type,
                               'filename':filename,
                               'content':content
                               }
                     return fields
                 def import_from_file(self, fp):
                     parms = {}
                     path, ext = os.path.splitext(fp.name)
                     if ext == '.json':
                         parms = json.load(fp)
                     if ext == '.dds':
                         lines = fp.readlines()
                         parms = dds_data.text_to_dict(lines)
                     if ext == '.racp':
                         if self.device.device_type.name == 'jars':
                             parms = RacpFile(fp).to_dict()
                             parms['filter_parms'] = json.loads(self.filter_parms)
                             return parms
                         parms = RCFile(fp).to_dict()
                     return parms
                 def status_device(self):
                     self.message = 'Function not implemented'
                     return False
                 def stop_device(self):
                     self.message = 'Function not implemented'
+                    print("BUENAS SEÑALES??? NO LO CREO2",flush=True)
                     return False
                 def start_device(self):
                     self.message = 'Function not implemented'
+                    print("BUENAS SEÑALES??? NO LO CREO",flush=True)
                     return False
                 def write_device(self, parms):
                     self.message = 'Function not implemented'
+                    print("BUENAS SEÑALES??? NO LO CREO3",flush=True)
                     return False
                 def write_device_mqtt(self, parms):
                     self.message = 'Function not implemented'
                     return False
                 def read_device(self):
                     self.message = 'Function not implemented'
                     return False
                 def get_absolute_url(self):
                     return reverse('url_%s_conf' % self.device.device_type.name, args=[str(self.id)])
                 def get_absolute_url_edit(self):
                     return reverse('url_edit_%s_conf' % self.device.device_type.name, args=[str(self.id)])
                 def get_absolute_url_delete(self):
                     return reverse('url_delete_dev_conf', args=[str(self.id)])
                 def get_absolute_url_import(self):
                     return reverse('url_import_dev_conf', args=[str(self.id)])
                 def get_absolute_url_export(self):
                     return reverse('url_export_dev_conf', args=[str(self.id)])
                 def get_absolute_url_write(self):
                     return reverse('url_write_dev_conf', args=[str(self.id)])
                 def get_absolute_url_write_mqtt(self):
                     return reverse('url_write_mqtt_dev_conf', args=[str(self.id)])
                 def get_absolute_url_read(self):
                     return reverse('url_read_dev_conf', args=[str(self.id)])
                 def get_absolute_url_start(self):
                     return reverse('url_start_dev_conf', args=[str(self.id)])
                 def get_absolute_url_stop(self):
                     return reverse('url_stop_dev_conf', args=[str(self.id)])
                 def get_absolute_url_stop_mqtt(self):
                     return reverse('url_stop_mqtt_dev_conf', args=[str(self.id)])
                 def get_absolute_url_status(self):
                     return reverse('url_status_dev_conf', args=[str(self.id)])

General Comments 0

Write
Preview

You need to be logged in to leave comments. Login now

No TODOs yet

	Site-wide shortcuts
/	Use quick search box
g h	Goto home page
g g	Goto my private gists page
g G	Goto my public gists page
g 0-9	Goto bookmarked items from 0-9
n r	New repository page
n g	New gist page

	Repositories
g s	Goto summary page
g c	Goto changelog page
g f	Goto files page
g F	Goto files page with file search activated
g p	Goto pull requests page
g o	Goto repository settings
g O	Goto repository access permissions settings
t s	Toggle sidebar on some pages