基于临界平衡假设的太阳风湍流各向异性的观测分析.doc

1、中国地质大学（北京）本科生毕业论文分类号密级中国地质大学（北京）本科毕业设计题目基于临界平衡假设的太阳风湍流各向异性的观测分析英文题目OBSERVATIONTHESOLARWINDTURBULENCEANISOTROPYBASEDONTHECRITICALBALANCEASSUMPTION学生姓名院（系）地球物理与信息技术专业地球物理学学号10108222指导教师职称讲师，研究员2012年6月中国地质大学（北京）本科生毕业论文摘要ELSASSER于1950年对磁流体方程引入新变量后GOLDREICH和SHIRIDA于1995年提出临界平衡理论假设。并将理论应用于太阳风湍流中得到功率谱各向异性特

2、征和功率谱的谱指标特征。本文的主要目的是做出速度场的功率谱，进一步检验临界平衡理论假设的正确性与适用性。文中采用了安装在WIND卫星上的磁场探测器MFI和三维等离子体探测器3DP在1995年1月31日所采集到的总磁场和等离子体速度场数据。并且采用基于小波变换的功率谱分析，做出速度场功率谱密度函数图和磁场功率谱密度函数图。其中横轴为局部磁场和太阳径向向外方向的夹角；纵轴为时间尺度。通过功率谱密度函数图的分析可以得出功率谱随角度变化的分布特征。通过计算机的计算结果表明当天的速度场功率谱显式出各向同性的特征，同时磁场的功率谱也同样显示出了各向同性的特征。这说明临界平衡理论假设并不是普遍适用的。从原始

3、数据看，由于速度扰动和速度化磁场扰动高度相关，笔者认为这可能是产生各向同性的原因之一。遗憾的是本文只针对1995年1月31日的数据进行了考察，无法同其他明显存在阿尔文波时间段的功率谱进行比较。同时当功率谱为各向同性时，功率谱谱指标特征也有待进一步探究。关键词临界平衡；太阳风湍流；功率谱；小波分析；各向异性；中国地质大学（北京）本科生毕业论文ABSTRACTGOLDREICHANDSHIRIDAIN1995,PROPOSEACRITICALBALANCEOFTHEORETICALASSUMPTIONSAFTERELSASSERINTERODUCEANEWVARIBLESINTHEMHDEQUAT

4、IONIN1950THETHEORYCANBEAPPLIEDTOTHESOLARWINDTURBULENCECHARACTERISTICSOFTHEANISOTROPYPOWERSPECTRUMANDCHARACTERISTICSOFTHEPOWERSPECTRUMSPECTRALINDEXTHEMAINPURPOSEOFTHISPAPERISTOMAKETHEPOWERSPECTRUMOFTHEVELOCITYFIELD,FURTHERTESTTHEACCURACYANDAPPLICABILITYOFTHECRITICALBALANCEOFTHEORETICALASSUMPTIONSTHIS

5、PAPERUSESAMAGNETICFIELDANDPLASMAVELOCITYFIELDDATAONJANUARY31,1995COLLECTEDBYMFIAND3DPWHICHINSTALLEDONTHEWINDSPACECRAFTANDBASEDONWAVELETTRANSFORMPOWERSPECTRUMANALYSIS,TOMAKETHEPOWERSPECTRALDENSITYOFTHEVELOCITYFIELDMAPSANDMAGNETICFIELDPOWERSPECTRALDENSITYMAPWHERETHEHORIZONTALAXISISTHEANGLEBETWEENTHELO

6、CALMAGNETICFIELDANDSOLARRADIALOUTWARDDIRECTION,THELONGITUDINALAXISISTHETIMESCALEBYTHEPOWERSPECTRALDENSITYMAPSCANBEEASILYOBSERVEDPOWERSPECTRUMCHANGESWITHTHEANGLEDISTRIBUTIONCHARACTERISTICSCOMPUTATIONALRESULTSSHOWTHATDAYTHEVELOCITYFIELDPOWERSPECTRUMSHOWSISOTROPICCHARACTERISTICSOFTHEMAGNETICFIELDPOWERS

