!*******************************************************
      
      program jupiter
!     Version: 00/00/03
      
!     etat_eos : equation d'etat saumon-chabrier
!     opa_f : opacite 
!     convp_jmj : convection
!     nuc_0 : reactions thermonucleaires
!     lim_0 : conditions limites externes atmosphere 
!     hopf : loi t(tau) utilisee
      
      implicit none
      integer irep
      character*20 c_etat,c_opa,c_lim
      character*30 type(0:5)
      data type/'ADIABATIQUE','NON-ADIABATIQUE','AD-PROTOPLANETE',      &
     &     'NA-PROTOPLANETE','A-POLYTROPE N=1','R-POLYTROPE N=3'/

      external etat_eos,etat_hhe,etat_div,                              &
     &     opaque,opa_f,opa_alex,opa_star,                              &
     &     convp_jmj,convp_ad,nuc_0,nuc_cte,                            &
     &     lim_0,lim_x,lim_xa3,lim_xa1000,hopf

!----------------------------------------------------------------------
      write(*,101)
 101  format(30('*'),' C E P A M ',30('*'),//,                          &
     &     'Indiquer l''equation d''etat a utiliser (defaut: etat_hhe)')
      read(*,'(a)')c_etat
      if (c_etat(:1).eq.' ') c_etat='etat_hhe'
      write(*,*)'Indiquer le programe d''opacite (defaut: opa_alex)'
      read(*,'(a)')c_opa
      if (c_opa(:1).eq.' ') c_opa='opa_alex'
      write(*,*)'Indiquer le programe de conditions limites (defaut:',  &
     &     ' lim_xa3)'
      read(*,'(a)')c_lim
      if (c_lim(:1).eq.' ') c_lim='lim_xa3'

      if (c_etat(:8).eq.'etat_hhe') then
       if (c_opa(:6).eq.'opaque') then
        if (c_lim(:5).eq.'lim_0') then
         call cepam(type(0),0,etat_hhe,opaque,convp_jmj,nuc_0,          &
     &      lim_0,hopf)
        elseif (c_lim(:7).eq.'lim_x  ') then
         call cepam(type(0),0,etat_hhe,opaque,convp_jmj,nuc_0,          &
     &      lim_x,hopf)
        elseif (c_lim(:7).eq.'lim_xa3') then
         call cepam(type(0),0,etat_hhe,opaque,convp_jmj,nuc_0,          &
     &      lim_xa3,hopf)
        elseif (c_lim(:10).eq.'lim_xa1000') then
         call cepam(type(0),0,etat_hhe,opaque,convp_jmj,nuc_0,          &
     &      lim_xa1000,hopf)
        endif
       elseif (c_opa(:5).eq.'opa_f') then
        if (c_lim(:5).eq.'lim_0') then
         call cepam(type(0),0,etat_hhe,opa_f,convp_jmj,nuc_0,           &
     &      lim_0,hopf)
        elseif (c_lim(:7).eq.'lim_x  ') then
         call cepam(type(0),0,etat_hhe,opa_f,convp_jmj,nuc_0,           &
     &      lim_x,hopf)
        elseif (c_lim(:7).eq.'lim_xa3') then
         call cepam(type(0),0,etat_hhe,opa_f,convp_jmj,nuc_0,           &
     &      lim_xa3,hopf)
        elseif (c_lim(:10).eq.'lim_xa1000') then
         call cepam(type(0),0,etat_hhe,opa_f,convp_jmj,nuc_0,           &
     &      lim_xa1000,hopf)
        endif
       elseif (c_opa(:8).eq.'opa_star') then
        if (c_lim(:5).eq.'lim_0') then
         call cepam(type(0),0,etat_hhe,opa_star,convp_jmj,nuc_0,        &
     &      lim_0,hopf)
        elseif (c_lim(:7).eq.'lim_x  ') then
         call cepam(type(0),0,etat_hhe,opa_star,convp_jmj,nuc_0,        &
     &      lim_x,hopf)
        elseif (c_lim(:7).eq.'lim_xa3') then
         call cepam(type(0),0,etat_hhe,opa_star,convp_jmj,nuc_0,        &
     &      lim_xa3,hopf)
        elseif (c_lim(:10).eq.'lim_xa1000') then
         call cepam(type(0),0,etat_hhe,opa_star,convp_jmj,nuc_0,        &
     &      lim_xa1000,hopf)
        endif
       elseif (c_opa(:8).eq.'opa_alex') then
        if (c_lim(:5).eq.'lim_0') then
         call cepam(type(0),0,etat_hhe,opa_alex,convp_jmj,nuc_0,        &
     &      lim_0,hopf)
        elseif (c_lim(:7).eq.'lim_x  ') then
         call cepam(type(0),0,etat_hhe,opa_alex,convp_jmj,nuc_0,        &
     &      lim_x,hopf)
        elseif (c_lim(:7).eq.'lim_xa3') then
         call cepam(type(0),0,etat_hhe,opa_alex,convp_jmj,nuc_0,        &
     &      lim_xa3,hopf)
        elseif (c_lim(:10).eq.'lim_xa1000') then
         call cepam(type(0),0,etat_hhe,opa_alex,convp_jmj,nuc_0,        &
     &      lim_xa1000,hopf)
        endif
       endif
      elseif (c_etat(:8).eq.'etat_div') then
       if (c_opa(:6).eq.'opaque') then
        if (c_lim(:5).eq.'lim_0') then
         call cepam(type(0),0,etat_div,opaque,convp_jmj,nuc_0,          &
     &      lim_0,hopf)
        elseif (c_lim(:7).eq.'lim_x  ') then
         call cepam(type(0),0,etat_div,opaque,convp_jmj,nuc_0,          &
     &      lim_x,hopf)
        elseif (c_lim(:7).eq.'lim_xa3') then
         call cepam(type(0),0,etat_div,opaque,convp_jmj,nuc_0,          &
     &      lim_xa3,hopf)
        elseif (c_lim(:10).eq.'lim_xa1000') then
         call cepam(type(0),0,etat_div,opaque,convp_jmj,nuc_0,          &
     &      lim_xa1000,hopf)
        endif
       elseif (c_opa(:5).eq.'opa_f') then
        if (c_lim(:5).eq.'lim_0') then
         call cepam(type(0),0,etat_div,opa_f,convp_jmj,nuc_0,           &
     &      lim_0,hopf)
        elseif (c_lim(:7).eq.'lim_x  ') then
         call cepam(type(0),0,etat_div,opa_f,convp_jmj,nuc_0,           &
     &      lim_x,hopf)
        elseif (c_lim(:7).eq.'lim_xa3') then
         call cepam(type(0),0,etat_div,opa_f,convp_jmj,nuc_0,           &
     &      lim_xa3,hopf)
        elseif (c_lim(:10).eq.'lim_xa1000') then
         call cepam(type(0),0,etat_div,opa_f,convp_jmj,nuc_0,           &
     &      lim_xa1000,hopf)
        endif
       elseif (c_opa(:8).eq.'opa_star') then
        if (c_lim(:5).eq.'lim_0') then
         call cepam(type(0),0,etat_div,opa_star,convp_jmj,nuc_0,        &
     &      lim_0,hopf)
        elseif (c_lim(:7).eq.'lim_x  ') then
         call cepam(type(0),0,etat_div,opa_star,convp_jmj,nuc_0,        &
     &      lim_x,hopf)
        elseif (c_lim(:7).eq.'lim_xa3') then
         call cepam(type(0),0,etat_div,opa_star,convp_jmj,nuc_0,        &
     &      lim_xa3,hopf)
        elseif (c_lim(:10).eq.'lim_xa1000') then
         call cepam(type(0),0,etat_div,opa_star,convp_jmj,nuc_0,        &
     &      lim_xa1000,hopf)
        endif
       elseif (c_opa(:8).eq.'opa_alex') then
        if (c_lim(:5).eq.'lim_0') then
         call cepam(type(0),0,etat_div,opa_alex,convp_jmj,nuc_0,        &
     &      lim_0,hopf)
        elseif (c_lim(:7).eq.'lim_x  ') then
         call cepam(type(0),0,etat_div,opa_alex,convp_jmj,nuc_0,        &
     &      lim_x,hopf)
        elseif (c_lim(:7).eq.'lim_xa3') then
         call cepam(type(0),0,etat_div,opa_alex,convp_jmj,nuc_0,        &
     &      lim_xa3,hopf)
        elseif (c_lim(:10).eq.'lim_xa1000') then
         call cepam(type(0),0,etat_div,opa_alex,convp_jmj,nuc_0,        &
     &      lim_xa1000,hopf)
        endif
       endif
      endif


!      elseif (irep.eq.4) then
!       call cepam(type(4),0,etat_gph,opaque,convp_ad,nuc_0,
!     &      lim_poly,hopf)
!      elseif (irep.eq.5) then
!       call cepam(type(5),0,etat_gph,opa_kramers,convp_jmj,nuc_cte,
!     &      lim_poly,hopf)
!      endif

      end