Linux 下面 fortran 的使用方法

Linux 下面 fortran 的使用方法一、安装gfortran

输入命令：sudo apt-get install gfortran 和密码

输入：y

二、在linux下面运行fortran

新建工作目录文件夹，如“hello”

在文件夹中创建“.f”文件，如“hello.f”

在“hello.f”里面写入：

program main

write(*,*)”hello”

stop

end

保存文件

每一行前面要键入

Tab键

在当前目录下打开终端，输入命令：gfortran hello.f

稍后会产生“a.out”的可执行文件

输入命令：./a.out 执行该文件

执行结果

CALYPSO教程

U s e r' s G u i d e CALYPSO version 1.3.0 October 28, 2011 Written by Yanming Ma, Yanchao Wang, Jian Lv and Li zhu State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, China Please contact with Yanchao Wang (wyc09@https://www.sodocs.net/doc/c112378713.html,), Jian Lv (jianlv10@https://www.sodocs.net/doc/c112378713.html,), Li Zhu (zhuli10@https://www.sodocs.net/doc/c112378713.html,) and Yanming Ma (mym@https://www.sodocs.net/doc/c112378713.html,) for any technical questions or any bugs. https://www.sodocs.net/doc/c112378713.html,/YanmingMa/mym.htm Copyright ? 2011 Jilin University. All Rights Reserved.

CONTENTS 1. INTRODUCTION (3) 2. COMPILATION (4) 3. EXECUTION OF THE PROGRAM (4) 4. INPUT FILE (5) 5. OUTPUT FILES (11) 6. THE RESULTS ANALYSIS (12) 7. SOME PUBLICATIONS FROM USE OF CALYPSO (13) 8. FUTURE DEVELOPMENT (13) 9. ACKNOWLEDGMENTS (13)

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程序实例

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的值'

fortran95教程4章

第2篇. 计算的叙述 算法的每一个步骤，都必须给予确切的定义。对于算法当中所 考虑的每一种情况，每一个有待执行的动作，都必须严格地和 不含混地加以规定。…对于以描述算法作为目的而设计出来 的，采用了形式的定义的程序设计语言，或者说计算机语言， 它的每一个语句都必须有非常确切的意义。 ---- D.E.Knuth[1]《The Art of Computer Programming》 本质上FORTRAN就是一门语言，一门人与计算机赖以进行有效交流的语言，在这个意义上和我们使用的中文，英文等没有本质差别。现在假设要来描述一种大家都陌生的语言，那么总是要分成两个方面来描述，即一方面要描述这门语言的表象和形态，也就是它使用哪些符号，哪些词汇，一般的句式如何，怎样才能完整叙述一个任务之类；另一方面需要说明这门语言的语义，也就是说这门语言是如何用来表达我们需要它表达的意思的。 第4章基本上就是描述FORTRAN作为一种语言的基本形态，也就是书写这种语言的书写规则。 接下来几章则逐步说明如何用FORTRAN来表达我们的要求，或者反过来说，FORTRAN提供了些什么表达方式，以便我们用来向计算机提出合理的任务： ● 表达基本数据； ● 表达数据的结构； ● 完整地描述数据； ● 构造表达式； ● 驱动计算的赋值； ● 计算过程的结构控制； 在整个第二篇，我们将领略到FORTRAN 95是如何能够做到精致地描述计算的，而把一个问题阐述清楚了，就意味着问题已经解决了一大半。

