收藏本站 劰载中...网站公告 | 吾爱海洋论坛交流QQ群:835383472

海洋数据分析 -海洋数据分析岛

[复制链接]
5 g8 n0 e0 d0 ~( g

要开始编程之前,首要要了解书写的规则以及能处理的数据类型。了解咱们手上的工具,才能更好得将其利用。

( }' Z+ o, }' a+ ?

书面格式

* O. q; S0 V; s+ M; |! S) y

代码的编写格式有两种,自由格式(Free Format)和固定格式(Fixed Format)。程序代码扩展名为.F / .FOR为固定格式,.F90为自由格式。推荐使用自由格式进行编码。

+ a; E V* z/ g2 N1 l/ f5 u

但是自由格式就是完全自由,可以随意编码吗?答案当然是否定的。与固定格式相比,它不会规定每一行的第几个字符有什么作用,但是它也有需要遵守的规则

每行可以编写132个字符叹号“!”后面的文本都是注释,没有(编程)意义代码不为区分大小写两行代码可以用符号“&”进行连接变量取名需要以字母开头,可以表现为字母和数字混合,且长度不超过32
自由格式的简单例子,引自高PPT

数据格式

/ z& p# Q% ^* m' V

与python不同,fortran程序中出现的每一个数据必须声明类型,且每一数只有一种类型。数据类型是指使用Fortran在计算机内存中记录文本、数值等数据的最小单位及方法。

* S. G/ [7 |5 n) i) E; O0 j9 N

数据类型:整数(integer)、浮点数(real)、复数(complex)、字符(character)、逻辑判断(logical)

整数需要提前声明 integer :: x,y,z数目大小不能太大,长整型占用32bits(4bytes)的空间,保存的范围为 −231+1∼231-2^{31}+1\sim2^{31} 之间,有的编译器还会提供短整型、更短整型整数和实数之间的转化,可以直接使用编译器提供的函数
real(9) = 9.0
实数/浮点数

1.精度问题。利用kind进行精度说明,值取4为单精度,值取8为双精度。实数的精度设置十分重要,将会影响最后的计算结果是否正确。

8 W; L5 s3 d. K$ L
real(kind=4) :: distance
9 I" M& i9 C' {$ _

2.两个实数相等问题。输出后,由于ab精度不同,他们的值也不同。

# h& t; C2 k1 I1 ~$ [( {7 j8 B
program test ; _' r0 }* z3 @ implicit none$ E+ K$ U0 [+ `- i real(KIND=4) :: a1 u: t2 W' ^" M" o real(KIND=8) :: b 3 z- C2 O( b7 f1 w2 E" @ a=0.123456789123456 ! `$ U% K; ]/ a, g! ]. H; h b=0.123456789123456 1 x5 s- e6 A4 @0 V" ^ write(*,*) a,b=,a,b , t. G/ |$ K8 t7 r" S end . Z- {& a0 R8 q1 I2 L !!!输出结果: a,b = 0.123456791 0.12345679104328156 3 j+ }* L) S6 N C" p
6 r! S$ |. r- _% N8 f

3.实数与整数之间的转换

, E! k! f; ]) b7 C; ^
int(7.8) = 7 !忽略小数部分) P, V$ ?" x; T5 y+ l4 g7 ]% b: M+ d nint(7.8) = 8 !四舍五入后的结果 2 P; W$ |) g e( A" _+ ?, W ceiling(7.8) = 8 !大于等于的整数 $ I9 T1 T. b5 F* H floor(7.8) = 7 !小于等于的整数
字符型
charcter(len=16) :: atmosphere 9 W7 p$ A3 W' a% N charcter(len=20) :: ocean6 O3 r! L/ F5 M3 N 8 c9 {1 K& f5 x- X* N atmosphere = 1234567890/ ?4 @$ j- J& \ ocean = 123456789123456 !字符赋值注意使用引号8 D4 `$ Z9 r# b' ?1 P 5 O7 Z4 {/ P7 S6 q7 u4 s! p9 r
逻辑型

对于较长的代码,灵活使用逻辑型变量,可以提高代码的多样性。

* E2 m4 S$ {, h$ E, E7 d" R
logical :: a,b ! }: G7 s6 D8 c a = .true.- h s7 D5 W) s6 H# ^# \1 o. G b = .false. 0 Q9 ]1 E! L) ]- D if(a)then & _- X. W, l% j; S 5 A; S( i$ F8 d: D4 R/ F endif3 |8 H: o }/ R+ V/ q/ \
复数型(不常用)
complex :: a 9 ^% O! x" f6 Y# e0 v a = (a1,a2)8 a0 J4 c$ }' D1 d
2 l v0 V/ X4 \0 ~4 O& Q- F

