fortran常见问题解决

楼主为了减少重复回答问题，特编此帖，并不定期添加和更新内容。错误难免，欢迎讨论，仅供参考。很多人问哪里可以找到Fortran编译器，有不少热心学友提供网址，特汇集在这里。虽然俺检验过这些链接，但是它们不一定总有效。Fortran编译器下载：CVF6.6/SoftDown.asp?ID=11937FTN95(License:Freeforpersonaluse)/Silverfrost-FTN95/3000-2069_4-10491440.html?tag=lst-0-10以下操作，如无特别说明，都是以CVF6.6为例。1.如何加大Stacksize?选Project=Settings=Link=Category:Output=StackallocationsReserve:这里填新值（默认为1M，若需要10M，则填10000000）2.如何用Fortran批量生成文件？设要生成4000个文件，文件名为AA1-AA4000，如何写循环生成文件,而不用写4000次write命令呢?用内部文件：character(len=80):filename,forminteger:idoi=1,4000selectcase(i)case(1:9)write(form,(i1)icase(10:99)write(form,(i2)icase(100:999)write(form,(i3)icase(1000:9999)write(form,(i4)iendselectwrite(filename,*)AA,trim(form),.TXTopen(10,file=filename)write(10,*)iclose(10)enddostopend3.如何用Fortran动态生成输出格式？设有一个数组data(100)，输出时，希望每行输出num个数，而num由用户输入，如何实现？用内部文件：character(len=80):formreal:data(100)integer:i,numdata=(/(i,i=1,100)/)/10.0read(*,*)numwrite(form,*)(,num,f10.3)write(*,form)datastopend4.MSPOWERSTATION4.0是不是很垃圾？POWERSTATION4.0是垃圾，其中Bug太多，多到不可用的地步！在这个主题里，换了CVF后问题就没了的人已有相当的数目。如果你用POWERSTATION4.0，遇到莫名其妙的错误，建议换CVF6.6，这是一个比较成熟的编译器。5.如何用F90/95生成随机数？注意：现在计算机产生的随机数都是伪随机数。random_number(x)产生一个0到1之间的随机数（x可以是向量），但是每次总是那几个数。用了random_seed()后，系统根据日期和时间随机地提供种子，使得随机数更随机了。programrandomimplicitnonereal:xcallrandom_seed()!系统根据日期和时间随机地提供种子callrandom_number(x)!每次的随机数就都不一样了write(*,*)xstopendprogramrandom6.函数/子程序超载的例子设要编一个两个变量值互换的子程序swap(a,b)，哑元a,b可能是实型数，整型数，数组，矩阵，字符串，派生类型等等。但是希望只用一个子程序接口swap(a,b)来实现。F90可以用类属接口来实现这种子程序超载：moduleUtilitiesimplicitnoneprivateI_Swap,R_Swap,RVec_Swap,RMat_Swap,Type_Swappublic:SwapinterfaceSwapmoduleprocedureI_Swap,R_Swap,RVec_Swap,RMat_Swap,Type_Swapendinterfacecontainssubroutinei_swap(a,b)!整型置换integer(ikind),intent(inout):a,binteger(ikind):t。!略endsubroutinei_swapsubroutiner_swap(a,b)!实型置换real(rkind),intent(inout):a,breal(rkind):tt=aa=bb=treturnendsubroutiner_swapsubroutineRVec_swap(a,b)!实型向量置换real(rkind),intent(inout):a(:),b(:)integer(ikind):idoi=1,size(a)callR_Swap(a(i),b(i)enddoreturnendsubroutineRVec_swapsubroutineRMat_swap(a,b)!实型矩阵置换。!略endsubroutineRMat_swapsubroutineType_swap(a,b)!派生类型置换。!略endsubroutineType_swapendmoduleUtilities7.如何在CVF中为多行代码加注释?见：/club/showbbs.asp?id=1243188.推荐好的代码风格根据F90子集语言ELF90和F的要求整理（部分）。“强迫用”的语言特性：+F90的自由格式的源代码。+implicitnone。+子过程的哑元都要有intent属性。+函数子程序的哑元必须指定为intent(in)。+所有子程序和函数都放在模块（module）中，然后引用（use）该模块；或者放在program中。+数组哑元要求是假定形状的，或者有固定的维数和大小。字符哑元要求是假定长度的。+对于recursivefunction（递归函数）语句，必须有result子句。+在所有派生类型（type）的定义语句中，必须用双冒号分隔符（:）。+主程序要求有program语句。+在程序单元的end语句中要求后跟程序单元的类型和名称。+在endtype语句中要求后跟类型的名称。+endprogram前必须有stop语句以表示停止执行。+子过程中必须有return语句，以表示返回。+subroutines()并且calls()，即必须有括号。“不得用”的语言特性：-allocatable、intent、pointer、save、dimension、parameter和target语句形式。(用属性形式代替。)-external语句形式。(用显式的接口代替。)-assign、赋值goto、交错return、continue、entry、和计算goto语句。-include文件。(用模块代替。)-data和blockdata。(在类型声明语句中进行初始化或赋值。)-common块。(将全局数据放在模块中，用模块代替。)-equivalence。(被认为是许多难以查找的编程错误的来源。)-doubleprecision语句。(用real语句声明双精度的实型数。)-语句函数。(用内部函数代替。)-专用固有函数。(用类属函数代替。)-假定大小数组。(用假定形状数组代替。)-don（其中n为语句标号）。(用do和enddo代替。)-非整数do变量和表达式。-同一行上多条语句。-逻辑型case表达式。-从if块外面分支到endif。-where语句形式。(用where结构形式。)-在open和inquire语句中的blank=说明符。-双字关键词之间要求有空格：inout，goto。不能写为inout，goto。9.将字符串改为大写的子程序subroutineUpCase(str)!=!changetouppercase!=character(len=*),intent(inout):strinteger(4):icha,LL,icvalinteger(4),parameter:diff=ichar(a)-ichar(A)LL=len_trim(str)doicha=1,LLicval=ichar(str(icha:icha)if(icval=ichar(a).and.icval=ichar(z)thenstr(icha:icha)=char(icval-diff)endifenddoreturnendsubroutineUpCase10.CVF中源代码信息浏览默认情况下，代码信息浏览是失效的（GotoDefinition/Reference都不起作用），你可以将其启用：1.Project/Settings/Fortran,选中GenerateSourceBrowseInformation。2.选BrowseInfo卡页，选中BuildBrowseinfofile，点击OK。3.Build程序。代码浏览器仅当成功Build程序后才可应用。4.Tools/SourceBrowser。5.在Browse对话框底部，找到Casesensitive复选框。Fortran是大小写不敏感的语言，去除Casesensitive的选中。点击OK。6.现在可以用Browse对话框来浏览函数调用关系以及变量声明和引用关系。7.光标置于一个变量名或函数名中，鼠标右键：GotoDefinition使光标跳到该变量声明或函数定义语句；GotoReference使光标跳到该变量的引用或函数引用语句。注：(Addison提供)IVF中现在不适用，因为此技术是属于MS的，不是INTEL的，所以，intel正在与MS就有关问题进行磋商，目前的进展比较顺利，可能要在9.1以后的版本中支持，目前还没有时间表。11.如何在CVF中检验代码是否符合F90/95标准？为了检验你的代码是否符合F90或F95标准，在CVF中选：Project/Settings.选Fortran页卡Category下拉框选：CompilationDiagnosticsFortranStandardChecking下拉框选：Fortran90或Fortran95这样，凡是CVF自家扩展的语法就都会有提示了。例：typetestinteger,allocatable:a(:)integer,pointer:b(:)endtypetest若选了标准检验，会提示：Warning:Allocatablefieldsofderivedtypesarenon-standardA也就是说按照F90/95标准，type中不允许用allocatable。臭石头雪球提供：VS2005+IVF9.1组合的方法差不多：项目属性，Fortran，Diagnostics，WarnForNon-standardFortran选择为Fortran90或Fortran9512.如何输出tab字符？char(9)。13.如何在CVF中设置把Tab转为空格Tab字符不在Fortran标准中，在不同环境下会有不同的显示，不推荐使用。采用下列方法，可将Tab键转为3个空格：依次选菜单Tools/Option，选Tab卡项，Filetype选Fortran，选Insertspaces，Tabsize选3，Indentsize选3。14.如何实现“代码自动补全”在CVF中：/club/post-238710.html在IVF中：/club/post-241489.html问题各位大侠，非常感谢你们的回复。置顶帖中关于批处理读取文件的帖子我已经看了。帖子里是批处理AA1-AA4000.txt文件，文件名是1-4000的数字序列。而我要读取的是时间序列的文件，即1月1日-12月31日，year+month+day的形式。因此置顶帖中的例子不适用我这个问题。文件名为（2006年365天的数据，共365个文件）pre.20060101.datpre.20060102.datpre.20060103.dat.pre.20060131.datpre.20060201.datpre.20060202.dat.pre.20061231.dat单个文件读取已经完成，但是如何批处理还没有解决。解决1：既然是有规律你就可以做循环来进行读取,每次循环读取一个文件.例如做3层循环,分别是年循环(如果就一年那就不用循环了),月循环,日循环.不难吧.解决2我的方法：先把文件名：pre.20060101.dat,pre.20060102.dat,.输入到一个文件中，然后依次读取，赋给字符型变量数组filenames(i),然后就可以循环打开所有文件。程序如下：programtestinteger:a(12)=31,28,31,30,31,30,31,31,30,31,30,31