Fortran自由格式和固定格式

1）Free Format（自由格式）

扩展名.F90

1行：一行为132列。可以有复数条程序语句，语句间用分号“;”分开。语句没有位置规定。一行可以是0~132个字符，空格有意义

语句最长不超过2640个字符

最多可有511个续行

2注释行：起始用注释符“!”号，此行其后的所有字符均作为注释（字符串中的!号除外，如 print *,’help!!!!’）。

3续行：当一个语句非常长以至于132列都书写不下时，允许有39个续行。在语句行最后加上续行符“&”号。如果字符串跨2行以上,则在续行的开始位置也要加&号。注意语句的有效字符是从“&”前和续行符“&”之后的位置算起。（一行程序代码的最后如果是符号&，代表下一行程序会和这一行连接。如果一行程序代码的开头是符号&，代表它会和上一行程序连接。如果把Fortran的关键字分为两行书写，则两行都要加续行符号。）

4空格：在语句名和变量名中间不能有空格。需要空格的地方必须有一个以上空格（GOTO和GO TO，ELSEIF和ELSE IF，END构造名(构造名有：DO, PROGRAM, FUNCTION, MODULE, SUBROUTINE)等有两种写法的除外），如关系运算符 ==，<= 不能写成= =，< =。这与F77有很大不同之处，因为F77的设计中将编译源程序的空格忽略。

（二）Fixed Format（固定格式）

扩展名：.F 或 .FOR

1标号区：第l—5列

可以写l至5位整数。也可以没有标号。标号区中的空格不起作用。如lOO与100 或1 0 0作用相同。标号应是无符号整数(无正负号或小数点)。标号大小顺序没有任何要求。假如第二行的标号为1000，第三行的标号可以是10，也可以是99999。

标号区内不得出现标号以外的内容，但注释行例外。注释行的内容可以写在标号区内，一行中第一列为C或*的，该行即被认为注释行，编译时对该行内容不作翻译，对程序运行

fortran调试经验

FORTRAN调试程序的时候注意的问题 调试程序的时候注意的问题。 程序编好，能够直接运行而且完全正确，基本不可能，这就有调试程序的问题。主要有一下几个方面： 其一，在每个子程序被调用的时候留个心眼，写个输出语句，表示程序已经运行到了这里。这样一个小提示会给调试带来巨大的方便，如果程序运行出错，至少你可以知道它是在运行到哪里出的错，这样，直接去检查那个程序就可以了。不必重头开始检查。 其二，注意对中间计算结果的输出。有时候，而且是很多的时候，程序编译成功，运行也没有问题，就是结果不对，这肯定是计算原理有问题，此时，输入一些重要步骤的中间结果，往往可以检查出问题所在。同时，就算查出了问题所在，也可以不删除这段输出中间计算结果的代码，有可能后面还会有用处，此时，在每行输出语句前加一个感叹号，把代码变成注释的绿体字就可以了。等到再次需要输出的时候，直接删除一个“！”比再写一遍输出代码，当然要简单的多。 其三，对WATCH功能的应用，FORTRAN提供的这个功能很实用，可以查很多问题，尤其是程序中间计算值，这个和上述的中间结果的输出有点相似。但两者的不同是前者可以进行中间结果的输出控制，就是只有符合了某个条件的才能被输出，这样可以便捷程序的调试，同时对中间结果输出后可以直接用STOP停止程序的运行，这样对于大型程序来说，节省了很多后面继续计算的时间——因为前面的结果已经不对了，后面的算也是白算。 其四，对中间结算结果输出形式的控制，一般来说，FORTRAN计算结果可以输出到文件里面和计算界面两个地方。对较大的计算结果，可以输出到文件里面，反之较少的结果可以直接输出到屏幕上，为了增强数据的可读性，最好进行有格式的数据输出，以利于相同性质的数据的比较。输出到屏幕上的结果直接用WRITE(6,*)就可以（无格式），对于输出到文件里面的数据，可以省些事情，直接用WRITE(X,*)就可以，其中X是一个任意的正整数，最好大于10，也不用事先对这个X设备进行说明，程序会将结果输出到一个FORT.X的文件里面，例如10，就是FORT.10，此时，用NOTEPAD或者ULTRA-EDIT都可以把它打开——FORT.10实质上就是一个.DAT的文件，你可以把它重命名。 3.对数据计算时的误差控制。 以前觉得小数点后的误差不是那回事，没有太在意，可经过实战，终于明白了小小的误差完全可以改变整个计算的结果。因此，如果程序能够输入结果而不正确时，除了寻找算法的问题，不要忽略的误差。一般认为，FORTRAN的REAL变量小数点后8位数字误差定义已经足够，而事实上,这个精度可能在一些情况下不满足，这个时候，需要用更精确的变量类型——REAL(8),同理，当要判断两个数是否相等的时候，一定要慎用相等判断（.EQ.）这个比较运算符，因为任何数据，别看着在现实中它们一定相等，在程序中就不一定了。一旦经过了计算，就不可避免的产生了舍入误差，对于整数和有限几位循环的有理数都问题不大，可一旦是一个无理数或者无限循环的小数，只有在判断了小数点后的每一位都相等的时候，程序才判断为相等成立。这个相等的标准是非常苛刻的，所以一般情况下，可行的方法是将