[1]Donald E. Knuth (高纳德),Stanford University的The Art of Computer Programming荣休教授，而The Art of Computer Programming（计算机程序设计技巧）正是他的伟大著作的名称。洋洋七大卷的《The Art of Computer Programming》是当今全世界每一个计算机科学家所膜拜的圣经。1974年在该书刚完成前面很少一部分时，就给他带来了计算机科学家们梦寐以求的图灵奖。 第4章FORTRAN 95语言的形貌 要说明一门语言的形态，必须回答以下问题： ●它使用哪些符号来表达信息？ ●它的词汇如何构成？ ●它的语句如何构成？ ●如何表达一个完整的任务？ 具体的对于一门计算机语言，把这几个问题更加明确地转换过来，就是： ●它使用键盘上的哪些符号，各个符号有哪些用途？ ●它的词汇如何由键盘字符构成？含有哪些固定的词汇？以及容许自由构成合法词汇 的规则是什么？ ●它具有哪些固定的语句格式？以及容许自由构成合法语句的规则是什么？ ●我们交待给计算机的任何任务，都必须明确说明任务的开始，执行步骤和完成，因 此一段完整的源代码应该具备什么样的形式？以及应该具备哪些要素？ 本章就是要回答这些问题。 4.1 FORTRAN语言所使用的字符 从最抽象的层面来看，人与计算机的交流只是信息的交流，而信息总是需要依靠某种信号来表示，对于人来说，最方便的就是字符。而对于计算机来说，自然就是键盘所能敲出的那些字符(信号)，因此下面就是要说明： ●FORTRAN 95能识别键盘上敲出的哪些字符？ ●每个字符对于FORTRAN 95来说又意味着什么？ 4.1.1 FORTRAN 95所使用的基本字符 按照FORTRAN 95标准的规定，一切FORTRAN 95的实现平台都必须使用下面表4-1所列出来的这个基本的字符集，或者说，这个字符集是所有遵循FORTRAN 95标准的编译

fortran安装

Intel Visual Fortran 2013 sp1 完整版下载及安装教程(windows) 10月16日版：更新支持win8.1、VS2013下载地址 【完整安装包】https://www.sodocs.net/doc/c112378713.html,/akdlm/irc_nas/3651/w_fcompxe_ novsshell_2013_sp1.1.139.exe 【联网安装包】https://www.sodocs.net/doc/c112378713.html,/akdlm/irc_nas/3651/w_fcompxe_ novsshell_online_2013_sp1.1.139.exe 【百度云副本】https://www.sodocs.net/doc/c112378713.html,/s/18E41z 更新内容（翻译自官网，不妥之处恳请指正）： 1.Intel? Fortran Compiler 更新至14.0.1 --增加/assume:std_value 编译选项 --增加/Q[a]xMIC-AVX512 编译选项 --增加/Qopt-gather-scatter-unroll=n 编译选项 2.Intel? Math Kernel Library 更新至11.1 Update 1 3.增加对Microsoft Windows 8.1*的支持 4.增加对日语的支持 5.修正问题 --https://www.sodocs.net/doc/c112378713.html,/en-us/articles/intel-composer-xe-2013-compilers-sp1-fixes-list --https://www.sodocs.net/doc/c112378713.html,/en-us/articles/intel-mkl-111-bug-fixes/ 安装说明： 1、下载安装包中的任意一个，前两个为官方地址，第三个为百度云副本（完整安装包） 2、运行w_fcompxe_novsshell_2013_sp1.0.103.exe 3、点击extract，解压安装包，解压完毕程序会自动运行安装向导 4、在第三步（Activation）中选择Choose alternative activation，点击下一步 5、选择Use a license file 6、选择你的授权许可文件（lib文件） 7、继续安装过程 备注： 1、安装Intel Visual Fortran Composer前，请确保已经安装了Microsoft Visual Studio 2010/2012/2013。若使用其他编辑器，需要手动配置编译环境。 Visual Studio官网：https://www.sodocs.net/doc/c112378713.html,/。 微软旗下的DreamPark计划为在校大学生免费提供VS2013Pro版，地址：https://https://www.sodocs.net/doc/c112378713.html, 2、安装过程中没有加以说明的步骤，请直接点击下一步。

fortran编译器操作

关于fortran语言的编译器使用问题 此文是一个简明教程，仅适合初学者来使用。老鸟们大牛们就不用在看了。本文档的目的是为了方便初学者快速的掌握基本的fortran编译器的使用。 一、关于cvf（compaq virtual fortran） CVF操作方便，限制不多。一般常用的有两种新建工程的方法： a)打开cvf，点击新建file----new，弹出选择对话框，简单的程序可以选择fortran console application。在左侧project name处填写工程名字，在location处填写（选择）工程保存路径。注意cvf中要求不允许出现汉语路径。操作顺序图如下 第一步： 第二步：

