WRF模式简易操作-中文指南

v1.0可编辑可修改WRF模式操作指南TheInstituteofAtmosphericPhysics，ChineseAcademyofSciencesNortheastInstituteofGeographyandAgroecology,ChineseAcademyofSciences中国科学院大气物理争论所中国科学院东北地理与农业生态争论所二○一七年三月二十日00名目\l“_TOC_250014“WRF模式简介 1WRF模式的 装 2\l“_TOC_250013“安装环境 2模式源程序… 2NetCDF函数库的安装… 2标准初始化〔SI〕的安装 6\l“_TOC_250012“WRF模式的安装… 9\l“_TOC_250011“WRF模式与T213 模 式 套… 17\l“_TOC_250010“嵌套方案 17\l“_TOC_250009“嵌套程序设计 17\l“_TOC_250008“编译嵌套程序 21嵌套的实现 22WRF模式系统 的 行… 29抱负大气方案… 29\l“_TOC_250007“真实大气方案 32WRF模式系统 作 卡… 47\l“_TOC_250006“源程序… 47真实大气方案 480\l“_TOC_250005“模式结果的显 示 理 61\l“_TOC_250004“Vis5D格式 61\l“_TOC_250003“MICAPS格式… 62\l“_TOC_250002“GrADS格式 65\l“_TOC_250001“附录1.WRF模式参数配置说明 68\l“_TOC_250000“附录2.T213场库 参 表 7810WRF/调试技术报告WRFWRF(WeatherResearchForecast)模式系统是由很多美国争论部门及大学的由国家自然科学基金和NOAA共同支持。现在，这项打算，得到了很多其他争论WRFWRF模式系统将成为改进从云尺度到天气尺度等不同尺度重要天气特征预不同地理位置的力量。它将很好的适应从抱负化的争论到业务预报等应用的需要，并具有便于进一步加强完善的敏捷性。WRFNCAR20001130200158WRF2001年116日，很快进展了模式的第三次公布，只是2002424第五次版公布模式系统，版本号为。原定于2002年10月份左右的第六次公布，20033202003112120045212004632006130为止最版本为。WRF安装环境硬件环境PCLINUX软件环境UNIX〔LINUX〕操作系统以上Fortran77〕CFortran77,Fortran90C〕MICAPSVIS5DGrADSRIP模式源程序WRFGRIBWRFNetCDF文件猎取WRFNetCDFWRFWRFNetCDFNetCDFNetCDFNetCDF安装环境不同平台要求也不一样，但必需有25兆的可用空间用来解压、编译和运行测试。由于WRF模式的程序承受了C，Fortran77和Fortran90编写，所以也要FORTRAN77，FORTRAN90C〔C＋＋〕编译器。netCDF将会编译并测试相应的语言接口。指定编译环境CCCFLAGS，FFLAGS，CXX，CXXFLAGS，F90F90FLAGS编译环境。环境变量也就是用户外壳程序〔shell〕shell〔例如设置脚本程序〔表1。1编译NetCDF的一些环境变量及其说明变量变量描述Notes假设未指定，设置脚本程序将试图找出适宜的C编译器，如：CCCcc，c89，xlc，orgcc。Fortran编译假设未指定，设置脚本程序将试图找出适宜的Fortran90orFC器Fortran77编译器。假设不需要FortranFC(假设存在)设置为““。Fortran90假设未指定，设置脚本程序将试图找出适宜的Fortran90F90 编译器(假设存在)CXX C++编译器

Fortran90FC假设不需要Fortran90F90““。C+C++接口时，可将CXX设置为““。CFLAGSCPPFLAGS

C选项C预处理参数选项Fortran编

例如“-0“或“-g“例如“-DNDEBUG“是无视声明局部的检测〔assertion〕FFLAGS 译器参数选例如“-O“或“-g“项Fortran90F90FLAGS 编译器参数例如“-O“或“-g“。假设未指定，则承受FFLAGS选项CXXFLAGS

C++编译器标志

例如“-O“或“-g“ARFLAGSNMFLAGSFPPM4FLAGSLIBSFLIBS

混合参数选项

在一些操作平台中可能用到这些参数中的一个或几个。Fortran90Fortran90NetCDFFortran90shell应用程序〔例如设置脚本程序。典型的环境变量有：USERHOMEPATH

steve/home/steve/bin:/usr/bin:/home/steve/bin执行命令行命令“envCshshellsetenvCCc89setenvCPPFLAGS“-0-DNDEBUG”如下设置：exportCC=c89exportCPPFLAGS=”-0-DNDEBUG”省略。否则，参数值只是设置为所要取值得一局部。变量可在命令行中定义和设置，并只在这一行中生效。例如：exportFOO=fooenv|grepFOOFOO=fooFOO=barenv|grepFOOFOO=barenv|grepFOOFOO=fooUNIXshellshell〔IBMSPAIX：exportCC=xlcCPPFLAGS=-DNDEBUG CFLAGS=-O FC=f77\FFLAGS=-O CXX=xlC F90=xlf90F90FLAGS=-qsuffix=f=f90IBMAIXSPNetCDFCC=/bin/xlcCPPFLAGS=-DNDEBUGCFLAGS=-OFC=/bin/xlfFFLAGS=-OF90=/bin/xlf90F90FLAGS=-qsuffix=f=f90 #Note:no“-O“optionCXX=/bin/xlCCXXFLAGS=-OAIXF90误信息为“1501-229Compilationendedduetolackofspace”的错误。通ulimit〔data、堆栈段〔stack〕以及内存段〔memory〕设置成〔unlimitedF9064CFortran-q64”外，还要设置以下环境变量：ARFLAGS=”-X64cru”NMFLAGS=”-X64”编译和安装NetCDF要生成编译netCDF所需的Makefile文件，必需先运行配置脚本程序netCDFsrc首先，选择安装的名目。用这个名目作为配置脚本程序参数“――prefix＝”的指定值。程序默认的安装名目〔prefix参数值〕执行设置脚本格式如下：./configure--prefix=名目名prefixnetCDF设置。检测结果将会显示在标准输出设备上〔如显示器。〔严峻的错误标志会导致下一步编译“make”程序失败。不过，记录文件中的编译器不能启动编译环境的程序接口。“makenetCDFnetCDF“libsrc/Fortran77Fortran90C＋＋C＋＋接口ncgenncdump。假设第一步的编译通过，则可以用命令“maketest”来检验netCDFC、FortranC＋＋的程序，以及应用程序“ncdump”ncgen***”开头的输出显示行中，过程中的编译器和连接器警告会被无视掉。定的名目中。在进展编译安装的时候，可以编写一个shell脚本程序。这样，在安装的时AIX比方取名为：makedistcleanexportCC=”/usr/bin/xlc”exportCPPFLAGS=-DNDEBUGexportCFLAGS=-OexportFC=”/usr/bin/xlf90”exportFFLAGS=”-O3-qfixed”exportF90FLAGS=”-O3-qfree-qsuffix=f=f90”exportF90=”/usr/bin/xlf90”exportCXX=”/usr/bin/xlC”exportCXXFLAGS=-Oconfigure--prefix=/u/wrf_xp/netcdfmakemaketestmakeinstall问题处理具体阅读本安装文档。netCDFSRC。执行命令“makedistclean假设存在“”文件，则将它删除。按上面描述的方法设置相应的环境变量。重执行配置脚本程序。执行命令“make执行命令“maketest执行命令“makeinstallCFortranFortran90旧，不能编译后缀名为“.F90”的程序，导致运行配置脚本时，无法通过对Fortran90Fortran90标准初始化〔SI〕的安装文件和资料的猎取WRFWRF/边值的局部。SI模式区域的定义和创立〔投影参数，范围大小，区域位置等。NetCDF〔如地形、土地利用、土壤类型、月叶面积指数、经度和纬度等。WRF〔包括水平方向和垂直方向〕上，为WRF模式供给初始场和边界条件。码格式的数据文件转换成“”程序能识别的格式。也可以不用“”程序，而是直接将其他数据格式的数据写成“”所能识别的格式即可。下载标准初始化〔SI〕的压缩〔土壤温度等：-10-10-1-10另外，还应当下载如下四个文件夹及其其中的数据文件：landuse_30s/ USGS24-categorylandusefractionstopo_30s/

