proj/lev0/apps/build_lev1_aia.c

/*-----------------------------------------------------------------------------
 * cvs/JSOC/proj/lev1/apps/build_lev1_aia.c
 *-----------------------------------------------------------------------------
 *
 * This is a module that runs with DRMS and processes lev0
 * filtergrams to lev1.
 * It is scheduled by build_lev1_mgr either by qsub to the cluster
 * or by fork to run on the local machine. It is called with 
 * 12 (NUMRECLEV1 defined in lev0lev1.h) lev0 images at a time.
 * See the -h nice_intro() for details of the calling options.
 */ 

#include <jsoc_main.h>
#include <cmdparams.h>
#include <drms.h>
#include <drms_names.h>
#include <stdio.h>
#include <stdlib.h>
#include <ctype.h>
#include <strings.h>
#include <sys/types.h>
#include <sys/time.h>
#include <sys/stat.h>   //for umask(2)
#include <unistd.h> //for alarm(2) among other things...
#include <printk.h>
#include <astro.h>
#include <fresize.h>
#include <gapfill.h>
//#include </home/jsoc/include/fftw3.h> 
#include "fftw3.h" 
#include "imgdecode.h"
#include "lev0lev1.h"
#include "quallev1.h"
#include "limb_fit.h"
#include "cosmic_ray.h"
#include "get_pointing_info.c"

//default in and out data series
//#define LEV0SERIESNAMEHMI "hmi.lev0e"
//#define LEV0SERIESNAMEAIA "aia.lev0e"
#define LEV0SERIESNAMEHMI "su_production.lev0f_hmi"
#define LEV0SERIESNAMEAIA "su_production.lev0f_aia"
//#define LEV1SERIESNAMEHMI "su_production.hmi_lev1e"   //temp test case
#define LEV1SERIESNAMEHMI "hmi.lev1"
#define LEV1SERIESNAMEAIA "su_production.aia_lev1e" //temp test case
//#define DSFFNAME "su_richard.flatfield"       //temp test case
//#define DSFFNAMEHMI "su_production.hmi_flatfield" //temp test case
#define DSFFNAMEHMI "hmi.flatfield"
//#define DSFFNAMEAIA "su_production.aia_flatfield" //temp test case
//#define DSFFNAMEAIA "aia_test.flatfield"  //temp test case
#define DSFFNAMEAIA "aia.flatfield"

#define LEV1LOG_BASEDIR "/usr/local/logs/lev1"
#define H1LOGFILE "/usr/local/logs/lev1/build_lev1.%s.log"
#define NUMTIMERS 8     //number of seperate timers avail
#define NOTSPECIFIED "***NOTSPECIFIED***"
#define LOGTEST 0
#define CAL_HCFTID 17       //image is cal mode 
#define STOP_FILE "/usr/local/logs/lev1/build_mgr_stop_aia"

#define kDpath    "dpath"
#define kDpathDef "/home/jsoc/cvs/Development/JSOC"

int compare_rptr(const void *a, const void *b);
static TIME SDO_to_DRMS_time(int sdo_s, int sdo_ss);

// List of default parameter values. 
ModuleArgs_t module_args[] = { 
  {ARG_STRING, "instru", NOTSPECIFIED, "instrument. either hmi or aia"},
  {ARG_STRING, "mode", NOTSPECIFIED, "either recnum or fsn"},
  {ARG_STRING, "dsin", NOTSPECIFIED, "dataset of lev0 filtergrams"},
  {ARG_STRING, "dsout", NOTSPECIFIED, "dataset of lev1 output"},
  {ARG_STRING, "dsff", NOTSPECIFIED, "dataset of darks flats bads"},
  {ARG_STRING, "dsaiabad", NOTSPECIFIED, "dataset of AIA bad pixels"},
  {ARG_STRING, "logfile", NOTSPECIFIED, "optional log file name. Will create one if not given"},
  {ARG_INTS, "brec", "0", "first lev0 rec# to process. 0=error must be given by build_lev1_mgr"},
  {ARG_INTS, "erec", "0", "last lev0 rec# to process. 0=error must be given by build_lev1_mgr"},
  {ARG_INTS, "bfsn", "0", "first lev0 fsn# to process. 0=error must be given by build_lev1_mgr"},
  {ARG_INTS, "efsn", "0", "last lev0 fsn# to process. 0=error must be given by build_lev1_mgr"},
  {ARG_INTS, "quicklook", "1", "1=quick look, 0 = definitive mode"},
  {ARG_FLAG, "v", "0", "verbose flag"},
  {ARG_FLAG, "h", "0", "help flag"},
  {ARG_FLAG, "r", "0", "restart flag"},
  {ARG_STRING, kDpath, kDpathDef, "path to the source code tree (the JSOC root directory)"},
  {ARG_END}
};

CmdParams_t cmdparams;
// Module name presented to DRMS. 
char *module_name = "build_lev1_aia";

FILE *h1logfp;      // fp for h1 ouput log for this run 
static struct stat stbuf;
//static IMG Image0, Image1;
//static IMG *Img0 = &Image0;
//static IMG *Img1 = &Image1;
static LEV0LEV1 lev0lev1;
static LEV0LEV1 *l0l1 = &lev0lev1;
//static CCSDS_Packet_t *Hk;
static DRMS_Record_t *rs;
static DRMS_Record_t *rs0, *rs1, *rsff, *rsbad_pix, *rec_bad_aia, *rt, *rresp;
static DRMS_Record_t *rptr;
static DRMS_Segment_t *segment;
static DRMS_Segment_t *segmentff;
static DRMS_Segment_t *darkseg;
static DRMS_Segment_t *badseg;
static DRMS_Segment_t *badoutpixseg;
static DRMS_Segment_t *spikeseg;
static DRMS_Array_t *segArray;
static DRMS_RecordSet_t *rset0, *rset1, *rsetff, *rsbad_aia, *rs_t=NULL,
                        *rs_resp=NULL;
static DRMS_Array_t *Array0;
static DRMS_Array_t *Arrayff;
static DRMS_Array_t *ArrayDark;
static DRMS_Array_t *ArrayBad;
static DRMS_Array_t *ArraySpike;
static TIME sdo_epoch;
static PTINFO *ptinfo = NULL;
static PTINFO ptdata;
static char bld_vers[16];
static char datestr[32];
static char open_dsname[256];
static char dsffname[128];
static char open_aiabad_dsname[256];
static char path[DRMS_MAXPATHLEN], bad_pix_path[DRMS_MAXPATHLEN];
static char bad_aia_path[DRMS_MAXPATHLEN];
static char rs1_path[DRMS_MAXPATHLEN];
static struct timeval first[NUMTIMERS], second[NUMTIMERS];
static char *orbseries = "sdo.fds_orbit_vectors";
//static char *orbseries = "sdo_ground.fds_orbit_vectors";

static int nspikes, respike, fid, aiftsid, *oldvalues, *spikelocs, *newvalues;
static int hcftid, aiagp6;
static short aifcps;
double aiascale = 1.0;

IORBIT_Info_t *IOinfo = NULL;
IORBIT_Info_t IOdata;
LIBASTRO_Error_t IOstatus;
unsigned int fsnarray[NUMRECLEV1];
unsigned int fsnx = 0;
//short data1[MAXPIXELS];
//int data1[MAXPIXELS];
float data1[MAXPIXELS];     //floats for HMI
float ftmp;
int data1A[MAXPIXELS];      //ints for AIA
int array_cosmic[16777216]; //4096x4096
double tgttimes[NUMRECLEV1];

long long brec, erec, bfsn, efsn; //begin and end lev0 rec/fsn. must be same data type
long long bnumx, enumx;       //has either the brec/erec of bfsn/efsn pair accord to mode
int verbose;
int hmiaiaflg = 0;      //0=hmi, 1=aia
int modeflg = 0;        //0=fsn, 1=recnum
int imagecnt = 0;       // num of images since last commit 
int restartflg = 0;     // set when build_lev1 is called for a restart
int abortflg = 0;
int sigalrmflg = 0;             // set on signal so prog will know 
int ignoresigalrmflg = 0;       // set after a close_image()
int quicklook;
//global quality flags
int flatmiss[NUMRECLEV1];
int orbmiss[NUMRECLEV1];
int asdmiss[NUMRECLEV1];
int mpdmiss[NUMRECLEV1];
int limbmiss[NUMRECLEV1];
int noimage[NUMRECLEV1];
int missflg[NUMRECLEV1];

char logname[128];
char argdsin[128], argdsout[128], arglogfile[128], arginstru[80];
char argbx[80], argex[80], argquick[80], argmode[80], argdsaiabad[80];
char timetag[32];
char tlmseriesname[128];    // e.g. hmi.tlm
char lev0seriesname[128];   // e.g. hmi.lev0
char *username;         // from getenv("USER") 
const char *logfile;            // optional log name passed in 
const char *instru;     // instument. hmi or aia 
const char *mode;       // given mode. recnum or fsn
const char *dsin;       // lev0 input dataset
const char *dsout;      // lev1 output dataset
const char *dsff;               // flat field dataset
const char *dsaiabad;           // AIA bad pixel dataset (kludge around bad ff)

typedef struct {
  float rsun_lf;
  float x0_lf;
  float y0_lf;
} LIMB_SOMETHING;