ABAQUS中Fortran子程序调用方法

第一种方法： / o/ J5 @6 U/ ^- o$ 1. 建立工作目录/ ]" 2. 将Abaqus安装目录\6.4-pr11\site下的aba_param_dp.inc或aba_param_sp.inc拷贝到工作目录，并改名为aba_param.inc； # ~/ |0 I0 E6 {, @4 X3 q: W3. 将编译的fortran程序拷贝到工作目录； 4. 将.obj文件拷贝到工作目录； 5. 建立好输入文件.inp； 6. 运行abaqusjob=inp_name user=fortran name即可。 第二种方法： 在Job模块里，创建工作，在EditJob对话框中选择General选项卡，在Usersubroutine file中点击Select 按钮，从弹出对话框中选择你要调用的子程序文件（后缀为.for或.f）。 , D8 i7 d/r c6 @" | 以下是网上摘录的资料，供参考：. |$ t/ }$W7 Y6 m4 h6 D6 j 用户进行二次开发时，要在命令行窗口执行下面的命令： 4 O. R+ ^,@( ? abaqus job=job_name user=sub_name ABAQUS会把用户的源程序编译成obj文件，然后临时生成一个静态库standardU.lib和动态库standardU.dll，还有其它一些临时文件，而它的主程序（如standard.exe和explicit.exe等）则没有任何改变，由此看来ABAQUS是通过加载上述2个库文件来实现对用户程序的连接，而一旦运行结束则删除所有的临时文件。这种运行机制与ANSYS、LS-DYNA、marc等都不同。 : j6 g' R-o( {0 [* N2 J3 X这些生成的临时文件要到文件夹C:\Documentsand Settings\Administrator\Local Settings\Temp\中才能找到，这也是6楼所说的藏了一些工作吧，大家不妨试一下。 1子程序格式(程序后缀是.f; .f90; .for;.obj??) 答：我试过，.for格是应该是不可以的，至少6.2和6.3版本应该是不行，其他的没用过，没有发言权。在Abaqus中，运行abaqusj=jobname user=username时，默认的用户子程序后缀名是.for（.f,.f90应该都不行的，手册上也有讲过），只有在username.for文件没有找到的情况下，才会去搜索username.obj，如果两者都没有，就会报错误信息。 如果username包括扩展名for或obj，那么就根据各自的扩展名ABAQUS会自动选择进行操作。 2CAE中如何调用？Command下如何调用？ 答：CAE中在creat job的jobmanager中的general中可以指定子程序； Command下用命令:abaqus j=jobnameuser=userfilename （无后缀）； 3若有多个子程序同时存在，如何处理 答：将其写在一个文件中即可，然后用一个总的子程序调用（具体参见手册） 4我对VF不是很熟，是否可以用VC，C＋＋编写子程序？ A: 若要在vf中调试，那么应该根据需要把SITE文件夹中的ABA_PARAM_DP.INC（双精度）或ABA_PARAM_SP.INC（单精度）拷到相应的位置，并改名为ABA_PARAM.INC即可。 据说6.4的将可以，6.3的你可以尝试着将VC，C＋＋程序编译为obj文件，没试过。在你的工作目录下应该已经存在ufield.obj和uvarm.obj这两个文件（这两个文件应该是你分别单独调试ufield.FOR和uvarm.FOR时自动编译生成的，你可以将他们删掉试试看），但是由于你的FOR文件中已经有了UV ARM 和UFIELD这两个subroutine，显然会造成重复定义，请查实。 用户子程序的使用 假设你的输入文件为：a.inp b.for 那么在ABAQUS Command 中的命令应该是这样的： abaqusjob=a user=b

fortran常见问题解决

楼主为了减少重复回答问题，特编此帖，并不定期添加和更新内容。 错误难免，欢迎讨论，仅供参考。 很多人问哪里可以找到Fortran编译器，有不少热心学友提供网址，特汇集在这里。虽然俺检验过这些链接，但是它们不一定总有效。 Fortran编译器下载： CVF？ FTN95(License:Freeforpersonaluse) 以下操作，如无特别说明，都是以为例。 1.如何加大Stacksize？ 选Project=>Settings=>Link=>Category:Output=>？ Stackallocations Reserve:这里填新值（默认为1M，若需要10M，则填） 2.如何用Fortran批量生成文件？ 设要生成4000个文件，文件名为AA1-AA4000，如何写循环生成文件,而不用写4000次write 命令呢？ 用内部文件： character(len=80)::filename,form integer::i doi=1,4000 selectcase(i) case(1:9) write(form,'(i1)')i case(10:99) write(form,'(i2)')i case(100:999) write(form,'(i3)')i case(1000:9999) write(form,'(i4)')i endselect write(filename,*)"AA",trim(form),".TXT" open(10,file=filename) write(10,*)i close(10)

enddo？ stop end 3.如何用Fortran动态生成输出格式？ 设有一个数组data(100)，输出时，希望每行输出num个数，而num由用户输入，如何实现？ 用内部文件： character(len=80)::form real::data(100) integer::i,num data=(/(i,i=1,100)/)/ read(*,*)num write(form,*)"(",num,"" write(*,form)data stop end 4.MS是不是很垃圾？ 是垃圾，其中Bug太多，多到不可用的地步！ 在这个主题里，换了CVF后问题就没了的人已有相当的数目。 如果你用，遇到莫名其妙的错误，建议换，这是一个比较成熟的编译器。 5.如何用F90/95生成随机数？ 注意： 现在计算机产生的随机数都是伪随机数。 random_number(x)产生一个0到1之间的随机数（x可以是向量），但是每次总是那几个数。用了random_seed()后，系统根据日期和时间随机地提供种子，使得随机数更随机了。 programrandom implicitnone real::x callrandom_seed()!系统根据日期和时间随机地提供种子 callrandom_number(x)!每次的随机数就都不一样了 write(*,*)x stop endprogramrandom 6.函数/子程序超载的例子

FORTRAN 90 程序编程规范

FORTRAN 90 程序编程规范 Fortran 90 编程规范，使程序代码高度组织化，更加易读、易懂、易于维护，程序更加高效。使编出的程序更易懂、易于维护。 1 语言选择 数值预报创新系统软件开发应避免使用Fortran77 的某些过时特征以Fortran 90不一致的特征。选择Fortran 90 作为开发语言，并采用Fortran 90 的新功能，如动态内存的分配(dynamic memory allocation)、递归（recursion ）, 模块(modules)、POINTER 、长变量名、自由格式等。 Fortran 77其中某些只是一些冗余的功能，这些功能已经过时，另外，还有一些在Fortran90 中被证明是不好的用法，建议不要使用。 2 Fortran 90 的新特性 2.1.1 建议使用的Fortran 90 新特性 建议使用Fortran 90 提供的模块（module ），并用Use ONLY 指定module 中哪些变量或派生类型定义可用于调用程序。 尽量使用数组下标三元组，这样可优化并减少所需的代码行数。为提高可读性，要在括号内表明数组的维数，例如： 1dArrayA(:) = 1dArrayB(:) + 1dArrayC(:) 2dArray(: , :) = scalar * Another2dArray(: , :) 当访问数组的子集时，例如在有限差分等式中，可以通过使用下标三元组实现。例如：2dArray(: , 2:len2) = scalar *（ & Another2dArray(:, 1:len2 -1) & - Another2dArray(:, 2:len2) & ) 对程序单元（program units ）命名，并使用End program ，End subroutine ，End interface ，End module 等结构再次指定“program unit ”的名称。 在逻辑表达式中使用>、 >=、 ==、 <、 <=、 /=，它们分别代 替.gt.、.ge.、.eq.、.lt.、.le.、.ne. 。新的表示方法更接近标准的数学符号 在变量定义中始终使用“：：”；始终用“DIMENSION ”定义数组形状；始终用（len=）的语法格式声明字符变量的长度。

