C COVSRT    SOURCE    PV        08/09/11    21:15:27     6150
      subroutine covsrt(covar,ncvm,ma)
c
c etant donne une matrice de covariance covar pour un jeu de mfit parmi ma
c parametres, et leur numerotation lista, reordonne ; met a zero les termes
c lies aux parametres fixes ; ncvm est la dimension de covar
c
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
      real*8 covar(ncvm,ncvm)

* met a zero les termes sous la diagonale
      do 12 j=1,ma-1
         do 11 i = j+1,ma
            covar(i,j) = 0.
 11      continue
 12   continue

      do 14 i=1,ma-1
         do 13 j = i+1,ma
            if(j.gt.i) then
              covar(j,i)=covar(i,j)
            else
              covar(i,j)=covar(i,j)
            endif
 13      continue
 14   continue

      swap = covar(1,1)
      do 15 j=1,ma
         covar(1,j) = covar(j,j)
         covar(j,j) = 0.
 15   continue

      covar(1,1) = swap
      do 16 j=2,ma
        covar(j,j)=covar(1,j)
 16   continue

      do 18 j=2,ma
       do 17 i=1,j-1
          covar(i,j) = covar(j,i)
 17    continue
 18   continue

      return
      end