int get_nspikes() { return nspikes; }
int get_respike(void) { return respike; }
int *get_spikelocs() { return spikelocs; }
int *get_oldvalues() { return oldvalues; }
int *get_newvalues() { return newvalues; }
void set_nspikes(int new_nspikes) { nspikes = new_nspikes; }
void set_spikelocs(int *new_spikelocs) { spikelocs = new_spikelocs; }
void set_oldvalues(int *new_oldvalues) { oldvalues = new_oldvalues; }
void set_newvalues(int *new_newvalues) { newvalues = new_newvalues; }

//Set the QUALITY keyword for lev1
void do_quallev1(DRMS_Record_t *rs0, DRMS_Record_t *rs1, int inx, unsigned int fsn)
{
  int quallev1 = 0;
  int rstatus, quallev0;;
  char *pchar;

  quallev1 = missflg[inx];
  if(flatmiss[inx]) quallev1 = quallev1 | Q_NOFLAT;
  if(orbmiss[inx]) quallev1 = quallev1 | Q_NOORB;
  if(limbmiss[inx]) quallev1 = quallev1 | Q_NOLIMB;
  if(asdmiss[inx]) {
    quallev1 = quallev1 | Q_NOASD;
    drms_setkey_string(rs1, "ASD_REC", DRMS_MISSING_STRING);
  }
  if(mpdmiss[inx]) quallev1 = quallev1 | Q_NOMPD;
  if(noimage[inx]) quallev1 = quallev1 | Q_MISSALL;
  if(ptinfo) {
    ptdata = ptinfo[inx];
    if(strcmp(ptdata.acs_mode, "SCIENCE")) { //not in science pointing mode
      quallev1 = quallev1 | Q_NOACS_SCI;
    }
    if(strcmp(ptdata.acs_eclp, "NO")) {      //in eclipse
      quallev1 = quallev1 | Q_ACS_ECLP;
    }
    if(strcmp(ptdata.acs_sunp, "YES")) {     //not in sun presence flag
      quallev1 = quallev1 | Q_ACS_SUNP;
    }
    if(strcmp(ptdata.acs_safe, "NO")) {     //safemode set
      quallev1 = quallev1 | Q_ACS_SAFE;
    }
  }
  pchar = drms_getkey_string(rs0, "IMG_TYPE", &rstatus);
  if(rstatus) {
    printk("ERROR: in drms_getkey_string(IMG_TYPE) fsn=%u\n", fsn);
  }
  else {
    if(strcmp(pchar, "LIGHT")) {    //dark image
      quallev1 = quallev1 | Q_IMG_TYPE;
    }
  }
  if(hmiaiaflg) {         //aia
    pchar = drms_getkey_string(rs0, "AISTATE", &rstatus);
    if(aiftsid >= 0xc000) quallev1 = quallev1 | Q_CAL_IMG;
    if((aifcps <= -20) ||(aifcps >=100)) quallev1 = quallev1 | Q_AIA_FOOR;
    if(aiagp6 != 0) quallev1 = quallev1 | Q_AIA_REGF;
  }
  else {              //hmi
    pchar = drms_getkey_string(rs0, "HWLTNSET", &rstatus);
    if((fid >= 1) && (fid <=9999)) {
      quallev1 = quallev1 | Q_CAL_IMG;      //cal image
    } 
    if(hcftid == CAL_HCFTID) quallev1 = quallev1 | Q_CALM_IMG; //hmi cal mode
  }
  if(rstatus) {
    printk("ERROR: in drms_getkey_string(HWLTNSET or AISTATE) fsn=%u\n", fsn);
  }
  else {
    if(!strcmp(pchar, "OPEN")) {    //ISS loop open
      quallev1 = quallev1 | Q_LOOP_OPEN;
    }
  }

  if(quicklook) {
    quallev1 = quallev1 | Q_NRT;
  }

  //for camera minimum value error
  quallev0 = drms_getkey_int(rs0, "QUALITY", &rstatus);
  if ( (quallev0 & 0x8000) != 0) {
   quallev1 = quallev1 | Q_CAM_ANOM1;
  }

  drms_setkey_int(rs1, "QUALITY", quallev1);
  drms_setkey_string(rs1, "BLD_VERS", bld_vers); //build vers to every record
}

int nice_intro ()
{
  int usage = cmdparams_get_int (&cmdparams, "h", NULL);
  if (usage)
    {
    printf ("Usage:\nbuild_lev1_aia [-vhr] "
    "mode=<recnum|fsn> instru=<hmi|aia> dsin=<lev0> dsout=<lev1>\n"
    "brec=<rec#>|bfsn=<fsn#> erec=<rec#>|efsn=<fsn#>\n"
    "quicklook=<0|1> [logfile=<file>]\n"
    "  -h: help - show this message then exit\n"
    "  -v: verbose\n"
    "  -r: restart. only used when we restart our selves periodically\n"
        "mode= recnum: brec and erec have the record # range to process \n"
        "      fsn: bfsn and efsn have the fsn # range to process\n"
        "      For safety, the mode and arg name used must be consistent\n"
    "instru= instrument. must be 'hmi' or 'aia'\n"
    "dsin= data set name of lev0 input\n"
    "      default hmi=hmi.lev0e   aia=aia.lev0e\n"
    "dsout= data set name of lev1 output\n"
    "      default hmi=su_production.hmi_lev1e   aia=su_production.aia_lev1e\n"
    "dsff= data set name of AIA flat field series\n"
        "      default aia_test.flatfield\n"
        "dsaiabad= AIA bad pixel series, default is to use BPL from dsff,\n"
               "but bugs need to be fixed\n"
    "brec= first lev0 rec# to process. 0=error must be given by build_lev1_mgr\n"
    "erec= last lev0 rec# to process. 0=error must be given by build_lev1_mgr\n"
    "bfsn= first fsn# to process. 0=error must be given by build_lev1_mgr\n"
    "efsn= last fsn# to process. 0=error must be given by build_lev1_mgr\n"
    "quicklook= 1 = quicklook mode, 0 = definitive mode\n"
    "logfile= optional log file name. If not given uses:\n"
        "         /usr/local/logs/lev1/build_lev1_aia.<time_stamp>.log\n");
    return(1);
    }
  verbose = cmdparams_get_int (&cmdparams, "v", NULL);
  restartflg = cmdparams_get_int (&cmdparams, "r", NULL);
  return (0);
}

TIME SDO_to_DRMS_time(int sdo_s, int sdo_ss)
{
static int firstcall = 1;
if (firstcall)
  {
  firstcall = 0;
  }
/* XXX fix build 3/18/2008, arta */
return(sdo_epoch + (TIME)sdo_s + (TIME)(sdo_ss & 0xFFFF)/65536.0);
}

// Setup global datestr[] like: 2008.07.14_08:29:31
char *do_datestr() {
  time_t tval;
  struct tm *t_ptr;

  tval = time(NULL);
  t_ptr = localtime(&tval);
  sprintf(datestr, "%d.%02d.%02d_%02d:%02d:%02d", 
      (t_ptr->tm_year+1900), (t_ptr->tm_mon+1),
      t_ptr->tm_mday, t_ptr->tm_hour, t_ptr->tm_min, t_ptr->tm_sec);
  return(datestr);
}

// Returns a time tag like  yyyy.mm.dd.hhmmss 
char *gettimetag()
{
  struct timeval tvalr;
  struct tm *t_ptr;

  gettimeofday(&tvalr, NULL);
  t_ptr = localtime((const time_t *)&tvalr);
  sprintf(timetag, "%04d.%02d.%02d.%02d%02d%02d",
        (t_ptr->tm_year+1900), (t_ptr->tm_mon+1), t_ptr->tm_mday, t_ptr->tm_hour, t_ptr->tm_min, t_ptr->tm_sec);
  return(timetag);
}


void BeginTimer(int n)
{
  gettimeofday (&first[n], NULL);
}

float EndTimer(int n)
{
  gettimeofday (&second[n], NULL);
  if (first[n].tv_usec > second[n].tv_usec) {
    second[n].tv_usec += 1000000;
    second[n].tv_sec--;
  }
  return (float) (second[n].tv_sec-first[n].tv_sec) +
    (float) (second[n].tv_usec-first[n].tv_usec)/1000000.0;
}

// Outputs the variable format message (re: printf) to the log file.
int h1log(const char *fmt, ...)
{
  va_list args;
  char string[32768];

  va_start(args, fmt);
  vsprintf(string, fmt, args);
  if(h1logfp) {
    fprintf(h1logfp, string);
    fflush(h1logfp);
  } 
  else {            // couldn't open log 
    printf(string);     // also print to stdout
    fflush(stdout);
  }
  va_end(args);
  return(0);
}

int send_mail(char *fmt, ...)
{
  va_list args;
  char string[1024], cmd[1024];

  va_start(args, fmt);
  vsprintf(string, fmt, args);
  sprintf(cmd, "echo \"%s\" | Mail -s \"build_lev1_aia mail\" lev0_user", string);
  system(cmd);
  va_end(args);
  return(0);
}

// Got a fatal error. 
void abortit(int stat)
{
  printk("***Abort in progress ...\n");
  printk("**Exit build_lev1_aia w/ status = %d\n", stat);
  if (h1logfp) fclose(h1logfp);
  exit(stat);
}