fortran课后习题答案

第一章 FORTRAN程序设计基础第15页 1、2 1．简述程序设计的步骤。 “程序设计”：反映了利用计算机解决问题的全过程，通常要经过以下四个基本步骤：(1)分析问题，确定数学模型或方法；(2)设计算法，画出流程图；(3)选择编程工具，编写程序；(4)调试程序，分析输出结果。 2. 什么是算法？它有何特征？如何描述算法？ 解决问题的方法和步骤称为算法。 算法的五个特征：(1) 有穷性。 (2) 确定性。 (3) 有效性。 (4) 要有数据输入。(5) 要有结果输出。 算法的描述有许多方法，常用的有：自然语言、一般流程图、N－S图等。 第二章顺序结构程序设计 第29页 1、2、3、4、5、6、7、8、9 1．简述符号常量与变量的区别？ 符号常量在程序运行过程中其值不能改变。变量在程序运行过程中其值可以改变。 2. 下列符号中为合法的FORTRAN 90标识符的有哪些? (1) A123B (2) M%10 (3) X_C2 (4) 5YZ (5) X+Y (6) F(X) (7) COS(X) (8) A.2 (9) ‘A’ONE (10) U.S.S.R. (11) min*2 (12) PRINT 3. 下列数据中哪一些是合法的FORTRAN常量？ (1) 9,87 (2) .0 (3) 25.82(4) -356231 (5) 3.57*E2 (6) 3.57E2.1 (7) 3.57E+2(8) 3,57E-2 4. 已知A=2,B=3,C=5(REAL)；且I=2,J=3(INTEGER),求下列表达式的值: (1) A*B+C 表达式的值: 11 (2) A*(B+C) 表达式的值: 16 (3) B/C*A 表达式的值: 1.2 (4) B/(C*A) 表达式的值: 0.3 (5) A/I/J 表达式的值: 0.33 (6) I/J/A 表达式的值: 0 (7) A*B**I/A**J*2 表达式的值: 4.5 (8) C+(B/A)**3/B*2. 表达式的值: 7.25 (9) A**B**I 表达式的值: 512 5. 将下列数学表达式写成相应的FORTRAN表达式： (1) 1E-2 （2）(-B+SQRT(B*B-4*A*C)/(2*A) (3) 1+X+X*X/2+X**3/2/3 (4) COS(ATAN((A**3+B**3)**(1.0/3)/(C*C+1))) (5) EXP(A*X**2+B*X+C) (6) COS(X*Y/SQRT(X*X+Y*Y))**3 6. 用FORTRAN语句完成下列操作： (1) 将变量I的值增加1。I=I+1 (2) I的立方加上J,并将结果保存到I中。 I=I**3+J (3) 将E和F中大者存储到G中。G=Max(E,F) (4) 将两位自然数N的个位与十位互换，得到一个新的数存储到M中(不考虑个位为0的情况) M=MOD(N,10)*10+N/10 第三章选择结构程序设计第43页 1、2、3、5、6、7、9 1．分析下列程序运行结果 (1) LOGICAL P INTEGER I,I1,I2,I3 P=.FALSE. READ*,I I1=MOD(I,10) I2=MOD(I/10,10) I3=I/100

FORTRAN程序设计复习题及答案

FORTRAN程序设计复习题 一、选择题 B (1)下列各FORTRAN表达式中合法的是 A) S+T*2P >= B) .NOT. (A*B+C) C) A2+B2/(C+D) <= D) (A+B).NOT．A*B．GT．(.NOT.只跟一个表达式) C (2)数学式(3/5)ex+y的FORTRAN表达式是 A) 3*EXP(X+Y)/5 B) 3*E* *(X+Y)/ C) (3/5)*EXP(X+Y)D) EXP(X+Y) D (3)下列FORTRAN77表达式中不合法的是 A) A．GT．B．EQV．C．GT．D B) A．AND．B．AND．C．AND．D C) .NOT．(X．LE．D) A．LT．B．LT．C．LT．D D(4)下列叙述中不正确的是 A) FORTRAN子程序可以单独编译 B) 对一个FORTRAN源程序进行编译和连接无误后可生成可执行文件 C) 即使编译和连接都正确无误，FORTRAN程序运行时仍可能出错 D) FORTRAN连接的主要任务是把函数库中的函数翻译成机器指令(正确描述：主要任务为连接目标文件) B (5)在下列FORTRAN77运算符中，优先级最高的是 A) .AND. B) .NOT. C) .OR. D) .EQ. B (6)FORTRAN表达式"6/5+9/2**3/2"的值为 A) 33 B) 1 C) 5 D) 3 A (7)下列FORTRAN77表达式中，合法的是： A) .AND.. B) 10.0 C) D) 提示：A）相当于 .AND.（.NOT.（）） D (8)关于编译一个FORTRAN源程序文件，下列说法中错误的是 A) 允许编译只有一个主程序而没有子程序的源文件 B) 允许编译有多个子程序的源文件 C) 允许编译只有一个子程序而没有主程序的源文件 D) 允许编译有多个主程序的源文件 C (9)在FORTRAN77源程序中，续行标志符必须放在 A) 第1列 B) 第1-6列C) 第6列D) 第5列 D (10)下列关于"SUBROUTIN E MAP(X，Y)"语句行的叙述中，不正确的是 A) 这是子程序的第一个语句 B) 字符串"MAP"是子程序名 C) 变量X是子程序的形参D) 子程序执行后，MAP将返回整型数据 提示：子程序无返回值，自定义函数才有） A (11)FORTRAN表达式"2/4+"的值是 A) B) 1 C) D) 0 提示：2/4默认等于整型，=》 D (12)FORTRAN表达式"MOD，"的值是 A) B)0.0 C) D) A (13下列FORTRAN运算符中，优先级最低的是 A）逻辑运算符.AND. B）算术运算符*

中国海洋大学Fortran程序考试典型题

1、找出100-150之间和400-450之间能被9整除的数 2、找出100-999之间的回文数 3、找出水仙花数

4、输出“*”图形 5、输入n个数，找出大于平均值的数和最小数

6、把一个数值型数组的相同数删除到只剩一个 7、形成一个5*5矩阵，对角线元素为“i”，其余为“j”program juzhen implicit none integer I,J integer,parameter::size=5 integer::a(size,size) forall(I=1:size,J=1:size,I>J) a(I,J)=j forall(I=1:size,J=1:size,I==J) a(I,J)=i forall(I=1:size,J=1:size,I

integer,allocatable::c(:) write(*,*) "输入数组A的数据个数" read(*,*) n write(*,*) "输入数组B的数据个数" read(*,*) m l=m+n allocate (a(n)) allocate (b(m)) allocate (c(l)) write(*,*) "从小到大输入A的元素" do i=1,n read(*,*) a(i) end do write(*,*) "从小到大输入B的元素" do i=1,m read(*,*) b(i) end do do i=1,n c(i)=a(i) end do i=1 do while(i=c(n))then n=n+1 c(n)=b(i) goto 10 else if(b(i)>c(j).and.b(i)

fortran常见错误

FAQ之常见错误 2014-02-02 13:45:35 来源：Fcode研讨团队评论：2点击：4419 本文从编译错误，链接错误，运行时错误，计算结果错误等四个方面介绍了常见的错误及解决思路。适合初学者阅读。 首先应该明确：错误有哪几种？我们当前遇到的是何种错误？ 阐述这些问题前，我们先讨论一下常规的应用程序开发的过程： 1>>编写代码，使用一个或多个源代码文件。 2>>对第一步的每一个源代码文件执行编译操作。得到一个或若干个目标代码。 3>>将目标代码，运行时库（Run-time Library）和其他使用到的函数库链接起来。得到一个可执行文件（EXE 或其他） 4>>编写程序的说明书，必要的（输入）数据文件 5>>将上述得到的结果发布给用户。（发布的方式可以是刻录成光盘，销售，放在网站上供别人下载，或者其他） 6>>用户得到程序后，运行，输入数据，得到计算结果。 对于很多 Fortran 程序员来说，可能用户就是自己，也可能仅仅是自己教研室的同事同学。所以第4，5，6步骤很多时候不明显。而如果使用集成开发环境（IDE）进行开发，第1，2，3步骤又可以一键完成。因此，很多初学者就认为，写程序就是：输入代码，运行，得到结果。这样的理解太狭义。 不管我们面对什么使用者来写代码，程序开发应该是上述的过程。我们的编译器，编译环境，也是为这个过程而设计的。 于是，我们将错误分为四种： 一. 编译错误（发生在第2步） 编译错误，一般是源代码书写格式不正确，不符合语法要求。 二. 链接错误（发生在第3步） 链接错误，一般是源代码结构不完整，运行时库或函数库使用不合理。 三. 运行时错误（发生在第6步） 运行时错误，一般是执行代码时，遇到了事先未料及的错误。比如内存不足了，磁盘空间不够了，输入文件格式不对了，输出文件写入失败了等等。 四. 计算结果不符合预期（程序代码不规范，或不符合你的设想） 计算结果不符合预期，可能性就很多了。语法与你的想法不一致，超出函数库的适用范围，执行流程控制不当等等。 这四种错误，其排查难度依次增大。也就是，编译错误最容易排查和修改，而计算结果不正确，最让人头疼。

fortran程序实例

1)实例3—求多个半径下的圆周长 ! z3.f90 --Fortran95 ! FUNCTIONS: ! z3 - Entry point of console application. !************************************************************************* ! PROGRAM: z3 ! PURPOSE: Entry point for the console application. !************************************************************************ program z3 ! 求多个半径下的圆周长 ! 主程序 ! PROGRAM Z3 PRINT *, 'R=',1.2,'C=',C(1.2) PRINT *, 'R=',3.4,'C=',C(3.4) PRINT *, 'R=',15.6,'C=',C(15.6) PRINT *, 'R=',567.3,'C=',C(567.3) END program z3 !子程序 FUNCTION C(R) PI=3.1415926 C=2*PI*R RETURN ! Body of z3 end

2)实例4—键盘与显示器输入/输出 a)Fortran 基本操作 b)程序指令 ! ZXZ_I_O.f90 ! FUNCTIONS: ! ZXZ_I_O - Entry point of console application. ! PROGRAM: ZXZ_I_O ! PURPOSE: Entry point for the console application. !***************输入、输出样式种种************************** program ZXZ_I_O implicit none !变量声明的位置 INTEGER(2) i; INTEGER(4) j; INTEGER(4) m; REAL n INTEGER A,B ! Variables PRINT*,'输入整数A'; READ*, A PRINT*,'输入整数B'; READ*, B B=A+B PRINT*,'B=A+B=',B WRITE(*,*) 'A*B=',A*B PRINT* ,'以上为计算机的计算结果，注意B的值'