7、PECTRUMALSOSHOWSISOTROPICCHARACTERISTICSTHISSHOWSTHATTHETHEORETICALASSUMPTIONSOFTHECRITICALBALANCEISNOTGENERALLYAPPLICABLEFROMTHEDATA,HIGHLYCORRELATEDTOTHEVELOCITYDISTURBANCEANDTHESPEEDOFTHEMAGNETICFIELDDISTURBANCE,ITHINKTHISMAYBEONEOFTHEREASONSTOPRODUCEISOTROPICUNFORTUNATELY,THISPAPERONLYEXAMINEDDA

8、TAFORJANUARY31,1995,CANNOTBECOMPAREDWITHOTHERPOWERSPECTRUMWHENTHEALFVENWAVESISAPPARENTEXISTENCEANDWHENTHEPOWERSPECTRUMISISOTROPIC,THECHARACTERISTICSOFTHEPOWERSPECTRUMSPECTRALINDEXALSONEEDSFURTHEREXPLORATIONKEYWORDSCRITICALBALANCESOLARWINDTURBULENCEPOWERSPECTRUMWAVELETANALYSISANISOTROPY中国地质大学（北京）本科生毕

9、业论文目录前言错误未定义书签。1背景介绍错误未定义书签。11太阳大气与太阳风错误未定义书签。12太阳风及湍流313阿尔文波414临界平衡515太阳风湍流功率谱的各向异性816卫星介绍102数据及数据处理原理1错误未定义书签。21数据来源及处理过程1322小波分析数据原理13221小波分析原理13222实际工作中小波分析的运用14223功率谱密度PSD1523功率谱随角度的分布163数据处理及处理结果错误未定义书签。831太阳风中的时间序列错误未定义书签。832磁场数据处理错误未定义书签。933速度场数据处理20结论21致谢22参考文献23附录一24附录二32附录三36中国地质大学（北京）本科生

10、毕业论文1前言1950年，ELSASSER对磁流体方程进行了改进，引入了ELSASSER新变量。新变量将磁流体的磁场和动量结合在了一起，不仅将磁流体的动量方程和磁感应方程结合在一起，同时新方程“显式”的描述了反向传播的阿尔文波的非线性相互作用。在此基础上GOLDREICH和SHIRIDA于1995年提出了临界平衡理论。他们假设磁流体中波动的传播效应和反向波动相互作用效应相当。并且将理论应用于太阳风中，得到了太阳风湍流各向异性的结论。理论也给出了垂直于磁场方向和平行于磁场方向的功率谱的谱指数分别是5/3和2。TIMOTHYSHORBURY于2008年的一篇文章中给出了垂直于磁场方向和平行于磁场方

11、向的功率谱的谱指数的观测值，结果与理论符合的很好。JJPODESTA于2009年的一篇文章中进一步细致的研究了局部磁场方向与太阳径向向外方向不同夹角不同时的功率谱谱指标特征，结果同样与理论符合的很好。在HEIDLPLOTS,TIME,COI,NOCLIP0YPADRETURNSTHEPADDEDTIMESERIESTHATWASACTUALLYUSEDINTHE中国地质大学（北京）本科生毕业论文25WAVELETTRANSFORMDAUGHTERIFINITIALLYSETTO1,THENRETURNTHEDAUGHTERWAVELETSTHISISACOMPLEXARRAYOFTHESAME

12、SIZEASWAVELETATEACHSCALETHEDAUGHTERWAVELETISLOCATEDINTHECENTEROFTHEARRAYSIGNIFOUTPUTSIGNIFICANCELEVELSASAFUNCTIONOFPERIODFFT_THEOROUTPUTTHEORETICALBACKGROUNDSPECTRUMSMOOTHEDBYTHEWAVELETFUNCTION,ASAFUNCTIONOFPERIODDEFUNCTINPUTSOCTTHEOFOCTAVESTOANALYZEOVERLARGESTSCALEWILLBES02OCTDEFAULTISLOG2N1VOICEOF