int rdout_mode_correct()
{
  return(0);
}

int orbit_calc()
{
  return(0);
}

//#include "aia_despike.c"
#include "do_flat.c"
#include "get_image_location.c"
#include "limb_fit_function.c"
#include "cosmic_ray.c"
#include "heightformation.c"

//Called with the range to do. The args are either rec#s or fsn#s accord to modeflg
int do_ingest(long long bbrec, long long eerec, const char *dpath)
{
  //FILE *fwt;
  Image_Location *p_imageloc;
  Image_Location imageloc[NUMRECLEV1];
  TIME t_obs0;
  TIME tobs[NUMRECLEV1];
  float percentd;
  float cdelt1, rsun, crpix1, crpix2, crota2;
  double rsun_lf, x0_lf, y0_lf;
  int rstatus, dstatus, lstatus, ncnt, fcnt, i, j, k, qualint, nobs;
  int hshiexp, hcamid, nbad, n_cosmic;
  int *spikedata, status, axes[2], nbytes;
  uint32_t missvals, totvals;
  long long recnum0, recnum1, recnumff;
  char recrange[128], lev0name[128], flatrec[128];
  //char tmpname[80];

  if(modeflg) sprintf(recrange, ":#%lld-#%lld", bbrec, eerec);
  else sprintf(recrange, "%lld-%lld", bbrec, eerec);
  sprintf(open_dsname, "%s[%s][? t_obs > 0 ?]", dsin, recrange);
  printk("open_dsname = %s\n", open_dsname);
  printk("#levnum recnum fsn\n");

    if(hmiaiaflg) t_obs0 = 0;
    rset0 = drms_open_records(drms_env, open_dsname, &rstatus); //open lev0
    if(!rset0 || (rset0->n == 0) || rstatus) {
      printk("Can't do drms_open_records(%s)\n", open_dsname);
      return(1);
    }
    drms_stage_records(rset0, 1, 0);
    ncnt = rset0->n;
    rptr = (DRMS_Record_t *)malloc(ncnt * sizeof(DRMS_Record_t));
    if(rptr == NULL) {
      printk("Can't malloc() for DRMS_Record_t sort\n");
      return(1);
    }
    //make opened records sequential in mem for sort
    for(i=0; i < ncnt; i++) {
      memcpy(&rptr[i], rset0->records[i], sizeof(DRMS_Record_t));
    }
    //Must sort to get in ascending t_obs order
    qsort(rptr, ncnt, sizeof(DRMS_Record_t), &compare_rptr);

    //this loop is for the benefit of the lev1view display to show
    //info on all lev0 records opened
    for(i=0; i < ncnt; i++) {
      flatmiss[i] = 0;      //init quality flags
      orbmiss[i] = 0;
      limbmiss[i] = 0;
      asdmiss[i] = 0;
      mpdmiss[i] = 0;
      noimage[i] = 0;
      missflg[i] = 0;
      rs0 = &rptr[i];
      recnum0 = rs0->recnum;
      //must get fsn in case got record from last lev1 record
      fsnx = drms_getkey_int(rs0, "FSN", &rstatus); 
      fsnarray[i] = fsnx;
      printk("*0 %u %u\n", recnum0, fsnx);
      //also set up call for get_pointing_info() and iorbit_getinfo()
      tobs[i] = drms_getkey_time(rs0, "t_obs", &rstatus);
      if(rstatus) {
        printk("Error on drms_getkey_time() fsn=%u. Use DRMS_MISSING_TIME\n", 
        fsnx);
       tobs[i] = DRMS_MISSING_TIME;
      }
    }
    if(rstatus = get_pointing_info(drms_env, tobs, ncnt, &ptinfo)) {
      printk("**ERROR: get_pointing_info() status = %d  fsn tobs ASD:\n", 
        rstatus);
      for(j=0; j < ncnt; j++) {  
        printk("%u %10.5f ", fsnarray[j], tobs[j]);
        if(ptinfo) {
          ptdata = ptinfo[j];
          printk("%s\n", ptdata.asd_rec);
        }
        asdmiss[j] = 1;     //set for QUALITY for ea image
      }
      //return(1);      //!!No,press on- changed 2/22/2011
    }
    if ((IOstatus = iorbit_getinfo(drms_env,
                       orbseries,
                       NULL,
                       IORBIT_Alg_Quadratic,
                       tobs,
                       ncnt,
                       kIORBIT_CacheAction_DontCache,
                       &IOinfo)) != kLIBASTRO_Success)
    {
      if(IOstatus == kLIBASTRO_InsufficientData) {
        printk("***ERROR in iorbit_getinfo: kLIBASTRO_InsufficientData\n");
      }
      else { 
       printk("***ERROR in iorbit_getinfo() status=%d\n", IOstatus);
      }
      for(j=0; j < ncnt; j++) {  //set qual bits
        orbmiss[j] = 1;
      }
      return(1);        //abort. new 2/22/2011
    }
    rset1 = drms_create_records(drms_env, ncnt, dsout, DRMS_PERMANENT,&dstatus);
    if(dstatus) {
      printk("**ERROR: Can't create records for %s\n", dsout);
      for(j=0; j < ncnt; j++) {  //set qual bits
        noimage[j] = 1;
      }
      return(1);    //new 2/22/2011
    }
    //Now fill in info for call to Carl's get_image_location()
    for(i=0; i < ncnt; i++) {
      rs0 = &rptr[i];
      imageloc[i].tobs = tobs[i];
      imageloc[i].camera = drms_getkey_int(rs0, "CAMERA", &rstatus);
      if(rstatus) {
        printk("ERROR: in drms_getkey_int(CAMERA) fsn=%u\n", fsnarray[i]);
      }
      imageloc[i].wavelength = drms_getkey_int(rs0, "WAVELNTH", &rstatus);
      if(rstatus) {
        printk("ERROR: in drms_getkey_int(WAVELNTH) fsn=%u\n", fsnarray[i]);
      }
      snprintf(imageloc[i].telescope, 10, "%s", 
        drms_getkey_string(rs0, "TELESCOP", &rstatus));
      if(rstatus) {
        printk("ERROR: in drms_getkey_string(TELESCOP) fsn=%u\n", fsnarray[i]);
      }
    }
    p_imageloc = imageloc;
    rstatus = get_image_location(drms_env, ncnt, &p_imageloc);
    if(rstatus) {       //error
      printk("ERROR: get_image_location() returns status=%d\n", rstatus);
      for(j=0; j < ncnt; j++) {  //set qual bits
        mpdmiss[i] = 0;
      }
      return(1);
    }

    for(i=0; i < ncnt; i++) {   //do for all the sorted lev0 records
      //StartTimer(2);  //!!TEMP
      rs0 = &rptr[i];
      recnum0 = rs0->recnum;
      fsnx = fsnarray[i]; 
      sprintf(lev0name, "%s[%u]", dsin, fsnx);
      if(drms_getkey_int(rs0, "QUALITY", 0) < 0) {
        printk("Bad QUALITY for %s, no lev1 made\n", lev0name);
        noimage[i] = 1;
        //continue; //make an image anyway so the lev1 fsn will exist
      }
      segment = drms_segment_lookupnum(rs0, 0);
      Array0 = drms_segment_read(segment, DRMS_TYPE_SHORT, &rstatus);
      if(!Array0) {
        printk("Can't do drms_segment_read() %s status=%d\n", 
            lev0name, rstatus);
        noimage[i] = 1;
        return(1);  //return until we learn
        continue;
      }
      l0l1->adata0 = (short *)Array0->data; //free at end
      l0l1->rs0 = rs0;
      l0l1->recnum0 = recnum0;
      l0l1->fsn = fsnx;
      if(hmiaiaflg) {           //aia
    
    // Lookup table #18 corruption 2017.02.07
    // FSN range 145600147 - 145616603, WAVELNTH 1700 only
    // Any pixel with value 16329 should be changed to 714
    
    int lutid = drms_getkey_int(rs0, "LUTID", &rstatus);
    if (fsnx >= 145600147 && fsnx <=145616603 && lutid == 18) {
        for (int i=0; i<4096*4096; ++i)
        if (l0l1->adata0[i] == 16329) l0l1->adata0[i] = 714;
    }

        l0l1->dat1.adata1A = &data1A[0];
        //l0l1->himgcfid = drms_getkey_int(rs0, "AIFDBID", &rstatus);
    l0l1->himgcfid = 90;    
        aiftsid = drms_getkey_int(rs0, "AIFTSID", &rstatus);
        aifcps = drms_getkey_short(rs0, "AIFCPS", &rstatus);
        aiagp6 = drms_getkey_int(rs0, "AIAGP6", &rstatus);
      }
      else {
        l0l1->dat1.adata1 = &data1[0];
        l0l1->himgcfid = drms_getkey_int(rs0, "HIMGCFID", &rstatus);
      }
      if(rstatus) {
        printk("Can't do drms_getkey_int(HIMGCFID) for fsn %u\n", fsnx);
        l0l1->himgcfid = 104; 
        //return(1);    //return until we learn
        //continue; //maybe cleanup and continue here
      }

      sprintf(open_dsname, "%s[%u]", dsout, fsnx);
      rs = rset1->records[i]; 
      drms_record_directory(rs, rs1_path, 0);
      if(!*rs1_path) {
        printk("***ERROR: No path to segment for %s\n", open_dsname);
        noimage[i] = 1;
        continue;
      }
      //printf("\npath to lev1 = %s\n", rs1_path);  //!!TEMP
      dstatus = drms_setkey_int(rs, "FSN", fsnx);
      dstatus = drms_setkey_string(rs, "LEV0SERIES", lev0name);
      if(!(segment = drms_segment_lookup(rs, "image_lev1"))) {
        printk("No drms_segment_lookup(rs, image_lev1) for %s\n", open_dsname);
        noimage[i] = 1;
        continue;
      }
      if(hmiaiaflg) {
        segArray = drms_array_create(DRMS_TYPE_INT,
                                       segment->info->naxis,
                                       segment->axis,
                                       &data1A,
                                       &dstatus);
      }
      else {
        segArray = drms_array_create(DRMS_TYPE_FLOAT,
                                       segment->info->naxis,
                                       segment->axis,
                                       &data1,
                                       &dstatus);
      }
      //transer the lev0 keywords
      rstatus = drms_copykeys(rs, rs0, 0, kDRMS_KeyClass_Explicit);
      if(rstatus != DRMS_SUCCESS) {
        printk("Error %d in drms_copykeys() for fsn %u\n", fsnx);
        return(1);      //new 2/22/2011
        continue;
      }
      qualint = drms_getkey_int(rs0, "QUALITY", &rstatus);
      drms_setkey_int(rs, "QUALLEV0", qualint);
      fid = drms_getkey_int(rs0, "FID", &rstatus);

      drms_setkey_time(rs, "T_OBS", tobs[i]);
      printk("t_obs for lev0 = %10.5f\n", tobs[i]);  
      drms_setkey_double(rs, "DATE", CURRENT_SYSTEM_TIME);
      if(ptinfo) {
        ptdata = ptinfo[i];
          drms_setkey_float(rs, "SAT_Y0", ptdata.sat_y0);
          drms_setkey_float(rs, "SAT_Z0", ptdata.sat_z0);
          drms_setkey_float(rs, "SAT_ROT", ptdata.sat_rot);
          drms_setkey_string(rs, "ACS_MODE", ptdata.acs_mode);
          drms_setkey_string(rs, "ACS_ECLP", ptdata.acs_eclp);
          drms_setkey_string(rs, "ACS_SUNP", ptdata.acs_sunp);
          drms_setkey_string(rs, "ACS_SAFE", ptdata.acs_safe);
          drms_setkey_string(rs, "ASD_REC", ptdata.asd_rec);
          drms_setkey_string(rs, "ACS_CGT", ptdata.acs_cgt);
      }
      if(IOinfo) {
        IOdata = IOinfo[i];
           drms_setkey_double(rs, "HAEX_OBS", IOdata.hciX);
           drms_setkey_double(rs, "HAEY_OBS", IOdata.hciY);
           drms_setkey_double(rs, "HAEZ_OBS", IOdata.hciZ);
           drms_setkey_double(rs, "GAEX_OBS", IOdata.gciX);
           drms_setkey_double(rs, "GAEY_OBS", IOdata.gciY);
           drms_setkey_double(rs, "GAEZ_OBS", IOdata.gciZ);
           //drms_setkey_float(rs, "DSUN_OBS", (float)IOdata.dsun_obs);
           drms_setkey_double(rs, "DSUN_OBS", IOdata.dsun_obs);
           drms_setkey_double(rs, "OBS_VR", IOdata.obs_vr);
           drms_setkey_double(rs, "OBS_VW", IOdata.obs_vw);
           drms_setkey_double(rs, "OBS_VN", IOdata.obs_vn);
           drms_setkey_double(rs, "RSUN_OBS", IOdata.rsun_obs);
           drms_setkey_float(rs, "CRLN_OBS", (float)IOdata.crln_obs);
           drms_setkey_float(rs, "CRLT_OBS", (float)IOdata.crlt_obs);
           drms_setkey_float(rs, "HGLT_OBS", (float)IOdata.crlt_obs);
           drms_setkey_float(rs, "HGLN_OBS", 0.0);
           drms_setkey_int(rs, "CAR_ROT", (int)IOdata.car_rot);
           drms_setkey_string(rs, "ORB_REC", IOdata.orb_rec);
      }
      drms_setkey_float(rs, "X0_MP", imageloc[i].x);
      drms_setkey_float(rs, "Y0_MP", imageloc[i].y);
      drms_setkey_float(rs, "INST_ROT", imageloc[i].instrot);
      //drms_setkey_float(rs, "INST_ROT", 180.0); //force this for now
      drms_setkey_float(rs, "IMSCL_MP", imageloc[i].imscale);
      drms_setkey_string(rs, "MPO_REC", imageloc[i].mpo_rec);

      int camera = drms_getkey_int(rs0, "CAMERA", &rstatus);
      if(rstatus) {
        printk("Can't do drms_getkey_int() for fsn %u\n", fsnx);
        noimage[i] = 1;
        goto TEMPSKIP;
        //return(1);
      }
      if ( 0 == drms_setkey_time(rs, "T_REC", tobs[i])) {
        int status, allstat = 0;
        double tr_step;
        long long tr_index;
        TIME t_rec, tr_epoch;
        tr_index = drms_getkey_longlong(rs, "T_REC_index", &status);
        allstat += status;
        tr_step =  drms_getkey_double(rs, "T_REC_step", &status);
        allstat += status;
        tr_epoch = drms_getkey_time(rs, "T_REC_epoch", &status);
        allstat += status;
        if (0 == allstat) {
          t_rec = tr_epoch + tr_index*tr_step;
          drms_setkey_time(rs, "T_REC", t_rec);
        }
      }

      if(hmiaiaflg && !quicklook) {
        float tempccd, tempgt, tempsmir, tempfpad;
        if(fabs(tobs[i] - t_obs0) > 300.0) {
          char *dstemp = "aia.temperature_summary_300s";
          char *selstr = "select max(t_start) from ";
          char *whrstr = "where t_start <= ";
          int nr;
          if(rs_t) {
            drms_close_records(rs_t, DRMS_FREE_RECORD);
            rs_t = NULL;
          }
          sprintf(open_dsname, "%s[? t_start=(%s %s %s %f) ?]",
                  dstemp, selstr, dstemp, whrstr, tobs[i]);
          //printf(" %s\n", open_dsname);
          rt = NULL;
          rs_t = drms_open_records(drms_env, open_dsname, &rstatus);
          if(rstatus) printk("Can not open temperature series.\n");
          else {
            nr = rs_t->n;
            if (nr>0) {
              if(nr != 1) printk("%d records != 1.\n", nr);
              rt = rs_t->records[0];
            }
          }
          t_obs0 = tobs[i];
        }
        if (rt) {
          switch (camera) {
            int st;
            case 1:
              tempccd  = drms_getkey_float(rt, "TC01_T1_CCD1_MEAN", &st);
              tempgt   = drms_getkey_float(rt, "T08_G1_1_MEAN", &st);
              tempsmir = drms_getkey_float(rt, "T02_T1_SMIR_MEAN", &st);
              tempfpad = drms_getkey_float(rt, "T01_T1_FADP_MEAN", &st);
              break;
            case 2:
              tempccd  = drms_getkey_float(rt, "TC03_T2_CCD1_MEAN", &st);
              tempgt   = drms_getkey_float(rt, "T21_G2_1_MEAN", &st);
              tempsmir = drms_getkey_float(rt, "T15_T2_SMIR_MEAN", &st);
              tempfpad = drms_getkey_float(rt, "T14_T2_FADP_MEAN", &st);
              break;
            case 3:
              tempccd  = drms_getkey_float(rt, "TC05_T3_CCD1_MEAN", &st);
              tempgt   = drms_getkey_float(rt, "T34_G3_1_MEAN", &st);
              tempsmir = drms_getkey_float(rt, "T28_T3_SMIR_MEAN", &st);
              tempfpad = drms_getkey_float(rt, "T27_T3_FADP_MEAN", &st);
              break;
            case 4:
              tempccd  = drms_getkey_float(rt, "TC07_T4_CCD1_MEAN", &st);
              tempgt   = drms_getkey_float(rt, "T47_G4_1_MEAN", &st);
              tempsmir = drms_getkey_float(rt, "T41_T4_SMIR_MEAN", &st);
              tempfpad = drms_getkey_float(rt, "T40_T4_FADP_MEAN", &st);
              break;
          }
          drms_setkey_float(rs, "TEMPCCD", tempccd);
          drms_setkey_float(rs, "TEMPGT", tempgt);
          drms_setkey_float(rs, "TEMPCEB", tempgt);
          drms_setkey_float(rs, "TEMPSMIR", tempsmir);
          drms_setkey_float(rs, "TEMPFPAD", tempfpad);
          drms_setkey_float(rs, "TEMPPMIR", tempfpad);
        }
      }
      if(hmiaiaflg) {
        int nr, st, ver_num;
        float dt, eperdn, dnperpht, eff_area, eff_wl, factor, p1, p2, p3;
        TIME t_start;
        char *dsresp = "aia.response";
        char *selstr = "select max(t_start) from ";
        char *whrstr = "where t_start <= ";
        char *wavstr = drms_getkey_string(rs0, "WAVE_STR", &rstatus);
        if(rs_resp) {
//          drms_close_records(rs_resp, DRMS_FREE_RECORD);
          rs_resp = NULL;
        }
//        sprintf(open_dsname, "%s[][%s][? t_start=(%s %s %s %f) ?]",
//                dsresp, wavstr, selstr, dsresp, whrstr, tobs[i]);
        //change from Rock on Jan 31, 2012
        sprintf(open_dsname, "%s[][%s][? t_start <= %10.5f and t_stop > %10.5f ?]",
                dsresp, wavstr, tobs[i], tobs[i]);
        rresp = NULL;
        rs_resp = drms_open_records(drms_env, open_dsname, &rstatus);
        if(rstatus) printk("Can not open aia.response series.\n");
        else {
          nr = rs_resp->n;
          if (nr>0) rresp = rs_resp->records[0];
        }
        if(rresp) {
          eperdn = drms_getkey_float(rresp, "EPERDN", &st);
          drms_setkey_float(rs, "DN_GAIN", eperdn);
          dnperpht = drms_getkey_float(rresp, "DNPERPHT", &st);
          drms_setkey_float(rs, "DNPERPHT", dnperpht);
          eff_wl = drms_getkey_float(rresp, "EFF_WVLN", &st);
          drms_setkey_float(rs, "EFF_WVLN", eff_wl);
          p1 = drms_getkey_float(rresp, "EFFA_P1", &st);
          p2 = drms_getkey_float(rresp, "EFFA_P2", &st);
          p3 = drms_getkey_float(rresp, "EFFA_P3", &st);
          t_start = drms_getkey_float(rresp, "T_START", &st);
          dt = (float) (tobs[i] - t_start)/86400.0;
          factor = ((p3*dt + p2)*dt + p1)*dt + 1.0;
          eff_area = drms_getkey_float(rresp, "EFF_AREA", &st);
          eff_area = eff_area*factor;
          drms_setkey_float(rs, "EFF_AREA", eff_area);
          ver_num = drms_getkey_int(rresp, "VER_NUM", &st);
          drms_setkey_int(rs, "DN_GN_V", ver_num);
          drms_setkey_int(rs, "EFF_AR_V", ver_num);
        }
      }

      //Now figure out what flat field to use
      if(drms_ismissing_time(tobs[i])) {
        printk("DRMS_MISSING_TIME for fsn=%u. Continue...\n", fsnx);
        noimage[i] = 1;
        goto TEMPSKIP;
      }
    if(!hmiaiaflg) {        //HMI
      if(quicklook) {
        sprintf(open_dsname, "%s[? t_start=(select max(t_start) from %s where t_start <= %10.5f and t_stop > %10.5f and CAMERA=%d) and CAMERA=%d ?]",
            dsffname, dsffname, tobs[i], tobs[i], camera, camera);
      }
      else {
        sprintf(open_dsname, "%s[? t_start <= %10.5f and t_stop > %10.5f and CAMERA=%d and flatfield_version >= 1 ?]",
        dsffname, tobs[i], tobs[i], camera);
      }
    }
    else {          //AIA
      char *wavstr = drms_getkey_string(rs0, "WAVE_STR", &rstatus);
      if(rstatus) {
        printk("Can't do drms_getkey_string() for WAVE_STR\n");
        return(1);
      }
      if(quicklook) {
        sprintf(open_dsname, "%s[? t_start=(select max(t_start) from %s where t_start <= %10.5f and t_stop > %10.5f and WAVE_STR='%s') and WAVE_STR='%s' ?]",
            dsffname, dsffname, tobs[i], tobs[i], wavstr, wavstr);
      }
      else {
        sprintf(open_dsname, "%s[? t_start <= %10.5f and t_stop > %10.5f and WAVE_STR='%s' ?]",
        dsffname, tobs[i], tobs[i], wavstr);
      }
    }
      //printf("!!TEMP Flat field query: %s\n", open_dsname); //!!TEMP
      rsetff = drms_open_records(drms_env, open_dsname, &rstatus); //open FF 
      if(!rsetff || (rsetff->n == 0) || rstatus) {
        printk("Can't do drms_open_records(%s)\n", open_dsname);
        flatmiss[i] = 1; noimage[i] = 1;
        goto TEMPSKIP;
        return(1);      //new 2/22/2011
      }
      fcnt = rsetff->n;
      if(fcnt > 1) {
        printk("More than one FF found for %s?\n", open_dsname);
        printk("Use last one of %d found\n", fcnt); 
        //return(1);        //!!TBD
      }
      //rsff = rsetff->records[0];
      rsff = rsetff->records[fcnt-1];
      recnumff = rsff->recnum;
      sprintf(flatrec, "%s[:#%lld]", dsffname, recnumff);
      if(dstatus = drms_setkey_string(rs, "FLAT_REC", flatrec )) {
        printk("**ERROR on drms_setkey_string() for %s\n", flatrec);
      }
      drms_record_directory(rsff, path, 1);
      if(!*path) {
        printk("***ERROR: No path to segment for %s\n", open_dsname);
        //goto TEMPSKIP;    //!!!TEMP until have good flatfield
        return(1);
      }
      //printf("\npath to FF = %s\n", path);    //!!TEMP
      segmentff = drms_segment_lookup(rsff, "flatfield");
      Arrayff = drms_segment_read(segmentff, DRMS_TYPE_FLOAT, &rstatus);
      if(!Arrayff) {
        printk("Can't do drms_segment_read() for Flat Field status=%d\n", 
            rstatus);
        return(1);
      }
      l0l1->adataff = (float *)Arrayff->data; 

      darkseg = drms_segment_lookup(rsff, "DARK");
      ArrayDark = drms_segment_read(darkseg, DRMS_TYPE_FLOAT, &rstatus);
      if(!ArrayDark) {
        printk("Can't do drms_segment_read() for DARK. status=%d\n", rstatus);
        return(1);
      }
      l0l1->adatadark = (float *)ArrayDark->data; //free at end

      if (strcmp(dsaiabad, NOTSPECIFIED)) {
        char *wavstr = drms_getkey_string(rs0, "WAVE_STR", &rstatus);
        sprintf(open_aiabad_dsname, "%s[%s][]", dsaiabad, wavstr);
        rsbad_aia = drms_open_records(drms_env, open_aiabad_dsname, &rstatus);
        if(!rsbad_aia || rsbad_aia->n == 0 || rstatus)
        printk("drms_open_records for rsbad_aia failed\n");
        fcnt = rsbad_aia->n;
        if(fcnt > 1) {
          printk("More than one bpl found for %s?\n", open_aiabad_dsname);
        }
        rec_bad_aia = rsbad_aia->records[fcnt-1];
        badseg = drms_segment_lookup(rec_bad_aia, "bad_pixel_list");
      } else {
        badseg = drms_segment_lookup(rsff, "BAD_PIXEL");
      }
      ArrayBad = drms_segment_read(badseg, DRMS_TYPE_INT, &rstatus);
      nbad = drms_array_size(ArrayBad)/sizeof(int);
      if(!ArrayBad) {
        printk("Can't do drms_segment_read() for BAD_PIXEL. status=%d\n", 
            rstatus);
        return(1);
      }
      l0l1->adatabad = (int *)ArrayBad->data; //free at end
      /************************************
      if(!(badoutpixseg = drms_segment_lookup(rs, "BAD_PIXEL"))) BRACKET
      if(!(badoutpixseg = drms_segment_lookup(rs, "bad_pixel"))) BRACKET
      ***************************************/
      if(!hmiaiaflg) {      //HMI. !!TBD make the aia & hmi jsd agree
        if(!(badoutpixseg = drms_segment_lookup(rs, "bad_pixel_list"))) {
          printk("No drms_segment_lookup(rs, bad_pixel_list) for lev1\n");
          return(1);
        }
      }
      else {
        if(!(badoutpixseg = drms_segment_lookup(rs, "bad_pixel"))) {
          printk("No drms_segment_lookup(rs, bad_pixel) for lev1\n");
          return(1);
        }
      }
      //move the segment write to after cosmic_rays() updates it
      /**************************************************************
      dstatus = drms_segment_write(badoutpixseg, ArrayBad, 0);
      if (dstatus) {
        printk("ERROR: drms_segment_write error=%d for lev1 bad_pixel_list\n");
        //noimage[i] = 1;
        return(1);
      }
      **************************************************************/
      l0l1->rs1 = rs;
      l0l1->rsff = rsff;
      l0l1->recnum1 = rs->recnum;  
      l0l1->darkflag = 0;
      hshiexp = drms_getkey_int(rs, "HSHIEXP", &rstatus);
      hcamid = drms_getkey_int(rs, "HCAMID", &rstatus);
      if(hmiaiaflg) {
        float sumdc=0.0;
        int idc, numdc=0;
        int aimgshce = drms_getkey_int(rs, "AIMGSHCE", &rstatus);
        if(aimgshce == 0) l0l1->darkflag = 1;
        if(rstatus = do_flat_aia(l0l1)) {
          printk("***ERROR in do_flat_aia() status=%d\n", rstatus);
          printf("***ERROR in do_flat_aia() status=%d\n", rstatus);
          flatmiss[i] = 1; noimage[i] = 1;
          //return(1);      //!!TBD what to do?
          goto FLATERR;
        }
        for (idc=2047; idc<4096*4096; idc+=4096) {
          if (data1A[idc] != DRMS_MISSING_INT) {
            sumdc = sumdc + data1A[idc];
            numdc++;
          }
          if (data1A[idc+1] != DRMS_MISSING_INT) {
            sumdc = sumdc + data1A[idc+1];
            numdc++;
          }
        }
        if (numdc) drms_setkey_float(rs, "DATACENT", sumdc/numdc);
      }
      else {
        if(hshiexp == 0) l0l1->darkflag = 1;
        //StartTimer(1);    //!!TEMP
        if(rstatus = do_flat(l0l1)) {
          printk("***ERROR in do_flat() status=%d\n", rstatus);
          printf("***ERROR in do_flat() status=%d\n", rstatus);
          flatmiss[i] = 1; noimage[i] = 1;
          //return(1);      //!!TBD what to do?
          goto FLATERR;
        }
        //ftmp = StopTimer(1);
        //printf( "\nTime sec for hmi do_flat() = %f\n\n", ftmp );
      } 
        //sprintf(tmpname, "/tmp/data_lev1.%u", fsnx);
        //fwt = fopen(tmpname, "w");
        //int wsize = 4096 * 4096;
        //fwrite(l0l1->adata1, sizeof(float), wsize, fwt);
        //fclose(fwt);
        drms_setkey_float(rs, "OSCNMEAN", l0l1->oscnmean);
        drms_setkey_float(rs, "OSCNRMS", l0l1->oscnrms);
        drms_setkey_int(rs, "DATAMIN", l0l1->datamin);
        drms_setkey_int(rs, "DATAMAX", l0l1->datamax);
        drms_setkey_int(rs, "DATAMEDN", l0l1->datamedn);
        drms_setkey_float(rs, "DATAMEAN", l0l1->datamean);
        drms_setkey_float(rs, "DATARMS", l0l1->data_rms);
        drms_setkey_float(rs, "DATASKEW", l0l1->dataskew);
        drms_setkey_float(rs, "DATAKURT", l0l1->datakurt);
        drms_setkey_int(rs, "DATAVALS", l0l1->datavals);
        drms_setkey_int(rs, "MISSVALS", l0l1->missvals);
        missvals = (uint32_t)l0l1->missvals;
        totvals = (uint32_t)l0l1->datavals + missvals;
        drms_setkey_int(rs, "TOTVALS", (int)totvals);
        percentd = (float)((100.0 * (float)l0l1->datavals)/(float)totvals);
        drms_setkey_float(rs, "PERCENTD", percentd);
        if(missvals > 0) missflg[i] = missflg[i] | Q_1_MISS0;
        if(missvals > (uint32_t)(totvals * 0.01))
           missflg[i] = missflg[i] | Q_1_MISS1;
        if(missvals > (uint32_t)(totvals * 0.05)) 
           missflg[i] = missflg[i] | Q_1_MISS2;
        if(missvals > (uint32_t)(totvals * 0.25)) 
           missflg[i] = missflg[i] | Q_1_MISS3;
        if(l0l1->datavals == 0) 
           missflg[i] = missflg[i] | Q_MISSALL; //high bit, no data
        if(hmiaiaflg && nspikes) {
          if (spikeseg = drms_segment_lookup(rs,"spikes") ) {
            nbytes = nspikes*sizeof(int);
            axes[0] = nspikes;
            axes[1] = 3;
            spikedata = (int *) malloc (3*nspikes*sizeof(int));
            ArraySpike = drms_array_create(DRMS_TYPE_INT, 2, axes,
                       (void *) spikedata, &status);
            memcpy((void *)spikedata, (void *)spikelocs, nbytes);
            memcpy((void *)(spikedata+nspikes), (void *)oldvalues, nbytes);
            memcpy((void *)(spikedata+2*nspikes), (void *)newvalues, nbytes);
            status = drms_segment_write(spikeseg, ArraySpike, 0);
            drms_free_array(ArraySpike);
          } // else { printf("spikes segment not found\n"); }
        }
      if(!hmiaiaflg) {          //only call for hmi
        //StartTimer(1);    //!!TEMP
        dstatus = cosmic_rays(rs, l0l1->dat1.adata1, l0l1->adatabad, nbad,
                array_cosmic, &n_cosmic, 4096, 4096);
        //ftmp = StopTimer(1);
        //printf( "\nTime sec for cosmic_rays() = %f\n\n", ftmp );
        if(dstatus) {
          printk("ERROR: cosmic_rays() error=%d\n", dstatus);
          noimage[i] = 1;
        }
        else {
          drms_setkey_int(rs, "NBADPERM", nbad);
          drms_setkey_int(rs, "NBADTOT", n_cosmic);
        }
        if(nbad != n_cosmic) {  //pretend the bad pix list is like spike
          nbytes = n_cosmic*sizeof(int);
          axes[0] = n_cosmic;
          spikedata = (int *)malloc(nbytes);
          ArraySpike = drms_array_create(DRMS_TYPE_INT, 1, axes,
                       (void *)spikedata, &status);
          memcpy((void *)spikedata, (void *)array_cosmic, nbytes);
          status = drms_segment_write(badoutpixseg, ArraySpike, 0);
          drms_free_array(ArraySpike);
        }
        else {
          dstatus = drms_segment_write(badoutpixseg, ArrayBad, 0);
        }
        if (dstatus) {
          printk("ERROR: drms_segment_write error=%d for lev1 bad_pixel_list\n",
             dstatus);
          noimage[i] = 1;
        }
      }
FLATERR:
      drms_close_records(rsetff, DRMS_FREE_RECORD);
      free(ArrayDark->data);
      free(Arrayff->data);
      free(ArrayBad->data);

TEMPSKIP:
      if(!hmiaiaflg) {
        segArray->type = DRMS_TYPE_FLOAT;
        segArray->bscale = 1.0; 
        segArray->bzero = 0.0; 
        dstatus = drms_segment_write(segment, segArray, 0);
        if (dstatus) {
          printk("ERROR: drms_segment_write error=%d for fsn=%u\n", dstatus,fsnx);
          noimage[i] = 1;
        }
        recnum1 = rs->recnum;
        printk("*1 %u %u\n", recnum1, fsnx);
        free(Array0->data);
      }
 
  x0_lf = DRMS_MISSING_DOUBLE;
  y0_lf = DRMS_MISSING_DOUBLE;
  rsun_lf = DRMS_MISSING_DOUBLE;
  //Don't send calmode or darks to limb_fit()
  //calmode: HCFTID=17
  //dark: HSHIEXP=0 and HCAMID=0 or 1
  int skiplimb = 0;
  hcftid = drms_getkey_int(rs, "HCFTID", &rstatus);
  if(hcftid == CAL_HCFTID) {    //don't call limb_fit()
    printk("Cal mode image fsn=%u\n", fsnx);
    skiplimb = 1;
  }
  else {
    if(hshiexp == 0) {
      if(hcamid == 0 || hcamid == 1) {
        printk("Dark image fsn=%u\n", fsnx);
        skiplimb = 1;
      }
    }
  }
  lstatus = 1;              //default bad limb_fit
  //WCS calculations for condition 1
  //(For WCS conditions see mail from Rock 10Aug2010 16:40)
  if(!skiplimb && !hmiaiaflg) {     //only call for HMI
    //StartTimer(1);    //!!TEMP
     lstatus = limb_fit(rs,l0l1->dat1.adata1,&rsun_lf,&x0_lf,&y0_lf,4096,4096,0);
    if(lstatus) {
      printk("ERROR: limb_fit() %d error for fsn=%u\n", lstatus, fsnx);
      //noimage[i] = 1;
      limbmiss[i] = 1;
    }
    else {      //limb_fit() returns good values
      //ftmp = StopTimer(1);
      //printf( "\nTime sec for limb_fit() = %f\n\n", ftmp );
      //printf("Calling WCS condition 1\n"); //!!TEMP
      drms_setkey_float(rs, "RSUN_LF", (float)rsun_lf);
      drms_setkey_float(rs, "R_SUN", (float)rsun_lf);
      drms_setkey_float(rs, "X0_LF", (float)x0_lf);
      drms_setkey_float(rs, "Y0_LF", (float)y0_lf);
      drms_setkey_float(rs, "CRVAL1", 0.0);
      drms_setkey_float(rs, "CRVAL2", 0.0);
      drms_setkey_string(rs, "CTYPE1", "HPLN-TAN");
      drms_setkey_string(rs, "CTYPE2", "HPLT-TAN");
      drms_setkey_string(rs, "CUNIT1", "arcsec");
      drms_setkey_string(rs, "CUNIT2", "arcsec");
      cdelt1 = (float)IOdata.rsun_obs/rsun_lf;
      drms_setkey_float(rs, "CDELT1", cdelt1);
      drms_setkey_float(rs, "CDELT2", cdelt1);
      crpix1 = (float)x0_lf + 1;  //set up for correction by heightformation()
      crpix2 = (float)y0_lf + 1;
      drms_setkey_float(rs, "CRPIX1", crpix1);
      drms_setkey_float(rs, "CRPIX2", crpix2);
      rsun = (float)rsun_lf;
      crota2 = imageloc[i].instrot + ptdata.sat_rot;
      drms_setkey_float(rs, "CROTA2", crota2);
      goto WCSEND;
    }
  }
  //WCS calculations for condition 2
  //(As a first pass, this is following Rock's notes and can be simplified)
  int cond2 = 0;
  if(!strcmp(ptdata.acs_mode, "SCIENCE") || !strcmp(ptdata.acs_mode, DRMS_MISSING_STRING)) {
    //if(quicklook) {
      if(!skiplimb && !hmiaiaflg) { //hmi, not dark or cal
        cond2 = 1;
      }
      if(hmiaiaflg) {
        cond2 = 1;
        if(hshiexp == 0) {
          if(hcamid == 0 || hcamid == 1) {  //it's aia dark
            cond2 = 0;
          }
        }
      }
      if(cond2) {           //this is WCS condition 2
/**********************************
        if(!hmiaiaflg) {
          cond2 = 0;
          if((rsun_lf == DRMS_MISSING_DOUBLE) || (x0_lf == DRMS_MISSING_DOUBLE) || (y0_lf == DRMS_MISSING_DOUBLE)) {
            cond2 = 1;
          }
        }
**************************************/
        if(cond2) {
          printf("Calling WCS condition 2\n"); 
          rsun = (float)IOdata.rsun_obs/imageloc[i].imscale;
          drms_setkey_float(rs, "R_SUN", rsun);
          drms_setkey_string(rs, "CTYPE1", "HPLN-TAN");
          drms_setkey_string(rs, "CTYPE2", "HPLT-TAN");
          drms_setkey_float(rs, "CRVAL1", 0.0);
          drms_setkey_float(rs, "CRVAL2", 0.0);
          drms_setkey_string(rs, "CUNIT1", "arcsec");
          drms_setkey_string(rs, "CUNIT2", "arcsec");
          cdelt1 = (float)imageloc[i].imscale;
          drms_setkey_float(rs, "CDELT1", cdelt1);
          drms_setkey_float(rs, "CDELT2", cdelt1);
          crpix1 = imageloc[i].x + 1;
          crpix2 = imageloc[i].y + 1;
          drms_setkey_float(rs, "CRPIX1", crpix1);
          drms_setkey_float(rs, "CRPIX2", crpix2);
          crota2 = imageloc[i].instrot + ptdata.sat_rot;
          drms_setkey_float(rs, "CROTA2", crota2);
          goto WCSEND;
        }
      }
    //}
  } 
  //WCS calculations for condition 3
  int cond3 = 0;
  if(strcmp(ptdata.acs_mode, "SCIENCE")) {  //not is sci pointing mode
    if(!hmiaiaflg && !skiplimb && lstatus) { //hmi,limb fit NG, not dark or cal
      cond3 = 1;
    }
    else {
      if(hmiaiaflg) {
        cond3 = 1;
        if(hshiexp == 0) {
          if(hcamid == 0 || hcamid == 1) {  //it's aia dark
            cond3 = 0;
          }
        }
      }
    }
  }
/******extraneous already lstatus means limb fit ng
  if(cond3) {
    if(!hmiaiaflg) {
      cond3 = 0;
      if((rsun_lf == DRMS_MISSING_DOUBLE) || (x0_lf == DRMS_MISSING_DOUBLE) || (y0_lf == DRMS_MISSING_DOUBLE)) {
        cond3 = 1;
      }
    }
  }
**************************************/
  if(cond3) {
    printf("Calling WCS condition 3\n"); 
    rsun = (float)IOdata.rsun_obs/imageloc[i].imscale;
    drms_setkey_float(rs, "R_SUN", rsun);
    drms_setkey_string(rs, "CTYPE1", "HPLN-TAN");
    drms_setkey_string(rs, "CTYPE2", "HPLT-TAN");
    drms_setkey_string(rs, "CUNIT1", "arcsec");
    drms_setkey_string(rs, "CUNIT2", "arcsec");
    drms_setkey_float(rs, "CRVAL1", 0.0);
    drms_setkey_float(rs, "CRVAL2", 0.0);
    cdelt1 = (float)imageloc[i].imscale;
    drms_setkey_float(rs, "CDELT1", cdelt1);
    drms_setkey_float(rs, "CDELT2", cdelt1);
    //below missing - SC_Y/Z_INRT_BIAS/IMSCL_MP + 1
    if(hmiaiaflg) { //aia
      crpix1 = imageloc[i].x + (ptdata.sat_y0 - imageloc[i].yinrtb)/imageloc[i].imscale + 1;
      crpix2 = imageloc[i].y + (ptdata.sat_z0 - imageloc[i].zinrtb)/imageloc[i].imscale + 1;
    }
    else {      //hmi
      crpix1 = imageloc[i].x - (ptdata.sat_y0 - imageloc[i].yinrtb)/imageloc[i].imscale + 1;
      crpix2 = imageloc[i].y - (ptdata.sat_z0 - imageloc[i].zinrtb)/imageloc[i].imscale + 1;
    }
    drms_setkey_float(rs, "CRPIX1", crpix1);
    drms_setkey_float(rs, "CRPIX2", crpix2);
    crota2 = imageloc[i].instrot + ptdata.sat_rot;
    drms_setkey_float(rs, "CROTA2", crota2);
    goto WCSEND;
  }
  //WCS calculations for condition 4
  int cond4 = 0;
  if(!hmiaiaflg && skiplimb) {  //hmi, dark or cal
    cond4 = 1;
  }
  else {
    if(hshiexp == 0) {
      if(hcamid == 0 || hcamid == 1) {  //it's aia dark
        cond4 = 1;
      }
    }
  }
  if(cond4) {
    printf("Calling WCS condition 4\n"); 
    rsun = DRMS_MISSING_FLOAT;
    drms_setkey_float(rs, "R_SUN", rsun);
    drms_setkey_string(rs, "CTYPE1", "RAW");
    drms_setkey_string(rs, "CTYPE2", "RAW");
    drms_setkey_string(rs, "CUNIT1", DRMS_MISSING_STRING);
    drms_setkey_string(rs, "CUNIT2", DRMS_MISSING_STRING);
    drms_setkey_float(rs, "CRVAL1", DRMS_MISSING_FLOAT);
    drms_setkey_float(rs, "CRVAL2", DRMS_MISSING_FLOAT);
    cdelt1 = DRMS_MISSING_FLOAT;
    drms_setkey_float(rs, "CDELT1", cdelt1);
    drms_setkey_float(rs, "CDELT2", cdelt1);
    crpix1 = DRMS_MISSING_FLOAT;
    crpix2 = DRMS_MISSING_FLOAT;
    drms_setkey_float(rs, "CRPIX1", crpix1);
    drms_setkey_float(rs, "CRPIX2", crpix2);
    crota2 = imageloc[i].instrot + ptdata.sat_rot;
    drms_setkey_float(rs, "CROTA2", crota2);
    goto WCSEND;
  }

WCSEND:
    if(!hmiaiaflg && !lstatus) {    //only do for hmi and good limb fit
      //Now call Sebastien's heightformation() fuction (email 08/09/10 17:50)
      //21Sep2010 change -crota2 to crota2 and HFCORRVR to 2
      if(!(dstatus = heightformation(fid, IOdata.obs_vr, &cdelt1, &rsun, &crpix1, &crpix2, crota2))) {
        drms_setkey_float(rs, "CDELT1", cdelt1);
        drms_setkey_float(rs, "CDELT2", cdelt1);
        drms_setkey_float(rs, "R_SUN", rsun);
        drms_setkey_float(rs, "CRPIX1", crpix1);
        drms_setkey_float(rs, "CRPIX2", crpix2);
        drms_setkey_int(rs, "HFCORRVR", 2);
      }
      else {
        drms_setkey_int(rs, "HFCORRVR", 0);
        printk("ERROR: heightformation() returned error for FID=%d\n", fid);
      }
    }

  //if(hmiaiaflg) {                       //aia
  //  int wl = drms_getkey_int(rs, "WAVELNTH", &rstatus);
  //}

  do_quallev1(rs0, rs, i, fsnx);
  //ftmp = StopTimer(2); //!!!TEMP
  //printf( "\nTime sec in inside loop for fsn=%u : %f\n\n", fsnx, ftmp );
    if(hmiaiaflg) {
      dstatus = drms_segment_writewithkeys(segment, segArray, 0);
      if (dstatus) {
        printk("ERROR: drms_segment_write error=%d for fsn=%u\n", dstatus,fsnx);
        noimage[i] = 1;
      }
      recnum1 = rs->recnum;
      printk("*1 %u %u\n", recnum1, fsnx);
      free(Array0->data);
//      if (rs_t) drms_close_records(rs_t, DRMS_FREE_RECORD);
      if (rs_resp) drms_close_records(rs_resp, DRMS_FREE_RECORD);
      rs_resp = NULL;
    }
  }             //END do for all the sorted lev0 records

  drms_close_records(rset0, DRMS_FREE_RECORD);   //close lev0 records
  drms_close_records(rset1, DRMS_INSERT_RECORD); //close lev1 records
  return(0);
}