第三步： 第四步：运行

至此，所有步骤完成。 b)直接打开cvf，点击新建文档（new）然后单击保存，弹出保存对话框，选择保存路径及文件夹， 文件名改为“工程名.f90”格式此处注意如果是fortran自由格式，请保存为*.f90或*.f95格式，如果是固定格式请保存为*.for或*.f格式。因为编译器是根据后缀的不同调用不同的语言编译器，否则将出错。保存文件的路径和文件夹即为该工程所在的路径和文件夹。同样不能有汉语。示意图如下：第一步 第二步、

第三步、 至此，所有的工程完成。

二、关于CVF的调试（debug） 在相应代码编辑框左侧发灰色的竖线部分，鼠标变为反三角方向的形状时即可右键鼠标，选择insert/remove BreakPoint选项，在代码左侧可以看到一个红色的标示（代码行尽量不要设在代码最后。可以设置到自认为可能发生错误的地方。或者尽量靠前设置），此时就可以按F5进入调试状态。可以在watch窗口查看各个变量，数组的值与内容。按F11进行单步运行。查看错误出处。具体的更细致的调试，请参看相关文献或书籍。这类书籍不少。 具体操作如下：

几个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

Fortran6.5 使用流程

简单说明一次编译程序的过程： 1.建立一个新的Project（File-new-选择Project选项卡-Utility Project,给定project的名称）， 确定生成Project。 2.生成一个新的程序文件（File-new-选择Files选项卡-Fortran Free Format Source File，给 定文件名），或是插入一个已有的程序文件（选项Project/Add to Project/Files）. 3.用Build菜单中的Excute选项来编译并运行程序，或是只单击Build选项来只做编译， 不运行程序。 4.要写新的程序可以另外建立一个新的Project,或是直接跟换Project中的文件。千万不要 把两个独立的程序文件放入同一个Project中，否则导致编译过程出现错误。 5.下次修改程序时，可以直接使用File/Open Workspace来打开*.dsw的Project工程文件。 一、安装 Fortran6.5安装包中包括以下文件（夹） 本人电脑-x86 64位-win7 进入X86文件夹，运行，中间需要输入序列号，记事本打开 获取，进过中间一些文件位置的选择，到百分之九十多的时候，需要选择是否安装运行环境，选是，最后安装完成。 二、运行 开始---Compaq Visual Fortran 6---Developer Studio

如果提示打开失败，则选中Developer Studio，右键----管理员获得所有权---就可以了 三、打开 出现如下对话框： 熟悉一些tips: 第一条： If your program hangs应用程序瘫痪while running outside the debugger在调试器外运行程序, attach贴上、附加the debugger调试器to your program程序by choosing Start Debug调试、除错form the Build menu, and then choosing Attach to Process. 如果你的程序在调试器外运行，导致应用程序瘫痪，通过选择Build---Start Dubug---attach to program，讲调试器连接到你的程序。 第二条： Visual可视化Fortran【公式翻译程式语言（formula translator）】uses the file extension文件扩展名. F90 to indicate表示Free Form source code自由形式的源代码, and the extensions .FOR or .F to indicate Fixed Form source code. Visual Fortran采用拓展名为.F90的文件来表示自由形式的源代码，扩展名为.FOR或.F 的文件来表示固定形式的源代码。 第三条：

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

Fortran95程序设计课后习题答案(word版方便)

