当前位置:文档之家› fortran常见问题解决

fortran常见问题解决

楼主为了减少重复回答问题,特编此帖,并不定期添加和更新内容。

错误难免,欢迎讨论,仅供参考。

很多人问哪里可以找到Fortran编译器,有不少热心学友提供网址,特汇集在这里。

虽然俺检验过这些链接,但是它们不一定总有效。

Fortran编译器下载:CVF?FTN95(License:Freeforpersonaluse)以下操作,如无特别说明,都是以为例。

1.如何加大Stacksize?选Project=>Settings=>Link=>Category:Output=>?StackallocationsReserve:这里填新值(默认为1M,若需要10M,则填)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),".TXT"open(10,file=filename)write(10,*)iclose(10)enddo?stopend3.如何用Fortran动态生成输出格式?设有一个数组data(100),输出时,希望每行输出num个数,而num由用户输入,如何实现?用内部文件:character(len=80)::formreal::data(100)integer::i,numdata=(/(i,i=1,100)/)/read(*,*)numwrite(form,*)"(",num,""write(*,form)datastopend4.MS是不是很垃圾?是垃圾,其中Bug太多,多到不可用的地步!在这个主题里,换了CVF后问题就没了的人已有相当的数目。

如果你用,遇到莫名其妙的错误,建议换,这是一个比较成熟的编译器。

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中为多行代码加注释?见:8.推荐好的代码风格根据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就有关问题进行磋商,目前的进展比较顺利,可能要在以后的版本中支持,目前还没有时间表。

11.如何在CVF中检验代码是否符合F90/95标准?为了检验你的代码是否符合F90或F95标准,在CVF中选:Project/Settings...选Fortran页卡Category下拉框选:CompilationDiagnosticsFortranStandardChecking下拉框选:Fortran90或Fortran95这样,凡是CVF自家扩展的语法就都会有提示了。

例:typetestinteger,allocatable::a(:)integer,pointer::b(:)endtypetest若选了标准检验,会提示:Warning:Allocatablefieldsofderivedtypesarenon-standard[A]也就是说按照F90/95标准,type中不允许用allocatable。

相关主题