fortran中批处理实现

********************************************* fortran中批处理命令的实现函数： 利用systemqq命令（需要调用DFLIB 数据库） ********************************************* 例1： USE DFLIB character*100 CMD LOGICAL(4) res CMD="dir/a-d/b/s "//trim(fPath)//" >"//trim(outPut) res=SYSTEMQQ(CMD) 例2： USE DFLIB LOGICAL(4) result result = SYSTEMQQ('copy e:\dir.txt e:\test\dir.txt') !将e:\dir.txt 复制到e:\test\dir.txt文件中。！****************实例3:复制文件************************* 例3： programmain_pro USE DFLIB implicit none integer,parameter::sta_num=123 character(5),dimension(sta_num)::sta_ID character(500)::filein,fileout character(5000)::cmd logical(4)::judge

integer::status,is open(1,file='山东.txt',status='old',action='read',iostat=status) read(1,*) do is=1,sta_num read(1,*) sta_ID(is) filein='Z:\data\降水逐小时数据-戴至修\precip_data\'//sta_ID(is)//'_precip.txt' open(2,file= filein,status='old',action='read',iostat=status) if(status/=0) goto 1000 fileout='Z:\data\降水逐小时数据-戴至修\山东省-降水数据\'//sta_ID(is)//'_precip.txt' cmd='copy '//filein//' '//fileout judge=SYSTEMQQ( cmd) 1000 continue enddo end program

