/* --------------------------------------------------------------------------- */
/*     Conversion from users' observation sequence to KMTNet Script Input      */
/*                                        Developed by Seung-Lee Kim at KASI   */
/* --------------------------------------------------------------------------- */

#include <stdio.h>
#include <stdlib.h>
#include <string.h>

#define MAXCOL 50
#define MAXFIL 50

char version[200] = "Version 1.0 made on March-06-2018\nVersion 1.1 updated on December-21-2018";

void Write_Script(FILE *fp, int loop, int num, char com[MAXCOL][100])
{
    fprintf(fp, "%-13s %11s %11s %4s %8s %12s %6s %7s %19s %5s",
        com[0], com[1], com[2], com[3], com[4], com[5], com[6], com[7], com[8], com[9]);
    if(loop<3) loop = 5;
    if((num%loop)==1) fprintf(fp, " #%d\n", num);
    else if((num%loop)==0) fprintf(fp, "\n\n");
    else fprintf(fp, "\n");
}

int main(int argc, char *argv[])
{
    int i,j, flagloop, nfil, nrun, copt, tol, cLabel, cRA, cDEC, cObject, cFilter, cExposure, cUTObs, maxCol, eol;
    float exptime[MAXFIL];
    char pID[100], filter[MAXFIL][100], imgtyp[100];
    char seqfile[100], setfile[100], outfile[100], field[300], col[MAXCOL][100], prt[MAXCOL][100], tmp[100];
    FILE *in, *out;

    puts(version);
    if(argc < 3) {  printf("Error : <command> %s objects_filename setting_filename flag_loop (option, default : 0)\n", argv[0]);  return 0;  }
    flagloop = nfil = nrun = copt = tol = cLabel = cRA = cDEC = cObject = cFilter = cExposure = cUTObs = maxCol = 0;

    sprintf(seqfile, "%s", argv[1]);
    sprintf(setfile, "%s", argv[2]); 
    sprintf(outfile, "sc%s", argv[1]); 
    if(argc > 3 ) flagloop = atoi(argv[3]);

    if( (in = fopen(setfile, "r") ) == NULL) {  puts("Script setting file read error");  return 0;  }
    fgets(field, 300, in); sscanf(field, "%s", pID);
    fgets(field, 300, in); sscanf(field, "%d", &nfil);
    if(nfil!=0) for(i=0; i<nfil; i++) {  fgets(field, 300, in); sscanf(field, "%s %f", filter[i], &exptime[i]);  }
    fgets(field, 300, in); sscanf(field, "%d", &cLabel);   if(cLabel>maxCol) maxCol = cLabel;
    fgets(field, 300, in); sscanf(field, "%d", &cRA);      if(cRA>maxCol) maxCol = cRA;
    fgets(field, 300, in); sscanf(field, "%d", &cDEC);     if(cDEC>maxCol) maxCol = cDEC;
    fgets(field, 300, in); sscanf(field, "%d", &copt);
    fgets(field, 300, in); sscanf(field, "%s", imgtyp);
    fgets(field, 300, in); sscanf(field, "%d", &cObject);  if(cObject>maxCol) maxCol = cObject;
    fgets(field, 300, in); sscanf(field, "%d", &cFilter);  if(cFilter>maxCol) maxCol = cFilter;
    fgets(field, 300, in); sscanf(field, "%d", &cExposure);if(cExposure>maxCol) maxCol = cExposure;
    fgets(field, 300, in); sscanf(field, "%d", &cUTObs);   if(cUTObs>maxCol) maxCol = cUTObs;
    fgets(field, 300, in); sscanf(field, "%d", &tol);
    fclose(in);
    if(maxCol>MAXCOL) {  puts("Column number in the setting file is larger than maximum\n");  return 0;  }

    if( (in = fopen(seqfile, "r") ) == NULL) {  puts("Users' sequence file read error");  return 0;  }
    if( (out = fopen(outfile, "w") ) == NULL) {  puts("Script file write error");  return 0;  }
    fprintf(out, "# KMTNet Observation script\n");
    fprintf(out, "#   LABEL         RA          DEC     COPT  IMGTYP   OBJECT_NAME FILTER EXPTIME      UTOBS          UTTOL #Comments\n");
    fprintf(out, "#------------ ----------- ----------- ---- -------- ------------ ------ ------- ------------------- ----- ---------\n\n");
    for(;;) {
        for(i=0; i<MAXCOL; i++) strcpy(col[i], "-");
        for(i=0; i<MAXCOL; i++) {                       /* read column */
            do {  eol = fgetc(in);  } while(eol==' ');
                 if(eol=='\n' || eol=='\r') break;
            else if(eol=='#' || eol=='=') {  fgets(field, 300, in); break;  }
            else {  fseek(in, (long)sizeof(char)*(-1), SEEK_CUR); fscanf(in, "%s", col[i+1]);  }
        }
        if(feof(in)) break;
        if(i==0) continue;  /* for blank line */

        sprintf(prt[1], "%s", col[cRA]);
        sprintf(prt[2], "%s", col[cDEC]);
        sprintf(prt[3], "%d", copt);
        sprintf(prt[4], "%s", imgtyp);
        sprintf(prt[5], "%s", col[cObject]);
        sprintf(prt[8], "%s", col[cUTObs]); 
        sprintf(prt[9], "%d", tol);
        if(!strcmp(pID, "NEO")) sprintf(prt[8], "%s-%s-%sT%c%c:%c%c:%c%c", 
            col[cUTObs-3],col[cUTObs-2],col[cUTObs-1], col[cUTObs][0],col[cUTObs][1],col[cUTObs][2],col[cUTObs][3],col[cUTObs][4],col[cUTObs][5]);

        if(nfil!=0) 
            for(j=0; j<nfil; j++) {
                sprintf(prt[0], "%s(%s)", col[cLabel], filter[j]);
                sprintf(prt[6], "%s",     filter[j]);
                sprintf(prt[7], "%.0f",   exptime[j]);
                Write_Script(out, nfil, nrun+1, prt);
                nrun++;
            }
        else {
            sprintf(prt[0], "%s(%s)", col[cLabel], col[cFilter]);
            sprintf(prt[6], "%s",     col[cFilter]);
            sprintf(prt[7], "%s",     col[cExposure]);
            if(!strcmp(pID, "BLG")) sprintf(prt[0], "%s(%s,%s)", col[cLabel], col[cFilter], col[7]);
            Write_Script(out, nfil, nrun+1, prt);
            nrun++;
        }
    }
    if(flagloop > 0) fprintf(out, "\n+ostart %d\n", flagloop); 

    fclose(out);
    fclose(in);

    return 1;
}