Changeset 1300 in Sophya for trunk/SophyaExt

Timestamp:

Nov 8, 2000, 11:02:58 AM (25 years ago)

Author:

ansari

Message:

classe FITS_AutoReader

Location:

trunk/SophyaExt/FitsIOServer

Files:

• fitsbntbllineRW.cc (modified) (2 diffs)
• fitsbntbllineRW.h (modified) (1 diff)
• fitsfile.cc (modified) (23 diffs)
• fitsfile.h (modified) (8 diffs)
• fitslocalmap.cc (modified) (2 diffs)
• fitsntuple.cc (modified) (2 diffs)
• fitsspherehealpix.cc (modified) (2 diffs)
• fitsspherehealpix.h (modified) (1 diff)
• fitstarray.cc (modified) (2 diffs)
• fitsxntuple.cc (modified) (3 diffs)

trunk/SophyaExt/FitsIOServer/fitsbntbllineRW.cc

@@ -21 +21 @@
   InitNull();
   inFits_  = new FitsInFile (inputfile);
-  inFits_->ReadFInit(hdunum);
+  inFits_->ReadHeader(hdunum);
 
   //  if (!fn->IsFitsTable())
@@ -52 +52 @@
   for (k=0; k<nbcols;k++)
     {
-      //      char ss= fn->ColTypeFromFits(k);
-      char ss= inFits_->ColTypeFromFits(k);
-      if (ss == 'D') DfitsCol.push_back(k);
-      else if (ss == 'E') FfitsCol.push_back(k);
-      else if (ss == 'I') IfitsCol.push_back(k);
-      else if (ss == 'S') SfitsCol.push_back(k);
+       FitsFile::FitsDataType ss= inFits_->ColTypeFromFits(k);
+      if (ss == FitsFile::FitsDataType_double) DfitsCol.push_back(k);
+      else if (ss == FitsFile::FitsDataType_float) FfitsCol.push_back(k);
+      else if (ss == FitsFile::FitsDataType_int) IfitsCol.push_back(k);
+      else if (ss == FitsFile::FitsDataType_char) SfitsCol.push_back(k);
       else {
         cout << " FITS_XNTuple: colonne fits " << k << " type= " << ss << endl;

trunk/SophyaExt/FitsIOServer/fitsbntbllineRW.h

@@ -48 +48 @@
 
  FitsInFile* inFits_;
-  long nextLineToBeRead_;
-  BnTblLine ligne_;
+ long nextLineToBeRead_;
+ BnTblLine ligne_;
 };
 

trunk/SophyaExt/FitsIOServer/fitsfile.cc

@@ -104 +104 @@
 void FitsIOHandler::Read(FitsInFile& is, int hdunum) 
 {
-  is.ReadFInit(hdunum);
+  is.ReadHeader(hdunum);
   ReadFromFits(is);
 }
@@ -224 +224 @@
 void FitsInFile::InitNull()
 {
-  
-  bitpix_ = 0;
+  imageDataType_ = FitsDataType_NULL;
   naxis_ = 0;
   nbData_ = 0;
@@ -236 +235 @@
   dvl_.Clear();
   
-  
 }
 
@@ -268 +266 @@
     {
       typeOfExtension = FitsExtensionType_IMAGE;
-      int bitpix;
-      GetImageParameters (fileptr, bitpix, naxis, naxisn);
-      if(bitpix == DOUBLE_IMG) dataType = FitsDataType_double;
-      else
-        if(bitpix == FLOAT_IMG) dataType = FitsDataType_float;
-        else
-          if(bitpix == LONG_IMG || bitpix == SHORT_IMG ) dataType = FitsDataType_int;
-          else
-            if (bitpix == BYTE_IMG) dataType = FitsDataType_char;
-            else 
-            {
-              cout << " bitpix= " << bitpix << endl;
-              throw PException(" FitsFile::GetBlockType : unsupprted FITS data type");
-            }
-      
+      GetImageParameters (fileptr, dataType, naxis, naxisn);
     }
   else
@@ -289 +273 @@
         int nrows = 0;
         vector<string> noms; 
-        vector<char> types;   
+        vector<FitsDataType> types;   
         vector<int>  taille_des_chaines;
         GetBinTabParameters(fileptr, naxis, nrows, naxisn, noms, types, taille_des_chaines); 
@@ -302 +286 @@
           {
             typeOfExtension = FitsExtensionType_BINARY_TBL;
-            if(types[0] == 'D') dataType = FitsDataType_double;
-            else
-              if(types[0] == 'E') dataType = FitsDataType_float;
-              else
-                if(types[0] == 'I' ) dataType = FitsDataType_int;
-                else
-                  if(types[0] == 'S' ) dataType = FitsDataType_char;
-                  else 
-                    {
-                      cout << " types[0]= " << types[0] << endl;
-                      throw PException(" FitsFile::GetBlockType : unsupprted FITS data type");
-            }
+            dataType = types[0];
           }
       }