13、VOICESINEACHOCTAVEDEFAULTIS8HIGHERGIVESBETTERSCALERESOLUTION,BUTISSLOWERTOPLOTEXAMPLEIDLNTIME256IDLYRANDOMNS,NTIMECREATEARANDOMTIMESERIESIDLDT025IDLTIMEFINDGENNTIMEDTCREATETHETIMEINDEXIDLIDLWAVEWAVELETY,DT,PERIODPERIOD,COICOI,/PAD,SIGNIFSIGNIFIDLNSCALEN_ELEMENTSPERIODIDLLOADCT,39IDLCONTOUR,ABSWAVE2,

14、TIME,PERIOD,XSTYLE1,XTITLETIME,YTITLEPERIOD,TITLENOISEWAVELET,YRANGEMAXPERIOD,MINPERIOD,LARGESMALLPERIOD/YTYPE,MAKEYAXISLOGARITHMICNLEVELS25,/FILLIDLIDLSIGNIFREBINTRANSPOSESIGNIF,NTIME,NSCALEIDLCONTOUR,ABSWAVE2/SIGNIF,TIME,PERIOD,/OVERPLOT,LEVEL10,C_ANNOT95IDLPLOTS,TIME,COI,NOCLIP0ANYTHING“BELOW“THI

15、SLINEISDUBIOUS中国地质大学（北京）本科生毕业论文26COPYRIGHTC19952004,CHRISTOPHERTORRENCEANDGILBERTPCOMPOTHISSOFTWAREMAYBEUSED,COPIED,ORREDISTRIBUTEDASLONGASITISNOTSOLDANDTHISCOPYRIGHTNOTICEISREPRODUCEDONEACHCOPYMADETHISROUTINEISPROVIDEDASISWITHOUTANYEXPRESSORIMPLIEDWARRANTIESWHATSOEVERNOTICEPLEASEACKNOWLEDGETHEUSEOF

16、THEABOVESOFTWAREINANYPUBLICATIONSWAVELETSOFTWAREWASPROVIDEDBYCTORRENCEANDGCOMPO,ANDISAVAILABLEATURLHTTP/PAOSCOLORADOEDU/RESEARCH/WAVELETS/REFERENCETORRENCE,CANDGPCOMPO,1998APRACTICALGUIDETOWAVELETANALYSISBULLAMERMETEORSOC,79,6178PLEASESENDACOPYOFSUCHPUBLICATIONSTOEITHERCTORRENCEORGCOMPODRCHRISTOPHER

17、TORRENCEDRGILBERTPCOMPORESEARCHSYSTEMS,INCCLIMATEDIAGNOSTICSCENTER4990PEARLEASTCIRCLE325BROADWAYR/CDC1BOULDER,CO80301,USABOULDER,CO803053328,USAEMAILCHRISATRSINCDOTCOMEMAILCOMPOATCOLORADODOTEDUFUNCTIONMORLET,MORLETK0,SCALE,K,PERIOD,COI,DOFMIN,CDELTA,PSI0IFK0EQ1THENK06DNN_ELEMENTSKEXPNTSCALEKK02/2DKG

18、T0DT2PI/NK1NORMSQRT2PISCALE/DTPI025TOTALENERGYNEQN7MORLETNORMEXPEXPNT100DMORLETMORLETEXPNTGT100AVOIDUNDERFLOWERRORSMORLETMORLETKGT0HEAVISIDESTEPFUNCTIONMORLETISCOMPLEXFOURIER_FACTOR4PI/K0SQRT2K02SCALEFOURIERSEC3HPERIODSCALEFOURIER_FACTORCOIFOURIER_FACTOR/SQRT2CONEOFINFLUENCESEC3GDOFMIN2DEGREESOFFREE

19、DOMWITHNOSMOOTHINGCDELTA1IFK0EQ6THENCDELTA0776RECONSTRUCTIONFACTORPSI0PI025PRINT,SCALE,N,SQRTTOTALABSMORLET2,/DOUBLE中国地质大学（北京）本科生毕业论文27RETURN,MORLETENDFUNCTIONPAUL,PAULM,SCALE,K,PERIOD,COI,DOFMIN,CDELTA,PSI0IFMEQ1THENM4DNN_ELEMENTSKEXPNTSCALEKKGT0DT2DPI/NK1NORMSQRT2PISCALE/DT2M/SQRTMFACTORIAL2M1PAUL

