program elc
c         Header(1)=Rlat
c         Header(2)=Rlon
c         Header(3)=IASI_Satellite_Zenith(J,I)
c         Header(4)=Land_Sea_Alti(J,I)
c         Header(5)=cloud_cover(J,I)
c         Header(6)=hh*1.e4+mn*1.e2+ss
c         Header(7)=float(indtm*10+indjn)
c         Header(8)=yyyymmdd
      real head(8),mots(339)
      real ang_sol,norb
      integer indi(339)
      character*200 filein
c header : rlat,rlon,angle_zen,
      call getarg(1,filein)
      open(1,file=filein,form='unformatted')
      read(1) nbmots,(indi(i),i=1,nbmots)
       print *,'nbmots',nbmots,(indi(i),i=1,nbmots)

! nbre de temperatures de brillance extraites = nbmots - nbre mots header - angle solaire - numero orbite
      nb_TBs=nbmots-8-2
      write(*,*) nb_TBs
      icpt=0
      ic=0
      imn=0
      imj=0
      itn=0
      itj=0
! ixx calcule le nombre de points dans 4 classes combinant les indices terre/mer/jour/nuit
* xx=mn pour mer/nuit
 10   read(1,end=100) head,(mots(i),i=1,nb_TBs),ang_sol,norb
	print * ,"head=",nb_TBs
      sectta=1./cosd(head(3))
      !write(*,1) head(1),head(6),head(8),mots(1),ang_sol,int(norb)
1     format(f8.2,2f10.0,2f8.2,i7)
      if(head(7).eq.0) imn=imn+1
      if(head(7).eq.1) imj=imj+1
      if(head(7).eq.10) itn=itn+1
      if(head(7).eq.11) itj=itj+1
      icpt=icpt+1
      ic=ic+1
      goto 10
 100  continue
 14    format(8f11.3,2i4)
 15    format(3f11.3,2i4)
       print *,icpt,imn,imj,itn,itj
      end