int compare_rptr(const void *a, const void *b)
{
  TIME t1, t2;
  int rstatus;
  DRMS_Record_t *x=(DRMS_Record_t *)a, *y=(DRMS_Record_t *)b;

  t1 = drms_getkey_time(x, "t_obs", &rstatus); 
  if(rstatus) t1 = DRMS_MISSING_TIME;   //treat error as missing t_obs
  t2 = drms_getkey_time(y, "t_obs", &rstatus); 
  if(rstatus) t2 = DRMS_MISSING_TIME;
  if(t1 < t2) return(-1);
  if(t1 > t2) return(1);
  return(0);
}

// Initial setup stuff called when main is first entered.
void setup()
{
  FILE *fin;
  char string[128], cwdbuf[128], idstr[256], lfile[128];
  int tpid;

  sdo_epoch = sscan_time("1958.01.01_00:00:00_TAI");
  do_datestr();
  printk_set(h1log, h1log); // set for printk calls 
  printk("%s\n", datestr);
  gethostname(idstr, 256);
  printf("Host: %s\n", idstr);
  printk("Host: %s\n", idstr);
  getcwd(cwdbuf, 126);
  sprintf(idstr, "Cwd: %s\nCall: ", cwdbuf);
  sprintf(string, "build_lev1_aia started as pid=%d ppid=%d user=%s\n", 
        getpid(), getppid(), username);
  strcat(idstr, string);
  printk("%s", idstr);
  printf("%s", idstr);
  if(restartflg) printk("-r ");
  sprintf(argmode, "mode=%s", mode);
  sprintf(arginstru, "instru=%s", instru);
  sprintf(argdsin, "dsin=%s", dsin);
  sprintf(argdsout, "dsout=%s", dsout);
  if(modeflg) {
    sprintf(argbx, "brec=%lld", brec);
    sprintf(argex, "erec=%lld", erec);
  }
  else {
    sprintf(argbx, "bfsn=%lld", bfsn);
    sprintf(argex, "efsn=%lld", efsn);
  }
  sprintf(argquick, "quicklook=%d", quicklook);
  sprintf(arglogfile, "logfile=%s", logname);
  printk("%s %s %s %s %s %s %s %s %s\n", 
    argmode, arginstru, argdsin, argdsout, argbx, argex, argquick, arglogfile, argdsaiabad);
  printf("%s %s %s %s %s %s %s %s %s\n", 
    argmode, arginstru, argdsin, argdsout, argbx, argex, argquick, arglogfile, argdsaiabad);
  if(!restartflg) {
    //printk("tlmseriesname=%s\nlev0seriesname=%s\n", 
    //      tlmseriesname, lev0seriesname);
  }
  sprintf(bld_vers, "%s", jsoc_version);
  sprintf(idstr,  "ps -ef | grep %s", LEV1VIEWERNAME);
  fin = popen(idstr, "r");
  while(fgets(string, sizeof string, fin)) {  //get ps line
    if(!(strstr(string, "perl"))) continue;
    sscanf(string, "%s %d", idstr, &tpid); /* get user name & process id */
    sprintf(lfile, "%s/build_lev1_aia_restart_%d.touch", LEV1LOG_BASEDIR, tpid);
    sprintf(idstr, "/bin/touch %s", lfile);
    printk("%s\n", idstr);
    system(idstr);
  }
  umask(002);           // allow group write 
  //Image.initialized = 0;  // init the two image structures 
  //ImageOld.initialized = 0;
  //Img = &Image;
  //ImgO = &ImageOld;
  //Img->initialized = 0;
  //ImgO->initialized = 0;
}