@@ -327 +300 @@
 
 
-void FitsInFile::ReadFInit(int hdunum) 
-{
-  InitNull();
+void FitsInFile::ReadHeader(int hdunum) 
+{
+  //  InitNull();
   int status = 0;
 
@@ -352 +325 @@
 
 
-void FitsInFile::GetImageParameters (fitsfile* fileptr,int& bitpix,int& naxis,vector<int>& naxisn)
+void FitsInFile::GetImageParameters (fitsfile* fileptr,FitsDataType& dataType,int& naxis,vector<int>& naxisn)
 {
   int hdunum=0;
@@ -359 +332 @@
 
   // bits per pixels
+  int bitpix=0;
   fits_read_key(fileptr,TINT,"BITPIX",&bitpix,NULL,&status);
   if( status )  printerror( status );
+  if(bitpix == DOUBLE_IMG) dataType = FitsDataType_double;
+  else if(bitpix == FLOAT_IMG) dataType = FitsDataType_float;
+  else if(bitpix == LONG_IMG || bitpix == SHORT_IMG ) dataType = FitsDataType_int;
+  else if (bitpix == BYTE_IMG) dataType = FitsDataType_char;
+  else 
+    {
+      cout << " bitpix= " << bitpix << endl;
+      throw PException(" FitsFile::GetImageParameters : unsupported FITS data type");
+    }
 
   // number of dimensions in the FITS array
@@ -401 +384 @@
 void FitsInFile::getHeader()
 {
+  // si hdunum_ > 1 lit le header correspondant
+  // si hdunum_ = 1 se positionne au (et lit le) premier header qui 
+  //                contient reellement un objet 
   int status=0;
   if (hdunum_ < 1)  throw PException(" attempt to read hdunum < 1");
+   InitNull();
   if (hdunum_ == 1)
     {
@@ -412 +399 @@
         {
           hdutype_ = IMAGE_HDU;
-          GetImageParameters (fptr_, bitpix_, naxis_, naxisn_);
+          GetImageParameters (fptr_, imageDataType_, naxis_, naxisn_);
           nbData_ =  1;
           int k;
@@ -436 +423 @@
       if(hdutype_ == IMAGE_HDU) 
         {
-          GetImageParameters (fptr_, bitpix_, naxis_, naxisn_);
+          GetImageParameters (fptr_, imageDataType_, naxis_, naxisn_);
           nbData_ =  1;
           int k;
@@ -509 +496 @@
   */
 
-char FitsInFile::ColTypeFromFits(int nocol) const
+FitsFile::FitsDataType FitsInFile::ColTypeFromFits(int nocol) const
 {
   if(hdutype_ != ASCII_TBL && hdutype_ != BINARY_TBL)
@@ -548 +535 @@
   for (k=0; k<=nocol; k++) 
     {
-      if (types_[k] == 'S') index++;
+      if (types_[k] == FitsDataType_char) index++;
     }
   return  taille_des_chaines_[index];
@@ -578 +565 @@
       switch (types_[ncol])
         { 
-        case 'D' :
+        case FitsDataType_double :
           fits_read_col(fptr_,TDOUBLE,ncol+1,NoLine+1,1,1,&dnull,&ddata[dcount++],&anull,&status);
             break;
-        case 'E' :
+        case FitsDataType_float :
           fits_read_col(fptr_,TFLOAT,ncol+1,NoLine+1,1,1,&fnull,&fdata[fcount++],&anull,&status);
           break;
-        case 'I' :
+        case  FitsDataType_int :
           fits_read_col(fptr_,TINT,ncol+1,NoLine+1,1,1,&inull,&idata[icount++],
                         &anull,&status);
           break;
-        case 'S' :
+        case FitsDataType_char :
           fits_read_col(fptr_,TSTRING,ncol+1,NoLine+1,1,1,cnull,&cdata[ccount++],&anull,&status);
           break;
@@ -621 +608 @@
       switch (types_[ncol])
         { 
-        case 'D' :
+        case FitsDataType_double :
           fits_read_col(fptr_,TDOUBLE,ncol+1,NoLine+1,1,1,&dnull,&ligne.ddata_[dcount++],&anull,&status);
-            break;
-        case 'E' :
+          break;
+        case FitsDataType_float :
           fits_read_col(fptr_,TFLOAT,ncol+1,NoLine+1,1,1,&fnull,&ligne.fdata_[fcount++],&anull,&status);
           break;
-        case 'I' :
+        case FitsDataType_int :
           fits_read_col(fptr_,TINT,ncol+1,NoLine+1,1,1,&inull,&ligne.idata_[icount++],
                         &anull,&status);
           break;
-        case 'S' :
+        case FitsDataType_char :
           char* chaine = new char[taille_des_chaines_[ccount]];
           fits_read_col(fptr_,TSTRING,ncol+1,NoLine+1,1,1,cnull,&chaine,&anull,&status);
@@ -789 +776 @@
       long repeat,width;
       fits_get_coltype(fptr_, NoCol+1,&DTYPE,&repeat,&width,&status);
-      if( DTYPE != TSTRING ) 
-        {
-          throw IOExc("FitsFile::GetBinTabFCol, tentative de lecture non float");
+      if( DTYPE != TSTRING &&  DTYPE != TBYTE) 
+        {
+          throw IOExc("FitsFile::GetBinTabFCol, tentative de lecture non string");
         }
       long nels=nentries;
@@ -821 +808 @@
     {
   
-      if(bitpix_ != DOUBLE_IMG) 
+      if(imageDataType_ != FitsDataType_double) 
         {
           cout << " The data type on fits file is not double...";
@@ -855 +842 @@
   if(hdutype_ == IMAGE_HDU) 
     {
-      if(bitpix_ != FLOAT_IMG) 
+      if(imageDataType_ != FitsDataType_float) 
         {
           cout << " The data type on fits file is not float ";
@@ -888 +875 @@
   if(hdutype_ == IMAGE_HDU) 
     {
-      if(bitpix_ != LONG_IMG) 
+      if(imageDataType_ != FitsDataType_int) 
         {
           cout << " The data type on fits file is not int ";
@@ -916 +903 @@
                                   vector<int>& repeat,
                                   vector<string>& noms, 
-                                  vector<char>& types,   
+                                  vector<FitsDataType>& types,   
                                   vector<int>&  taille_des_chaines)
 {
@@ -961 +948 @@
   if( status ) printerror( status,"erreur lecture des noms de colonne");
   int rept=0;
-  for(ii = 0; ii <  nbcols; ii++)
-    {
-      int DTYPE;
-      long width;
-      long repete = 0;
-      fits_get_coltype(fileptr,ii+1,&DTYPE,&repete,&width,&status);
-      if( status ) printerror( status,"erreur lecture type de colonne");
-      rept = repete;
-      noms.push_back(string(ttype[ii]));
-      switch (DTYPE)
-        {
-        case TDOUBLE :
-          types.push_back('D');
-          break;
-        case TFLOAT :
-          types.push_back('E');
-          break;
-        case TLONG :
-          types.push_back('I');
-          break;
-        case TINT :
-          types.push_back('I');
-          break;
-        case TSHORT :
-          types.push_back('I');
-          break;
-        case TSTRING :
-          types.push_back('S');
-          taille_des_chaines.push_back(width);
-          rept/=width;
-          break;
-        default :
-          cout << " field " << ii+1 << " DTYPE= " << DTYPE << endl;
-          throw IOExc("FitsFile:: unknown type of field");
-        }
-      repeat.push_back(rept);
+  if(hdutype == ASCII_TBL)
+    {
+      for(ii = 0; ii <  nbcols; ii++)
+        {
+          int DTYPE;
+          long width;
+          long repete = 0;
+          fits_get_coltype(fileptr,ii+1,&DTYPE,&repete,&width,&status);
+          if( status ) printerror( status,"erreur lecture type de colonne");
+          rept = repete;
+          noms.push_back(string(ttype[ii]));
+          switch (DTYPE)
+            {
+            case TDOUBLE :
+              types.push_back(FitsDataType_double);
+              break;
+            case TFLOAT :
+              types.push_back(FitsDataType_float);
+              break;
+            case TLONG :
+              types.push_back(FitsDataType_int);
+              break;
+            case TSHORT :
+              types.push_back(FitsDataType_int);
+              break;
+            case TSTRING :
+              types.push_back(FitsDataType_char);
+              taille_des_chaines.push_back(width);
+              rept/=width;
+              break;
+            default :
+              cout << " field " << ii+1 << " DTYPE= " << DTYPE << endl;
+              throw IOExc("FitsFile::GetBinTabParameters, unsupported data type of field, for ASCII table");
+            }
+          repeat.push_back(rept);
+        }
+    }
+  else 
+    {
+      for(ii = 0; ii <  nbcols; ii++)
+        {
+          int DTYPE;
+          long width;
+          long repete = 0;
+          fits_get_coltype(fileptr,ii+1,&DTYPE,&repete,&width,&status);
+          if( status ) printerror( status,"erreur lecture type de colonne");
+          rept = repete;
+          noms.push_back(string(ttype[ii]));
+          switch (DTYPE)
+            {
+            case TDOUBLE :
+              types.push_back(FitsDataType_double);
+              break;
+            case TFLOAT :
+              types.push_back(FitsDataType_float);
+              break;
+            case TLONG :
+              types.push_back(FitsDataType_int);
+              break;
+            case TINT :
+              types.push_back(FitsDataType_int);
+              break;
+            case TSHORT :
+              types.push_back(FitsDataType_int);
+              break;
+            case TSTRING :
+              types.push_back(FitsDataType_char);
+              taille_des_chaines.push_back(width);
+              rept/=width;
+              break;
+            case TBYTE :
+              types.push_back(FitsDataType_char);
+              taille_des_chaines.push_back(width);
+              rept/=width;
+              break;
+            default :
+              cout << " field " << ii+1 << " DTYPE= " << DTYPE << endl;
+              throw IOExc("FitsFile::GetBinTabParameters, unsupported data type of field, for BINTABLE");
+            }
+          repeat.push_back(rept);
+        }
     }
   for (ii=0; ii < nbcols; ii++) delete [] ttype[ii];
@@ -1292 +1324 @@
 \param <taille_des_chaines> vector containing the number of characters of  data  for each char* typed column, with order of appearance in 'fieldType'
 */
-void FitsOutFile::makeHeaderBntblOnFits( string fieldType, vector<string> Noms, int nentries, int tfields, DVList* ptr_dvl, string extname, vector<int> taille_des_chaines) 
+void FitsOutFile::makeHeaderBntblOnFits(string fieldType, vector<string> Noms, int nentries, int tfields, DVList* ptr_dvl, string extname, vector<int> taille_des_chaines) 
 { 
   int k;
   int status = 0;
   long nrows;
-  // verifications de coherence
+  // verifications de coherences
 
   if (fieldType.length() != tfields)

trunk/SophyaExt/FitsIOServer/fitsfile.h

@@ -60 +60 @@
    };
    enum FitsDataType {
+     FitsDataType_NULL,
      FitsDataType_double,
      FitsDataType_float,
      FitsDataType_int,
      FitsDataType_char,
-     FitsDataType_ASCII
+     FitsDataType_ASCII,
    };
 
@@ -72 +73 @@
    inline  int   statusF() const { return fits_status_;}
    inline void firstImageOnPrimaryHeader(bool choice) {imageOnPrimary_=choice;}
+   inline int currentHeaderIndex() {return hdunum_;}
  
 
@@ -81 +83 @@
    inline void  InitNull() {fptr_ = NULL; hdutype_= 0; hdunum_ = 0; 
    fits_status_ = 0; imageOnPrimary_ = true;}
-
    fitsfile *fptr_;     /**<  pointer to the FITS file, defined in fitsio.h */
    int hdutype_;        /**<  image or bintable ? */
@@ -102 +103 @@
    static int  NbBlocks(char flnm[]);
    static void GetBlockType(char flnm[], int hdunum, FitsExtensionType& typeOfExtension, int& naxis, vector<int>& naxisn, FitsDataType& dataType, DVList& dvl  );
-   void        ReadFInit(int hdunum);
+   void        ReadHeader(int hdunum);
  
   /*! \return a reference on a DVList containing the keywords from FITS file */
@@ -139 +140 @@
 inline int nbOfImageData() const { return nbData_; }
 
+/*! \return data type of the current IMAGE extension */
+inline FitsFile::FitsDataType ImageType() const {return imageDataType_;}
 
 
@@ -155 +158 @@
                                   vector<int>& repeat,
                                   vector<string>& noms, 
-                                  vector<char>& types,   
+                                  vector<FitsDataType>& types,   
                                   vector<int>&  taille_des_chaines);
- char   ColTypeFromFits(int nocol) const;
+ FitsDataType   ColTypeFromFits(int nocol) const;
  string ColNameFromFits(int nocol) const;
  int    ColStringLengthFromFits(int nocol) const;
@@ -188 +191 @@
 void getHeader();
 static void KeywordsIntoDVList(fitsfile* fileptr, DVList& dvl, int hdunum); 
-static  void GetImageParameters (fitsfile* fileptr,int& bitpix,int& naxis,vector<int>& naxisn);
-
- int bitpix_;          /**< fits-Image parameter */
+static  void GetImageParameters (fitsfile* fileptr,FitsDataType& dataType,int& naxis,vector<int>& naxisn);
+
+ FitsDataType imageDataType_;          /**< fits-Image parameter (bitpix)*/
  int naxis_;           /**< fits-Image parameter */
  vector<int> naxisn_;  /**< fits-Image parameters : sizes of dimensions */
@@ -198 +201 @@
  int nbcols_;          /**< Bintable parameter */
  vector<string> noms_; /**< Bintable parameter: column names */
- vector<char> types_;  /**< Bintable parameters: types of columns (D: double, E: float, I: integers,  A: char*) */
+ vector<FitsDataType> types_;  /**< Bintable parameters: types of columns (D: double, E: float, I: integers,  A: char*) */
  DVList dvl_;          /**< DVList for transferring keywords */
  vector<int>  taille_des_chaines_; /**< Bintable parameters:   length of the char* variables */

trunk/SophyaExt/FitsIOServer/fitslocalmap.cc

@@ -159 +159 @@
     }
   dvl["Content"]= "LocalMap";
+  dvl.SetComment("Content", "name of SOPHYA object");
   // On ecrit les dataBlocks 
   vector<string> Noms;   
@@ -182 +183 @@
 #pragma define_template FITS_LocalMap<r_8>
 #pragma define_template FITS_LocalMap<r_4>
+#pragma define_template FITS_LocalMap<int_4>
 #endif
 #if defined(ANSI_TEMPLATES) || defined(GNU_TEMPLATES)
 template class FITS_LocalMap<r_8>;
 template class FITS_LocalMap<r_4>;
+template class FITS_LocalMap<int_4>;
 #endif

trunk/SophyaExt/FitsIOServer/fitsntuple.cc

@@ -89 +89 @@
   for (k=0; k<nbcols;k++)
     {
-      char ss= is.ColTypeFromFits(k);
+      FitsFile::FitsDataType ss= is.ColTypeFromFits(k);
       string type;
-      if (ss != 'E') 
+      if (ss !=  FitsFile::FitsDataType_float) 
         {
-          if (ss == 'D') type= string("double");
+          if (ss == FitsFile::FitsDataType_double) type= string("double");
           else
-            if (ss == 'I') type= string("integer");
+            if (ss == FitsFile::FitsDataType_int) type= string("integer");
             else
-              if (ss == 'A') type = string("char*");
+              if (ss == FitsFile::FitsDataType_char) type = string("char*");
               else
                 type = string("unknown");
@@ -159 +159 @@
   // get names and values from the join DVList object
   DVList dvl= dobj_->Info();
+  dvl["Content"]= "NTuple";
+  dvl.SetComment("Content", "name of SOPHYA object");
   // extension name
   string extname("NTuple_Binary_tbl"); 

trunk/SophyaExt/FitsIOServer/fitsspherehealpix.cc

@@ -97 +97 @@
   dvl.SetComment("LASTPIX", "Last pixel # (0 based)");
   dvl["Content"]= "SphereHEALPix";
-   
+  dvl.SetComment("Content", "name of SOPHYA object");
+  
   // On ecrit les dataBlocks 
   vector<string> Noms;   
@@ -273 +274 @@
 #pragma define_template FITS_SphereHEALPix<r_8>
 #pragma define_template FITS_SphereHEALPix<r_4>
+#pragma define_template FITS_SphereHEALPix<int_4>
 #endif
 #if defined(ANSI_TEMPLATES) || defined(GNU_TEMPLATES)
 template class FITS_SphereHEALPix<r_8>;
 template class FITS_SphereHEALPix<r_4>;
+template class FITS_SphereHEALPix<int_4>;
 #endif

trunk/SophyaExt/FitsIOServer/fitsspherehealpix.h

@@ -8 +8 @@
 #include "spherehealpix.h"
 #include "anydataobj.h"
-#include "ppersist.h"
 #include "fitsfile.h"
 

trunk/SophyaExt/FitsIOServer/fitstarray.cc

@@ -134 +134 @@
     }
   cout << " nombre total d'elements a copier " << nbels << endl;
-  os.makeHeaderImageOnFits(type, nbdim, naxisn, &dobj_->Info());
+  DVList dvl( dobj_->Info() );  
+  dvl["Content"]= "TArray";
+  dvl.SetComment("Content", "name of SOPHYA object");
+  os.makeHeaderImageOnFits(type, nbdim, naxisn, &dvl);
   if (!dobj_->IsPacked()) 
     {
@@ -156 +159 @@
 #pragma define_template FITS_TArray<r_8>
 #pragma define_template FITS_TArray<r_4>
+#pragma define_template FITS_TArray<int_4>
 #endif
 #if defined(ANSI_TEMPLATES) || defined(GNU_TEMPLATES)
 template class FITS_TArray<r_8>;
 template class FITS_TArray<r_4>;
+template class FITS_TArray<int_4>;
 #endif

trunk/SophyaExt/FitsIOServer/fitsxntuple.cc

@@ -80 +80 @@
   for (k=0; k<nbcols;k++)
     {
-      char ss= is.ColTypeFromFits(k);
-      if (ss == 'D') DfitsCol.push_back(k);
-      else if (ss == 'E') FfitsCol.push_back(k);
-      else if (ss == 'I') IfitsCol.push_back(k);
-      else if (ss == 'S') SfitsCol.push_back(k);
+      FitsFile::FitsDataType ss= is.ColTypeFromFits(k);
+      if (ss == FitsFile::FitsDataType_double) DfitsCol.push_back(k);
+      else if (ss == FitsFile::FitsDataType_float) FfitsCol.push_back(k);
+      else if (ss == FitsFile::FitsDataType_int) IfitsCol.push_back(k);
+      else if (ss == FitsFile::FitsDataType_char) SfitsCol.push_back(k);
       else {
         cout << " FITS_XNTuple: colonne fits " << k << " type= " << ss << endl;
@@ -126 +126 @@
       if (ownobj_)
         {
-      (*dobj_)= XNTuple(DfitsCol.size(), FfitsCol.size(), IfitsCol.size(), SfitsCol.size(),ColName);
+          (*dobj_)= XNTuple(DfitsCol.size(), FfitsCol.size(), IfitsCol.size(), SfitsCol.size(),ColName);
         }
       else
         {
-        if (DfitsCol.size() != dobj_->NDVar() || FfitsCol.size() != dobj_->NFVar() ||  IfitsCol.size() != dobj_->NIVar() || SfitsCol.size() != dobj_->NSVar())
-
-          throw  SzMismatchError("FITS_XNTuple : structure incorrecte du ntuple");
+          if (DfitsCol.size() != dobj_->NDVar() || FfitsCol.size() != dobj_->NFVar() ||  IfitsCol.size() != dobj_->NIVar() || SfitsCol.size() != dobj_->NSVar())
+            {
+              cout << " WARNING : FITS_XNTuple : XNTuple reconfigured " << endl;
+              (*dobj_)= XNTuple(DfitsCol.size(), FfitsCol.size(), IfitsCol.size(), SfitsCol.size(),ColName);
+            }
         }
     }
@@ -227 +229 @@
   vector<int> StringSizes(dobj_->NSVar());
   for (k=0; k< StringSizes.size(); k++) StringSizes[k]=dobj_->mStrSz;
+  dvl["Content"]= "XNTuple";
+  dvl.SetComment("Content", "name of SOPHYA object");
   os.makeHeaderBntblOnFits(types, Noms, nrows, ncols, &dvl, extname,StringSizes);