20、NORMSCALEKMEXPEXPNT100DEXPNTGT100PAULPAULKGT0FOURIER_FACTOR4PI/2M1PERIODSCALEFOURIER_FACTORCOIFOURIER_FACTORSQRT2DOFMIN2DEGREESOFFREEDOMWITHNOSMOOTHINGCDELTA1IFMEQ4THENCDELTA1132RECONSTRUCTIONFACTORPSI02MFACTORIALM/SQRTPIFACTORIAL2MPRINT,SCALE,N,NORM,SQRTTOTALPAUL2,/DOUBLESQRTNRETURN,PAULENDFUNCTION

21、DOG,DOGM,SCALE,K,PERIOD,COI,DOFMIN,CDELTA,PSI0IFMEQ1THENM2NN_ELEMENTSKEXPNTSCALEK2/2DDT2DPI/NK1NORMSQRT2PISCALE/DTSQRT1D/GAMMAM05IDCOMPLEX0,1GAUSSNORMIMSCALEKMEXPEXPNT100DEXPNTGT100FOURIER_FACTOR2PISQRT2/2M1PERIODSCALEFOURIER_FACTORCOIFOURIER_FACTOR/SQRT2DOFMIN1DEGREESOFFREEDOMWITHNOSMOOTHINGCDELTA1

22、PSI01IFMEQ2THENBEGINCDELTA3541RECONSTRUCTIONFACTORPSI00867325ENDIF中国地质大学（北京）本科生毕业论文28IFMEQ6THENBEGINCDELTA1966RECONSTRUCTIONFACTORPSI0088406ENDIFPRINT,SCALE,N,NORM,SQRTTOTALABSGAUSS2,/DOUBLESQRTNRETURN,GAUSSENDWAVELETFUNCTIONWAVELET,Y1,DT,REQUIREDINPUTSS0S0,DJDJ,JJ,OPTIONALINPUTSPADPAD,MOTHERMOTHER,

23、PARAMPARAM,VERBOSEVERBOSE,NO_WAVENO_WAVE,RECONRECON,LAG1LAG1,SIGLVLSIGLVL,DOFDOF,GLOBALGLOBAL,OPTIONALINPUTSSCALESCALE,PERIODPERIOD,YPADYPAD,OPTIONALOUTPUTSDAUGHTERDAUGHTER,COICOI,SIGNIFSIGNIF,FFT_THEORFFT_THEOR,OCTOCT,VOICEVOICEDEFUNCTINPUTSON_ERROR,2RCHECK_MATH0,1NN_ELEMENTSY1N1NBASE2FIXALOGN/ALOG

24、204999POWEROF2NEARESTTONCHECKKEYWORDSEQN10IFN_ELEMENTSMOTHERLT1THENMOTHERMORLETIFN_ELEMENTSPARAMLT1THENPARAM1IFN_ELEMENTSSIGLVLLT1THENSIGLVL095IFN_ELEMENTSLAG1LT1THENLAG100LAG1LAG10VERBOSEKEYWORD_SETVERBOSEDO_DAUGHTERKEYWORD_SETDAUGHTERDO_WAVENOTKEYWORD_SETNO_WAVERECONKEYWORD_SETRECONIFKEYWORD_SETGL

25、OBALTHENMESSAGE,PLEASEUSEWAVE_SIGNIFFORGLOBALSIGNIFICANCETESTS中国地质大学（北京）本科生毕业论文29CONSTRUCTTIMESERIESTOANALYZE,PADIFNECESSARYYPADY1TOTALY1/NREMOVEMEANIFKEYWORD_SETPADTHENBEGINPADWITHEXTRAZEROES,UPTOPOWEROF2YPADYPAD,FLTARR2LBASE21NNN_ELEMENTSYPADENDIFCONSTRUCTSCALEARRAYOFSCALESSCALEDINDGENNADJARRAYOFJ

