计算材料学1.pptVIP

计算材料学Computational materials science 张晖 材料科学与工程学院 Email: huizhang@ahut.edu.cn; Mobile 本堂课内容： 计算材料学概论 本课程研究对象和内容 学习要求 Fortran语言语法结构简单回顾 1.2.1编程实例 C------求两种平均値------ PROGRAM Example_1_1 REAL a, b, av1, av2 READ (*,*) a, b av1 = (a + b)/2 av2 = sqrt(a*b) WRITE(*,*) av1, av2 END 程序中第1行是注释行，对程序起说明作用。F77注释行是以“C”或“*”作为该行第一个字符的，F90可在任意一行末以“!”开始作为注释符。第2行是主程序名，第3行是变量类型定义，第4行是输入语句，第5－6行是赋值部分，第7行是打印输出语句，最后是程序结束。注意F77的固定书写格式和F90的自由格式。F90中用“；”将两行并为一行。 常量和变量类型 a)??? 常量 常量是指其值始终不变的一些量。整型、实型、双精度和复型常量是算数型常量，也为常数。整型：实型：复? 型：逻辑型：字符型：数? 组： b)??? 变量 变量是指在程序运行期间其值是可以变化的量。系统为程序中的每一个变量开辟一个存储单元，用来存放变量的值。 常量是分为类型的，而变量是用来存放常量的，因此变量也相应地区分为整型变量INTEGER、实型变量REAL、双精度变量DOUBLE PRECISION、复型变量COMPLEX、逻辑型变量LOGICAL、字符型变量CHARACTER。在程序中应当说明哪些变量是整型变量，哪些变量是实型变量。变量在内存中所占的字节数和数据存储形式与相应类型的常数相同。例如，实型变量一般占4个字节，按指数形式存放。在程序中规定变量的类型可以用以下几种方法。 隐含约定：Fortran规定，凡以字母I，J，K，L，M，N六个字母开头的变量名，如无另外说明则为整型变量。以其它字母开头的变量为实型变量。可以将这个隐含约定称为“I—N规则” ，表示用I到N之间的字母开头的变量为整型。例如：I，J，IMAX，NUMBER，LINE，JOB，Kl为整型变量，而A，Bl，COUNT，AMOUNT，TOTAL，BOOK为实型变量。 内在函数 1.4.4赋值语句 f)??? 执行时的类型转换 一个算术赋值语句中的被赋值的变量(V)和表达式(e)的类型可以相同，也可以不相同。F77作如下规定： ????????? 如果变量V与表达式e的类型相同，则直接进行赋值。 ????????? 如果类型不同，则应先进行表达式的求值，然后将该表达式的值转换为被赋值变量的类型。如： I＝3.5*2.1 表达式3.5*2.1的值为7.35，实型。而变量I为整型，因此I的值等于7。又如: T＝3*5/7 表达式的值为2，整型。由于T为实型，故系统先将整数2转换成实数2.0，再赋给变量T，T的值为2.0。 当赋值号两侧的类型不同时，往往会产生程序设计者事先预想不到的结果。所以在编写程序时，应尽可能使两侧保持同类型。 关系表达式 a) 关系运算符 关系运算符就是关系比较符。Fortran中有六个关系运算符: 选择结构 2.3.1 块IF构造 a) 块IF的组成和执行 块IF的一般形式可写成： IF(逻辑表达式) THEN???? ← 块IF语句 块1??????????????? ← then块 ELSE 块2??????????????? ← else块 END IF 例：打印学生考试成绩，大于等于80分的为“A”，大于等于60分而小于80分的为“B”，小于60分的为“C”： READ *,GRADE IF(GRADE.GE.8O) PRlNT *，”A” IF(GRADE.GE.6O .AND. GRADE.LT.80) PRlNT *，”B” IF(GRADE.LT.6O) PRlNT *，”C” 程 序 举 例 例4-3 输入三个数，把它们 按从大到小的顺序输出。 循环结构 3.1.2有循环变量的DO构造 a)??? DO语句和循环次数 DO构造可分为不带循环变量与带循环变量两种形式。 当需要执行的循环次数为已知时，用DO语句实现循环比较方便。它的一般形式为： DO [[标号][，]] 循环变量＝初值式，终值式[，增量式] 例：循环读入学生的学号和成绩30次并打