; Programme de detection des signaux "Pulsadis" EJP
; Ne nécessite qu'une entrée secteur (ici AN8) sur un PIC16F690
; L'affichage du signal se fait sur un ecran LCD 2x16
;
; Contact: pulsadis at free . fr
;
; Mars 2015 - V 1.1
;
; Code GCBASIC
;
;
;Chip Settings
#chip 16F690, 8
#config FOSC_INTOSC, WDTE_OFF

;Defines (Constants)                 'parametrage ecran LCD affichage
#define LCD_IO 4
#define LCD_RS PORTC.5
#define LCD_NO_RW
#define LCD_Enable PORTC.4
#define LCD_DB4 PORTC.0
#define LCD_DB5 PORTC.1
#define LCD_DB6 PORTC.2
#define LCD_DB7 PORTC.3
#define ADSpeed InternalClock
#define MaxjEjp 22
#define BoutonEjp    porta.0         'Bouton ajustement nbre jours EJP sur A0


Dir PORTC.6 In
Dir BoutonEjp In

InitTimer1 OSC,PS1_1                 'Initialisation timer 1

dim Jour,Demain,Jour_Left,Jour_Save as string*3

dim cumul1,cumul2,cumul3,cumul4 as word
Dim TMR1Load alias TMR1H,TMR1L as word  'Permet de reinitialiser Timer1
dim treshold1,target as word
dim vect(40)
vs5=0
vs15=0

target=62700                         'Variable-clé de Detect175Hz !
treshold1=30
treshold2=15
EPRead 0,NjEJP                     'Lect des jours restant EJP dans EEPROM
if NjEjp > MaxjEjp then NjEjp = 0
jour="---"
demain="???"

Function Detect175Hz                    'Détection du signal 175Hz
; on initialise la variable timer1 à la valeur de target et on attend l'overflow
; du timer lorsqu'il atteind sa valeur max soit 65535
; Si target=62700 il s'écoulera 65535-62700=2835 soit à 8Mhz 2835/2=1417,5 msec
; La valeur 62700 a été déterminée par optimisation (voir blog)
; Elle devra impérativement être revue si le code de la fonction ou si la
; fréquence du processeur sont modifiés
  cumul1=0
  cumul2=0
  cumul3=0
  cumul4=0
  index4=0
  for index = 1 to 28                   '28 relevés sur 40 ms
    index4=index4+1
    TMR1ON = 0                          'Arrêt timer1
    TMR1Load = target                   'Initialisation de timer1
    TMR1IF = 0                          'Flag d'interruption timer1 initialisé
    TMR1ON = 1                          'Démarrage timer1
    readAN8=ReadAD(an8)                 'Lecture signal analogique sur AN8
    select case index4
      case 1: cumul1=cumul1+ReadAN8
      case 2: cumul2=cumul2+ReadAN8
      case 3: cumul3=cumul3+ReadAN8
      case 4: cumul4=cumul4+ReadAN8
    end select
    if index4 > 3 then
      index4=0
    end if
    wait until TMR1IF                   'on attend l'overflow du timer1
  next
  if (cumul1 > cumul3) then
    cumul1 = cumul1 - cumul3
  else
    cumul1 = cumul3 - cumul1
  end if
  if (cumul2 > cumul4) then
    cumul2 = cumul2 - cumul4
  else
    cumul2 = cumul4 - cumul2
  end if
  cumul1=cumul1+cumul2
  if (cumul1 > treshold1) then
    Detect175Hz = 1
  else
    Detect175Hz = 0
  end if
end function

Sub Startbit                       'On attend la trame de démarrage
  Do
    if Detect175Hz = 1 then        'Si signal 175Hz sur 40 ms détecté ...
      nbtrame = 1+testbit(24)      '... on decode les 24 x 40ms suivantes
      returnsub = 1
      if (nbtrame > treshold2 ) then
        savetrame=nbtrame
        returnsub = 2              'signal OK :-)
        exit sub
      end if
    end if
    Jour_Save = Jour
    Check_Bouton                    'On vérifie état du bouton
  Loop