// Module main function. 
int DoIt(void)
{
  long long numofrecs, frec, lrec;
  int numrec, numofchunks, i;
  char line[80];
  const char *dpath = NULL;

  if (nice_intro())
    return (0);
  if(!(username = (char *)getenv("USER"))) username = "nouser"; 
  instru = cmdparams_get_str(&cmdparams, "instru", NULL);
  if(strcmp(instru, "hmi") && strcmp(instru, "aia")) {
    printf("Error: instru=%s must be given as 'hmi' or 'aia'\n", instru);
    printk("Error: instru=%s must be given as 'hmi' or 'aia'\n", instru);
    return(0);
  }
  mode = cmdparams_get_str(&cmdparams, "mode", NULL);
  if(strcmp(mode, "recnum") && strcmp(mode, "fsn")) {
    printf("Error: mode= must be given as 'recnum' or 'fsn'\n");
    return(0);
  }
  if(!strcmp(instru, "aia")) hmiaiaflg = 1;
  if(!strcmp(mode, "recnum")) modeflg = 1;
  dsin = cmdparams_get_str(&cmdparams, "dsin", NULL);
  dsout = cmdparams_get_str(&cmdparams, "dsout", NULL);
  dsff = cmdparams_get_str(&cmdparams, "dsff", NULL);
  dsaiabad = cmdparams_get_str(&cmdparams, "dsaiabad", NULL);
  if (strcmp(dsaiabad, NOTSPECIFIED)) sprintf(argdsaiabad, "dsaiabad=%s", dsaiabad);
  else sprintf(argdsaiabad, " ");
  brec = cmdparams_get_int(&cmdparams, "brec", NULL);
  erec = cmdparams_get_int(&cmdparams, "erec", NULL);
  bfsn = cmdparams_get_int(&cmdparams, "bfsn", NULL);
  efsn = cmdparams_get_int(&cmdparams, "efsn", NULL);
  quicklook = cmdparams_get_int(&cmdparams, "quicklook", NULL);

  dpath = cmdparams_get_str(&cmdparams, kDpath, NULL);

  //quicklook = 1; //!!TEMP for test
  if(modeflg) {     //recnum mode
    if(brec == 0 || erec == 0) {
      fprintf(stderr, "brec and erec must be given for recnum mode. 0 not allowed\n");
      return(0);
    }
    if(brec > erec) {
      fprintf(stderr, "brec must be <= erec\n");
      return(0);
    }
    bnumx = brec;
    enumx = erec;
  }
  else {        //fsn mode
    if(bfsn == 0 || efsn == 0) {
      fprintf(stderr, "bfsn and efsn must be given for fsn mode. 0 not allowed\n");
      return(0);
    }
    if(bfsn > efsn) {
      fprintf(stderr, "bfsn must be <= efsn\n");
      return(0);
    }
    bnumx = bfsn;
    enumx = efsn;
  }
  logfile = cmdparams_get_str(&cmdparams, "logfile", NULL);
  if (strcmp(dsin, NOTSPECIFIED) == 0) {
    if(hmiaiaflg == 0) dsin = LEV0SERIESNAMEHMI;
    else dsin = LEV0SERIESNAMEAIA;
  }
  if (strcmp(dsout, NOTSPECIFIED) == 0) {
    if(hmiaiaflg == 0) dsout = LEV1SERIESNAMEHMI;
    else dsout = LEV1SERIESNAMEAIA;
  }
  if(hmiaiaflg) {
    if (strcmp(dsff, NOTSPECIFIED)) sprintf(dsffname, "%s", dsff);
    else sprintf(dsffname, "%s", DSFFNAMEAIA);
    if(strstr(dsin, "hmi") || strstr(dsout, "hmi")) {
      printf("Warning: You said instru=aia but have 'hmi' in ds name?\n");
      printf("Do you want to abort this [y/n]? ");
      if(fgets(line, 4, stdin) == NULL) { return(0); }
      if(strcmp(line, "n")) { return(0); }
    }
  }
  else {
    sprintf(dsffname, "%s", DSFFNAMEHMI);
    if(strstr(dsin, "aia") || strstr(dsout, "aia")) {
      printf("Warning: You said instru=hmi but have 'aia' in ds name?\n");
      printf("Do you want to abort this [y/n]? ");
      if(fgets(line, 4, stdin) == NULL) { return(0); }
      if(strcmp(line, "n")) { return(0); }
    }
  }
  if (strcmp(logfile, NOTSPECIFIED) == 0) {
    sprintf(logname, H1LOGFILE, gettimetag());
  }
  else {
    sprintf(logname, "%s", logfile);
  }
  if(restartflg || LOGTEST) {
    if((h1logfp=fopen(logname, "a")) == NULL)
      fprintf(stderr, "**Can't open for append the log file %s\n", logname);
  }
  else {
    if((h1logfp=fopen(logname, "w")) == NULL)
      fprintf(stderr, "**Can't open the log file %s\n", logname);
  }
  setup();
  numofrecs = (enumx - bnumx) + 1;
  numrec = NUMRECLEV1;     //# of records to do at a time
  numofchunks = numofrecs/numrec;
  if((numofrecs % numrec) != 0) numofchunks++; //extra loop for partial chunk
  lrec = bnumx-1;
  for(i = 0; i < numofchunks; i++) {
    frec = lrec+1; lrec = (frec + numrec)-1;
    if(lrec > enumx) lrec=enumx;
    if(do_ingest(frec, lrec, dpath)) {  //do a chunk to get files from the lev0
      printf("build_lev1_aia abort\nSee log: %s\n", logname); 
      send_mail("build_lev1_aia abort\nSee log: %s\n", logname); 
      return(0);
    }
    if(modeflg) {               //only do for recnum mode
      if(stat(STOP_FILE, &stbuf) == 0) {
        printf("Stop file %s seen. Exit build_lev1_aia.\n", STOP_FILE);
        break;
      }
    }
  }
  printf("build_lev1_aia done last fsn=%u\n", fsnx); 
  return(0);
}

---