26、VALUESSCALE2D0SCALES0ARRAYOFSCALES2JEQN9PERIODFLTARRNA,/NOZEROEMPTYPERIODARRAYFILLEDINBELOWWAVECOMPLEXARRN,NA,/NOZEROEMPTYWAVELETARRAYIFDO_DAUGHTERTHENDAUGHTERWAVEEMPTYDAUGHTERARRAYCONSTRUCTWAVENUMBERARRAYUSEDINTRANSFORMEQN5KDINDGENN/212PI/DOUBLENDTK0D,K,REVERSEK0N1/21COMPUTEFFTOFTHEPADDEDTIMESERIES

27、YFFTFFTYPAD,1,/DOUBLEEQN3IFVERBOSETHENBEGINVERBOSEPRINTPRINT,MOTHERPRINT,POINTS,N1,S0,S0,DJ,DJ,J,FIXJIFN1NENTHENPRINT,PADDEDWITH,NN1,ZEROESPRINT,J,SCALE,PERIOD,VARIANCE,MATHFLAG,FORMAT/,A3,3A11,A10ENDIFVERBOSEIFN_ELEMENTSFFT_THEOREQNTHENFFT_THEOR_KFFT_THEORELSEFFT_THEOR_K1LAG12/12LAG1COSKDTLAG12EQN1

28、6FFT_THEORFLTARRNALOOPTHRUEACHSCALEFORA10,NA1DOBEGINSCALEPSI_FFTCALL_FUNCTIONMOTHER,PARAM,SCALEA1,K,PERIOD1,COI,DOFMIN,CDELTA,PSI0IFDO_WAVETHENWAVE,A1FFTYFFTPSI_FFT,1,/DOUBLEWAVELETTRANSFORMEQN4PERIODA1PERIOD1SAVEPERIODFFT_THEORA1TOTALABSPSI_FFT2FFT_THEOR_K/NIFDO_DAUGHTERTHEN中国地质大学（北京）本科生毕业论文30DAUGH

29、TER,A1FFTPSI_FFT,1,/DOUBLESAVEDAUGHTERIFVERBOSETHENPRINT,A1,SCALEA1,PERIODA1,TOTALABSWAVE,A12,CHECK_MATH0,FORMATI3,3F113,I6ENDFORSCALECOICOIFINDGENN11/2,REVERSEFINDGENN1/2DTCOISEC3GIFDO_DAUGHTERTHENSHIFTSODAUGHTERSAREINMIDDLEOFARRAYDAUGHTERDAUGHTERNN1/2,DAUGHTER0N1/21,SIGNIFICANCELEVELSSEC4SDEVMOMEN

30、TY11FFT_THEORSDEVFFT_THEORINCLUDETIMESERIESVARIANCEDOFDOFMINSIGNIFFFT_THEORCHISQR_CVF1SIGLVL,DOF/DOFEQN18IFRECONTHENBEGINRECONSTRUCTIONEQN11IFCDELTAEQ1THENBEGINY11MESSAGE,/INFO,CDELTAUNDEFINED,CANNOTRECONSTRUCTWITHTHISWAVELETENDIFELSEBEGINY1DJSQRTDT/CDELTAPSI0FLOATWAVE1/SQRTSCALEY1Y10N11ENDELSEENDIF

31、RETURN,WAVE0N11,GETRIDOFPADDINGBEFORERETURNINGEND中国地质大学（北京）本科生毕业论文31附录二单位转换程序PROGET_BVECTOR_IN_ALFVENUNIT_VECT_WIND_19950131SUB_DIR_DATE19950131WIND_DATA_DIRWIND_DATA_DIRWIND_FIGURE_DIRWIND_FIGURE_DIRSETENV,WIND_DATA_DIRSETENV,WIND_FIGURE_DIRSTEP1STEP1DIR_RESTOREGETENVWIND_DATA_DIRSUB_DIR_DATEFILE_R

