#include <string.h>
#include <stdio.h>
#include "fitsio.h"

int main() {
    fitsfile *fptr;
    int status = 0, i;

    fits_open_file(&fptr, "WFPC2ASSNu5780205bx.fits", READONLY, &status);
    long naxes[2];
    fits_get_img_size(fptr, 3, naxes, &status);

    fitsfile *ofptr;
    fits_create_file(&ofptr, "o_nasa.fits", &status);
    fits_copy_header(fptr, ofptr, &status);

    long fp[2] = {1, 1};
    long nelements = naxes[0];

    float arr[nelements];
    for (i = 0; i < nelements; i++) {
        arr[i] = 100;
    }

    int ii, jj, kk;
    for (ii = 1; ii <= naxes[0]; ii++) {
        for (jj = 1; jj <= naxes[1]; jj++) {
            fits_write_pix(ofptr, TFLOAT, fp, nelements, arr, &status); // this is not working
        }
    }
    fits_close_file(ofptr, &status);
    fits_close_file(fptr, &status);

    return status;
}