fortran问题

imsl7.0中用use linear_operators这句话是会出错的，当时intel论坛上也有人问，后来intel 给出了个X64 imsl的补丁，但是32位的没有。 要使用包含在linear_operators这个库中的函数时，要用use+原函数 例如：上面那个例子把use linear_operators改成use operation_xt即可 另外，imsl7.0引用函数和以前版本不一样的 补充一点： linear_operators这个文件是这样子的（一看就知道原因了）： modulelinear_operators usecond_int usedet_int usediag_int usediagonals_int usefft_int useifft_int useeye_int uselin_eig_self_int uselin_sol_self_int usenorm_int useoperation_i useoperation_ix useoperation_t useoperation_h useoperation_tx useoperation_hx useoperation_x useoperation_xi useoperation_xt useoperation_xh useorth_int userand_int userank_int usesvd_int useunit_int useeig_int usechol_int useisnan_int end module 1. 如何加大Stack size? 选Project => Settings => Link => Category: Output =>

Fortran 运行中给出的系统错误及解决方法

. Fortran 运行中给出的系统错误及解决方法 以下均为linker tools errors and warnings Linker Tools Error LNK1000 unknown error; consult documentation for technical support options Note the circumstances of the error, try to isolate the problem and create a reproducible test case, then contact technical support. Linker Tools Error LNK1101 incorrect MSPDBxx.DLL version; recheck installation of this product The version of MSPDBxx.DLL available on your system does not match the version required by this tool. Linker Tools Error LNK1102 out of memory There was not enough memory for the tool to run. Probably the paging file exceeded available disk space. If a shortage of disk space is not the cause, note the circumstances of the error, try to isolate the problem and create a reproducible test case, then request technical support. Linker Tools Error LNK1103 debugging information corrupt; recompile module Probably the compilation was terminated before a valid object file was created.Recompile the given object file. If recompiling does not correct the problem,note the circumstances of the error, try to isolate the problem and create a reproducible test case, then consult technical support. Linker Tools Error LNK1104 cannot open file "filename" The tool could not open the given file. One of the following may be a cause: l There was not enough disk space. l The file does not exist. l The filename or its path was incorrectly specified. l The specified drive is invalid. l The file does not have the appropriate permissions. l The path for filename expands to more than 260 characters. l If the given file is named LNKn, which is a filename generated by the linker for a temporary file, then the directory specified in the TMP environment variable may not exist, or more than one directory is specified for the TMP environment variable. (Only one directory path should be specified for the TMP environment variable.) l If the error occurs on the executable filename, an earlier version of the executable may still be running. You will need to terminate the executable before linking it. In Windows NT (including Windows 2000) or Windows 95,you can use the utility PVIEW to look for and kill instances of the application. l If the error message occurs for a library name, and you recently ported the .MAK file from a previous Microsoft Visual C++ development

