int dsel(struct dirent *dp) {
#else
int dsel(const struct dirent *dp) {
#endif
if (dp->d_name[0]=='.') return 0;
return 1;
}
int test_file_epoch(char *fname,double stime,double etime) {
long yr,mo,dy,hr,mt,sc,ms;
double fstime=0,fetime=0;
int max=0;
CDFid id;
CDFstatus status;
double *dptr;
char *varlist[]={"Epoch",NULL};
status=CDFopen(fname,&id);
if (status !=CDF_OK) return -1;
max=RCDFMaxRecR(id,varlist[0]);
if (max>0) {
status=RCDFReadR(id,0,varlist,data);
dptr=(double *) data[0].data;
EPOCHbreakdown(dptr[0],&yr,&mo,&dy,&hr,&mt,&sc,&ms);
fstime=TimeYMDHMSToEpoch(yr,mo,dy,hr,mt,sc+(ms/1000.0));
status=RCDFReadR(id,max-1,varlist,data);
dptr=(double *) data[0].data;
EPOCHbreakdown(dptr[0],&yr,&mo,&dy,&hr,&mt,&sc,&ms);
fetime=TimeYMDHMSToEpoch(yr,mo,dy,hr,mt,sc+(ms/1000.0));
}
max=RCDFMaxRecZ(id,varlist[0]);
if (max>0) {
status=RCDFReadZ(id,0,varlist,data);
dptr=(double *) data[0].data;
EPOCHbreakdown(dptr[0],&yr,&mo,&dy,&hr,&mt,&sc,&ms);
fstime=TimeYMDHMSToEpoch(yr,mo,dy,hr,mt,sc+(ms/1000.0));
status=RCDFReadZ(id,max-1,varlist,data);
dptr=(double *) data[0].data;
EPOCHbreakdown(dptr[0],&yr,&mo,&dy,&hr,&mt,&sc,&ms);
fetime=TimeYMDHMSToEpoch(yr,mo,dy,hr,mt,sc+(ms/1000.0));
}
CDFclose(&id);
if (stime>fetime) return 0;
if (etime<fstime) return 0;
return 1;
}
int impmag_pos(CDFid id,struct posdata *ptr,double stime,double etime) {
long yr,mo,dy,hr,mt,sc,ms;
double rtime;
int i,cnt=0;
char *varlist[]={"Epoch","SC_pos_sm","SC_pos_se",0,NULL};
struct RCDFData *data=NULL;
int max;
int rmax;
int status;
double *dptr;
float *fptr;
strcpy(ptr->sat,"i8");
strcpy(ptr->ins,"mag");
cnt=ptr->cnt;
rmax=RCDFMaxRecR(id,varlist[0]);
max=cnt+rmax;
if (max==0) return -1;
data=RCDFMake(3);
if (ptr->time==NULL) ptr->time=malloc(sizeof(double)*max);
else ptr->time=realloc(ptr->time,sizeof(double)*max);
if (ptr->PGSM==NULL) ptr->PGSM=malloc(3*sizeof(float)*max);
else ptr->PGSM=realloc(ptr->PGSM,3*sizeof(float)*max);
if (ptr->PGSE==NULL) ptr->PGSE=malloc(3*sizeof(float)*max);
else ptr->PGSE=realloc(ptr->PGSE,3*sizeof(float)*max);
for (i=0;i<rmax;i++) {
status=RCDFReadR(id,i,varlist,data);
if (status==0) {
dptr=(double *) data[0].data;
EPOCHbreakdown(dptr[0],&yr,&mo,&dy,&hr,&mt,&sc,&ms);
rtime=TimeYMDHMSToEpoch(yr,mo,dy,hr,mt,sc+(ms/1000.0));
if (rtime<stime) continue;
if (rtime>etime) break;
ptr->time[cnt]=rtime;
fptr=(float *) data[1].data;
ptr->PGSM[cnt*3]=fptr[0]/6370.0;
ptr->PGSM[cnt*3+1]=fptr[1]/6370.0;
ptr->PGSM[cnt*3+2]=fptr[2]/6370.0;
fptr=(float *) data[2].data;
ptr->PGSE[cnt*3]=fptr[0]/6370.0;
ptr->PGSE[cnt*3+1]=fptr[1]/6370.0;
ptr->PGSE[cnt*3+2]=fptr[1]/6370.0;
cnt++;
}
}
RCDFFree(data,3);
ptr->cnt=cnt;
return 0;
}
int acemfi_imf(CDFid id,struct imfdata *ptr,double stime,double etime,
int cnv) {
long yr,mo,dy,hr,mt,sc,ms;
double rtime;
int i,cnt=0;
char *varlist[]={"Epoch","BGSEc","BGSM",NULL};
struct RCDFData *data=NULL;
int max;
int rmax;
int status;
double *dptr;
float *fptr;
double gsex,gsey,gsez;
double gsmx,gsmy,gsmz;
if (cnv==1) varlist[2]=NULL;
strcpy(ptr->sat,"ac");
strcpy(ptr->ins,"mfi");
cnt=ptr->cnt;
rmax=RCDFMaxRecR(id,varlist[0]);
max=cnt+rmax;
if (rmax==0) return -1;
data=RCDFMake(3);
if (ptr->time==NULL) ptr->time=malloc(sizeof(double)*max);
else ptr->time=realloc(ptr->time,sizeof(double)*max);
if (ptr->BGSMc==NULL) ptr->BGSMc=malloc(3*sizeof(float)*max);
else ptr->BGSMc=realloc(ptr->BGSMc,3*sizeof(float)*max);
if (ptr->BGSEc==NULL) ptr->BGSEc=malloc(3*sizeof(float)*max);
else ptr->BGSEc=realloc(ptr->BGSEc,3*sizeof(float)*max);
for (i=0;i<rmax;i++) {
status=RCDFReadR(id,i,varlist,data);
if (status==0) {
dptr=(double *) data[0].data;
EPOCHbreakdown(dptr[0],&yr,&mo,&dy,&hr,&mt,&sc,&ms);
rtime=TimeYMDHMSToEpoch(yr,mo,dy,hr,mt,sc+(ms/1000.0));
if (rtime<stime) continue;
if (rtime>etime) break;
ptr->time[cnt]=rtime;
fptr=(float *) data[1].data;
ptr->BGSEc[cnt*3]=fptr[0];
ptr->BGSEc[cnt*3+1]=fptr[1];
ptr->BGSEc[cnt*3+2]=fptr[2];
if (cnv==1) {
gsex=fptr[0];
gsey=fptr[1];
gsez=fptr[2];
GeoPackRecalc(yr,mo,dy,hr,mt,(sc+1.0*ms/1000.0));
GeoPackGseGsm(gsex,gsey,gsez,&gsmx,&gsmy,&gsmz);
ptr->BGSMc[cnt*3]=gsmx;
ptr->BGSMc[cnt*3+1]=gsmy;
ptr->BGSMc[cnt*3+2]=gsmz;
} else {
fptr=(float *) data[2].data;
ptr->BGSMc[cnt*3]=fptr[0];
ptr->BGSMc[cnt*3+1]=fptr[1];
ptr->BGSMc[cnt*3+2]=fptr[2];
}
cnt++;
}
}
RCDFFree(data,3);
ptr->cnt=cnt;
return 0;
}
int windswe_plasma(CDFid id,struct plasmadata *ptr,double stime,double etime) {
long yr,mo,dy,hr,mt,sc,ms;
double rtime;
int i,cnt=0;
char *varlist[]={"Epoch","V_GSM","V_GSE","THERMAL_SPD","Np",NULL};
struct RCDFData *data=NULL;
int max;
int rmax;
int status;
double *dptr;
float *fptr;
float Mp=1.6e-27;
float Vx;
int pcnt=0;
strcpy(ptr->sat,"we");
strcpy(ptr->ins,"swe");
cnt=ptr->cnt;
rmax=RCDFMaxRecR(id,varlist[0]);
max=cnt+rmax;
if (rmax==0) return -1;
data=RCDFMake(5);
if (ptr->time==NULL) ptr->time=malloc(sizeof(double)*max);
else ptr->time=realloc(ptr->time,sizeof(double)*max);
if (ptr->VGSM==NULL) ptr->VGSM=malloc(3*sizeof(float)*max);
else ptr->VGSM=realloc(ptr->VGSM,3*sizeof(float)*max);
if (ptr->VGSE==NULL) ptr->VGSE=malloc(3*sizeof(float)*max);
else ptr->VGSE=realloc(ptr->VGSE,3*sizeof(float)*max);
if (ptr->vth==NULL) ptr->vth=malloc(sizeof(float)*max);
else ptr->vth=realloc(ptr->vth,sizeof(float)*max);
if (ptr->den==NULL) ptr->den=malloc(sizeof(float)*max);
else ptr->den=realloc(ptr->den,sizeof(float)*max);
if (ptr->pre==NULL) ptr->pre=malloc(sizeof(float)*max);
else ptr->pre=realloc(ptr->pre,sizeof(float)*max);
for (i=0;i<rmax;i++) {
status=RCDFReadR(id,i,varlist,data);
if (status==0) {
dptr=(double *) data[0].data;
EPOCHbreakdown(dptr[0],&yr,&mo,&dy,&hr,&mt,&sc,&ms);
rtime=TimeYMDHMSToEpoch(yr,mo,dy,hr,mt,sc+(ms/1000.0));
if (rtime<stime) continue;
if (rtime>etime) break;
ptr->time[cnt]=rtime;
fptr=(float *) data[1].data;
ptr->VGSM[cnt*3]=fptr[0];
ptr->VGSM[cnt*3+1]=fptr[1];
ptr->VGSM[cnt*3+2]=fptr[2];
fptr=(float *) data[2].data;
ptr->VGSE[cnt*3]=fptr[0];
ptr->VGSE[cnt*3+1]=fptr[1];
ptr->VGSE[cnt*3+2]=fptr[2];
fptr=(float *) data[3].data;
ptr->vth[cnt]=fptr[0];
fptr=(float *) data[4].data;
ptr->den[cnt]=fptr[0];
Vx=ptr->VGSM[3*cnt];
if ((ptr->den[cnt]>0) && (Vx > -10000)) {
ptr->pre[cnt]=(ptr->den[cnt]*Mp*1e6)*
fabs(Vx*1e3)*fabs(Vx*1e3)*1e9;
pcnt++;
} else ptr->pre[cnt]=ptr->den[cnt];
cnt++;
}
}
if (pcnt==0) memset(ptr->pre,0,sizeof(float)*cnt);
ptr->cnt=cnt;
return 0;
}