WMO/FAO16-category2-layersoiltypesgreenfrac_10m/pctland_10m/topo_10m/soiltemp_1deg/landuse_30s/soiltype_bot_30s/soiltype_top_30s/topo_30s/解压后的各名目中所得到的数据文件大多数都是10°×10°的数据块文件。用WRFSI要做相应的设定。安装环境以上的版本〕make〔以上的版本。由于Fortran90C〔C＋＋〕编译器。SINetCDFNetCDFSISISISI的编译和运分别为：SRCROOT这个名目在编译和安装SISI是在用户解压tarwrfsi_SRCROOTSRCROOTperl10SRCROOTINSTALLROOTSISI可执行程序〔包括可执行的二进制文件和可执行的perlINSTALLROOT：bin、dataetcbindata；etcperl件。INSTALLROOT和SRCROOT可以一样，这样的话，INSTALLROOT名目下的三个SRCROOTINSTALLROOTINSTALLROOTSRCROOTINSTALLROOT5MOAD_DATAROOTofAllDomain/SRCROOTMOAD_DATAROOT是格外有用的。它可以使用户在一次编译后，处理不同的模式区域。INSTALLROOT/dataINSTALLROOT/datacdlstatic可。EXT_DATAROOT这个名目在运行SI时使用。它是用来存放一些所谓的“外部”数据。数据extprdlogextprdSIANALPATH)分开，使用两个不同的路径。编译和安装SIsrc/include/〔如FC〕,cpp〔CPPFLAGS-DFORTRANDOUBLEUNDERSCORE,而有些计算机〔比方AlphaUnix工作站〕则要求为-DFORTRANUNDERSCORE。用户还需为选项。perl正确的版本进展安装。netCDFperlNETCDF库函数的正确名目，例如：setenvNETCDF/usr/local/然后，输入以下安装命令：perlnetCDF和Perl的安装名目其他的一些参数选项如下：perl \--srcroot=/u/denglt/wrfsi_src\--installroot=/u/denglt/wrfsi\--dataroot=/u/denglt/wrfsi/data\--path_to_netcdf=/usr/local/netcdf\--path_to_perl=/usr/local/perl\--machine=ibm>logWRFSI，datarootSI用户没有指定dataroot的名目，则会在安装名目〔$installroot〕中创立一个“/data”perl\”前，最好留一个两行不留空格地合并起来，造成语法错误。的名目，并且创立有如下四个可执行文件：-rwxr-xr-x1wrfhelpusers458752Apr1715:36-rwxr-xr-x1wrfhelpusers647168Apr1715:36-rwxr-xr-x1wrfhelpusers614400Apr1715:37-rwxr-xr-x1wrfhelpusers63120Apr1715:37siscan-rwxr-xr-x1wrfhelpusers630221Apr1715:37以上才有)-rwxr-xr-x1wrfhelpusers188416Apr1715:37GRIB〔投影WRFsiscan”是一个调试程序，它可以查看“”和“”所产生的NetCDFwrfstatic_d01。在进展编译安装的时候，可以编写一个shell脚本程序。这样，在安装的时AIX比方取名为：#!/bin/shmy_homeroot=/u/wrf_xp/denglt/wrfsysmy_dataroot=/dgpfs/fs2/wrf_xp/ddeng/wrfsysexportNETCDF=“/u/wrf_xp/denglt/netcdf“exportSOURCE_ROOT=“$my_homeroot/sisrc“exportINSTALLROOT=“$my_homeroot/wrfsi“exportDATAROOT=“$my_dataroot/moad“exportEXT_DATAROOT=“$my_dataroot/moad/extroot“exportTEMPLATES=“$my_dataroot/moad/temp“exportGEOG_DATAROOT=“$my_dataroot/geog“--install_ui=y问题处理假设失败，则可以用makerealclean来去除失败的安装程序。 进入到src名目下的各个子名目〔如grib_prep、hinterp、vinterpgrid〕下，输入make找出编译安装时的错误所在。WRF安装环境有FORTRAN77，FORTRAN90和标准C〔或C＋＋〕编译器。SI的数据格式承受了NetCDFNetCDFWRFnetCDFWRF/netCDFnetCDF用的系统名目〔/usr/local/，则必需指明netCDF函数库以及引用netCDF。安装程序的猎取WRFWRF〔WRFWRFWRF局部则是模式的积分计算。编译环境配置WRFconfigure〔shell：s〕compile〔CshshellcsharchPerlarch/UNIX“makearchAIX块：############################################################ARCH#AIX(single-threaded)FC=xlf_rCC=xlc_r......############################################################ARCH#AIX(single-threaded) (NONETCDF)FC=xlf_rCC=xlc_r......############################################################ARCHAIXSM(OpenMP)#OMP=-qsmp=noauto-qthreadedOMPCPP=-D_OPENMP......############################################################ARCHAIXDM (RSL-IO,IBM-MPI)#DMPARALLEL=1OMP=SFC =xlf_r......############################################################ARCHAIXDM (RSL-IO,IBM-MPI) (NONETCDF)#DMPARALLEL = 1......############################################################ARCHAIXDM/SM (RSL-IO,IBM-MPI,OpenMP)#DMPARALLEL=1OMP=-qsmp=noauto-qthreaded......OpenMPMPIRSLMPI译、带RSLHDF5MPIRSLMPIOpenMPMCELMPIRSLESMFRSLRSLWRF译选项、编译规章等。进展WRF的编译环境配置时，进入WRF模式的顶级名目中，键入./configure〔或更高版本Perl那么设置文件“arch/”中的模板将会被扫描检测，提取相应的配置模块信息，并且会列出一个选择菜单共用户选择：$configurecheckingforperl5...nocheckingforperl...found/usr/bin/perl(perl)WilluseNETCDFindir:/usr/local/netcdfPleaseselectfromamongthefollowingsupportedplatforms.1. AIX(single-threaded,nonesting)2. AIXSM(OpenMP,nonesting)3. AIXDM-Parallel(RSL_LITE,IBM-MPI,Allowsnesting)4. AIXDM-Parallel(RSL,IBM-MPI,allowsnesting)5. AIXDM-Parallel(RSL,IBM-MPI,allowsnesting)(PARALLELHDF5)6. AIXDM-Parallel(RSL_LITE,IBM-MPI,Allowsnesting)(PARALLELHDF5)AIXDM-Parallel/SM-Parallel(notrecommended) (RSL,IBM-MPI,OpenMP,allowsnesting)AIXDM-Parallel (RSL,IBM-MPI,MCEL) May2003,EXPERIMENTALAIXDM-ParallelESMF (RSL,IBM-MPI,ESMFcoupling,nonesting)EXPERIMENTALAIX (single-threaded,nestingusingRSLwithoutMPI)AIX (OpenMP,nestingusingRSLwithoutMPI)Enterselection[1-11]:1～11shell〔如按“Control-C。NETCDFnetCDF假设没有设置环境变量，而设置脚本也没有找到netCDF，那么程序将会要求用netCDF状态来进展编译；不过，由于当前版本的程序只能用netCDF的输出/输入接口，所以必定会造成程序的运行失败。模式名目中还供给了用于卸载的去除程序clean〔Csh类shell脚本〔*.o〕以及编译时产生的一些中间文-acleane命令，则会将模式恢复到模式的安装最原始的状态，而文件“”中的一些修改信息都会丧失〔runtest执行文件删除掉。只是在编译器、编译选项或者编译规章需要变动时，才有必要对“”文件进展修改。在WRFWRF模式承受动态安排数组的方式来开销较大。在IBMSP上编译时，默认的内存设置可能满足不了模式计算时的内可以将“”的原始设置：FCBASEOPTS = -qarch=auto-qfree=f90 -qzerosize-qnosave-qmaxmem=-1-qspill=20000$(FCDEBUG)$(OMP)......LDFLAGS ......wrfio_nf:( cd ../external/io_netcdf ; makeNETCDFPATH=/u/denglt/netcdf-qarch=auto -qfree=f90