32、ESTOREB_RTN_3S_ARRTIMESAVFILE_ARRAYFILE_SEARCHDIR_RESTOREFILE_RESTORE,COUNTNUM_FILESFORI_FILE0,NUM_FILES1DOBEGINPRINT,I_FILE,FILE,I_FILE,FILE_ARRAYI_FILEENDFORI_SELECT0READ,I_SELECT,I_SELECTFILE_RESTOREFILE_ARRAYI_SELECTRESTORE,FILE_RESTORE,/VERBOSEDATA_DESCRIPGOTFROM“CONVERT_B_GSE_TO_B_RTN_WIND_200

33、107PRO“SAVE,FILENAMEDIR_SAVEFILE_SAVE,DATA_DESCRIP,JULDAY_3S_VECT,BXYZ_GSE_3S_ARR,B_RTN_3S_ARRJULDAY_VECT_BJULDAY_3S_VECTDIR_RESTOREGETENVWIND_DATA_DIRSUB_DIR_DATEMFIFILE_RESTOREV_RTN_3S_ARRTIMESAVFILE_ARRAYFILE_SEARCHDIR_RESTOREFILE_RESTORE,COUNTNUM_FILES中国地质大学（北京）本科生毕业论文32FORI_FILE0,NUM_FILES1DOBE

34、GINPRINT,I_FILE,FILE,I_FILE,FILE_ARRAYI_FILEENDFORI_SELECT0READ,I_SELECT,I_SELECTFILE_RESTOREFILE_ARRAYI_SELECTRESTORE,FILE_RESTORE,/VERBOSEDATA_DESCRIPGOTFROM“CONVERT_V_GSE_TO_V_RTN_WIND_200107PRO“SAVE,FILENAMEDIR_SAVEFILE_SAVE,DATA_DESCRIP,JULDAY_3S_VECT,VXYZ_GSE_3S_ARR,V_RTN_3S_ARRJULDAY_VECT_VJU

35、LDAY_3S_VECTSTRPOS_TMPSTRPOSFILE_RESTORE,TIMESTRPOS_TMP_V2STRPOSFILE_RESTORE,SAVTIME_STRSTRMIDFILE_RESTORE,STRPOS_TMP,STRPOS_TMP_V2STRPOS_TMPDIR_RESTOREGETENVWIND_DATA_DIRSUB_DIR_DATE3DPFILE_RESTOREWI_PM_3DP_SAVFILE_ARRAYFILE_SEARCHDIR_RESTOREFILE_RESTORE,COUNTNUM_FILESFORI_FILE0,NUM_FILES1DOBEGINPR

36、INT,I_FILE,FILE,I_FILE,FILE_ARRAYI_FILEENDFORI_SELECT0READ,I_SELECT,I_SELECTFILE_RESTOREFILE_ARRAYI_SELECTRESTORE,FILE_RESTOREDATA_DESCRIPGOTFROM“READ_PM_3DP_WIND_CDF_19950131PRO“SAVE,FILENAMEDIR_SAVEFILE_SAVE,DATA_DESCRIP,JULDAY_VECT,VXYZ_GSE_ARR,NUMDENS_VECT,TEMPER_VECTSTEP2STEP2GETNEWB/V_R/T/N_VE

37、CT_NEWJULDAY_MINMAXMINJULDAY_VECT_B,MINJULDAY_VECT_VJULDAY_MAXMINMAXJULDAY_VECT_B,MAXJULDAY_VECT_V中国地质大学（北京）本科生毕业论文33NUM_TIMES_BN_ELEMENTSJULDAY_VECT_BNUM_TIMES_VN_ELEMENTSJULDAY_VECT_VDJULDAY_PIX_BMEANJULDAY_VECT_B1NUM_TIMES_B1JULDAY_VECT_B0NUM_TIMES_B2DJULDAY_PIX_VMEANJULDAY_VECT_V1NUM_TIMES_V1JUL