end Sub

Sub Check_Bouton
  if BoutonEjp = On Then
    wait 10 ms                       'Délai de stabilisation (debounce time)
    ButtonCount = 0
    Do While BoutonEjp = On
      wait 10 ms
      ButtonCount += 1
    Loop
    if ButtonCount > 40 then         'Appui long: Règlage puis affichage
      NjEjp = NjEjp + 1
      if NjEjp > MaxjEjp then NjEjp = 0
      EPWrite 0, NjEjp
    end if                            'Appui court: Affichage seulement
    jour="EJP"
    Affiche_EJP                      'On affiche les jours du Bouton
    wait 1 s
    Jour = Jour_save
    Affiche_EJP                     'On restaure l'état initial de l'écran
  end if
end sub

Sub Decode_frame
  wait 1248 ms                      'On attend 1248 ms (1250ms en théorie)
  for i1=1 to 40                    'On décode les 40 signaux Pulsadis...
    wait 1500 ms                    '... espacés de 1,5 sec
    vect(i1)=testbit(25)            'on stocke le nombre de trames détectées
  next
  returnsub = 3
end sub

Function testbit(ntime)             'appel ntime fois fonction Detect175Hz
  nbtrame=0
  repeat ntime
    nbtrame = nbtrame + Detect175Hz
  end repeat
  testbit = nbtrame                 '...et retourne le nbre trames détectées
End function

Sub Process_frame                   'Traitement du signal pour affichage

  for index = 1 to 40               'Vect est transformé en 0 ou 1 selon seuil
    if (vect(index) > treshold2 ) then
      vect(index) = 1
    else
      vect(index)=0
    end if
  next
  ; démarrage effctif de la journée EJP après 30' de préavis
  ; ce flag restera actif jusqu'à l'apparition du flag fin de journée EJP
  if vs5 & vs15 & (1-vect(5)) & (1-vect(15)) then
    ejpstart = 1
  end if
  ; fin de la journée EJP
  if (1-vect(5)) & vect(15) then
    ejpstart = 0
  end if
  jour = "---"
  if ejpstart then
    jour = "EJP"
    njEjp = njEjp + 1
    if njEjp > maxjEjp then NjEjp = 0
    EPWrite 0, NjEjp
  end if
; signal de préavis 30' avant début journée EJP
  flash_ejp = 0
  if vs5 & (1-vs15) & vect(5) & vect(15) then
    jour = "EJP"
    flash_ejp = 1
  end if
  demain = "---"
  if vect(5) & (1-vect(15)) then
    demain = "EJP"
  end if
  if vs5 & (1-vs15) & (1-vect(5)) & vect(15) then
    demain = "EJP"
  end if
; Save bits 5 & 15 pour traitement suivant
  vs5=vect(5)
  vs15=vect(15)

  Affiche_EJP

  returnsub = 1
end sub

sub Affiche_EJP              'Affichage du statut EJP sur LCD 2x16
  Select case jour
    case "---": Jour_left = "     "
    case "EJP":
      Jour_left = Pad0(NjEjp)+"/"
      Jour_left = Jour_left+Pad0(MaxjEjp-NjEjp)
  end Select
  cls
  locate 0,8
  print Jour_Left
  locate 1,2
  print "demain: "
  locate 1,10
  print Demain
  if flash_ejp then
    repeat 1620                      'flash pendant 27 mn ... (27 x 60sec)
      locate 0,2
      print jour
      wait 500 ms
      locate 0,2
      print "   "
      wait 500 ms
    end repeat
  else
    locate 0,2
    print jour
  end if
end sub

function pad0(in var) as string*2
  if var < 10 then
    pad0=" "+str(var)
  else
    pad0=str(var)
  end if
end Function

; Main program
;
cls
wait 1 s
Affiche_EJP
returnsub = 1
do
  Select case returnsub
    case 1: Startbit
    case 2: Decode_frame
    case 3: Process_frame
  end Select
loop