FC=“$(SFC)-qzerosize

$(FCDEBUG)-qnosave“ ;/bin/cp ../../inc)wrfio_int:(cd../external/io_int;\make CC=$(CC) FC=“$(SFC) $(FCDEBUG) -qarch=auto-qfree=f90 -qzerosize-qnosave“alldiffwrf):$(RM)$@$(CPP)-I../inc$(CPPFLAGS)$*.F>$*.f$(FC)-c-qfree=f90 -qspillsize=32767-I.$(MODULE_DIRS)$*.f修改为：FCBASEOPTS = -qarch=auto-qfree=f90-bmaxdata:0x-bmaxstack:0x-qzerosize-qnosave-qmaxmem=-1-qspill=20000$(FCDEBUG)$(OMP)......LDFLAGS=-bmaxdata:0x-bmaxstack:0x-qmaxmem=-1......wrfio_nf:(cd../external/io_netcdf; makeNETCDFPATH=/u/denglt/netcdf-qarch=auto -qfree=f90

FC=“$(SFC)-bmaxdata:0x

$(FCDEBUG)-bmaxstack:0x-qmaxmem=-1-qzerosize-qnosave“;/bin/cp ../../inc)wrfio_int:(cd../external/io_int;\make CC=$(CC) FC=“$(SFC) $(FCDEBUG) -qarch=auto-qfree=f90 -bmaxdata:0x -bmaxstack:0x -qmaxmem=-1-qzerosize-qnosave“alldiffwrf):$(RM)$@$(CPP)-I../inc$(CPPFLAGS)$*.F>$*.f$(FC)-c-qfree=f90-bmaxdata:0x-bmaxstack:0x-qmaxmem=-1-qspillsize=32767-I.$(MODULE_DIRS)$*.fcleanarch信息。“-O2archIBMSP2\\\\\\当文件中的ARCHFLAG标志参数包含-DDM_PARALLEL参数值，并且连接好相-DWRF_RSL_IO_DRSL。CFLAGS-DWRF_RSL_IO，-DMAXDOM_MAKE=$(MAX_DOMAINSDMAXPROC_MAKE=$(MAX_PROC)DMAXPROC_MAKE=$(MAX_PROC)，-I../external/RSL/RSL。在设置文件arch/的用于设置分布式内存并行编译的OpenMP编译代码在配置完编译环境变量之后，就可以输入编译命令“compile”进展程序代码编译了。首先，输入“compileUsage:compilewrf compilewrfinrundirtestcases(seeREADME_test_cases):compileem_b_wave (欧拉质量坐标斜压波方案)compileem_grav2d_x (欧拉质量坐标重力波方案)compileem_hill2d_x 质量坐标二维钟形山体方案)compileem_quarter_ss (欧拉质量坐标三维超级单体云体方案)compileem_real (欧拉质量坐标真实大气方案)compileem_squall2d_xcompileem_squall2d_y

