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 imppla_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_GSE","Np",NULL};
struct RCDFData *data=NULL;
int max;
int rmax;
int status;
double gsex,gsey,gsez;
double gsmx,gsmy,gsmz;
double *dptr;
float *fptr;
float Mp=1.6e-27;
float Vx;
int pcnt=0;
cnt=ptr->cnt;
rmax=RCDFMaxRecZ(id,varlist[0]);
max=rmax+cnt;
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->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=RCDFReadZ(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;
if (fptr[0]>9999) {
ptr->VGSE[cnt*3]=1e31;
ptr->VGSE[cnt*3+1]=1e31;
ptr->VGSE[cnt*3+2]=1e31;
ptr->VGSM[cnt*3]=1e31;
ptr->VGSM[cnt*3+1]=1e31;
ptr->VGSM[cnt*3+2]=1e31;
ptr->vth[cnt]=1e31;
ptr->den[cnt]=1e31;
ptr->pre[cnt]=1e31;
cnt++;
continue;
}
ptr->VGSE[cnt*3]=fptr[0];
ptr->VGSE[cnt*3+1]=fptr[1];
ptr->VGSE[cnt*3+2]=fptr[2];
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->VGSM[cnt*3]=gsmx;
ptr->VGSM[cnt*3+1]=gsmy;
ptr->VGSM[cnt*3+2]=gsmz;
ptr->vth[cnt]=0;
fptr=(float *) data[2].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);
RCDFFree(data,3);
ptr->cnt=cnt;
return 0;
}