Fortran的数学表达

2 d' Y1 r2 w+ D a

数学运算符的运算优先级,由高到低排列如下

( ) 括号里优先计算* * 乘幂* 乘法 / 除法+ 加法 - 减法

需要注意:乘号容易忘记写

+ E0 N3 ]* f3 X. t2 S3 S* t& z

Fortran的输入输出

- H4 j8 }4 T6 p4 S7 F/ {9 V

输入

4 A' v' [$ i2 D3 K1 f

建议使用格式化进行输入输出,这样不仅可以使得显示更加美观,在某些情况下,设置恰当的输入格式才能得到正确的数据。

方式1
integer :: a) W$ h. l* S4 u$ d( {6 Q: G7 ~ wirte(*,100) a5 T% w3 V0 Z2 v! n) x 100 format(I4) 8 J5 f$ m" p% ]& }( u% s
, s$ l+ e: ]& @" Q, X, { V

其中代号为100的这一行,作为上一行设置的输出格式

- }/ f* p1 ]# N% s0 b$ t1 u

当不设置输出格式时,我们发现数字前面会出现一段额外的空格,当设置的输出格式为(I4)时,也是表示用4列来显示整数a的数值。

方式2,推荐
write(*,"(1X,I5)")& ^8 n0 g* i- P 7 z+ [% N; f7 C2 p% r- E4 m write(*,"((,i2,,,i2,))") i,j8 X/ b# M$ R. m b5 v !输出:( 1, 2)$ O5 L" M- K3 L" \' C, { % `, p2 m. c6 A! S write(*,(5f7.2) list( . h: \. u3 J0 d !每行5个的形式输出list数组4 ]$ Y7 w7 L' N5 ]: S% p+ g$ p4 [ 2 n3 U$ h. B7 E" Z; N F9 w T3 A write(*,(mf7.2) list(9 p* D6 v/ j. o# R m>=size(list()) e2 G& |8 G# \ !在一行内全部输出 J, j# Y; k' y- ~! U% J- ~
方式3
character(len=10) :: fmtstring) h! j+ ?' O! A% [2 w& V fmtstring = (I3)! E/ V9 t$ {" X9 K$ a' p% @ write(*,fmtstring) 33 z1 v- K% `8 [( |& l5 W
: | | L3 z" ^3 I2 R( ^

声明长度为10的字符串,内容设置为输出一个整数的格式。

5 P: i$ w! V! S' o7 z. u; M

常用的格式化

, D3 H6 s, D7 d/ S* [# u% C& |

一般来说,最常用的几种格式为“I, F, E, A, X”,不同格式之间可用用,隔开,或者不用

Iw[.m], egI4,I0是特殊描述符,此时w=整数的长度以w个字符宽度输出整数,至少输出m个数字,当数值宽度大于w时输出*。Fw.d, egF9.3以w个字符宽度输出浮点数,小数字符占d个宽,数值大雨输出宽度时出现*Ew.d[Ee], egE15.7,E9.2E3w个宽度输出浮点数,小数部分占d个,指数部分占e个Aw以w个字符宽来输出字符串,不足w个字符的部分会自动补上空格nX输出位置向右移动n位/换行输出Tc把输出的位置移动到本行的第c个字节

输出

2 c S6 u9 L/ b0 T. X# X4 L
write(*,*) "nice to meet you!"
write,作用是显示后面双引号中的字符串。其中括号里第一个*表示输出位置使用默认值,即屏幕,如果写成unit=6表示的也是屏幕;第二个*表示不特别设置输出格式。write每执行一次命令后,会在佛那个换到下一行准备下一次的输出其中“”是为了输出字符,若想要输出“”需要使用两层双引号双引号和单引号都可以用来输出字符,但是若想要输出双引号时,就需要连续使用两个双引号
print *, "nice to meet you too!"
* C: }3 f& N( ~8 d

这里的*表示不限定输出个格式,且该命令只能正对屏幕来使用。因此建议使用write来作为输出工具

. b# ]; Y0 Q$ m6 p+ F; _, |# Z & T {+ p4 N0 V) O) { g9 i# g | ?8 Q5 v3 R* j r7 L 2 t) n4 L+ H! Q. h1 c2 n. Q6 I 7 M: V) b5 b! }+ \
回复

举报 使用道具

相关帖子

全部回帖
暂无回帖,快来参与回复吧
懒得打字?点击右侧快捷回复 【吾爱海洋论坛发文有奖】
您需要登录后才可以回帖 登录 | 立即注册
太美
活跃在2026-4-10
快速回复 返回顶部 返回列表