第四章 1.program main implicit none write(*,*) "Have a good time." write(*,*) "That's not bad." write(*,*) '"Mary" isn''t my name.' end program 2.program main real, parameter :: PI=3 implicit none.14159 real radius write(*,*) "请输入半径长" read(*,*) radius write(*,"(' 面积='f8. 3)") radius*radius*PI end program 3.program main implicit none real grades write(*,*) "请输入成绩" read(*,*) grades write(*,"(' 调整后成绩为'f8.3)") SQRT(grades)*10.0 end program 4.integer a,b real ra,rb a=2 b=3 ra=2.0 rb=3.0 write(*,*) b/a ! 输出1, 因为使用整数计算, 小数部分会无条件舍去write(*,*) rb/ra ! 输出1.5 5.p rogram main implicit none type distance real meter, inch, cm end type type(distance) :: d write(*,*) "请输入长度:" read(*,*) d%meter d%cm = d%meter*100 d%inch = d%cm/2.54 write(*,"(f8.3'米='f8.3'厘米='f8.3'英寸')") d%meter, d%cm, d%inch end program 第五章 1.program main implicit none integer money real tax write(*,*) "请输入月收入" read(*,*) money if ( money<1000 ) then tax = 0.03 else if ( money<5000) then tax = 0.1 else tax = 0.15 end if write(*,"(' 税金为'I8)") nint(money*tax) end program 2.program main implicit none integer day character(len=20) :: tv write(*,*) "请输入星期几" read(*,*) day select case(day) case(1,4) tv = "新闻" case(2,5) tv = "电视剧" case(3,6) tv = "卡通" case(7) tv = "电影" case default write(*,*) "错误的输入" stop end select write(*,*) tv end program 3.program main implicit none integer age, money real tax write(*,*) "请输入年龄" read(*,*) age write(*,*) "请输入月收入" read(*,*) money if ( age<50 ) then if ( money<1000 ) then tax = 0.03 else if ( money<5000 )then tax = 0.10 else tax = 0.15 end if else if ( money<1000 ) then tax = 0.5 else if ( money<5000 )then tax = 0.7 else tax = 0.10 end if end if write(*,"(' 税金为'I8)") nint(money*tax) end program 4.program main implicit none integer year, days logical mod_4, mod_100, mod_400 write(*,*) "请输入年份" read(*,*) year mod_4 = ( MOD(year,4) == 0 ) mod_100 = ( MOD(year,100) == 0 ) mod_400 = ( MOD(year,400) == 0 ) if ( (mod_4 .NEQV. mod_100) .or. mod_400 ) then days = 366 else days = 365 end if write(*,"('这一年有'I3'天')") days stop end program 第六章 1.program main implicit none integer i do i=1,5 write(*,*) "Fortran" end do stop end program 2.program main implicit none integer i,sum sum = 0 do i=1,99,2 sum = sum+i end do write(*,*) sum stop end program 3.program main implicit none integer, parameter :: answer = 45 integer, parameter :: max = 5 integer weight, i do i=1,max write(*,*) "请输入体重" read(*,*) weight if ( weight==answer ) exit end do if ( i<=max ) then write(*,*) "猜对了" else write(*,*) "猜错了" end if stop end program 4.program main implicit none integer, parameter :: max=10 integer i real item real ans ans = 1.0 item = 1.0 do i=2,max item = item/real(i) ans = ans+item

Fortran程序总结