38、DAY_VECT_V0NUM_TIMES_V2DJULDAY_PIXMAXDJULDAY_PIX_B,DJULDAY_PIX_VNUM_TIMESFLOORJULDAY_MAXJULDAY_MIN/DJULDAY_PIX1JULDAY_VECTJULDAY_MINDINDGENNUM_TIMESDJULDAY_PIXB_RTN_3S_ARRTRANSPOSEB_RTN_3S_ARRV_RTN_3S_ARRTRANSPOSEV_RTN_3S_ARRB_R_VECT_OLDB_RTN_3S_ARR,0B_T_VECT_OLDB_RTN_3S_ARR,1B_N_VECT_OLDB_RTN_3S_AR

39、R,2V_R_VECT_OLDV_RTN_3S_ARR,0V_T_VECT_OLDV_RTN_3S_ARR,1V_N_VECT_OLDV_RTN_3S_ARR,2V_R_VECT_OLDREFORMVXYZ_GSE_3S_ARR0,V_T_VECT_OLDREFORMVXYZ_GSE_3S_ARR1,V_N_VECT_OLDREFORMVXYZ_GSE_3S_ARR2,NUMDENS_VECT_OLDNUMDENS_VECTB_R_VECT_NEWINTERPOLB_R_VECT_OLD,JULDAY_VECT_B,JULDAY_VECTB_T_VECT_NEWINTERPOLB_T_VECT

40、_OLD,JULDAY_VECT_B,JULDAY_VECTB_N_VECT_NEWINTERPOLB_N_VECT_OLD,JULDAY_VECT_B,JULDAY_VECTV_R_VECT_NEWINTERPOLV_R_VECT_OLD,JULDAY_VECT_V,JULDAY_VECTV_T_VECT_NEWINTERPOLV_T_VECT_OLD,JULDAY_VECT_V,JULDAY_VECTV_N_VECT_NEWINTERPOLV_N_VECT_OLD,JULDAY_VECT_V,JULDAY_VECTNUMDENS_VECT_NEWINTERPOLNUMDENS_VECT_O

41、LD,JULDAY_VECT_V,JULDAY_VECTSTEP3STEP3GETNEWB_R/T/N_VECT_ALFUNITIS_INSTANT_OR_CONSTANT_N2PRINT,USEINSTANTORCONSTANTDENSITY1OR2,IS_INSTANT_OR_CONSTANT_NIS_CONTINUEREAD,IS_CONTINUE,IS_CONTINUE中国地质大学（北京）本科生毕业论文34IFIS_INSTANT_OR_CONSTANT_NEQ2THENBEGINAVER_NUMDENSMEANNUMDENS_VECT_NEWNUMDENS_VECT_NEW_V2AV

42、ER_NUMDENSFLTARRNUM_TIMESENDIFIFIS_INSTANT_OR_CONSTANT_NEQ1THENBEGINNUMDENS_VECT_NEW_V2NUMDENS_VECT_NEWENDIFB_R_VECT_ALFUNITGET_ALFVENSPEEDB_R_VECT_NEW,NUMDENS_VECT_NEW_V2B_T_VECT_ALFUNITGET_ALFVENSPEEDB_T_VECT_NEW,NUMDENS_VECT_NEW_V2B_N_VECT_ALFUNITGET_ALFVENSPEEDB_N_VECT_NEW,NUMDENS_VECT_NEW_V2VA_

43、R_VECT_NEWB_R_VECT_ALFUNITVA_T_VECT_NEWB_T_VECT_ALFUNITVA_N_VECT_NEWB_N_VECT_ALFUNITDIR_SAVEGETENVWIND_DATA_DIRSUB_DIR_DATEFILE_SAVEVPROREAD_H0_MFI_WIND_CDF_19950131SUB_DIR_DATE19950131STEP1STEP1SETENV,WIND_DATA_DIRDIR_READGETENVWIND_DATA_DIRSUB_DIR_DATEMFIFILE_READWI_H0_MFI_CDFFILE_ARRAYFILE_SEARCH

