JSOC_Documentation: proj/lev0/apps/hmi_time

00001 /* block of code to do the HMI shutter time calculations
00002    and set some keywords.  This block sets:
00003     T_OBS
00004     EXPTIME
00005     EXPSDEV
00006     DATE__OBS
00007 
00008  Call: HMI_compute_exposure_times(DRMS_Record_t *rec, HK_Keyword_t *isp, int flg)
00009 
00010  rec is lev0 record just before ISP keywords are added.
00011  isp is lev0 ISP linked-list from the ISP HK telem decode.
00012 
00013 */
00014 
00015 #include <drms_keyword.h>
00016 #include "packets.h"
00017 #include "printk.h"
00018 
00019 extern int INVALtime;
00020 
00021 static int HK_getkey_int(HK_Keyword_t *isp, char *key)
00022 {
00023 while (isp)
00024   {
00025   if (strcmp(key,isp->fitsname)==0)
00026      return((int)isp->raw_value);
00027   isp = isp->next;
00028   }
00029 return(DRMS_MISSING_INT);
00030 }
00031 
00032 static TIME SDO_to_DRMS_time(int sdo_s, int sdo_ss)
00033 {
00034 static int firstcall = 1;
00035 static TIME sdo_epoch;
00036 if (firstcall)
00037   {
00038   firstcall = 0;
00039   sdo_epoch = sscan_time("1958.01.01_00:00:00_TAI");
00040   }
00041 /* XXX fix build 3/18/2008, arta */
00042 return(sdo_epoch + (TIME)sdo_s + (TIME)(sdo_ss & 0xFFFF)/65536.0);
00043 }
00044 
00045 static void sprint_time_ISO (char *tstring, TIME t)
00046 {
00047 sprint_time(tstring,t,"UTC",0);
00048 tstring[4] = tstring[7] = '-';
00049 tstring[10] = 'T';
00050 tstring[19] = '\0';
00051 } 
00052 
00053 /* Two ways to calculate counter rollovers. Take your pick.
00054 static int nrollct(int cmdt, int clt)
00055 {
00056   if (cmdt <  50332) return 0;
00057   if (cmdt <  83886) return ((clt > 0x800000) ? 0 : 1);
00058   if (cmdt < 117441) return 1;
00059   if (cmdt < 150995) return ((clt > 0x800000) ? 1 : 2);
00060   if (cmdt < 184550) return 2;
00061   if (cmdt < 218104) return ((clt > 0x800000) ? 2 : 3);
00062   if (cmdt < 251659) return 3;
00063   return ((clt > 0x800000) ? 3 : 4);
00064 }
00065 */
00066 static int nrollct(int cmdt, int clt)
00067 {
00068   if (cmdt < 50332) return 0;
00069   int nr = (cmdt - 16777)/67109;
00070   int rem = cmdt%67109;
00071   if ((rem < 16777) || (rem > 50332)) return nr + (clt > 0x800000 ? 0 : 1);
00072   return nr;
00073 }
00074 
00075 void HMI_compute_exposure_times(DRMS_Record_t *rec, HK_Keyword_t *isp, int flg)
00076 {
00077 static int waveltab[10] = { 335, 131, 211, 193, 1600, 1700, 4500, 171, 304, 94 };
00078 
00079 static char date_obs[100]; 
00080 //static char datestr[100];
00081 TIME t_obs,MJD_epoch = -3727641600.000; /* 1858.11.17_00:00:00_UT  */
00082 double exptime, expsdev;
00083 float int_time;
00084 
00085 if (flg) {  /* AIA */
00086   char *aia_instru[] = { "AIA_1", "AIA_2", "AIA_3", "AIA_4" };
00087   int aimgshce = HK_getkey_int(isp, "AIMGSHCE"); /* AIA_IMG_SH_CMDED_EXPOSURE */
00088   int aimgots  = HK_getkey_int(isp, "AIMGOTS");  /* AIA_IMG_OBT_TIME_SH_SEC  */
00089   if((aimgots == 0) && (aimgshce != 0)) INVALtime = 1;//used by do_quallev0()
00090   else INVALtime = 0;
00091   int aimgotss = HK_getkey_int(isp, "AIMGOTSS"); /* AIA_IMG_OBT_TIME_SH_SS   */
00092   int aimshcbc = HK_getkey_int(isp, "AIMSHCBC"); /* AIA_IMG_SH_CLOSE_BOT_CENTR*/
00093   aimshcbc += 0x1000000*nrollct(aimgshce, aimshcbc);
00094   int aimshcbe = HK_getkey_int(isp, "AIMSHCBE"); /* AIA_IMG_SH_CLOSE_BOT_EDGE */
00095   aimshcbe += 0x1000000*nrollct(aimgshce, aimshcbe);
00096   int aimshctc = HK_getkey_int(isp, "AIMSHCTC"); /* AIA_IMG_SH_CLOSE_TOP_CENTR*/
00097   aimshctc += 0x1000000*nrollct(aimgshce, aimshctc);
00098   int aimshcte = HK_getkey_int(isp, "AIMSHCTE"); /* AIA_IMG_SH_CLOSE_TOP_EDGE */
00099   aimshcte += 0x1000000*nrollct(aimgshce, aimshcte);
00100   int aimshobc = HK_getkey_int(isp, "AIMSHOBC"); /* AIA_IMG_SH_OPEN_BOT_CENTR*/
00101   int aimshobe = HK_getkey_int(isp, "AIMSHOBE"); /* AIA_IMG_SH_OPEN_BOT_EDGE */
00102   int aimshotc = HK_getkey_int(isp, "AIMSHOTC"); /* AIA_IMG_SH_OPEN_TOP_CENTR*/
00103   int aimshote = HK_getkey_int(isp, "AIMSHOTE"); /* AIA_IMG_SH_OPEN_TOP_EDGE */
00104   int aicfgdl3 = HK_getkey_int(isp, "AICFGDL3");
00105   int aicfgdl4 = HK_getkey_int(isp, "AICFGDL4");
00106   int aiagp9   = HK_getkey_int(isp, "AIAGP9");
00107   int aiagp10   = HK_getkey_int(isp, "AIAGP10");
00108   
00109   if ((aicfgdl3 == DRMS_MISSING_SHORT) || (aicfgdl4 == DRMS_MISSING_SHORT)) {
00110     int_time = DRMS_MISSING_FLOAT; /* no info to calculate integration time */
00111   } else if (aiagp9 && aiagp10) {  /* non wrapped cfg delays are available */
00112     int_time = (float)(aiagp10 - aiagp9)/128.0; /* this value already seconds */
00113   } else { /* have to use potentially wrapped around configuration delays */
00114     int nroll = 256*(int)((aimgshce + 273)/2000);
00115     int_time = ((float)(aicfgdl4 - aicfgdl3 + nroll)*1000.0)/128.0;
00116     if((int)int_time < aimgshce) int_time = int_time + 2000.0;
00117     int_time = int_time/1000.0;     //make seconds
00118   }
00119   drms_setkey_float(rec, "INT_TIME", int_time);
00120 
00121   double shebc = (aimshcbc - aimshobc)*4.0e-6;
00122   double shebe = (aimshcbe - aimshobe)*4.0e-6;
00123   double shetc = (aimshctc - aimshotc)*4.0e-6;
00124   double shete = (aimshcte - aimshote)*4.0e-6;
00125   double offset = ( aimshobc + aimshcbc + aimshobe + aimshcbe +
00126                     aimshotc + aimshctc + aimshote + aimshcte )/8.0e6;
00127   if(aimgshce == 0) {           //use pkt time for t_obs
00128     drms_setkey_string(rec, "IMG_TYPE", "DARK");
00129     int axsec = HK_getkey_int(isp, "ATCSISP");
00130     int axssec = HK_getkey_int(isp, "ATCSSISP");
00131     if((axsec == DRMS_MISSING_INT) || (axssec == DRMS_MISSING_INT)) {
00132       axsec = HK_getkey_int(isp, "ATCS211");
00133       axssec = HK_getkey_int(isp, "ATCSS211");
00134     }
00135     if((axsec == DRMS_MISSING_INT) || (axssec == DRMS_MISSING_INT)) {
00136       axsec = HK_getkey_int(isp, "ATCS239");
00137       axssec = HK_getkey_int(isp, "ATCSS239");
00138     }
00139     if((axsec != DRMS_MISSING_INT) && (axssec != DRMS_MISSING_INT)) {
00140       axssec = axssec >> 16;
00141       t_obs = SDO_to_DRMS_time(axsec, axssec);
00142     }
00143   }
00144   else {
00145     drms_setkey_string(rec, "IMG_TYPE", "LIGHT");
00146     t_obs = SDO_to_DRMS_time(aimgots, aimgotss) + offset;
00147   }
00148   exptime = (shebc + shebe + shetc + shete)/4.0;
00149   if (aimgshce < 73) exptime = exptime*7.0/20.0; //from jps 2Jun2010
00150   expsdev = sqrt((shebc*shebc + shebe*shebe + shetc*shetc +
00151                  shete*shete)/4.0 - exptime*exptime);
00152   int asqhdr = HK_getkey_int(isp, "ASQHDR");
00153   if (DRMS_MISSING_INT == asqhdr) asqhdr = HK_getkey_int(isp, "A827A");
00154   int ahtelid = (asqhdr >> 30) & 0x3;
00155   int ahtlfsn = asqhdr & 0x3FFFFFF;
00156   int asqtnum = HK_getkey_int(isp, "ASQTNUM");
00157   drms_setkey_string(rec, "INSTRUME", aia_instru[ahtelid]);
00158 //  CAMERA set elsewhere, as is FSN
00159 //  drms_setkey_int(rec, "CAMERA", ahtelid + 1);
00160   int wavel = HK_getkey_int(isp, "AIAWVLEN");    /* AIA_IMG_WAVELENGTH */
00161   if (wavel == DRMS_MISSING_INT) wavel = HK_getkey_int(isp, "A856C");
00162   if (wavel != DRMS_MISSING_INT) {
00163     if (wavel < 10) drms_setkey_int(rec, "WAVELNTH", waveltab[wavel]);
00164     else drms_setkey_int(rec, "WAVELNTH", DRMS_MISSING_INT);
00165   }
00166 
00167   int aecmode = HK_getkey_int(isp, "AECMODE");   /* AIA_IMG_AEC_MODE         */
00168   if (aecmode == DRMS_MISSING_INT) aecmode = HK_getkey_int(isp, "A8328");
00169   if (aecmode == DRMS_MISSING_INT) aecmode = HK_getkey_int(isp, "A8215-00");
00170   if (aecmode != DRMS_MISSING_INT) {
00171     if (aecmode) drms_setkey_string(rec, "A8215_00", "On");
00172     else drms_setkey_string(rec, "A8215_00", "Off");
00173   }
00174 } else {    /* HMI */
00175   int hshiexp = HK_getkey_int(isp, "HSHIEXP"); /* HMI_FSW_IMG_CMDED_EXPOSURE */
00176   int hobitsec = HK_getkey_int(isp, "HOBITSEC"); /* HMI_OBT_IMG_TIME_SHM_SEC */
00177   if((hobitsec == 0) && (hshiexp != 0)) INVALtime = 1;//used by do_quallev0()
00178   else INVALtime = 0;
00179   int hobitss  = HK_getkey_int(isp, "HOBITSS");  /* HMI_OBT_IMG_TIME_SHM_SS  */
00180   int hshmiclb = HK_getkey_int(isp, "HSHMICLB"); /* HMI_SHM_IMG_CLOSE_BOTTOM */
00181   int hshmiclm = HK_getkey_int(isp, "HSHMICLM"); /* HMI_SHM_IMG_CLOSE_MIDDLE */
00182   int hshmiclt = HK_getkey_int(isp, "HSHMICLT"); /* HMI_SHM_IMG_CLOSE_TOP    */
00183   int hshmiopb = HK_getkey_int(isp, "HSHMIOPB"); /* HMI_SHM_IMG_OPEN_BOTTOM  */
00184   int hshmiopm = HK_getkey_int(isp, "HSHMIOPM"); /* HMI_SHM_IMG_OPEN_MIDDLE  */
00185   int hshmiopt = HK_getkey_int(isp, "HSHMIOPT"); /* HMI_SHM_IMG_OPEN_TOP     */
00186   double sheb = (hshmiclb - hshmiopb)*1.0e-6;
00187   double shem = (hshmiclm - hshmiopm)*1.0e-6;
00188   double shet = (hshmiclt - hshmiopt)*1.0e-6;
00189   double offset = ( hshmiclb + hshmiopb + hshmiclm + hshmiopm + hshmiclt + hshmiopt )/6.0e6;
00190   
00191   if(hshiexp == 0) {            //use pkt time for t_obs
00192     drms_setkey_string(rec, "IMG_TYPE", "DARK");
00193     int hxsec = HK_getkey_int(isp, "HTCSISP");
00194     int hxssec = HK_getkey_int(isp, "HTCSSISP");
00195     if((hxsec == DRMS_MISSING_INT) || (hxssec == DRMS_MISSING_INT)) {
00196       hxsec = HK_getkey_int(isp, "HTCS1BD");
00197       hxssec = HK_getkey_int(isp, "HTCSS1BD");
00198     }
00199     if((hxsec == DRMS_MISSING_INT) || (hxssec == DRMS_MISSING_INT)) {
00200       hxsec = HK_getkey_int(isp, "HTCS1DB");
00201       hxssec = HK_getkey_int(isp, "HTCSS1DB");
00202     }
00203     if((hxsec != DRMS_MISSING_INT) && (hxssec != DRMS_MISSING_INT)) {
00204       hxssec = hxssec >> 16;
00205       hxsec = hxsec + 1;        //add to dark 26July2010
00206       t_obs = SDO_to_DRMS_time(hxsec, hxssec);
00207     }
00208   }
00209   else {
00210     drms_setkey_string(rec, "IMG_TYPE", "LIGHT");
00211     t_obs = SDO_to_DRMS_time(hobitsec,hobitss) + offset;
00212   }
00213   exptime = (sheb + shem + shet)/3.0;
00214   expsdev = sqrt((sheb*sheb + shem*shem + shet*shet)/3.0 - exptime*exptime);
00215   
00216   int camid = HK_getkey_int(isp, "HCAMID");
00217   camid = camid & 0x01;
00218   camid++;
00219   drms_setkey_int(rec, "CAMERA", camid);
00220   if(camid == 1) drms_setkey_string(rec, "INSTRUME", "HMI_SIDE1");
00221   else drms_setkey_string(rec, "INSTRUME", "HMI_FRONT2");
00222 }
00223   
00224 //sprint_time_ISO(date_obs, t_obs - exptime/2.0);
00225 TIME date__obs = t_obs - exptime/2.0;
00226 TIME mjd = date__obs - MJD_epoch;
00227 double  mjd_day = floor(mjd / 86400.0);
00228 double  mjd_time = mjd - 86400.0 * mjd_day;
00229 //sprint_time(datestr, CURRENT_SYSTEM_TIME, "ISO", 0);
00230 //sprint_time(date_obs, date__obs, "ISO", 0);
00231   
00232 drms_setkey_double(rec, "T_OBS", t_obs);
00233 drms_setkey_double(rec, "EXPTIME", exptime);
00234 drms_setkey_float(rec, "EXPSDEV", expsdev);
00235 //drms_setkey_string(rec, "DATE__OBS", date_obs);
00236 drms_setkey_double(rec, "DATE__OBS", date__obs);
00237 drms_setkey_double(rec, "MJD", mjd_day);
00238 drms_setkey_double(rec, "TIME", mjd_time);
00239 //NOTE: DATE now set in close_image() in ingest_lev0.c
00240 //drms_setkey_string(rec, "DATE", datestr);
00241 //drms_setkey_double(rec, "DATE", CURRENT_SYSTEM_TIME);
00242 }
00243
proj/lev0/apps/hmi_time_setting.c