1.行的书写(行的长度、分行、续行) 一行可以是0~132个字符，空格有意义， 语句最长不超过2640个字符 一行可以有多个语句，用“;”分隔 一个语句可分行写，读行标记为&（放在尾部），但如为关键字，首尾均加&。 最多可有511个续行。 2.说明语句必须出现在可执行语句之前，格式说明语句（FORMAT语句）除外。 3.注释标志符： 自由格式：！固定格式：C * 语句分隔符：分号；（仅自由格式可以使用） 续行符：自由格式：& 申明标号：1到5位无符号整数 空格：关键字、变量、常量内部不能用空格，但相邻两者之间须用空格 4.信息处理的分类：数值处理和信息处理 现代计算机工作原理：程序存储和程序控制（冯·诺依曼） 1、运算器——算术运算、逻辑运算 2、控制器——根据指令控制计算机工作 运算器、控制器和寄存器称为中央处理器CPU 3、I/O设备——提供数据传输服务 4、总线——数据传输的公共通道 1.机器语言：二进制代码形式，可以被计算机直接执行，不可移植 2.汇编语言：用助记符来代替机器指令，容易记忆，不可移植 3.高级语言：接近自然语言（英语）的程序设计语言，又称算法语言，易学、易用、易 维护，可移植性好例：FORTRAN、BASIC、PASCAL、C、LISP、PROLOG等 5.FORTRAN90程序是一种分块结构，由若干个程序单元块组成： 主程序、外部子程序、模块、块数据单元无论是主程序单元，还是子程序单元，都是独立的程序单位，应该独立编写，它们的形式相似。 非语句行：注释语句：！后的所有字符都被编译器忽略。可独占一行，也可在其它语句之后，空行为注释行（固定格式用C和*） 6.常量的定义：常量是程序执行过程中不能变化的量。 基本数据类型有五种：整型、实型、复型、字符型和逻辑型 前三种属于数值类型，后两种为非数值类型，主要用于文字处理和控制。 FORTRAN 90通过KIND值确定整数的存储开销、最大位数和取值范围，如表所示。缺省KIND 值则取默认值4。 KIND值字节数二进位数取值范围 1 1 8 -128~127 2 2 16 -32768~32767 4 4 32 -2147483648~2147483647 格式：SELECTEND_INT_KIND(N) ： 功能：返回范围在-10N~10N的整数的kind-Value 在FORTRAN90中，也可以表示二进制、八进制、十六进制形式的无符号整数。其形式如下：二进制数：B’101101’ 或B”101101” 八进制数：O’76210’ 十六进制数：Z’1FA2’

平面四边形四节点等参单元Fortran源程序

C ************************************************ C * FINITE ELEMENT PROGRAM * C * FOR Two DIMENSIONAL ELASticity PROBLEM * C * WITH 4 NODE * C ************************************************ PROGRAM ELASTICITY character*32 dat,cch DIMENSION SK(80000),COOR(2,300),AE(4,11),MEL(5,200), & WG(4),JR(2,300),MA(600),R(600),iew(30),STRE(3,200) COMMON /CMN1/ NP,NE,NM,NR COMMON /CMN2/ N,MX,NH COMMON /CMN3/ RF(8),SKE(8,8),NN(8) WRITE(*,*)'PLEASE ENTER INPUT FILE NAME' READ(*,'(A)')DAT OPEN(4,FILE=dat,STATUS='OLD') OPEN(7,FILE='OUT',STATUS='UNKNOWN') READ(4,*)NP,NE,NM,NR WRITE(7,'(A,I6)')'NUMBER OF NODE---------------------NP=',np WRITE(7,'(A,I6)')'NUMBER OF ELEMENT------------------NE=',ne WRITE(7,'(A,I6)')'NUMBER OF MATERIAL-----------------NM=',nm WRITE(7,'(A,I6)')'NUMBER OF surporting---------------NC=',Nr CALL INPUT (JR,COOR,AE,MEL) CALL CBAND (MA,JR,MEL) DO I=1,NH SK(I)= enddo CALL SK0(SK,MEL,COOR,JR,MA,AE) do I=1,N R(I)= enddo pause 'aaa' stop READ(4,*)NCP,NBE,iz WRITE(*,'(5i8)')NCP,NBE,iz WRITE(7,'(5i8)')NCP,NBE,iz IF CONCR(NCP,R,JR) IF CALL BODYR(NBE,R,MEL,COOR,JR,AE) IF do jj=1,iz READ (4,*)Js,nse,(WG(I),I=1,4) read(4,*)(iew(m),m=1,nse) CALL FACER(iew,NSE,R,MEL,COOR,JR,WG) enddo endif CALL DECOP (SK,MA)