44、DIR_READFILE_READ,COUNTNUM_FILESFORI_FILE0,NUM_FILES1DOBEGINPRINT,I_FILE,FILE,I_FILE,FILE_ARRAYI_FILEENDFORI_SELECT0READ,I_SELECT,I_SELECTFILE_READFILE_ARRAYI_SELECTFILE_READSTRMIDFILE_READ,STRLENDIR_READ,STRLENFILE_READSTRLENDIR_READCD,CURRENTWK_DIRCD,DIR_READCDF_IDCDF_OPENFILE_READCDF_CONTROL,CDF_

45、ID,VARIABLEEPOCH3,GET_VAR_INFOINFO_EPOCH3ANUM_RECORDSINFO_EPOCHMAXALLOCRECNUM_RECORDSINFO_EPOCH3MAXREC1LCDF_VARGET,CDF_ID,EPOCH3,EPOCH_3S_VECT,REC_COUNTNUM_RECORDSCDF_VARGET,CDF_ID,B3GSE,BXYZ_GSE_3S_ARR,REC_COUNTNUM_RECORDSCDF_CONTROL,CDF_ID,VARIABLEEPOCH,GET_VAR_INFOINFO_EPOCHNUM_RECORDSINFO_EPOCHM

46、AXREC1LCDF_VARGET,CDF_ID,EPOCH,EPOCH_1MIN_VECT,REC_COUNTNUM_RECORDSCDF_VARGET,CDF_ID,BGSE,BXYZ_GSE_1MIN_ARR,REC_COUNTNUM_RECORDSCDF_VARGET,CDF_ID,PGSE,XYZ_GSE_1MIN_ARR,REC_COUNTNUM_RECORDS中国地质大学（北京）本科生毕业论文36CDF_CLOSE,CDF_IDCD,WK_DIREPOCH_3S_VECTREFORMEPOCH_3S_VECTEPOCH_BEGEPOCH_3S_VECT0CDF_EPOCH,EPO

47、CH_BEG,YEAR_BEG,MON_BEG,DAY_BEG,HOUR_BEG,MIN_BEG,SEC_BEG,MILLI_BEG,/BREAKDOWN_EPOCHJULDAY_BEGJULDAYMON_BEG,DAY_BEG,YEAR_BEG,HOUR_BEG,MIN_BEG,SEC_BEGMILLI_BEG1E3JULDAY_3S_VECTJULDAY_BEGEPOCH_3S_VECTEPOCH_BEG/1E3246060EPOCH_1MIN_VECTREFORMEPOCH_1MIN_VECTEPOCH_BEGEPOCH_1MIN_VECT0CDF_EPOCH,EPOCH_BEG,YEA

48、R_BEG,MON_BEG,DAY_BEG,HOUR_BEG,MIN_BEG,SEC_BEG,MILLI_BEG,/BREAKDOWN_EPOCHJULDAY_BEGJULDAYMON_BEG,DAY_BEG,YEAR_BEG,HOUR_BEG,MIN_BEG,SEC_BEGMILLI_BEG1E3JULDAY_1MIN_VECTJULDAY_BEGEPOCH_1MIN_VECTEPOCH_BEG/1E3246060STEP2STEP2BADVAL1E31BX_GSE_3S_VECTREFORMBXYZ_GSE_3S_ARR0,SUB_NANWHEREBX_GSE_3S_VECTEQBADVA

49、LORABSBX_GSE_3S_VECTGE200IFSUB_NAN0NE1THENBEGINSUB_VALWHEREBX_GSE_3S_VECTNEBADVALANDABSBX_GSE_3S_VECTLT200JULDAY_3S_VECT_V2JULDAY_3S_VECTSUB_VALBX_GSE_3S_VECT_V2BX_GSE_3S_VECTSUB_VALBX_GSE_3S_VECTINTERPOLBX_GSE_3S_VECT_V2,JULDAY_3S_VECT_V2,JULDAY_3S_VECTENDIFBY_GSE_3S_VECTREFORMBXYZ_GSE_3S_ARR1,SUB_NANWHEREBY_GSE_3S_VECTEQBADVALORABSBY_GSE_3S_VECTGE200IFSUB_NAN0NE1THENBEGIN中国地质大学（北京）本科生毕业论文37SUB_VALWHEREBY_GSE_3S_VECTNEBADVALAND