#include <stdio.h>                                                              
#define MAX_ROW 50                                                              
#define MAX_COL 100                                                             
#define SIZE_ARRAY (MAX_ROW * MAX_COL)                                          
                                                                                
main()                                                                          
{                                                                               
   int            result,             i,                 j            ;         
   int            rank,               inrank                          ;         
   long int       shape??( 2 ??),     inshape??( 2 ??)                ;         
   float          pressure  ??( MAX_ROW ??)??( MAX_COL ??)            ;         
   float          inpressure??( MAX_ROW ??)??( MAX_COL ??)            ;         
   float          xscales   ??( MAX_ROW ??)                           ;         
   float          inxscales ??( MAX_ROW ??)                           ;         
   float          yscales   ??( MAX_COL ??)                           ;         
   float          inyscales ??( MAX_COL ??)                           ;         
   float          maxpressure,        inmaxpressure                   ;         
   float          minpressure,        inminpressure                   ;         
   char           *datalabel,         indatalabel??( 256 ??)          ;         
   char           *dataunit,          indataunit ??( 256 ??)          ;         
   char           *datafmt,           indatafmt  ??( 256 ??)          ;         
   char           *dimlabels??( 2 ??),indimlabels??( 2 ??)??( 256 ??) ;         
   char           *dimunits ??( 2 ??),indimunits ??( 2 ??)??( 256 ??) ;         
   char           *dimfmts  ??( 2 ??),indimfmts  ??( 2 ??)??( 256 ??) ;         
   char           indummy   ??( 256 ??)                               ;         
                                                                                
   void           jesarray();                                                   
                                                                                
   rank               = 2           ;                                           
   shape??( 0 ??)     = MAX_ROW     ;                                           
   shape??( 1 ??)     = MAX_COL     ;                                           
                                                                                
   datalabel          = "Pressure 1";                                           
   dataunit           = "Pascals"   ;                                           
   datafmt            = "E15.9"     ;                                           
   dimlabels??( 0 ??) = "x"         ;                                           
   dimlabels??( 1 ??) = "y"         ;                                           
   dimunits ??( 0 ??)=                                                          
   dimunits ??( 1 ??) = "cm"        ;                                           
   dimfmts  ??( 0 ??) =                                                         
   dimfmts  ??( 1 ??) = "F10.2"     ;                                           
                                                                                
                                                                                
                                                                                
   result = dsgdims( "dd:input", &inrank, inshape, 2 )      ;                   
      printf( "return code from dsgdims   %i\n", result )  ;                    
      printf( "inrank, %i  \n",inrank  );                                       
      printf( "inshape, %i %i \n",inshape??(0??),inshape??(1??) );              
      printf( "rank is %i \n", inrank );                                        
      printf( "shape is %i by %i \n", inshape??( 0 ??), inshape??( 1 ??) );     
                                                                                
   result = dsgdata( "dd:input", 2, inshape, inpressure )   ;                   
      printf( "return code from dsgdata   %i\n", result )  ;                    
      for ( i=0; i<MAX_ROW-20; i++ )                                            
         {                                                                      
            for ( j=0; j<MAX_COL-20; j++ )                                      
                printf( " %10.5f ", inpressure??( i ??)??( j ??) );             
            printf( "\n" );                                                     
         }                                                                      
                                                                                
   result = dsgdast( indatalabel, indataunit, indatafmt, "indummy" )    ;       
      printf( "return code from dsgdast   %i\n", result )    ;                  
      printf( "label is %s  \n", indatalabel );                                 
      printf( "unit is %s  \n", indataunit );                                   
      printf( "fmt is %s  \n", indatafmt );                                     
                                                                                
   result = dsgdist( 1, indimlabels??( 0 ??),                                   
                               indimunits ??( 0 ??),                            
                               indimfmts  ??( 0 ??)             )    ;          
      printf( "return code from dsgdist    %i\n", result )    ;                 
                                                                                
   result = dsgdist( 2, indimlabels??( 1 ??),                                   
                               indimunits ??( 1 ??),                            
                               indimfmts  ??( 1 ??)             )    ;          
      printf( "return code from dsgdist    %i\n", result )    ;                 
      printf( "dim lablels are %s and %s \n",                                   
              indimlabels??( 0 ??), indimlabels??( 1 ??) );                     
      printf( "dim units are %s and %s \n",                                     
              indimunits??( 0 ??), indimunits??( 1 ??) );                       
      printf( "dim formats are %s and %s \n",                                   
              indimfmts??( 0 ??), indimfmts??( 1 ??) );                         
                                                                                
   result = dsgdisc( 1, inshape??( 0 ??), inxscales )          ;                
      printf( "return code from dsgdisc %i\n", result )    ;                    
      for ( i=0; i<MAX_ROW-10; i++ )                                            
         printf( "%10.5f", inxscales??( i ??) );                                
                                                                                
      printf( "\n" );                                                           
                                                                                
   result = dsgdisc( 2, inshape??( 1 ??), inyscales )          ;                
      printf( "return code from dsgdisc %i\n", result )    ;                    
      for ( i=0; i<MAX_COL-10; i++ )                                            
          printf( "%8.3f", inyscales??( i ??) );                                
                                                                                
      printf( "\n" );                                                           
                                                                                
   result = dsgmaxm( &inmaxpressure, &inminpressure )              ;            
      printf( "return code from dsgmaxm   %i\n", result )    ;                  
      printf("maxpressure is %10.3f \n", inmaxpressure);                        
      printf("minpressure is %10.3f \n", inminpressure);                        
                                                                                
}