Computes the Sun positions for the *t* time vector.
        
        *t* have to be in absolute minutes (0 at 00:00 01 Jan). The and 
        in Sun positions calculated are in solar time, that is, maximun 
        solar zenit exactly at midday.
        
        The generated information is stored in:
            * **.azimuth_zenit** (*ndarray*)
            * **.true_time** (*datetime*): local time
        
        :param t: Absolute minutes vector.
        :type t: ndarray (dtype=int)
        :param true_time: If True, a datetime vector with the true local
            time will be stored at ``.true_time``
        :type true_time: bool
        :returns: Equivalent times in absolute minutes in year.
        :rtype: ndarray (dtype=int)
        
        :returns: None
        
        .. seealso:: :func:`to_minutes` to easily genetare valid input 
            t.

def get_sunpos(self, t, true_time=False):
        """
        Computes the Sun positions for the *t* time vector.
        
        *t* have to be in absolute minutes (0 at 00:00 01 Jan). The and 
        in Sun positions calculated are in solar time, that is, maximun 
        solar zenit exactly at midday.
        
        The generated information is stored in:
            * **.azimuth_zenit** (*ndarray*)
            * **.true_time** (*datetime*): local time
        
        :param t: Absolute minutes vector.
        :type t: ndarray (dtype=int)
        :param true_time: If True, a datetime vector with the true local
            time will be stored at ``.true_time``
        :type true_time: bool
        :returns: Equivalent times in absolute minutes in year.
        :rtype: ndarray (dtype=int)
        
        :returns: None
        
        .. seealso:: :func:`to_minutes` to easily genetare valid input 
            t.
        """
        import numpy as np
        lat = self.arg_latitude
        long = self.arg_longitude
        alphamin = self.arg_zenitmin
        
        # Solar calculations
        day = np.modf(t/1440)[0]
        fractional_year = 2*np.pi/(365*24*60)*(-24*60+t)
        declination = 0.006918 - \
                      0.399912*np.cos(fractional_year) + \
                      0.070257*np.sin(fractional_year) - \
                      0.006758*np.cos(2*fractional_year) + \
                      0.000907*np.sin(2*fractional_year) - \
                      0.002697*np.cos(3*fractional_year) + \
                      0.00148*np.sin(3*fractional_year)
        
        hour_angle = np.tile(np.arange(-np.pi, np.pi, 2*np.pi/(24*60), 
                                       dtype='float'), 365)[t]
        solar_zenit = np.arcsin(np.sin(lat)*np.sin(declination) + \
                      np.cos(lat)*np.cos(declination)*np.cos(hour_angle))
        solar_zenit[solar_zenit<=0+alphamin] = np.nan
        #### Avoiding numpy warning
        aux = (np.sin(solar_zenit)*np.sin(lat) - np.sin(declination))/ \
              (np.cos(solar_zenit)*np.cos(lat))
        not_nan = np.logical_not(np.isnan(aux))
        aux_1 = aux[not_nan]
        aux_1[aux_1>=1] = np.nan
        aux[not_nan] = aux_1
        ####
        solar_azimuth = np.arccos(aux)
        solar_azimuth[day==0.5] = 0
        solar_azimuth[day<0.5] *= -1
        self.azimuth_zenit = np.vstack((solar_azimuth, solar_zenit)).T
                
        # True time
        if true_time:
            import datetime as dt
            long = np.rad2deg(long)
            instant_0 = dt.datetime(1,1,1,0,0,0) # Simulator time 
            
            # Real time
            equation_time = 229.18*(0.000075+0.001868*np.cos(fractional_year) - \
                            0.032077*np.sin(fractional_year) - \
                            0.014615*np.cos(2*fractional_year) - \
                            0.040849*np.sin(2*fractional_year))
        
            time_offset = equation_time + 4*long + 60*self.arg_UTC
            true_solar_time = t + time_offset
            delta_true_date_objs = np.array([dt.timedelta(minutes=i) 
                                             for i in true_solar_time])
            self.true_time = instant_0 + delta_true_date_objs