compileexp_real (真实大气试验方案)compilenmm_real (欧拉高度坐标真实大气方案)compile-h 〔帮助信息〕WRFem〔比方和〕会有所不同。除了“nmm_real”和“em_real”两个编译方案是真实大气方案外，其他方案都是用抱负数据进展大气模拟试验的方案。编译时，输入带参数的编译命令“compile$compilecase_name比方说要编译质量动力框架下的二维飚线方案时，则可键入命令：$compileem_squall2d_x然后，编译程序启用注册机制：在“tools/directory”cinc文件〔.f。留意，其中名目“frame/”中的一个源代码文件在下一次编译时强制重运行注册机制，则将文件“frame/”移除掉即可。main/且在名目“test/em_squall2d_x”中有他们的连接文件。对于每一个不同抱负大气试验方案，都会不同，这是由于它们建立方案时是编译了不同的模块module_initializedyn_em/”或比方边界环境、热力湍流〔thermalbubble、模式层数以及地形等。假设模式动力框架不同，那么“”也会不同。clean来进展。当“clean”将删除源程序名目中的目标文件〔*.o〕以及编译时产生的一些中间文件，文件“”会保存。当指定“-a”为“clean”的参数时，则会连同文件“configure”以及一些在“external”名目中的库函〔留意，去除程序并不runtest。时，编译程序会在“external”名目中创立一些库函数。些麻烦。F90编译器允许在向外部子程序传递实参时，产生f90型伪形F90〔无需分割成子的。-DDEREF_KLUDGE指示注册机制在构建调用子程序的实参列表时创立一个额WRFT213嵌套方案WRFT213模式作简洁地嵌套。WRFWRF具体实现的方法有两种，一种是从T213场库猎取要素场，然后修改WRFSIT213WRFT213MM5TERRAIN、REGRIDINTERPF等程序，生成MM5MM52WRFWRF但由于MM5WRF模式没有资料同化方案的缺乏。嵌套程序设计T213IBMSPT213*〔度T213gdbrck参数。具体各接口参数说明如下：callgdbrck(ftype,idate,iftime,ielemcd,level,ii,jj,kk,record,data,index_dic,iret)场类型，:暴雨系统为”BJLP”,t213系统为”BJGL’idate:integer idate(4),日期数组，年〔四位、月、日、时次iftime:integer,时效ielemcdinteger,levelintegerlevel(n+1),层次数组，第一元为调用的总层数n,2到n+1level(3),datalevel/2,500,700/iiinteger，纬度格点数jjinteger,kk:integer,nii*jj。T213640*321datareal*8dimensiondata(ii,jj,kk)index_dic:character*80,字符串，模式字典文件全名iret:integer，返回码0：正常6：场没生成或记录不存在9：其他错值得留意的是，由于场库资料中整数和浮点数都是以64位格式存储的，也8-qintsize=8和-qrealsize=8IBMSP在SP:/gpfs/fs2/nwp/gdbt213目录下用index_dic(1:28)=”/u/nwp/gdbt213/dic/index_dic”/asm2/nwpgdb/t213gdb。

目录下用 index_dic(1:32)=”/u/nwp/gdbt213读取T213场库资料的程序设计T213〔10。parse_tableVtableVtableGRIBCode|LevelCode|Level1|Level2|REGRIDName|REGRIDUnits|REGRIDDescription|-----+------+------+------+----------+----------+------------------------------------------+80|1|0||SST|K|SeaSurfaceTemperature|1|1|0||PSFC|Pa|SurfacePressure|2|102|0||PMSL|Pa|Sea-levelPressure|7|1 | 0| |HGT |m |Height |TerrainHeight211|1|0||SKINTEMP|K|Skintemperature(canuseforSSTalso)|229|105|2||T|K|Temperature|At2m230|105|2||RH|%|RelativeHumidity|At2m202|105|10||U|ms{-1}|U|At10m203|105|10||V|ms{-1}|V|At10m7|100|||HGT|m|Height|11|100|*||T|K|Temperature|33|100|*||U|ms{-1}|U|34|100|*||V|ms{-1}|V|52|100|*||RH|%|RelativeHumidity|-----+------+------+------+----------+----------+---------------------------------------------+Vtableparse_table标志就是第一列包含中画线“-”的其次行。也就是说，程序将无视第一次消灭重要，重要的是第三次和第四次必需同时配套消灭。T213场库 调用子程序parse_table，读取变量列表文件Vtable读取要从T213读取参数列表文件读取时间参数及输出格式参数初始化计算时间时间循环：当计算时间在起止时间之内文件翻开标志赋值flopen＝.true.1≤loopj≤变量总数SST场数据，并按要求 是

要读取的变量为地形否要读取的变量为SST否

是翻开地形文件:读取地形数据、callproc_sst 读取场库气象场

关闭地形文件否气象场是否为高斯格点场是callGAUSS2LATLON11否flopen为真是=1 =输出类型＝=2创立文本格式 创立AVN第4类

AVN3calllandsea算区域范围中的水陆标志否

