int MDSPlusData::getECE_Center(int shotnumber)
{
int dtype_long = DTYPE_LONG;
string treename = "\\ECE_CENT";
#if 1
int status = MdsOpen("electron", &shotnumber);
if ( !status_ok(status) )
{
fprintf(stderr,"Error opening electron tree for shot %d\n",shotnumber);
return -1;
};
#endif
int len = 0;
int ece_center;
int null = 0;
/* create a descriptor for this signal */
int sigdesc = descr(&dtype_long,&ece_center,&null,&len);
/* retrieve signal */
status = MdsValue((char*)treename.c_str(), &sigdesc, &null, &len );
if ( !status_ok(status) )
{
/* error */
fprintf(stderr,"Error retrieving signal**ECE_CENT\n");
return -1;
};
MdsClose((char*)string("electron").c_str(), &shotnumber);
return (ece_center);
}
string MDSPlusData::getShotTime()
{
/* use get shot time */
int dtype_str = DTYPE_CSTRING;
string ftreename = "\\T0_STR";
#if 1
int status = MdsOpen("operation", &mshot);
if ( !status_ok(status) )
{
fprintf(stderr,"Error opening rdata tree for shot %d\n",mshot);
return "";
};
#endif
//int len = 26;
//char vCstring[26]; //format "0000/00/00 00:00:00.000000";
int len = 19;
char vCstring[20]; //format "0000/00/00 00:00:00.000000";
int null = 0;
/* create a descriptor for this signal */
int sigdesc = descr(&dtype_str,vCstring,&null,&len);
/* retrieve signal */
status = MdsValue((char*)ftreename.c_str(), &sigdesc, &null, &len );
if ( !status_ok(status) )
{
/* error */
fprintf(stderr,"Error retrieving signal**1\n");
return "";
};
MdsClose((char*)ftreename.c_str(), &mshot);
vCstring[20]='\0';
return string(vCstring);
}
int mdsplus_Tree_create(int storeMode, int shotNumber)
{
char buf[50];
int null = 0;
int status;
int dtype_float = DTYPE_FLOAT;
int dtypeLong = DTYPE_LONG;
int tstat, len;
int idesc = descr(&dtypeLong, &tstat, &null);
if(storeMode == 0){
sprintf(buf, "TCL(\"set tree %s/shot=-1\")",TREE);
status = MdsValue(buf, &idesc, &null, &len);
if (!status_ok(status))
{
fprintf(stderr,"Error TCL Command Set tree -1 \n");
return -1;
}
sprintf(buf, "TCL(\"create pulse %ld\")",shotNumber);
status = MdsValue(buf, &idesc, &null, &len);
if (!status_ok(status))
{
fprintf(stderr,"Error TCL Command Create Pulse \n");
return -1;
}
sprintf(buf, "TCL(\"close\")",TREE);
status = MdsValue(buf, &idesc, &null, &len);
}
return 1;
}
int MDSPlusData::getShotData(double *timebase, double *shotdata, int size)
{
/* use get_signal_length to get size of signal */
int dtype_float = DTYPE_DOUBLE;
string ftreename = string("\\")+mtreename;
int status = MdsOpen((char*)mnodename.c_str(), &mshot);
if ( !status_ok(status) )
{
fprintf(stderr,"Error opening rdata tree for shot %d\n",mshot);
return -1;
};
int null = 0;
/* create a descriptor for this signal */
int sigdesc = descr(&dtype_float, shotdata, &size, &null);
int len =0;
/* retrieve signal */
status = MdsValue((char*)ftreename.c_str(), &sigdesc, &null, &len );
if ( !status_ok(status) )
{
/* error */
fprintf(stderr,"Error retrieving signal**2\n");
return -1;
}
/* create a descriptor for the timebase */
int timedesc = descr(&dtype_float, timebase, &size, &null);
/* retrieve timebase of signal */
string strtree_dim= string("DIM_OF(")+ftreename+string(")");
status = MdsValue((char*)strtree_dim.c_str(), &timedesc, &null, 0);
if ( !status_ok(status) )
{
/* error */
fprintf(stderr,"Error retrieving timebase\n");
return -1;
};
MdsClose((char*)mnodename.c_str(), &mshot);
#if 0
// possible function end.
for ( int i = 0 ; i < size; i++ )
{
printf("%i X:%f Y:%f\n", i, timebase[i], shotdata[i]);
};
#endif
//double mxVal = *max_element(&mshotdata[0], &mshotdata[msize]);
//double mnVal = *min_element(&mshotdata[0], &mshotdata[msize]);
//printf("MAX:%f,Min:%f,Diff:%f\n", mxVal, mnVal,mxVal-mnVal);
return 0;
}
double MDSPlusData::getNBIPower(int shotnumber)
{
int dtype_long = DTYPE_DOUBLE;
string treename = "\\NB11_PNB";
#if 1
int status = MdsOpen("heating", &shotnumber);
if ( !status_ok(status) )
{
fprintf(stderr,"Error opening heating tree for shot %d\n",shotnumber);
return -1.0;
};
#endif
int len = 0;
double nb11pnb = 0.0, nb12pnb = 0.0, nbipower = 0.0;
int null = 0;
/* create a descriptor for this signal */
int sigdesc = descr(&dtype_long,&nb11pnb,&null,&len);
/* retrieve signal */
status = MdsValue((char*)treename.c_str(), &sigdesc, &null, &len );
if ( !status_ok(status) )
{
/* error */
fprintf(stderr,"Error retrieving signal**NB11_PNB\n");
nb11pnb = 0.0;
};
treename = "\\NB12_PNB";
/* create a descriptor for this signal */
sigdesc = descr(&dtype_long,&nb12pnb,&null,&len);
/* retrieve signal */
status = MdsValue((char*)treename.c_str(), &sigdesc, &null, &len );
if ( !status_ok(status) )
{
/* error */
fprintf(stderr,"Error retrieving signal**NB12_PNB\n");
nb12pnb = 0.0;
};
nbipower = (nb11pnb+nb12pnb);
printf("NB11_PNB[%f] + NB12_PNB[%f] = [%f]\n", nb11pnb, nb12pnb, nbipower);
MdsClose((char*)string("heating").c_str(), &shotnumber);
return (nbipower);
};
static int controller_ready(unsigned int drive, unsigned int head)
{
int retry = 100;
do {
if (controller_busy() & BUSY_STAT)
return 0;
outb_p(0xA0 | (drive<<4) | head, HD_CURRENT);
if (status_ok())
return 1;
} while (--retry);
return 0;
}
int MDSPlusData::getNodeLen()
{
/* local vars */
int dtype_long = DTYPE_LONG;
char buf[1024];
int size;
int null = 0;
int idesc = descr(&dtype_long, &size, &null);
/* init buffer */
memset(buf,0,sizeof(buf));
string ftreename = string("\\")+mtreename;
/* put SIZE() TDI function around tree name */
snprintf(buf,sizeof(buf)-1,"SIZE(%s)",ftreename.c_str());
//printf("node:%s, shot:%d\n", mnodename.c_str(), mshot);
int status = MdsOpen((char*)mnodename.c_str(), &mshot);
if ( !status_ok(status) )
{
fprintf(stderr,"Error opening rdata tree for [node:%s, tree:%s, shot:%d]\n",mnodename.c_str(), mtreename.c_str(), mshot);
return -1;
};
//printf("buf:%s, shot:%d\n", buf, mshot);
/* use MdsValue to get the tree length */
status = MdsValue(buf, &idesc, &null, 0);
if ( !status_ok(status) )
{
fprintf(stderr,"Unable to get length of %s.\n",ftreename.c_str());
return -1;
};
/* return tree length */
MdsClose((char*)mnodename.c_str(), &mshot);
return size;
}
void pars_whohere(void) {
strncpy(status->nick, packet, NICK_MAXLEN);
make_imhere();
status_ok(PARS_SHERE);
}
void pars_beepreq(void) { /* FIXME: who needs this??? */
status_ok(PARS_SBEEPREQ);
}
void mdsPlusPut_Task(int param)
{
unsigned long i;
int socket;
int null=0;
int dtypeFloat = DTYPE_FLOAT;
int dataDesc,timeDesc;
int status=1, j;
float tdT;
float toffT;
int tnoRd;
int tshotID;
char ttagName[20];
char ttreeID[20];
char tserverID[20];
float *timeBase;
DBADDR *paddr;
paddr = (DBADDR *)dbCalloc(1, sizeof(struct dbAddr));
if(genSubDebug > 0)
printf("Start MdsPut. ##############.\n");
status = dbNameToAddr(mpActivePV, paddr);
enum16Val = 1;
status = dbPutField(paddr, DBR_ENUM, &enum16Val, 1);
status = mdsPlusCreatePulse();
/* loop numNode times */
for(j = 0;j < numNode; j++) {
if(genSubDebug > 5)
printf("genSub: mdsPlusPut_Task() started. Task Position=%d, Flag=%ld\n",j,pmdsPutData[j].putFlag);
if(pmdsPutData[j].putFlag == 0)
goto endloop;
epicsThreadSleep(0.01);
tnoRd=(int)(pmdsPutData[j].noRd);
tshotID=(int)(pmdsPutData[j].shotID);
strcpy(ttagName,pmdsPutData[j].tagName);
strcpy(ttreeID,pmdsPutData[j].treeID);
strcpy(tserverID,pmdsPutData[j].serverID);
tdT=pmdsPutData[j].dT;
toffT=pmdsPutData[j].offT;
if(genSubDebug > 3)
printf("genSub: mdsPlusPut_Task() MdsConnect Ready. ServerID=[%s] TreeID=[%s] tagName=[%s] dt=[%f] off=[%f] noRd=[%d]\n",tserverID, ttreeID, ttagName, tdT, toffT, tnoRd);
/* Connect server and open tree*/
socket=MdsConnect(tserverID);
if(socket == -1) {
printf("genSub: Error connecting to mdsip server[%s].\n",tserverID);
break;
}
if(genSubDebug > 5)
printf("genSub: mdsPlusPut_Task() MdsConnected[%s]\n",tserverID);
status=MdsOpen(ttreeID, &tshotID);
if( !status_ok(status) ) {
printf("genSub: Error opening tree [%s] for shot [%d].\n",
ttreeID,tshotID );
break;
}
if(genSubDebug > 5)
printf("genSub: mdsPlusPut_Task() MdsOpened [%s] shot number [%d].\n",
ttreeID,tshotID);
/* put data */
timeBase = (float *)malloc(tnoRd*sizeof(float));
for(i=0;i<tnoRd;i++)
*(timeBase + i) = ((float)i)*tdT + toffT;
dataDesc=descr(&dtypeFloat,pmdsPutData[j].dataArray,&tnoRd, &null);
timeDesc=descr(&dtypeFloat,timeBase, &tnoRd, &null);
status=MdsPut(ttagName,"BUILD_SIGNAL($1,,$2)",&dataDesc,&timeDesc,&null);
if( !status_ok(status) ) {
printf("genSub: Error writing signal.\n");
break;
}
if(genSubDebug > 5)
printf("genSub: mdsPlusPut_Task() MdsPutted to tag [%s]. shot number [%d], noRd=[%d].\n",
ttagName,tshotID,tnoRd);
if(genSubDebug > 10)
for(i=0;i<tnoRd;i++) {
printf("timeBase=%f,data=%f\n",
*((float *)timeBase+i),*((float *)pmdsPutData[j].dataArray+i));
}
free(timeBase);
status=MdsClose(ttreeID, &tshotID);
if( !status_ok(status) ) {
printf("genSub: Error closing tree for shot [%d].\n",tshotID );
break;
}
if(genSubDebug > 5)
printf("genSub: mdsPlusPut_Task() MdsClosed [%s] shot number [%d]\n",ttreeID,tshotID);
endloop:
if(genSubDebug > 5)
printf("genSub: mdsPlusPut_Task() Data discarded for taskPos[%d].\n",j);
}
status = mdsPlusPrepNext();
/* end of mdsput */
startMdsPut = 0;
status = dbNameToAddr(mpActivePV, paddr);
enum16Val = 0;
status = dbPutField(paddr, DBR_ENUM, &enum16Val, 1);
free(paddr);
}
long mdsPlusPrepNext()
{
int socket;
int null=0;
int dtype_long = DTYPE_LONG;
int mdsstat, len;
int idesc;
int status=1,i;
int bufsize=40;
char buf[bufsize];
char treeID[20];
char serverID[20];
int shotID = 0;
DBADDR *paddr;
long options, nRequest;
paddr = (DBADDR *)dbCalloc(1, sizeof(struct dbAddr));
/* Get server information */
sprintf(buf, "icrf:daq:mptree:i.VAL");
status = dbNameToAddr(buf, paddr);
options = 0;
nRequest = 1;
status = dbGet(paddr, DBR_STRING, buf, &options, &nRequest, NULL);
i=sscanf(buf,"%s %s",serverID, treeID);
/* Connect server */
socket=MdsConnect(serverID);
if(socket == -1) {
printf("genSub: Error connecting to mdsip server[%s].\n",serverID);
return(-1);
}
if(genSubDebug > 5)
printf("genSub: mdsPlusPrepNext() MdsConnected[%s]\n",serverID);
/* Get pulse id */
sprintf(buf, "icrf:pulseid.VAL");
status = dbNameToAddr(buf, paddr);
options = 0;
nRequest = 1;
status = dbGetField(paddr, DBR_LONG, &shotID, &options, &nRequest, NULL);
/* open MDSPlus*/
status=MdsOpen(treeID, &shotID);
if( !status_ok(status) ) {
printf("genSub: Error opening tree [%s] for shot [%d].\n",
treeID,shotID );
return(-1);
}
if(genSubDebug > 5)
printf("genSub: mdsPlusPrepNext() MdsOpened [%s] shot number [%d].\n",
treeID,shotID);
/* generate MDSPlus event*/
sprintf(buf, "TCL(\"SETEVENT cmpl\")");
idesc = descr(&dtype_long, &mdsstat, &null);
status = MdsValue(buf,&idesc,&null,&len);
if( !status_ok(status) ) {
printf("genSub: Error generating event for tree [%s].\n", treeID);
return(-1);
}
if(genSubDebug > 0)
printf("genSub: mdsPlusPrepNext() MdsValue %s\n", buf);
/* Put next pulse number */
sprintf(buf, "icrf:pulseid.VAL");
status = dbNameToAddr(buf, paddr);
options = 0;
nRequest = 1;
shotID++;
status = dbPutField(paddr, DBR_LONG, &shotID, 1);
status=MdsClose(treeID, &shotID);
free(paddr);
return(0);
}
long mdsPlusCreatePulse()
{
int socket;
int null=0;
int dtype_long = DTYPE_LONG;
int dtype_cstring = DTYPE_CSTRING;
int dtype_float = DTYPE_FLOAT;
int mdsstat, len;
int idesc;
int sdesc;
int fdesc;
int status=1,i;
int bufsize=40;
char buf[bufsize];
char treeID[20];
char serverID[20];
int shotID = 0;
int tmpShotID = 0;
int nextShotID;
int tok_pulse;
float coef=1.0;
DBADDR *paddr;
long options, nRequest;
paddr = (DBADDR *)dbCalloc(1, sizeof(struct dbAddr));
/* Get server information */
sprintf(buf, "icrf:daq:mptree:i.VAL");
status = dbNameToAddr(buf, paddr);
options = 0;
nRequest = 1;
status = dbGetField(paddr, DBR_STRING, buf, &options, &nRequest, NULL);
i=sscanf(buf,"%s %s",serverID, treeID);
/* Connect server */
/* socket=MdsConnect(serverID);
if(socket == -1) {
printf("genSub: Error connecting to mdsip server[%s].\n",serverID);
return(-1);
}
if(genSubDebug > 5)
printf("genSub: mdsPlusCreatePulse() MdsConnected[%s]\n",serverID);
*/
/* Get pulse id */
/* sprintf(buf, "icrf:pulseid.VAL");*/
/* sprintf(buf, "icrf:daq:mcont:pulseid.VAL");
status = dbNameToAddr(buf, paddr);
options = 0;
nRequest = 1;
status = dbGetField(paddr, DBR_LONG, &nextShotID, &options, &nRequest, NULL);
*/ /* open MDSPlus */
/* should be opened with shot -1 */
/* shotID = -1;
status=MdsOpen(treeID, &shotID);
if( !status_ok(status) ) {
printf("genSub: Error opening tree [%s] for shot [%d].\n",
treeID,shotID );
return(-1);
}
if(genSubDebug > 5)
printf("genSub: mdsPlusCreatePulse() MdsOpened [%s] shot number [%d].\n",
treeID,shotID);
*/ /* create new pulse */
/* shot increment */
/* sprintf(buf, "TCL(\"SET CURRENT %s /INCREMENT\")", treeID);
idesc = descr(&dtype_long, &mdsstat, &null);
status = MdsValue(buf,&idesc,&null,&len);
if(genSubDebug > 0)
printf("genSub: mdsPlusCreatePulse() MdsValue %s\n", buf);
*/
/* get increased current shot id. increased shot ID should be same to the nextShotID */
/* sprintf(buf, "TCL(\"SHOW CURRENT %s\",_output)", treeID);
idesc = descr(&dtype_long, &mdsstat, &null);
status = MdsValue(buf,&idesc,&null,&len);
if( !status_ok(status) ) {
printf("genSub: Error with %s.", buf);
return (-1);
}
if(genSubDebug > 3)
printf("genSub: mdsPlusCreatePulse() MdsValue %s\n",buf);
if( status_ok(mdsstat) ) {
sdesc = descr(&dtype_cstring, buf, &null, &bufsize);
status = MdsValue("_output",&sdesc, &null, &len);
if( !status_ok(status) ) {
printf("genSub: Error getting output with %s.",buf);
}
sscanf(buf,"Current shot is %d", &tmpShotID);
if(tmpShotID != nextShotID ) {
printf("genSub: Shot ID is incorrect.\n");
return(-1);
}
if(genSubDebug > 0)
printf("genSub: mdsPlusCreatePulse() MdsValue SHOW CURRENT gets shot number [%d].\n",
tmpShotID);
*/ /* create current pulse */
/* sprintf(buf, "TCL(\"CREATE PULSE %d\")", tmpShotID);
status = MdsValue(buf,&idesc,&null,&len);
if(genSubDebug > 3)
printf("genSub: mdsPlusCreatePulse() MdsValue %s\n",buf);
}
*/ /* close and reopen with increased shot id*/
/* status=MdsClose(treeID, &shotID);*/
/* shotID=nextShotID;*/
/* Connect server */
socket=MdsConnect(serverID);
if(socket == -1) {
printf("genSub: Error connecting to mdsip server[%s].\n",serverID);
return(-1);
}
if(genSubDebug > 5)
printf("genSub: mdsPlusCreatePulse() MdsConnected[%s]\n",serverID);
/* open MDSPlus */
shotID = 0;
status=MdsOpen(treeID, &shotID);
if( !status_ok(status) ) {
printf("genSub: Error opening tree [%s] for shot [%d].\n",
treeID,shotID );
return(-1);
}
if(genSubDebug > 5)
printf("genSub: mdsPlusCreatePulse() MdsOpened [%s] shot number [%d].\n",
treeID,shotID);
/* write data at current(increased) shot id*/
/* Get tokamak shot ID */
/* sprintf(buf, "icrf:shotid.VAL");
status = dbNameToAddr(buf, paddr);
options = 0;
nRequest = 1;
status = dbGetField(paddr, DBR_LONG, &tok_pulse, &options, &nRequest, NULL);
*/ /* write tokamak shot id */
/* idesc = descr(&dtype_long, &tok_pulse, &null);
status = MdsPut("\\TOK_SHOT", "$", &idesc, &null);
if(genSubDebug > 3)
printf("genSub: mdsPlusCreatePulse() MdsPut [\\TOK_SHOT=%d]\n",tok_pulse);
*/
/* Get current time string */
/* sprintf(buf, "icrf:ioc:time.VAL");
status = dbNameToAddr(buf, paddr);
options = 0;
nRequest = 1;
status = dbGetField(paddr, DBR_STRING, buf, &options, &nRequest, NULL);
*/ /* write time */
/* sdesc = descr(&dtype_cstring, buf, &null, &bufsize);
status = MdsPut("\\PULSE_TIME", "$", &sdesc, &null);
if(genSubDebug > 3)
printf("genSub: mdsPlusCreatePulse() MdsPut [\\PULSE_TIME=%s]\n",buf);
*/ /* write coefficient */
sprintf(buf, "icrf:daq:rfd:coef0.VAL");
status = dbNameToAddr(buf, paddr);
status = dbGetField(paddr, DBR_FLOAT, &coef, &options, &nRequest, NULL);
fdesc = descr(&dtype_float, &coef, &null);
status = MdsPut("\\ICRF_0C", "$", &fdesc, &null);
sprintf(buf, "icrf:daq:rfd:coef1.VAL");
status = dbNameToAddr(buf, paddr);
status = dbGetField(paddr, DBR_FLOAT, &coef, &options, &nRequest, NULL);
fdesc = descr(&dtype_float, &coef, &null);
status = MdsPut("\\ICRF_1C", "$", &fdesc, &null);
sprintf(buf, "icrf:daq:rfd:coef2.VAL");
status = dbNameToAddr(buf, paddr);
status = dbGetField(paddr, DBR_FLOAT, &coef, &options, &nRequest, NULL);
fdesc = descr(&dtype_float, &coef, &null);
status = MdsPut("\\ICRF_2C", "$", &fdesc, &null);
sprintf(buf, "icrf:daq:rfd:coef3.VAL");
status = dbNameToAddr(buf, paddr);
status = dbGetField(paddr, DBR_FLOAT, &coef, &options, &nRequest, NULL);
fdesc = descr(&dtype_float, &coef, &null);
status = MdsPut("\\ICRF_3C", "$", &fdesc, &null);
status=MdsClose(treeID, &shotID);
free(paddr);
return(0);
}
long mdsPlusFindInit()
{
int socket;
int null=0;
int dtype_long = DTYPE_LONG;
int dtype_cstring = DTYPE_CSTRING;
int mbufsize = 1024;
char mbuf[1024];
int tstat, len;
int idesc = descr(&dtype_long, &tstat, &null);
int sdesc = descr(&dtype_cstring, mbuf, &null, &mbufsize);
int status=1,i;
char buf[40];
char treeID[20];
char serverID[20];
int initShotID = 0;
long currentShotID = 0;
if(mdsPlusManInit_first == 0) return(0);
mdsPlusManInit_first = 0;
sprintf(buf,"172.17.101.180:8000 adata");
i=sscanf(buf,"%s %s",serverID, treeID);
if(genSubDebug > 5)
printf("genSub: mdsPlusManInit() Trying MdsConnect[%s,%s]\n",serverID,treeID);
socket=MdsConnect(serverID);
if(socket == -1) {
printf("genSub: Error connecting to mdsip server[%s].\n",serverID);
} else {
if(genSubDebug > 5)
printf("genSub: mdsPlusManInit() MdsConnected[%s]\n",serverID);
status=MdsOpen(treeID, &initShotID);
if( !status_ok(status) ) {
printf("genSub: Error opening tree [%s] for shot [%d].\n",
treeID,initShotID );
} else {
if(genSubDebug > 5)
printf("genSub: mdsPlusManInit() MdsOpened [%s] shot number [%d].\n",
treeID,initShotID);
/* put data */
sprintf(mbuf, "TCL(\"SHOW CURRENT %s\",_output)", treeID);
status = MdsValue(mbuf,&idesc,&null,&len);
if( !status_ok(status) ) {
printf("genSub: Error with %s.\n", mbuf);
}
if(genSubDebug > 5)
printf("genSub: mdsPlusManInit() MdsValue with %s ok.\n",mbuf);
if( status_ok(tstat) ) {
status = MdsValue("_output",&sdesc, &null, &len);
if( !status_ok(status) ) {
printf("genSub: Error getting output with SHOW CURRENT.");
}
sscanf(mbuf,"Current shot is %ld", ¤tShotID);
if(genSubDebug > 0)
printf("genSub: mdsPlusManInit() MdsValue SHOW CURRENT gets shot number [%ld].\n",
currentShotID);
shotId = currentShotID + 1;
epicsThreadCreate("mdsPlusShotID",
epicsThreadPriorityLow,
epicsThreadGetStackSize(epicsThreadStackSmall),
(EPICSTHREADFUNC) mdsPlusShotID,
0);
}
}
status=MdsClose(treeID, &initShotID);
}
return(currentShotID+1);
}
int mdsplus_data_put_sif(char* fileNamePath, int shotNumber, int storeMode, float blipTime, float daqTime)
{
char channelTagName[5][15] = {"\\HR_VSS01:FOO", "\\HR_VSS02:FOO", "\\HR_VSS03:FOO", "\\HR_VSS04:FOO", "\\HR_VSS05:FOO"};
char channelNodeName1[7][22] = {"\\HR_INFO:IMAGE_NO", "\\HR_INFO:IMAGE_RIGHT", "\\HR_INFO:IMAGE_TOP", "\\HR_INFO:SUBIMAGE_NO", "\\HR_INFO:T_EXPOSURE", "\\HR_INFO:T_INTEGRATE", "\\HR_INFO:T_KINETIC"};
/* char channelNodeName2[6][22] = {"\\HR_VSS01:POS_BOTTOM", "\\HR_VSS02:POS_LEFT", "\\HR_VSS03:RIGHT", "\\HR_VSS04:TOP", "\\HR_VSS05:HOR", "\\HR_VSS05:VER"}; */
char channelNodeName2[6][22] = {":POS_BOTTOM", ":POS_LEFT", ":RIGHT", ":TOP", ":HOR", ":VER"};
char channelWaveName[5][17] = {"\\HR_WAVE:FOO", "\\HR_WAVE:X_CAL0", "\\HR_WAVE:X_CAL1", "\\HR_WAVE:X_CAL2", "\\HR_WAVE:X_CAL3"};
char proterty[91][40] = {ATSIF_PROP_TYPE,ATSIF_PROP_ACTIVE,ATSIF_PROP_VERSION,ATSIF_PROP_TIME,ATSIF_PROP_FORMATTED_TIME,ATSIF_PROP_FILENAME,
ATSIF_PROP_TEMPERATURE,ATSIF_PROP_UNSTABILIZEDTEMPERATURE,ATSIF_PROP_HEAD,ATSIF_PROP_HEADMODEL,ATSIF_PROP_STORETYPE,
ATSIF_PROP_DATATYPE,ATSIF_PROP_SIDISPLACEMENT,ATSIF_PROP_SINUMBERSUBFRAMES,ATSIF_PROP_PIXELREADOUTTIME,ATSIF_PROP_TRACKHEIGHT,
ATSIF_PROP_READPATTERN,ATSIF_PROP_READPATTERN_FULLNAME,ATSIF_PROP_SHUTTERDELAY,ATSIF_PROP_CENTREROW,ATSIF_PROP_ROWOFFSET,
ATSIF_PROP_OPERATION,ATSIF_PROP_MODE,ATSIF_PROP_MODE_FULLNAME,ATSIF_PROP_TRIGGERSOURCE,ATSIF_PROP_TRIGGERSOURCE_FULLNAME,
ATSIF_PROP_TRIGGERLEVEL,ATSIF_PROP_EXPOSURETIME,ATSIF_PROP_DELAY,ATSIF_PROP_INTEGRATIONCYCLETIME,ATSIF_PROP_NUMBERINTEGRATIONS,
ATSIF_PROP_KINETICCYCLETIME,ATSIF_PROP_FLIPX,ATSIF_PROP_FLIPY,ATSIF_PROP_CLOCK,ATSIF_PROP_ACLOCK,ATSIF_PROP_IOC,
ATSIF_PROP_FREQUENCY,ATSIF_PROP_NUMBERPULSES,ATSIF_PROP_FRAMETRANSFERACQMODE,ATSIF_PROP_BASELINECLAMP,ATSIF_PROP_PRESCAN,
ATSIF_PROP_EMREALGAIN,ATSIF_PROP_BASELINEOFFSET,ATSIF_PROP_SWVERSION,ATSIF_PROP_SWVERSIONEX,ATSIF_PROP_MCP,ATSIF_PROP_GAIN,
ATSIF_PROP_VERTICALCLOCKAMP,ATSIF_PROP_VERTICALSHIFTSPEED,ATSIF_PROP_OUTPUTAMPLIFIER,ATSIF_PROP_PREAMPLIFIERGAIN,ATSIF_PROP_SERIAL,
ATSIF_PROP_DETECTORFORMATX,ATSIF_PROP_DETECTORFORMATZ,ATSIF_PROP_NUMBERIMAGES,ATSIF_PROP_NUMBERSUBIMAGES,ATSIF_PROP_SUBIMAGE_HBIN,
ATSIF_PROP_SUBIMAGE_VBIN,ATSIF_PROP_SUBIMAGE_LEFT,ATSIF_PROP_SUBIMAGE_RIGHT,ATSIF_PROP_SUBIMAGE_TOP,ATSIF_PROP_SUBIMAGE_BOTTOM,
ATSIF_PROP_BASELINE,ATSIF_PROP_CCD_LEFT,ATSIF_PROP_CCD_RIGHT,ATSIF_PROP_CCD_TOP,ATSIF_PROP_CCD_BOTTOM,ATSIF_PROP_SENSITIVITY,ATSIF_PROP_DETECTIONWAVELENGTH,
ATSIF_PROP_COUNTCONVERTMODE,ATSIF_PROP_ISCOUNTCONVERT,ATSIF_PROP_X_AXIS_TYPE,ATSIF_PROP_X_AXIS_UNIT,ATSIF_PROP_Y_AXIS_TYPE,ATSIF_PROP_Y_AXIS_UNIT,
ATSIF_PROP_Z_AXIS_TYPE,ATSIF_PROP_Z_AXIS_UNIT,ATSIF_PROP_USERTEXT,ATSIF_PROP_ISPHOTONCOUNTINGENABLED,ATSIF_PROP_NUMBERTHRESHOLDS,ATSIF_PROP_THRESHOLD1,
ATSIF_PROP_THRESHOLD2,ATSIF_PROP_THRESHOLD3,ATSIF_PROP_THRESHOLD4,ATSIF_PROP_AVERAGINGFILTERMODE,ATSIF_PROP_AVERAGINGFACTOR,ATSIF_PROP_FRAMECOUNT,
ATSIF_PROP_NOISEFILTER,ATSIF_PROP_THRESHOLD,ATSIF_PROP_TIME_STAMP};
AT_U32 atu32_ret, atu32_noFrames, atu32_frameSize, atu32_noSubImages;
float * imageBuffer;
double calibValue;
char * sz_propertyValue;
AT_U32 wave_mode;
AT_U32 atu32_left,atu32_bottom,atu32_right,atu32_top,atu32_hBin, atu32_vBin;
double at32_delay, at32_exposeTime, at32_integCycleTime, at32_kineticCycleTime;
double x_cal0;
double x_cal1;
double x_cal2;
double x_cal3;
int imageNumber; /* imageNumber, image count at time 30720/5120 = 6 int => unsigned long int */
int number_subImage; /* image count, subImage = 5 */
int right_pixel; /* pixel number : 1024 int => unsigned long int */
int image_length; /* 5120 = SubImage * 1024(frameSize) */
int image_top; /* Image Top ex.255 */
int image_bottom; /*Image Bottom ex.1 */
AT_U32 total_length;
double cycle_time; /* same value delay & kinetic_cycle & integration cycle time */
double exposure_time; /* exposure time ex.0.01700 */
double integration_cycle; /* Integration cycle time ex.0.21200 */
char buf[50];
char bufNode[50];
int null = 0;
int status;
int socket;
int j, i, m, iwave;
unsigned long int k, jum;
int dtype_float = DTYPE_FLOAT;
int dtype_double = DTYPE_DOUBLE;
int timeCount=0; /* imageNumber * right_pixel */
int dtypeLong = DTYPE_LONG;
double *timeArray; /* number_image = imageNumber = Image.no_images */
float *dataArray;
double *data;
float *image_buffer1;
int tstat;
int dataDesc;
int timeDesc;
int wavedataDesc;
int idesc = descr(&dtypeLong, &tstat, &null);
int kkk;
char tmpBuf[50];
/* This function is used to select if the entire SIF file should be read or just the header section */
atu32_ret = ATSIF_SetFileAccessMode(ATSIF_ReadAll);
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not set File access Mode. Error: %u\n", atu32_ret);
}
else {
/* This function is used to open a SIF file where the file name and path are contained
in the character array _sz_filename */
atu32_ret = ATSIF_ReadFromFile(fileNamePath);
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not open File : %s.\nError: %u\n", fileNamePath,atu32_ret);
}
else {
atu32_ret = ATSIF_GetNumberFrames(ATSIF_Signal, &atu32_noFrames);
/* This function is used to retrieve the number of frames in the SIF file. */
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not Get Number Frames. Error: %u\n", atu32_ret);
} else {
printf("Image contains %u frames.\n", atu32_noFrames);
atu32_ret = ATSIF_GetFrameSize(ATSIF_Signal, &atu32_frameSize);
/* This function is used to retrieve the number of pixels in each frame in the SIF file. */
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not Get Frame Size. Error: %u\n", atu32_ret);
} else {
printf("Each frame contains %u pixels.\n", atu32_frameSize);
atu32_ret = ATSIF_GetNumberSubImages(ATSIF_Signal, &atu32_noSubImages);
/* This function is used to retrieve the number of sub-images in each frame in the SIF file. */
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not Get Number Sub Images. Error: %u\n", atu32_ret);
} else {
printf("Each frame contains %u sub images.\n", atu32_noSubImages);
}
}
/* ################# Property Type and Value ##################*/
sz_propertyValue = (char*)malloc((sizeof(char)*MAX_PATH));
/* We need This function. Because conversion and used the value from string to int value */
atu32_ret = ATSIF_GetPropertyValue(ATSIF_Signal, ATSIF_PROP_EXPOSURETIME, sz_propertyValue, MAX_PATH);
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not get Property Value.\n");
}
else {
at32_exposeTime = atof(sz_propertyValue);
printf("Property Value ExposureTime : %s, %0.4f\n", sz_propertyValue,at32_exposeTime);
}
atu32_ret = ATSIF_GetPropertyValue(ATSIF_Signal, ATSIF_PROP_DELAY, sz_propertyValue, MAX_PATH);
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not get Property Value.\n");
}
else {
at32_delay = atof(sz_propertyValue);
printf("Property Value Delay : %s\n", sz_propertyValue);
}
atu32_ret = ATSIF_GetPropertyValue(ATSIF_Signal, ATSIF_PROP_INTEGRATIONCYCLETIME, sz_propertyValue, MAX_PATH);
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not get Property Value.\n");
}
else {
at32_integCycleTime = atof(sz_propertyValue);
printf("Property Value IntegrationCycleTime : %s, %0.4f\n", sz_propertyValue,at32_integCycleTime );
}
atu32_ret = ATSIF_GetPropertyValue(ATSIF_Signal, ATSIF_PROP_KINETICCYCLETIME, sz_propertyValue, MAX_PATH);
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not get Property Value.\n");
}
else {
at32_kineticCycleTime = atof(sz_propertyValue);
printf("Property Value KineticCycleTime : %s, %0.4f \n", sz_propertyValue, at32_kineticCycleTime);
}
atu32_ret = ATSIF_GetPropertyValue(ATSIF_Signal, ATSIF_PROP_MODE, sz_propertyValue, MAX_PATH);
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not get Property Value.\n");
}
else {
printf("Property Value Mode : %s\n", sz_propertyValue);
}
for(i = 0; i < 91; i++) {
atu32_ret = ATSIF_GetPropertyValue(ATSIF_Signal, proterty[i], sz_propertyValue, MAX_PATH);
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not get Property Value : %s.\n",proterty[i]);
}
else {
printf("Property Value %s : %s\n",proterty[i], sz_propertyValue);
}
}
for(i = 0; i < (int)atu32_noFrames; i++) {
sprintf(tmpBuf,"TimeStamp %d",i);
atu32_ret = ATSIF_GetPropertyValue(ATSIF_Signal, tmpBuf, sz_propertyValue, MAX_PATH);
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not get Property Value : %s.\n",tmpBuf);
}
else {
printf("Property Value %s : %s\n",tmpBuf, sz_propertyValue);
}
}
/* ################# Property Type and Value ##################*/
for(i = 0; i < (int)atu32_noSubImages; i++) {
printf("SubImage %u Properties:\n", (i + 1));
/* This function is used to retrieve the information about each sub-image in the SIF file. */
atu32_ret = ATSIF_GetSubImageInfo(ATSIF_Signal,
i,
&atu32_left,&atu32_bottom,
&atu32_right,&atu32_top,
&atu32_hBin,&atu32_vBin);
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not Get Sub Image Info. Error: %u\n", atu32_ret);
}
else {
printf("\tleft\t: %u\tbottom\t: %u\n", atu32_left, atu32_bottom);
printf("\tright\t: %u\ttop\t: %u\n", atu32_right, atu32_top);
printf("\thBin\t: %u\tvBin\t: %u\n", atu32_hBin, atu32_vBin);
}
}
/* I just check the data of frame in SIF file */
imageBuffer = (float*)malloc((sizeof(float))* atu32_frameSize);
atu32_ret = ATSIF_GetFrame(ATSIF_Signal,0, imageBuffer, atu32_frameSize);
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not Get Frame. Error: %u\n", atu32_ret);
} else {
printf("The first 20 pixel values are : \n");
for(i = 0; i < 20; i++) {
printf("%f\n", imageBuffer[i]);
}
}
free(imageBuffer);
}
}
}
#if 1
/* Wavelength = x_cal(0) + x_cal(1).X + x_cal(2).X^2 + x_cal(3).X^3 */
/* X is the pixel number on in left hand column. double to float type changed 2009.05.15 */
/* AT_U32 atu32_ret, atu32_noFrames, atu32_frameSize, atu32_noSubImages; */
/* AT_U32 atu32_left,atu32_bottom,atu32_right,atu32_top,atu32_hBin, atu32_vBin; */
/* AT_32 at32_delay, at32_exposeTime, at32_integCycleTime, at32_kineticCycleTime; */
/* time count number = total_length/image_length */
cycle_time = at32_kineticCycleTime; /* same value delay & kinetic_cycle & integration cycle time */
exposure_time = at32_exposeTime; /* exposure time ex.0.01700 */
integration_cycle = at32_integCycleTime; /* Integration cycle time ex.0.21200 */
/* This function is used to retrieve the information about each sub-image in the SIF file. */
atu32_ret = ATSIF_GetSubImageInfo(ATSIF_Signal, 0, &atu32_left,&atu32_bottom, &atu32_right,&atu32_top, &atu32_hBin,&atu32_vBin);
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not Get Sub Image Info. Error: %u\n", atu32_ret);
}
else {
printf("\tleft\t: %u\tbottom\t: %u\n", atu32_left, atu32_bottom);
printf("\tright\t: %u\ttop\t: %u\n", atu32_right, atu32_top);
printf("\thBin\t: %u\tvBin\t: %u\n", atu32_hBin, atu32_vBin);
}
image_top = atu32_top; /* Image Top ex.255 */
image_bottom = atu32_bottom; /*Image Bottom ex.1 */
image_length = atu32_frameSize * atu32_noSubImages; /* 5120 = SubImage * 1024(frameSize) */
total_length = atu32_noFrames * atu32_frameSize * atu32_noSubImages;
image_buffer1 = (float*)malloc(4*(size_t)total_length);
atu32_ret = ATSIF_GetAllFrames(ATSIF_Signal, image_buffer1, total_length);
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not Get Frame. Error: %u\n", atu32_ret);
} else {
printf(" Success All frames data read. : \n");
}
imageNumber = atu32_noFrames; /* imageNumber, image count at time 30720/5120 = 6 int => unsigned long int */
number_subImage = atu32_noSubImages; /* image count, subImage = 5 */
right_pixel = atu32_frameSize; /* pixel number : 1024 int => unsigned long int */
printf("total_length : %d \n",total_length);
printf("image_length : %d \n",image_length);
printf("number_image : %d \n",imageNumber);
printf("imageNumber : %d \n",imageNumber);
printf("number_subImage : %d \n",number_subImage);
printf("left_pixel : %d \n",atu32_left);
printf("right_pixel : %d \n",right_pixel);
printf("cycle_time : %f \n",cycle_time);
printf("exposure_time : %f \n",exposure_time);
printf("integration_cycle : %f \n",integration_cycle);
printf("SIF MDSplus Data Put Part. \n");
timeCount = imageNumber*right_pixel;
timeArray = (double*)malloc((sizeof(double))*imageNumber*right_pixel); /* 1D plot Data Time - OK */
dataArray = (float*)malloc(4*(size_t)timeCount); /* 1D plot Data */
data = (double *)malloc(sizeof(double)*right_pixel); /* waveLength and Pixel Number */
iwave = 0;
for(i=(int)atu32_frameSize; i >= 1; i--) {
atu32_ret = ATSIF_GetPixelCalibration(ATSIF_Signal, ATSIF_CalibX, i, &calibValue);
/* This function is used to retrieve wavelength in each frame in the SIF file. */
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not Get ATSIF_GetPixelCalibration. Error: %u\n", atu32_ret);
}else {
/* printf("Wavelength pixel:%d , X: %f ", i, calibValue); */
data[iwave] = calibValue;
iwave ++;
}
}
printf("\n");
x_cal0 = data[0];
x_cal1 = data[1] - data[0];
x_cal2 = -0.000000;
x_cal3 = -0.000000;
if(x_cal1 < 0) {
iwave = 0;
for(i=1; i <= (int)atu32_frameSize; i++) {
atu32_ret = ATSIF_GetPixelCalibration(ATSIF_Signal, ATSIF_CalibX, i, &calibValue);
/* This function is used to retrieve wavelength in each frame in the SIF file. */
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not Get ATSIF_GetPixelCalibration. Error: %u\n", atu32_ret);
}else {
/* printf("Wavelength pixel:%d , X: %f ", i, calibValue); */
data[iwave] = calibValue;
iwave ++;
}
}
x_cal0 = data[0];
x_cal1 = data[1] - data[0];
x_cal2 = -0.000000;
x_cal3 = -0.000000;
}
printf("x_cal0 : %f \n",x_cal0);
printf("x_cal1 : %f \n",x_cal1);
printf("x_cal2 : %f \n",x_cal2);
printf("x_cal3 : %f \n",x_cal3);
if(storeMode == 1){
socket = MdsConnect(MDSIP);
if ( socket == -1)
{
fprintf(stderr,"Error connecting to Server.\n");
free(data);
free(dataArray);
free(timeArray);
return 0;
}
fprintf(stdout, "\nMdsplus >>> MdsConnect(\"%s\")...OK\n",MDSIP);
}
#if 0
if(storeMode == 0){
sprintf(buf, "TCL(\"set tree %s/shot=-1\")",TREE);
status = MdsValue(buf, &idesc, &null, &len);
if (!status_ok(status))
{
fprintf(stderr,"Error TCL Command Set tree -1 \n");
return -1;
}
sprintf(buf, "TCL(\"create pulse %ld\")",shotNumber);
status = MdsValue(buf, &idesc, &null, &len);
if (!status_ok(status))
{
fprintf(stderr,"Error TCL Command Create Pulse \n");
return -1;
}
sprintf(buf, "TCL(\"close\")",TREE);
status = MdsValue(buf, &idesc, &null, &len);
}
#endif
status = MdsOpen(TREE, &shotNumber);
if (!status_ok(status))
{
fprintf(stderr,"Error openning tree for shot %l. \n",shotNumber);
return -1;
}
fprintf(stdout, "\nMdsplus >>> MdsOpen(\"%s %d\")...OK\n",TREE,shotNumber);
/* MdsPut Visible Spectrometer Data Infomation */
sprintf(buf, "FS_FLOAT(%d)", imageNumber);
status = MdsPut(channelNodeName1[0], buf,&null);
sprintf(buf, "FS_FLOAT(%d)", right_pixel);
status = MdsPut(channelNodeName1[1], buf,&null);
sprintf(buf, "FS_FLOAT(%d)", image_top);
status = MdsPut(channelNodeName1[2], buf,&null);
sprintf(buf, "FS_FLOAT(%d)", number_subImage);
status = MdsPut(channelNodeName1[3], buf,&null);
sprintf(buf, "FS_FLOAT(%f)", exposure_time);
status = MdsPut(channelNodeName1[4], buf,&null);
sprintf(buf, "FS_FLOAT(%f)", integration_cycle);
status = MdsPut(channelNodeName1[5], buf,&null);
sprintf(buf, "FS_FLOAT(%f)", cycle_time);
status = MdsPut(channelNodeName1[6], buf,&null);
wavedataDesc = descr(&dtype_double, data, &right_pixel, &null);
atu32_ret = ATSIF_GetPropertyValue(ATSIF_Signal, ATSIF_PROP_X_AXIS_TYPE, sz_propertyValue, MAX_PATH);
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not get Property Value : %s.\n",ATSIF_PROP_X_AXIS_TYPE);
} else {
if(strcmp(sz_propertyValue,"Wavelength") != 0){
wave_mode = 0; /*pixel mode */
} else {
wave_mode = 1; /*wavelength mode */
}
printf("Property Value %s : %s\n",ATSIF_PROP_X_AXIS_TYPE, sz_propertyValue);
}
if (wave_mode==1)
{
fprintf(stdout, "DAQ WAVEMode OK : %s \n",MODEWAVE);
#if 0
fprintf(stdout,"1 Wave Data Cal 0 : %f \n", data[0]);
fprintf(stdout,"1 Wave Data Cal 1 : %f \n", data[1]);
for (i=0; i<right_pixel; i++)
{
data[i] = x_cal0 + x_cal1*(i) + x_cal2*((i)^2) + x_cal3*((i)^3);
}
fprintf(stdout,"2 Wave Data Cal 0 : %f \n", data[0]);
fprintf(stdout,"2 Wave Data Cal 1 : %f \n", data[1]);
#endif
status = MdsPut(channelWaveName[0], "BUILD_SIGNAL($1,,)", &wavedataDesc, &null);
sprintf(buf, "FS_FLOAT(%f)", x_cal0);
status = MdsPut(channelWaveName[1], buf,&null);
sprintf(buf, "FS_FLOAT(%f)", x_cal1);
status = MdsPut(channelWaveName[2], buf,&null);
sprintf(buf, "FS_FLOAT(%f)", x_cal2);
status = MdsPut(channelWaveName[3], buf,&null);
sprintf(buf, "FS_FLOAT(%f)", x_cal3);
status = MdsPut(channelWaveName[4], buf,&null);
}
else if (wave_mode!=1)
{
fprintf(stdout, "DAQ PIXMode OK : %s \n",MODEPIX);
for (i=0; i<right_pixel; i++)
{
data[i] = i+1;
}
}
else
{
fprintf(stdout, "DAQ PIXMode No Good : %s \n",MODEPIX);
for (i=0; i<right_pixel; i++)
{
data[i] = i+1;
}
fprintf(stderr, "Error mode is not waveLength and Pixel number \n");
}
fprintf(stdout, "CA_DAQ Start Time : %f & BLIP Time : %f \n",daqTime, blipTime);
timeDesc = descr(&dtype_double, timeArray, &imageNumber, &null);
for (i=0; i<imageNumber; i++)
{
timeArray[i] = (double)(daqTime-blipTime)+ i*cycle_time;
printf("timeArray : %f, frame nu : %d ",timeArray[i],i);
}
printf("\n");
/* 1D Data Time Plot image length = imagenumber * right_pixel OK - 2008.02.11
timeDesc = descr(&dtype_float, timeArray, &timeCount, &null);
fprintf(out,"int: timeCount %d\n",timeCount);
for (i=0; i<timeCount; i++)
{
timeArray[i] = (float)i/10;
fprintf(out,"int2: timeArray2 : %f, i2 : %d\n",timeArray[i],i);
}
dataDesc = descr(&dtype_float, dataArray, &timeCount, &null); //descr(TYPE, Data1D, 1Dcount,0)
*/
kkk=1;
dataDesc = descr(&dtype_float, dataArray, &right_pixel, &imageNumber, &kkk, &null); /* descr(TYPE, Data3D, nx,ny,nz,0) ?? */
if (wave_mode==1)
{
fprintf(stdout, "DAQ WAVEMode OK : %s \n",MODEWAVE);
for (m=0,k=0; m<number_subImage; m++)
{
for (jum=0,j=0; j<imageNumber; j++)
{
for (i=0; i<right_pixel; i++, k++)
{
dataArray[j*right_pixel + i] = image_buffer1[k];
/* printf("dataArray Put : %f, image_buffer Data : %f, k = %d\n", dataArray[j*right_pixel + i],image_buffer1[k],k); */
}
if(imageNumber>1)
{
jum=(j+1)*image_length + (m*right_pixel); /* right_pixel = 1024, image_length = 5120 */
k= jum;
}
}
k = (m+1)*right_pixel ;
status = MdsPut(channelTagName[m], "BUILD_SIGNAL($1,,$2)", &dataDesc, &timeDesc, &null); /* Data Put 2D and Time TEST ING */
/* MDSput The VSS Data Infomations about Image */
atu32_ret = ATSIF_GetSubImageInfo(ATSIF_Signal,
m,
&atu32_left,&atu32_bottom,
&atu32_right,&atu32_top,
&atu32_hBin,&atu32_vBin);
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not Get Sub Image Info. Error: %u\n", atu32_ret);
}
sprintf(buf, "FS_FLOAT(%d)", atu32_bottom);
sprintf(bufNode, "\\HR_VSS0%d%s",m+1,channelNodeName2[0]);
status = MdsPut(bufNode, buf,&null);
sprintf(buf, "FS_FLOAT(%d)", atu32_left);
sprintf(bufNode, "\\HR_VSS0%d%s",m+1,channelNodeName2[1]);
status = MdsPut(bufNode, buf,&null);
sprintf(buf, "FS_FLOAT(%d)", atu32_right);
sprintf(bufNode, "\\HR_VSS0%d%s",m+1,channelNodeName2[2]);
status = MdsPut(bufNode, buf,&null);
sprintf(buf, "FS_FLOAT(%d)", atu32_top);
sprintf(bufNode, "\\HR_VSS0%d%s",m+1,channelNodeName2[3]);
status = MdsPut(bufNode, buf,&null);
sprintf(buf, "FS_FLOAT(%d)", atu32_hBin);
sprintf(bufNode, "\\HR_VSS0%d%s",m+1,channelNodeName2[4]);
status = MdsPut(bufNode, buf,&null);
sprintf(buf, "FS_FLOAT(%d)", atu32_vBin);
sprintf(bufNode, "\\HR_VSS0%d%s",m+1,channelNodeName2[5]);
status = MdsPut(bufNode, buf,&null);
printf("check Node name for each channel VSS Name: %s\n",bufNode);
}
}
else if (wave_mode!=1)
{
fprintf(stdout, "DAQ PIXMode OK : %s \n",MODEPIX);
for (m=0,k=0; m<number_subImage; m++)
{
for (jum=0,j=0; j<imageNumber; j++)
{
for (i=0; i<right_pixel; i++, k++)
{
/* if(k==0){
fprintf(stdout, "MDSplus Put the Data image_lastBuffData : %f Num k : %d \n", image_buffer1[k], k);
}
*/
dataArray[j*right_pixel + i] = image_buffer1[k];
/* printf("dataArray Put : %f, image_buffer Data : %f, k = %d\n", dataArray[j*right_pixel + i],image_buffer1[k],k); */
/* if(k==1023 || k==1024){
fprintf(stdout, "MDSplus Put the Data image_lastBuffData 1023 or 1024 : %f Num k : %d \n", image_buffer1[k], k);
}
*/
}
if(imageNumber>1)
{
/* jum=jum + (m*right_pixel) + (image_length); //right_pixel = 1024, image_length = 5120 */
jum=(j+1)*image_length + (m*right_pixel); /* right_pixel = 1024, image_length = 5120 */
k= jum;
}
}
k = (m+1)*right_pixel ;
/*
// status = MdsPut(channelTagName[m], "BUILD_SIGNAL($1,,)", &dataDesc, &null); //Data Put OK 2008.01.28
// status = MdsPut(channelTagName[m], "BUILD_SIGNAL($1,,$2)", &dataDesc, &timeDesc, &null); //Data Put 1D and Time OK 2008.02.11
// status = MdsPut(channelTagName[m], "BUILD_SIGNAL($1,,$2)", &dataDesc, &timeDesc, &null); //Data Put 2D and Time TEST ING
*/
status = MdsPut(channelTagName[m], "BUILD_SIGNAL($1,,BUILD_WITH_UNITS($2,'S')", &dataDesc, &timeDesc, &null);
/* MDSput The VSS Data Infomations about Image */
atu32_ret = ATSIF_GetSubImageInfo(ATSIF_Signal,
m,
&atu32_left,&atu32_bottom,
&atu32_right,&atu32_top,
&atu32_hBin,&atu32_vBin);
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not Get Sub Image Info. Error: %u\n", atu32_ret);
}
sprintf(buf, "FS_FLOAT(%d)", atu32_bottom);
sprintf(bufNode, "\\HR_VSS0%d%s",m+1,channelNodeName2[0]);
status = MdsPut(bufNode, buf,&null);
sprintf(buf, "FS_FLOAT(%d)", atu32_left);
sprintf(bufNode, "\\HR_VSS0%d%s",m+1,channelNodeName2[1]);
status = MdsPut(bufNode, buf,&null);
sprintf(buf, "FS_FLOAT(%d)", atu32_right);
sprintf(bufNode, "\\HR_VSS0%d%s",m+1,channelNodeName2[2]);
status = MdsPut(bufNode, buf,&null);
sprintf(buf, "FS_FLOAT(%d)", atu32_top);
sprintf(bufNode, "\\HR_VSS0%d%s",m+1,channelNodeName2[3]);
status = MdsPut(bufNode, buf,&null);
sprintf(buf, "FS_FLOAT(%d)", atu32_hBin);
sprintf(bufNode, "\\HR_VSS0%d%s",m+1,channelNodeName2[4]);
status = MdsPut(bufNode, buf,&null);
sprintf(buf, "FS_FLOAT(%d)", atu32_vBin);
sprintf(bufNode, "\\HR_VSS0%d%s",m+1,channelNodeName2[5]);
status = MdsPut(bufNode, buf,&null);
printf("check Node name for each channel VSS Name: %s\n",bufNode);
/*
// status = MdsPut(channelTagName[m], "BUILD_SIGNAL(BUILD_WITH_UNITS($1,'Counts'),BUILD_WITH_UNITS($2,'wave'),BUILD_DIM(BUILD_WINDOW(0,2,10),BUILD_SLOPE($3)))", &dataDesc, &data, &timeDesc, &null);
// status = MdsPut(channelTagName[m], "BUILD_SIGNAL(BUILD_WITH_UNITS($1,'Counts'),BUILD_WITH_UNITS($2,'wave'),BUILD_WITH_UNITS($3,'sec'))", &dataDesc, &data, &timeDesc, &null);
// status = MdsPut(channelTagName[m], "BUILD_SIGNAL(BUILD_WITH_UNITS($1,'Counts'),BUILD_WITH_UNITS($2,'wave'),BUILD_WITH_UNITS($3,'sec'))", &dataDesc, &data2, &timeDesc, &null);
// status = MdsPut(channelTagName[m], "BUILD_SIGNAL(BUILD_WITH_UNITS($1,'Counts'),BUILD_WITH_UNITS($2,'wave'),BUILD_WITH_UNITS($3,'sec'))", &dataDesc, 1024, &timeDesc, &null);
// status = MdsPut(channelTagName[m], "BUILD_SIGNAL(BUILD_WITH_UNITS($1,'Counts'),BUILD_WITH_UNITS($2,'Sec'),BUILD_WITH_UNITS($3,'wave'))", &dataDesc, &timeDesc, &data2, &null);
// status = MdsPut(channelTagName[m], "BUILD_SIGNAL(BUILD_WITH_UNITS($1,'Counts'),BUILD_WITH_UNITS($2,'Sec'),BUILD_WITH_UNITS($3,'wave'))", &dataDesc, &timeDesc, &right_pixel, &null);
// status = MdsPut(channelTagName[m], "BUILD_SIGNAL(BUILD_WITH_UNITS($1,'Counts'),BUILD_WITH_UNITS($2,'Sec'),BUILD_WITH_UNITS($3,'wave'))", &dataDesc, &timeDesc, &right_pixel, &null);
// Image Test Put mdsput,"\TOP.X_IMAGE:FAST_CCD1","BUILD_SIGNAL($,*,MAKE_DIM(MAKE_WINDOW(0,425,0.0),MAKE_SLOPE(0.001)),MAKE_DIM(MAKE_WINDOW(0,401,0.0),MAKE_SLOPE(0.001)),MAKE_DIM(MAKE_WINDOW(0,$,0.0),MAKE_SLOPE(MAKE_WITH_UNITS($,'s'))))",jpgimage,filecount_1,dtime
// fprintf(stdout, "MDSplus Put the Data channelTagName2 : %s int Size : %f \n", channelTagName[m],data2);
// free(dataArray);
*/
}
}
else
{
fprintf(stderr, "Error mode is not waveLength and Pixel number \n");
}
free(data);
free(timeArray);
free(dataArray);
/* MDSplus Close Tree */
status = MdsClose(TREE, &shotNumber);
if (!status_ok(status))
{
fprintf(stderr,"Error closing tree for shot %l. \n",shotNumber);
return -1;
}
fprintf(stdout, "\nMdsplus >>> MdsClose(\"%d\")...OK\n",shotNumber);
atu32_ret = ATSIF_CloseFile();
if(atu32_ret != ATSIF_SUCCESS) {
printf("Could not get Property Value.\n");
}
#endif
return 1;
}