计算材料学Fortran程序汇总

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 ============================================================ ! Monte Carlo Simulation of Two Dimensional Diffusion INTEGER X,Y,XY(1:1000,1:1000) REAL XYM(1:1000) ! X:INSTANTANEOUS POSITION OF ATOM ! XY(J,I)：X*Y ,J:第几天实验,I:第几步跳跃 ! XYM(I): THE MEAN OF XY WRITE(*,*) "实验天数JMAX，实验次数IMAX" READ(*,*) JMAX,IMAX ISEED=RTC() DO J=1,JMAX !第几天实验 X=0 !!! Y=0 !!! DO I=1,IMAX !第几步跳跃 RN=RAN(ISEED) IF(RN.LT.0.25)THEN x=x y=y-1 END IF IF(RN.LT.0.5.AND.RN.GE.0.25)THEN x=x y=y+1 END IF IF(RN.LT.0.75.AND.RN.GE.0.5)THEN

FORTRAN经典入门程序20例

对于FORTRAN的初学者。这些例子可作为小练习。 1.例题：计算工资问题。每小时工资为RATE,如果工作超过40小时，加班呢部分工资是正常时间工资的1.5倍。 C Payroll with overtime program payroll real rate, hours, pay read (*,*) rate, hours if (hours>40) then pay=40*rate+(hours-40)*1.5*rate else pay=hours*rate END IF print *,"rate=", rate print *, "hours=", hours print *,"pay=",pay end 2.学生成绩问题。大于80为A级。大于60小于80为B级。小于60为C级。IF的嵌套。注意空格可以看清楚else if ,end if,print的内容. PROGRAM GRADESTUDENT READ *,GRADE IF (GRADE .GE. 80) THEN PRINT *,GRADE,"GRADE=>A" ELSE IF (GRADE.LT.60) THEN PRINT*,GRADE,"GRADE=>C" ELSE PRINT*,GRADE,"GRADE=>B" END IF END IF END 3.三个数按从小到大排序。 PROGRAM MAXMIN REAL A,B,C,T READ *,A,B,C IF (A.GT.B) THEN T=A A=B B=T ELSE END IF IF (B.GT.C) THEN T=B B=C

PROGRAM EQUATION READ (*,*) X IF (X .GE. 0.0) Y=3*X+6 IF (X .LT. 0.0) Y=-X**2+2*X-8 PRINT *,"X=",X, "Y=",Y END 6.CONTINUE语句。为了使循环的起止范围清晰，使用终端语句 例：求5！ program ex1 FACT=1 do 10, N=1,5,1 ##DO 语句标号，循环不变量=初值，终值，步长FACT=FACT*N 10continue print*,FACT end 求1！+2！+3！+···+N！ READ *,N FACT=1

优秀的fortran程序编程规范

Programming Guidelines for PARAMESH Software Development(NOTE: This document is heavily based upon the Introduction This document describes the programming guidelines to be used by software developers wishing to contribute software to the PARAMESH, parallel, adaptive mesh refinement software. We welcome people to contribute software and/or bug fixes to the PARAMESH AMR software. Software to be added to PARAMESH can come in 2 forms: ?Improvements to the basic PARAMESH kernal software found in the mpi_source, source and hearders directories. ?Software the addes additional functionality to PARAMESH. This type of software should be added as separate entities within the utilities directory. Complete applications should not be added as part of PARAMESH. PARAMESH is only meant to be a tool which supports parallel adaptive mesh applications and any software which supports this goal will be considered for acceptance into PARAMESH. For instance, a solver for the poisson equation that works with PARAMESH would be acceptable, but an application that solves the equation of gas dynamics would not. The PARAMESH software is slowly being evolved to be consistent with this document. Any new software which is contributed should follow these guidlines. If not, it will be rejected. This document deals mainly with Fortran 90, since most new PARAMESH software will probably be written in that language. [Throughout this document, the term "Fortran" should be understood to mean Fortran 90.] Since we expect C and C++ also to be used, a separate document dealing with them will be developed in the future. In the meantime, this document can serve as a general guideline for developing code to be used with PARAMESH in those programming languages. The guidlines in this document should be adhered to by ANY software which will be released as part of the PARAMESH package of source code. This includes software 'utilities' (stored in the paramesh/utilities directory) which add functionality to PARAMESH for different algorithms. It also should be applied to any new code developed and added to the main source code for PARAMESH in the paramesh/source, paramesh/mpi_source, or paramesh/headers directories.