A complete and efficient use of the GILDAS format from a program is beyond the scope of this chapter : refer to the programmer's guide for that. There is however a simple and easy-to-use subroutine that accepts a Fortran array and writes a GILDAS image. The price for simplicity is that only minimal header information is written: many useful stuff like World coordinates, Blanking value, and so on, will be missing. With this restriction, you can simply write a Fortran array in GILDAS format using the routine GDF_IMAGE.
SUBROUTINE GDF_IMAGE(NAME,NX,NY,NZ,NT,Z,ERROR) Where the arguments are NAME C*(*) the GILDAS file name NX I the number of pixels in the first dimension NY I the number of pixels in the second dimension NZ I the number of pixels in the third dimension NT I the number of pixels in the fourth dimension Z R*4 the array of dimensions NZ,NY,NZ,NT ERROR L a logical error flag (Output)
If you need more complete information, refer to the programmers guide. You may also later modify the file header using command HEADER in GRAPHIC, or through the SIC variables defined by commands DEFINE HEADER or DEFINE IMAGE.