noscf : obtention d'orbitales optimisees par iterations sur la
        matrice densite.

genere les donnees et appelle les programmes :

- motra
- molcost ou molcsd
- casdet
- casdi
- locnats

puis itere jusqu'a convergence

Utilisation :
cp "Orbitales" INPORB   (avec schmud : cp LOCORB INPORB)
export CurrDir WorkDir Project
export PATH=$PATH:"lieu_de_cost"/bin
noscf donnees   (ET NON i"noscf < donnees !!!)
orbitales de sortie : LOCORB2



Donnees :
 &nofil PREFIX=
 /
 &no 
 data=' '         'AUTO' lit quelques donnees d'un fichier de schmudorb
                  les donnees de schmudorb pilotent ainsi une grande partie de 
                  la suite du calcul (defaut ' ')
 FERMI=0          par sym, nb d'orb occupees
                  (lu de schmudorb si data='AUTO')
 MOLCAS=5         (ou 4)
 frozen=0         donne par sym
                  (lu de schmudorb si data='AUTO')
 delete=0         donne par sym
                  (lu de schmudorb si data='AUTO')
 maxit=50         nb max de cycles motra...locnats
 r=' > '          (les fichiers de sortie sont dans la
                  directory de travail (sinon r=' # ')
 nmat=1           nb de fichiers de mat densite 
 netat=           nb d'etats par fichier (nmat donnees)
 metat=           numeros des etats demandes (Somme(netat)) donnees
 coef=1.d0        coefs des etats (Somme(netat)) donnees
 stop=' '         ('MOLCOST' pour arret apres molcsd, donne la numerotation
                  des orbitales)
 seuil=1.d-4      seuil de convergence du processus (sur les elements des
                  blocs exterieurs de la matrice densit)    
 /
 &casdi           (nmat fois)
 syspin='0',
 nessai=0,
 iprec=6
 noac=0,          (lu de schmudorb si data='AUTO')
 numac=           (lu de schmudorb si data='AUTO')
 NELAC=0,         (Attention, non lu ! A donner imperativement)
 gener='CAS+S',
 is0=1
 ms2=0
 /

les variables en majuscules cont celles obligatoires
les valeurs donnees sont celles par defaut


Donnees supplementaires:

noscf genere, a partir de &nofil, &no et &casdi les donnees pour
motra, molcost, casdet, casdi et locnats.

Il peut arriver qu'on veuille donner a un de ces programmes un donnee 
supplementaire. Par exemple, si on veut donner la donnee twvec a &dav
pour le 2eme calcul casdi (alors nmat >=2).
On ajoute a la fin du fichier :
DAV2
 twvec=0.01
FIN

explication : DAV pour &dav, 2 pour 2eme casdi
              la (les) donnee(s)
              FIN  (de l'ajout a &dav)
la ligne est ajoutee avant la fin de namelist "/"

Rem : si on veut ajouter un vecteur d'essai en donnee
DAV1 
 nessai=1
 /
 &ess file='ESSOUT0' 
FIN
le fichier de donnees devient apres rajout des 3 lignes :
 &dav
 ...
 ...
 nessai=1                  \
 /                          >  3 lignes supplementaires
 &ess file='ESSOUT0'       /
 /       (deja present, donc pas a rajouter)