设置输出区域范围是变量为地形call设置地形场相关callcall=1 =3输出类型＝=2call call callout_mm5完毕：变量循环完毕：时间循环完毕311.读取T21332在变量列表Vtable文件中，这些变量信息设计程表格形式，其中第一列是变量的GRIB码，这些GRIB码是T213场库中规定的编码〔详见附录。当要使用WRFSI的程序时，还必需将T213场库的GRIB编码转换成AVNGRIB1001051021面与陆面的地球外表层。第三、四列是相应于具体的次类型的垂直层的值。比方垂直层类型为100，则第三列就是去各个等压面的值，如1000、950、850等等，则第三列用通配符“*”表示，而第四列为空；假设垂直层类型为105，则第三列则101010，1021，0，第四列为空。第五列为变量名称，依据GRIB编码的标准进展变量命名。第六列为变量单42在表格中，每一列都是用竖线“|”分隔。第一列的起始处不要加分隔符“|这些信息会被无视掉。场库的起止时间、输出区域的范围以及输出文件的格式。例如：&filetimespecSTART_YEAR=2003START_MONTH=04START_DAY=07START_HOUR=00START_MINUTE=00START_SECOND=00END_YEAR=2003END_MONTH=04END_DAY=10END_HOUR=00END_MINUTE=00END_SECOND=00INTERVAL=21600START_X=1START_Y=1DELTA_X=320DELTA_Y=/160!OUT_TYPEISAOUTPUTTYPE,THEVALIDVALUEIS:!1:PREGRIDFORMAT!2:HINTERPFORMAT! 3:MM5FORMAT&typespecOUT_TYPE=2/子程序build_hdate“yyyy-mm-dd_hh:mm:ss”的时间串。以起始时间为开头，在场库中找符合条件的时次，进展场库资料的读取，直到用户定义的调用终止时间为止。GAUSS2LATLON〔这个子程序是在陈启英、姚明明供给的转换程序的根底上修改而成。最终模式的输出是依据用于定义的输出格式来输出的。在文件中定义的v1.0可编辑可修改v1.0可编辑可修改34OUT_TYPE就是定义输出格式的量。当OUT_PUT＝1时，输出为SI程序中能识别格式，固然，的程序源代码要稍加修改，这在后面将会讲到。输出时调用子程序out_pregrid，输出文件名为“avn_*；当OUT_TYPE＝2时，输出为SI程序中能识别处理的格式，调用输出子程序 out_hinterp，输出文件名为AVN_FILE:*；当OUT_TYPE＝3时，输出为MM5的前处理程序能识别处理的格式。输出时调用子程序out_mm5来输出，输出文件名为“MM5_FILE:*。嵌套的实现11方案一首先运行程序读取场库资料，将场库资料输出为能够识别处理的格式。WRFSISI，它会调用程序。GRIBWRFSI./src/grib_prep/名目Fortran翻开文件：源代码：if(iuarr(iunit).thenifthencallcopen(iunit,nunit,gribflnm,1,ierr,iprint)else

gribflnm,1,ierr,-1)endififthencalldeallogribierr=2returncall copen(iunit, nunit, endifiuarr(iunit)=nunitendif修改为：if(iuarr(iunit).thenopen(iunit, file=gribflnm,status=”old”,iostat=ierr)nunit=iunitifthenif(iprint==1)print*,”Openfile源代码：callget_sec1(ksec1)callget_sec2(ksec2)

errorcode:”,ierrcalldeallogribierr=2returnendifiuarr(iunit)=nunitendifread(iunit,*,end=9112)(ksec1(i),i=1,24)read(iunit,*,end=9112)(ksec2(i),callginfo)

get_gridinfo(infogrid,

i=1,4)read(iunit,*,end=9112)(ginfo(i),i=1,24)读变量场数据：源程序：callgribdata(DataArray,map%nx*map%ny)修改为：read(iunit,*,end=9112)(DataArray(i),i=1,ii)应的修改：源代码：callcclose(iuarr(nunit1),iprint,ierr)修改为：close(iuarr(nunit1))做完修改后，重编译安装WRFSI源程序或者单独编译程序。36执行脚本文件，它会依次调用和程序。 第一次运行或者重计算的区域时，则在运行程序前还会运行程序。运行成功后，会生成WRF所需要的预处理数据文件real_input*。方案二首先运行程序，读取场库资料，将场库资料输出为能够识别处理的格式。维(nx,ny)实型数组。而前两个记录数据包括的记录元是固定不变的，第三成员及其说明如下：记录1：versionfield,units,decription,level假设proj_flag=0startlon,deltalat,deltalon假设proj_flag=1startlon,dx,dy,truelat1假设proj_flag=3记录3：startloc,startlat,startlon,dx,dy,orient,truelat1,truelat2假设proj_flag=53：startloc,startlat,startlon,dx,dy,orient,truelat11：INTEGER*4 ::version2：CHARACTER(LEN=24) ::hdate。REAL*4 ::xfcst数据块的预报时效，单位为小时。CHARACTER(LEN=32) ::source数据块的描述信息。CHARACTER(LEN=9) ::fieldVtable中的变量名称，而且必需是模式能够识别的名称。CHARACTER(LEN=25) ::units而模式本身不会依据此单位来实行自动转换。CHARACTER(LEN=46) ::description数据块中气象要素的具体描述信息。REAL*4 ::level帕。INTEGER*4 ::nx二维数据块的x方向的范围大小。INTEGER*4 ::ny二维数据块的x方向的范围大小。INTEGER*4 ::proj_flag地图投影的标志变量。其中：0：圆柱等距投影(经/纬格点)1:麦卡托投影3：兰勃特等角投影5：极射赤面投影(PolarStereographic)3：CHARACTER(LEN=8) ::startloc为：‘SWCORNER’:经纬度相对起始点为左下角‘CENTER’: 相对起始点为投影中心REAL*4 ::startlat起始点的纬度。REAL*4 ::startlon起始点的经度。REAL*4 ::deltalat经纬格点的纬度增量REAL*4 ::deltalon经纬格点的经度增量REAL*4 ::dx。REAL*4 ::dy。REAL*4 ::orient形的经度值。REAL*4 ::truelat1投影不变形的第一纬度值。REAL*4 ::truelat2值。WRFSInamelist。 执行脚本文件，它会依次调用和程序。 第一次运行或者重计算的区域时，则在运行程序前还会运行程序。运行成功后，会生成WRF所需要的预处理数据文件real_input*。方案三MM5的格式。维(nx,ny)实型数组。而前两个记录数据包括的记录元是固定不变的，第三成员及其说明如下：记录1：versionunits,decription,level假设proj_flag=0deltalat,deltalon假设proj_flag=1dx,dy,truelat1假设proj_flag=3dx,dy,orient,truelat1,truelat2假设proj_flag=53：startlat,startlon,dx,dy,orient,truelat11：INTEGER*4 ::version数据的版本号，此版本号为3。2：2：40CHARACTER(LEN=24) ::hdate。REAL*4 ::xfcst数据块的预报时效，单位为小时。CHARACTER(LEN=9) ::fieldVtable中的变量名称，而且必需是模式能够识别的名称。CHARACTER(LEN=25) ::units而模式本身不会依据此单位来实行自动转换。CHARACTER(LEN=46) ::description数据块中气象要素的具体描述信息。REAL*4 ::level帕。INTEGER*4 ::nx二维数据块的x方向的范围大小。INTEGER*4 ::ny二维数据块的x方向的范围大小。INTEGER*4 ::proj_flag地图投影的标志变量。其中：0：圆柱等距投影(经/纬格点)1:麦卡托投影3：兰勃特等角投影5：极射赤面投影(PolarStereographic)3：3：41REAL*4 ::startlat起始点的纬度。REAL*4 ::startlon起始点的经度。REAL*4 ::deltalat经纬格点的纬度增量REAL*4 ::deltalon经纬格点的经度增量REAL*4 ::dx。REAL*4 ::dy。REAL*4 ::orient形的经度值。REAL*4 ::truelat1投影不变形的第一纬度值。REAL*4 ::truelat2值。四版本的数据比第三版本的数据多了两个字符型记录元：source和startloc。运行MM5前处理程序，生成MM5的初始输入场和边界条件。 WRFreal_input*。42T213场库数据T213场库数据readt21MM5REGRIDTERRAININTERPFMM52WRFWRF2.WRFT213编译嵌套程序编译嵌套程序时，目前必需用64位形式编译，或者用编译选项“-qintsize=8”和“-qrealsize=8”来编译。另外，场库连接的库函数以及路径必需在编译连接时指定。IBMSP·SP:/gpfs/fs2/nwp/gdbt213-L/u/nwp/gdbt213/lib-lpubindex_dic(1:28)=”/u/nwp/gdbt213/dic/index_dic”·场库在盘阵库/asm2/nwpgdb/t213gdb-L/u/nwp/gdbt213/lib-lpub_bxindex_dic(1:32)=”/u/nwp/gdbt213/dic/index_nfs_dic’编译时，依据以下步骤进展：进入到t213/mod名目下，输入make命令编译嵌套程序所需的模块程序；进入到t213/util名目下，输入make命令编译一些关心子程序；编译成功后，可产生一个可执行程序。WRF模式系统的运行抱负大气方案运行流程WRF数据初始化程序WRF数据初始化程序运行WRF模式主程序wrfout后处理图3.WRF模式系统抱负大气方案的运行流程运行步骤以进入相应方案的名目下，如欧拉质量坐标x-z平面二维飚线方案对应的名目为./test/em_squall2d_x。 设置该名目下的参数配置文件；具体参数说明参见《WRF模式技术报告》中1时间步长，物理参数以前其他的一些参数。一些参数〔比方时间步长，积分〕的转变，可以不需要重执行“”程序。然而，当转变计算范围，边界条件，以及物理参数时，则必需重执行“”程序。 1个文件：wrfbdy_d01 wrfinput_d01：WRF 以下形式的结果文件：wrfout_d01_000000另外，还有一些程序运行时记录运行状态的中间记录文件。WRF目前，已经初步地将WRFVIS5D、GrADSMICAPSWRF运行个例〔比方时间步长，积分时间和输〕的转变，可以不需要重执行“”程序。然而，当转变计算范围，边界条件，以及物理参数时，则必需重执行“”程序。编辑原始程序代码，比方编辑dyn_eh/，已修改山体或圆形障碍物，或修改外围环境时，应当在WRFcompileeh_squall2d_xx-zKessler粘性系数固定为300m2/s。并利用y方向的周期性边条件，从而使用三维的WRFy理量的计算输出。经过以上步骤地运行，得到结果文件有：-rw-r--r--1dengltnumerica7806948Apr1406:33wrfbdy-rw-r--r--1dengltnumerica1630620Apr1406:33wrfinput-rw-r--r--1dengltnumericaApr1414:43wrfout_01_000000-rw-r--r--1dengltnumerica7381448Apr1414:05wrfrst_01_000600-rw-r--r--1dengltnumerica7381448Apr1414:43wrfrst_01_001200间文件，它的作用是使模式可以分各个时间段进展计算。wrfout_01_000000文件是模式的输出结果文件，从其中可以得到模式的各物理量的计算结果。图4是欧拉质量坐标二维(x-z)飚线方案模拟的水平风速的分布图。从模拟的图中，可以分析得出模式模拟的结果根本正确。v1.0可编辑可修改4.欧拉质量坐标二维(x-z)飚线方案模拟水平风速U4748真实大气方案运行流程数据预处理数据预处理运行WRF数据初始化程序运行WRF数据初始化程序运行WRF模式主程序wrfout后处理图5.WRF模式系统真实大气方案的运行流程运行步骤数据预处理WRFT213T213WRF运行步骤：第一步：设置变量列表文件“Vtable”T213气象场的相关信息。其格式如下：GRIB码|层次代号|第一层|其次层|变量名称|变量单位|变量描述|----------+-------------+-----------+-----------+--------------+-------------+---------------------||11|100|* ||11|100|* ||T33|100|* ||U||ms{-1} |U34 |100 | * | | V |ms{-1} |V|2 |102 | 0 | | PMSL |Pa |Sea-levelP |......-----------+-------------+-----------+-----------+---------------+------------+-------------------+Vtableparse_tableparse_table标志就是第一列包含中画线“-”的其次行。也就是说，程序将无视第一次消灭重要，重要的是第三次和第四次必需同时配套消灭。在变量列表Vtable文件中，这些变量信息设计程表格形式，其中第一列是GRIBGRIBT213〔2。当要使WRFSIT213GRIBAVN其次列是气象要是变量所处的垂直层的GRIB编码。其中100指气压层，105指四列是相应于具体的次类型的垂直层的值。比方垂直层类型为100，则第三列就1000、950、850*”表示，101010，而第四列为空；假设10210，第四列为空。第五列为变量名称，依据GRIBGRIB编码的单位标42|”分隔。第一列的起始处不要加分隔符“|必需以分隔符完毕。后面可以跟一些说明文字信息，但这些信息会被无视掉。其次步：设置参数配置文件“”场库的起止时间、输出区域的范围以及输出文件的格式。例如：&filetimespecSTART_YEAR=2003START_MONTH=04START_DAY=07START_HOUR=00START_MINUTE=00START_SECOND=00END_YEAR=2003END_MONTH=04END_DAY=10END_HOUR=00END_MINUTE=00END_SECOND=00INTERVAL=21600START_X=1START_Y=1DELTA_X=320DELTA_Y=160/!OUT_TYPEISAOUTPUTTYPE,THEVALIDVALUEIS:! 1:PREGRIDFORMAT! 2:HINTERPFORMAT! 3:MM5FORMAT&typespecOUT_TYPE=2/时间间隔和区域大小的。其中，START_YEARSTART_MONTHSTART_DAY，START_HOUR，START_MINUTE，START_SECOND分别是表示起始年、月、日、时、START_XSTART_YX〔纬向〕Y〔经向〕的起始格点；DELTA_XDELTA_YX〔纬向〕和Y〔经向〕的总格点数。&tyepspec局部是定义猎取资料的输出类型。当OUT_PUT＝1时，输出为SI程序中修改后的能识别处理的格式。输出时调用子程序out_pregriavn_*OUT_TYPE＝2SIout_hinterp，AVN_FILE:*当OUT_TYPE＝3MM5out_mm5MM5_FILE:*第三步：运行嵌套程序SIavn_*”或“AVN_FILE:*MM5程序所能识别的数据文件“MM5_FILE:*第四步：运行标准初始化程序SI上面的步骤中第一步到第三步是将WRF模式本地化，与T213模式嵌套时实行SI。不过，这两种方案在运行操作上还是存在一些差异的，运行时需要加SIMM5MM5WRFWRFWRF跳过前面三个步骤，而从第本次步骤开头运行。转到$MOAD_DATAROOT/static/名目中，设置参数配置文件。SI的参数配置文件中，一共有十个列表子块，但其中有三个子块目前没有使laps_analysis_control.其他子块列表信息说明如下：project_id:此，此子块是为了用户在描述元数据便利而设的。SIMULATION_NAME：字符串型，指定本次运行试验的名称USER_DESC：字符串型，描述运行此配置程序的用户信息filetimespec：隔可以修改，其他局部都会在运行命令〔〕时被修改，但当用户间。start_year(例如2001)：四位数字表示的起始年份。start_month(例如04)：两位数字(01-12)表示的起始月份。start_day(例如20)：两位数字(01-31)表示的起始天数。start_hour(例如12)：两位数字(00-23)表示的起始小时数。start_second(例如00)：两位数字(00-59)表示的起始秒数。end_year(例如2001)：四位数字表示的终止年份。end_month(例如04)：两位数字(01-12)表示的终止月份。end_day(例如21)：两位数字(01-31)表示的终止天数。end_hour(例如00)：两位数字(00-23)表示的终止小时数。end_minute(例如00)：两位数字(00-59)表示的终止分钟数。end_second(例如00)：两位数字(00-59)表示的终止秒数。interval_seconds43200)：前处理程序的两次输出时间之间的时间间隔，以秒为单位。hgridspec:模式的水平区域和格点定义信息，包括xdim,ydim,map_proj_name,moad_known_lat,moad_kown_lon,moad_stand_lats,moad_stand_lons，moad_delta_x,moad_delta_y(moad_delta_x=moad_delta_y)WRF则还需对$MOAD_DATAROOT/cdl/文件做相应的修改。整型数，嵌套区域数目。假设不使用嵌套，则此1。XDIM/YDIM：〔母区域〕东西方向和南北方向的格点数。1，即有嵌套时，在中也仅设置一个母区域的Xxdimydim动计算〔5。Arakawa-C3〔U，V〕U〔XDIM，YDIM〕=〔4，4：+V+V+V+UTUTUTU+V+V+V+UTUTUTU+V+V+V+UTUTUTU+V+V+V+其中，加号〔＋〕格点是在参数列表文件namelistWRFUWRF的输出中包含这些格点值〔VVWRFWRFSI设置了区域范围为〔XDIM，YDIM，则模式会输出：质量变量 U动量变量〔XDIM，YDIM－1〕V〔XDIM－1，YDIM〕PARENT_ID：整型数，嵌套区域的上一级区域的编号。留意，PARENT_ID1。其次列必需等于1。总列数必需等于NUM_DOMAINS典型值为2和3，但是也可以更大一些。该参数的列数必需等于NUM_DOMAINS。DOMAIN_ORIGIN_[LLI/LLJ/URI/URJ]：这四个参数定义每个子区域在其母区域中西南角〔[LLC]〕和东北角〔右上[URC]〕的i/jMOADSW1,1；NExdim,ydimxdim,ydim值可由下式计算〔SI：xdim=(uri-lli)*ratio_to_parent+1,ydim=(urj-llj)*ratio_to_parent+1.字符串型，指定地图投影的类型。有效的类型名称包括：“polar”->极射赤面投影“lambert”->兰勃特等角投影〔相割或者相切〕“mercator”->麦卡托投影MOAD_KNOWN_LAT/MOAD_KNOWN_LON：实型数，网格中的指定点的纬度/经度。对于目前来说，这个指定点必需是网格的中心点。90～90－180～180MOAD_STAND_LATS：2在－90～90个值必需是+/-90南极。值。格中心点的纬度。其次个值无定义。MOAD_STAND_LONS：实型数，投影的不变形经度值，取值范围为－180～180，单位为度。通常，将它取为网格中心点的经度。MOAD_DELTA_X/MOAD_DELTA_Y：浮点型数，指定网格在东西方向/南北方向上的格距，单位为米。目前，这两个值必需一样。嵌套子区域的格距由SI和GUI自动计算(见文件src/lib/functiongrid_spacing_wrf_m(nest))：grid_spacing_wrf_m=moad_delta_xi=nest:dowhile grid_spacing_wrf_m=grid_spacing_wrf_m/ratio_to_parent(i)i=parent_id(i)enddoSILAVWT_PARM_WRF：参见README_toptwvl_silavwtinsrc/gridTOPTWVL_PARM_WRF：参见README_toptwvl_silavwtinsrc/gridsfcfiles:指定静态数据〔地形高度，土地利用等〕存放的正确名目。TOPO_30S：USGS-derived30PCTLAND_10M：10〔陆地在整个格点中占有的百分比〕率数据的路径LANDUSE_30S：30USGS24FAO3016类型数据的路径FAO3016类型数GREENFRAC：指定植被分类数据的路径。径。指定月平均的气候反照率资料的路径〔归一化到normalizedtolocalzenith。SSTEMP：海表温度气候值。目前没有使用。interp_control:ptop_Pa,hinterp_method,output_coord,levels.output_coord=“ETAP”用于质量坐标，而output_coord=”ZETA”用于高度坐标。levelsETAP～0，而在ZETA0～xx19000ptop_Pa=5000。NUM_ACTIVE_SUBNESTS:假设ChgridspecNUM_DOMAINS1〔MOAD域进展水平和垂直插值来产生初始条件。假设仅处理MOAD，则NUM_ACTIVE_SUBNESTS必需设为0。NUM_ACTIVE_SUBNESTS设定ACTIVE_SUBNESTS0NUM_DOMAINS-1(含)之间.ACTIVE_SUBNESTS:NUM_ACTIVE_SUBNESTS1，该参数为NUM_ACTIVE_SUBNESTSID>=2以MOAD区时(NUM_ACTIVE_SUBNEST0)，NUM_ACTIVE_SUBNESTSWRF区域来处理。PTOP_PA:5000Pa。HINTERP_METHOD:整型数，指定大气物理量的水平插值方法，有效的取值为：0: 取邻近值，最近的格点值，不推举使用。1: 42: 16LSM_HINTERP_METHOD:整型数，指定做过水陆叠置分析的气象值方法常用于陆面过程需要的气象场有效取值同上，但推举值用0或者1。留意：假设想使用背景资料的土地利用和土地分类，该参数必需设为 0，而且必需由输入数据集中猎取“VEGCAT“和“SOILCAT“，其中VEGCAT为主要的土壤利用分类〔USGS为24〕2SOILCATFAO2维数组。NUM_INIT_TIMES:整型数，用来掌握初始场(INIT_ROOT)和陆面过程气象场(LSM_ROOT)1。此参nudging”预备的。程序会首先利用INIT_ROOT/LSM_ROOT中指定的时次(NUM_INIT_TIMES)数据，然LBC_ROOT中取剩余的时次数据。假设将此参数设为0，那么全部的数据都会从LBC_ROOTCONSTANTS_FULL_NAME1式使用不同模式的数据场来做初始场/陆面过程气象场和边界条件。INIT_ROOT:字符串型，指定初始场数据文件名前缀，在指定了NUM_INIT_TIMES时使用。脚本程序会在 ANALPATH路径(见SI_PATHS子块)中查找带有此前缀以及相应时间串的数据文件此参数只有当NUM_INIT_TIMES>0时有用。 LBC_ROOT:字符串型，指定边界场数据文件名前缀。脚本程序会在LBCPATHSI_PATHS子块)中查找带有此前缀以及相应时间串的数据文件。 LSM_ROOT:字符串型，指定陆面过程气象场数据文件名前缀，，NUM_INIT_TIMES>0。脚本程序会在LSMPATHSI_PATHS子块)中查找带有此前缀以及相应时间串的数据文件。此参数是用来供给来自不同模式(与INIT_ROOTNOAHLSM CONSTANTS_FULL_NAME:字符串列表，指定CONSTANTS_PATHLSM_ROOT/INIT_ROOT/LBC_ROOT VERBOSE_LOG:规律型，调试信息开关。设为真时，能够打印很多具体的运行信息。4，指定数据输出的垂直坐标类型。有效的取值为：WRF”NMMH”:LEVELSNMMWRFNCEPNMM形式的输出数据。WRF不再支持该坐标，选此不能正常运行！实型数列表，WRF”ZETA”，那么这OUTPUT_COORD“ETAP“，那么这些值则从到。2〔放弃〕si_paths:GRIB路径名。SI〔版GRIB〔T213〕的存放路径时，假设指定的输入模式为AVN或T213，则只需设置ETA，则需设置“path_to_ETA”路径。path_to_AVNEXTDATAROOTSI，则需设置如下路径，但大多数状况下他们都会设成一个一样的路径〔$EXT_DATAROOT/extprd。 NUM_INIT_TIMES>0INIT_ROOT LBC_ROOTLSMPATH: 字符串类型，指定0～NUM_INIT_TIMES时次带前缀:LSM_ROOT的文件所在的名目。 CONSTANTS_PATH: 字符串类型，指定全部时次数带前缀CONSTANTS_FULL_NAMEanalysis_control:确认是否是正确设置为“wrfsi〔放弃〕设置好计算数据生成的间隔，文件存放路径以及计算区域配置等信息。$MOAD_DATAROOT/cdl/名目下的作相应的修改。②转到$INSTALLROOT/etc/名目下，运行文件INSTALLROOT和MOAD_DATAROOT：export INSTALLROOT=/u/denglt/wrfsiexport 脚本文件要调用两个另外的脚本文件：和。ETAAVNGRIBT213一种能被程序识别的中间格式，放入到$EXTDAROOT\extprd名目中“si_paths”中指定)。：调用和和程序。文件需要四个命令行输入参数：yyyymmdd，hh〕――预报长度〔单位为小时〕――输入数据源〔ETAAVN〕例如：1212AVNwrfsitestEXTDATAROOT./wrok/：工作名目./log/: 运行文件的记录文件wrfsiextprd-rw-rw-r--1wrfhlpusersApr1716:37AVN_FILE:2001-06-11_12-rw-rw-r--1wrfhlpusersApr1716:38AVN_FILE:2001-06-12_00在$MOAD_DATAROOT/static/名目下创立以下几个文件：-rw-rw-r--1wrfhlpusers124Apr1719:44-rw-rw-r--1wrfhlpusers59704Apr1719:44-rw-rw-r--1wrfhlpusers2974548Apr1719:44-rw-rw-r--1wrfhlpusers29856Apr1719:44-rw-rw-r--1wrfhlpusers119474Apr1719:44在$MOAD_DATAROOT/名目中创立的文件夹还有：./silog/:运行脚本的记录文件./siprd/:最终结果文件。假设程序运行成功，则在./siprd/名目中可以看到如下的文件：-rw-rw-r--1wrfhelpusers6424828Apr1719:311wrfhelpusers6424828Apr1719:321wrfhelpusers368Apr1719:311wrfhelpusers5983520Apr1719:321wrfhelpusers5983520Apr1719:331wrfhelpusers368Apr1719:33其中，前缀“real_input_em”WRFSIOUTPUT_COORD=‘ZETA“real_input_ehWRFSI结果文件。WRF的初始化程序WRF./test/em_real你可以找到和的连接文件以及文件。设置参数配置文件“”WRF小，时间步长，物理参数以前其他的一些参数。一些参数〔比方时间步长，积分时间和输出时间，一些变化参数等〕的转变，可以不需要重执行“”程序。然而，当转变计算范围，边界条件，以及物理参数时，则必需重执行“”程序。运行WRF模式的数据初始化程序。CPU1需将输入数据拷贝和连接到当前名目下才能正确运行。而SP有限，将数据拷贝到当前名目下是不适宜的。因此，我将程序略微修改了一IN_ROOT”,它可以通过IN_ROOT：wrfinput_d01：WRFWRF模式的主程序可串行运行，也可并行运行。同样，由于运行时读取和存放数据的默认名目为当前名目，运行时需将输入数据拷贝和连接到当前名目下才能正确运行。因此，我将的源程序略微修改了一下，添加了