aspenfortran子程序的编译、连接、运行

Aspen plus fortran 子程序的编译、连接与运行 Aspen Plus可以用来模拟反应精馏，但提供的动力学方程形式比较死板，很多类型的动力学方程都不能直接输入，不得不外接一个fortran 子程序。本来用到fortran 的机会就不多，所以会的人也不多。何况编写好fortran 子程序后，还要用fortran 编译器来编译生成obj 文件，而且每次运行的时候都得连接一次，比较麻烦。我的电脑以前是装了编译器这些软件 的，但最近出了问题，只能编译fortran 子程序，运行aspen外接fortran 子程序的时候就会出错。恰好这几天要帮师兄写个fortran 子程序，写完后需要调试，但又运行不了，重装 呢也没把握，想起网上有兄弟说过可以用asplink命令生成一个dll文件，这样就不用运行 aspen的时候每次都连接编译器了。于是找来资料学习学习，居然成功了，欣喜之余，分享给大家。第一步当然是编写动力学子程序了（此处省略若干字） 第二步编译写好的动力学子程序（以网友soaty的动力学子程序为例） 打开Aspen Plus Simulation Engine JS M... 打开后是这样的 接下来改变文件夹路径，转到存放动力学子程序的那个文件夹 比如我的文件路径是C:\ProgramData\Aspe nTech\Aspen Plus \IPP 于是输入cd C:\ProgramData\AspenTech\Aspen Plus \IPP

