static int DoCamMulti(char *routine, char *name, int a, int f, int count, void *data, int mem, short *iosb)
{
int serverid = RemoteServerId();
int status = 0;
int writeData;
if (serverid)
{
struct descrip data_d = {8,1,{0,0,0,0,0,0,0},0};
struct descrip ans_d = {0,0,{0,0,0,0,0,0,0},0};
char cmd[512];
writeData = (!(f &0x08)) && (f > 8);
sprintf(cmd,"CamMulti('%s','%s',%d,%d,%d,%s,%d,_iosb)",routine,name,a,f,count,writeData ? "_data=$" : "_data",mem);
if (writeData)
{
data_d.dtype = mem < 24 ? DTYPE_SHORT : DTYPE_LONG;
data_d.dims[0] = count;
data_d.ptr = data;
status = MdsValue(serverid,cmd,&data_d,&ans_d,0);
}
else
{
status = MdsValue(serverid,cmd,&ans_d,0);
}
if (status & 1 && ans_d.dtype == DTYPE_LONG && ans_d.ptr)
{
memcpy(&status,ans_d.ptr,4);
free(ans_d.ptr);
ans_d.ptr = 0;
if (data && f < 8)
getdata(serverid,data);
getiosb(serverid,iosb);
}
}
return status;
}
int main( int argc, void **argv)
{
int status;
struct descrip ans;
float val = 9876;
struct descrip vald = {DTYPE_FLOAT,0};
long sock = ConnectToMds((argc > 1) ? argv[1] : "lost.pfc.mit.edu:9000");
if (sock != -1)
{
printf("status from MdsOpen = %d\n",MdsOpen(sock,"main",-1));
ans.ptr = 0;
if (MdsValue(sock,"f_float(member)",&ans,0) & 1)
{
printf("%g\n",*(float *)ans.ptr);
val = *(float *)ans.ptr;
val = val + (float)1.;
}
else
printf("%s\n",ans.ptr);
if (ans.ptr)
{
free(ans.ptr);
ans.ptr = 0;
}
vald.ptr = (void *)&val;
status = MdsPut(sock,"member","$",&vald,0);
if (!(status & 1)) printf("Error during put %d\n",status);
if (MdsValue(sock,"42.0",&ans,0) & 1)
printf("%g\n",*(float *)ans.ptr);
else
printf("%s\n",ans.ptr);
if (ans.ptr) free(ans.ptr);
}
return 1;
}
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;
}
STATIC_ROUTINE int MdsLoginVMS(SOCKET sock, char *username, char *password)
{
struct descrip loginget_arg;
struct descrip loginpwd_arg;
struct descrip pwd_arg;
struct descrip ansarg;
int status;
ansarg.ptr = 0;
status = MdsValue(sock, LOGINREQUEST, MakeDescrip(&loginget_arg,DTYPE_CSTRING,0,0,LOGINGETP1), &ansarg, NULL);
if (status & 1)
{
if (ansarg.ptr && ansarg.dtype == DTYPE_CHAR)
{
char alg = *(char *)ansarg.ptr;
free(ansarg.ptr);
status = MdsValue(sock, LOGINREQUEST, MakeDescrip(&loginget_arg,DTYPE_CSTRING,0,0,LOGINGETP2), &ansarg, NULL);
if (status & 1)
{
if (ansarg.ptr && ansarg.dtype == DTYPE_SHORT)
{
struct dscr { short length; char dtype; char class; void *pointer;};
int hash[2];
int i;
struct dscr hashd = {8, 9, 1, 0};
struct dscr pwdd = {0,14,1,0};
struct dscr userd = {0,14,1,0};
short salt = *(short *)ansarg.ptr;
free(ansarg.ptr);
ansarg.ptr = 0;
hashd.pointer = hash;
pwdd.length = strlen(password);
pwdd.pointer = strcpy((char *)malloc(pwdd.length+1),password);
userd.length = strlen(username);
userd.pointer = strcpy((char *)malloc(userd.length+1),username);
for (i=0;i<pwdd.length;i++)
((char *)pwdd.pointer)[i] = __toupper(((char *)pwdd.pointer)[i]);
for (i=0;i<userd.length;i++)
((char *)userd.pointer)[i] = __toupper(((char *)userd.pointer)[i]);
status = Lgihpwd(&hashd,&pwdd,alg,salt,&userd);
if (status & 1)
{
int two = 2;
status = MdsValue(sock, LOGINREQUEST, MakeDescrip(&loginpwd_arg,DTYPE_CSTRING,0,0,LOGINPWD),
MakeDescrip(&pwd_arg,DTYPE_LONG,1,&two,hash), &ansarg, NULL);
}
else
printf("Login error: Error checking password\n");
}
else
{
if (ansarg.ptr)
{
free(ansarg.ptr);
ansarg.ptr = 0;
}
printf("Login error: Unrecognized response from server\n");
status = 0;
}
}
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 doMdsValue(int argc, char **argv)
{
if (getenv("MDS_SOCK")){
int sock = atoi(getenv("MDS_SOCK"));
fprintf(stderr, "MDS_SOCK is %s\n", getenv("MDS_SOCK"));
switch(argc){
case 2: {
char* expr = argv[1];
struct descrip ans = {DTYPE_FLOAT, 0};
int status;
fprintf(stderr, "MdsValue sock:%d expr %s\n",
sock, expr);
status = MdsValue(sock, expr, &ans, 0);
if (!(status & 1)){
fprintf(stderr,
"Error during value %d\n",status);
}else{
fprintf(stderr, "Success %g\n",
*(float*)ans.ptr);
printf("%g\n", *(float*)ans.ptr);
}
return 0;
} default:
fprintf(stderr, "ERROR: args field value\n");
}
}else{
fprintf(stderr, "MDS_SOCK environment var is missing\n");
}
return -1;
}
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 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);
}
int MdsClose(SOCKET sock)
{
struct descrip ansarg;
STATIC_CONSTANT char *expression = "TreeClose()";
int status = MdsValue(sock, expression, &ansarg, NULL);
if ((status & 1) && (ansarg.dtype == DTYPE_LONG)) status = *(int *)ansarg.ptr;
if (ansarg.ptr) free(ansarg.ptr);
return status;
}
static void checkConnection(struct IpDescriptor *ipDescr)
{
struct descrip resDsc;
int status = 0;
if(ipDescr->id != -1)
status = MdsValue(ipDescr->id, "1", &resDsc, NULL);
if(!(status & 1))
ipDescr->id = ConnectToMds(ipDescr->ipAddr);
}
static void getdata(int serverid, void *data)
{
int status;
struct descrip ans_d = {0,0,{0,0,0,0,0,0,0},0};
status = MdsValue(serverid,"_data",&ans_d,0);
if (status & 1 && (ans_d.dtype == DTYPE_USHORT || ans_d.dtype == DTYPE_LONG) && ans_d.ptr)
memcpy(data,ans_d.ptr,((ans_d.dtype == DTYPE_USHORT) ? 2 : 4) * ans_d.dims[0]);
if (ans_d.ptr)
free(ans_d.ptr);
}
int MdsSetDefault(SOCKET sock, char *node)
{
struct descrip nodearg;
struct descrip ansarg;
STATIC_CONSTANT char *expression = "TreeSetDefault($)";
int status = MdsValue(sock, expression, MakeDescrip(&nodearg,DTYPE_CSTRING,0,0,node), &ansarg, NULL);
if ((status & 1) && (ansarg.dtype == DTYPE_LONG)) status = *(int *)ansarg.ptr;
if (ansarg.ptr) free(ansarg.ptr);
return status;
}
int MdsOpen(SOCKET sock, char *tree, int shot)
{
struct descrip treearg;
struct descrip shotarg;
struct descrip ansarg;
STATIC_CONSTANT char *expression = "TreeOpen($,$)";
int status = MdsValue(sock, expression, MakeDescrip((struct descrip *)&treearg,DTYPE_CSTRING,0,0,tree),
MakeDescrip((struct descrip *)&shotarg,DTYPE_LONG,0,0,&shot),
(struct descrip *)&ansarg, (struct descrip *)NULL);
if ((status & 1) && (ansarg.dtype == DTYPE_LONG)) status = *(int *)ansarg.ptr;
if (ansarg.ptr) free(ansarg.ptr);
return status;
}
static void getiosb(int serverid, short *iosb)
{
int status;
struct descrip ans_d = {0,0,{0,0,0,0,0,0,0},0};
status = MdsValue(serverid,"_iosb",&ans_d,0);
if (status & 1 && ans_d.dtype == DTYPE_USHORT && ans_d.ndims == 1 && ans_d.dims[0] == 4)
{
memcpy(RemCamLastIosb,ans_d.ptr,8);
if (iosb)
memcpy(iosb,ans_d.ptr,8);
}
if (ans_d.ptr)
free(ans_d.ptr);
}
static char * mdsplus_translate_status(int socket,int status)
{
char expression[100];
struct descrip ans;
int tmpstat;
sprintf(expression,"getmsg(%d)",status);
tmpstat=MdsValue(socket,expression,&ans,0);
if (tmpstat & 1 && ans.ptr)
return (char *)ans.ptr;
else
{
char *msg = malloc(100);
sprintf(msg,"Unable to determine error string, error status=%d",status);
return msg;
}
}
static void EventRegisterRemote(char *eventName)
{
int i, status;
struct descrip evNameDsc, thisIpDsc, resDsc;
struct ClientEventDescriptor *currEventDescr, *prevEventDescr;
char thisIp[512];
strncpy(thisIp, getenv("REVENT_THIS_MDSIP"), 511);
configure();
if(!thisIp || !*thisIp || numClientIpDescriptors == 0)
return; //Missing information for managing rempte events
MakeDescrip(&thisIpDsc, DTYPE_CSTRING, 0, 0, thisIp);
MakeDescrip(&evNameDsc, DTYPE_CSTRING, 0, 0, eventName);
lock();
prevEventDescr = currEventDescr = clientEventDescrHead;
while(currEventDescr && strcmp(currEventDescr->eventName, eventName))
{
prevEventDescr = currEventDescr;
currEventDescr = currEventDescr->nxt;
}
if(!currEventDescr)
{
currEventDescr = (struct ClientEventDescriptor *)malloc(sizeof(struct ClientEventDescriptor));
currEventDescr->eventName = malloc(strlen(eventName) + 1);
strcpy(currEventDescr->eventName, eventName);
currEventDescr->refCount = 1;
for(i = 0; i < numClientIpDescriptors; i++)
{
status = MdsValue(clientIpDescriptors[i].id, "RtEventsShr->EventRegisterExecute($1,$2)", &evNameDsc, &thisIpDsc, &resDsc, NULL);
if(!(status & 1))
printf("Error sending event: %s\n", MdsGetMsg(status));
}
if(prevEventDescr)
prevEventDescr->nxt = currEventDescr;
else
clientEventDescrHead = currEventDescr;
}
else //already found: someone already registered fior this event on this machine
currEventDescr->refCount++;
unlock();
}
static void eventHandler(char *evName, char *buf, int size, void *arg)
{
int i, status;
struct descrip evNameDsc, bufDsc, resDsc;
struct descrip bufSizeDsc;
struct ServerEventDescriptor *currEventDescr;
printf("event handler %s\n", evName);
lock();
currEventDescr= (struct ServerEventDescriptor *)arg;
if(currEventDescr->busy)
return; //Do not propagate events received from outside
unlock();
MakeDescrip(&evNameDsc,DTYPE_CSTRING,0,0,evName);
MakeDescrip(&bufDsc, DTYPE_UCHAR, 1, &size, buf);
MakeDescrip(&bufSizeDsc, DTYPE_LONG, 0,0, (char *)&size);
for(i = 0; i < currEventDescr->numIp; i++)
if(currEventDescr->clients[i]->id != -1)
status = MdsValue(currEventDescr->clients[i]->id, "RtEventsShr->EventTriggerExecute($1,$2,$3)", &evNameDsc, &bufDsc, &bufSizeDsc, &resDsc, NULL);
}
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;
}
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);
}
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);
}
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;
}