回车后就转到这个文件夹了，这时候就可以进行编译了：编译命令：aspcomp 输入命令aspcomp ,回车 到文件夹一看你会发现已经生成了obj文件了,

几个fortran程序

计算圆周率 REAL R,R1,R2,PI ISEED=RTC() N0=0 N=300000 DO I=1,N R1=RAN(ISEED) R2=RAN(ISEED) R=SQRT(R1*R1+R2*R2) IF(R<1.0)N0=N0+1 END DO PI=4.0*N0/N WRITE(*,*)PI END 一）蒙特卡洛计算生日问题 假设有N个人在一起，各自的生日为365天之一，根据概率理论，与很多人的直觉相反，只需23个人便有大于50％的几率人群中至少有2个人生日相同。 INTEGER M(1:10000), NUMBER1(0:364), NUMBER2 REAL X,Y ISEED=RTC() DO J=1, 10000 NUMBER1=0 X=RAN(ISEED) NUMBER1(0)=INT(365*X+1) JJJ=1 DO I=1,364 Y=RAN(ISEED) NUMBER2=INT(365*Y+1) ETR=COUNT(NUMBER1.EQ.NUMBER2) IF (ETR= =1) THEN EXIT ELSE JJJ=JJJ+1 M(J)=JJJ NUMBER1(I)=NUMBER2 END IF END DO END DO DO I=1,10000 IF(M(I).LE.23) SUM=SUM+1 END DO PRINT *,SUM/10000 END 二）MONTE CARLO SIMULATION OF ONE DIMENSIONAL DIFFUSION 蒙特卡罗计算一维扩散问题

INTEGER X,XX(1:1000,1:1000) REAL XXM(1:1000) ! X:INSTANTANEOUS POSITION OF ATOM ! XX(J,I)：X*X ,J:第几天实验,I:第几步跳跃 ! XXM(I): THE MEAN OF XX WRITE(*,*) "实验天数JMAX，实验次数IMAX" READ(*,*) JMAX,IMAX ISEED=RTC() DO J=1,JMAX !第几天实验 X=0 !!! DO I=1,IMAX !第几步跳跃 RN=RAN(ISEED) IF(RN<0.5)THEN X=X+1 ELSE X=X-1 END IF XX(J,I)=X*X END DO END DO OPEN(1,FILE="C:\DIF1.DAT") DO I=1,IMAX XXM=0.0 XXM(I)=1.0*SUM(XX(1:JMAX,I))/JMAX !! WRITE(1,*) I, XXM(I) END DO CLOSE(1) END 三维的！ 三）通过该程序了解FORTRAN语言如何画图（通过像素画图）USE MSFLIB INTEGER XR，YR !在的区域中画一个圆 PARAMETER XR=400,YR=400 INTEGER R, S(1:XR,1:YR) X0=XR/2 ! 圆心位置X0，YO Y0=YR/2 R=MIN(X0-10,Y0-10) !圆半径 S=0 !像素的初始状态（颜色） DO I=1,XR DO J=1,YR IF((I-X0)**2+(J-Y0)**2<=R**2)S(I,J)=10 IER=SETCOLOR(S(I,J)) IER=SETPIXEL(I,J) END DO END DO END

fortran程序30个

程序1 PROGRAM chaper1272 INTEGER :: A(10)=(/5,7,4,8,12,2,10,3,9,11/) INTEGER :: sum=0,ave !打开一数据文件,设置一个有格式直接存取文件,将10个数分2个记录写入文件。 !数据文件生成2个记录,每个记录行长度相同,记录长度为25 OPEN(1,FILE='input21.dat',FORM='FORMATTED',ACCESS='D IRECT',RECL=25) WRITE(1,"(5I5)",REC=1)(A(I)+10,I=1,5) !按格式说明将头5个数写入第1个记录 WRITE(1,"(5I5)",REC=2)(A(I)+10,I=6,10)!按格式说明将后5个数写入第2个记录 READ(1,"(5I5)",REC=2)(A(I),I=6,10) !按格式说明从第2个记录中读取后5个数 READ(1,"(5I5)",REC=1)(A(I),I=1,5) !按格式说明从第1个记录中读取头5个数 DO I=1,10 sum=sum+A(I) ENDDO ave=sum/10 !打开一个最大记录长度为22的有格式顺序存取文件

OPEN(2,FILE='input22.dat',FORM='FORMATTED',ACCESS='D IRECT',RECL=22) WRITE(2,"('10个数之和为：',I5)",REC=1) sum !输出1记录行,记录长度为22 WRITE(2,"('10个数平均值为：',I5)",REC=2)ave !输出1记录行,记录长度为22 WRITE(2,"(A)",REC=3) '程序运行正常结束。' !输出一个记录行,记录长度为22 END 程序2 PROGRAM average REAL sum,ave INTEGER n OPEN(1,file='score.dat') PRINT*,'正在统计平均成绩,请等待。' sum=0.0;n=0 DO READ(1,*,END=100) S sum=sum+s n=n+1 ENDDO 100 ave=sum/n