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

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

[复制链接]
$ A* L8 R- z. G2 Y) n5 {

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

9 [4 o. ]5 x/ c! l

书面格式

: Q# `9 s1 w& Y B; p; N+ t$ M: K

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

/ {4 L. |: j0 f- D* `

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

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

数据格式

5 ` Q6 q' D; O# B! {. `; ^# S0 D9 r

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

% ~; h% s" w1 s# a: {" v

数据类型:整数(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为双精度。实数的精度设置十分重要,将会影响最后的计算结果是否正确。

0 I# l9 T3 G% M v7 N- O0 X
real(kind=4) :: distance
+ |/ L+ }1 B5 `, @/ `

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

' Z. c+ \* D0 V4 n- C5 v
program test 4 z. a: C5 E: p! S# i" z implicit none + r, U4 P' x- N/ A real(KIND=4) :: a 1 o. p3 S! x: f4 D real(KIND=8) :: b/ G. V* j& W, A) C4 r/ ^" D* a a=0.123456789123456! K5 J' {# A$ P# i4 z9 B: j b=0.123456789123456) a! f# k0 _% E6 M6 O( y0 Q) U write(*,*) a,b=,a,b/ m `( W# i! }% k# V* _/ O& L end/ U7 M7 w9 i" U5 u1 `. o$ D !!!输出结果: a,b = 0.123456791 0.12345679104328156; ?- b( h, B- u/ b& ^5 a4 Y
" k# J ^( S- c" l, f- I! _, H+ h

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

X( X. b) ?) Q0 I( i3 w/ t
int(7.8) = 7 !忽略小数部分 8 C$ `: d' l: M$ _0 ^8 b nint(7.8) = 8 !四舍五入后的结果* h! N5 M; t m* l5 c ceiling(7.8) = 8 !大于等于的整数 2 O. K. G5 q9 u3 s! L& {3 J floor(7.8) = 7 !小于等于的整数
字符型
charcter(len=16) :: atmosphere' b/ x& Z( D6 j: ^ charcter(len=20) :: ocean : z7 M! _ B4 ?3 c) e) \, _/ }1 { b- }. h. B& G v atmosphere = 12345678905 b" V, s) ?* l) {3 D$ P: I4 A7 ~ D ocean = 123456789123456 !字符赋值注意使用引号 $ y+ k+ {% M; W4 s" _ 9 R+ U: P6 A8 p: p- D
逻辑型

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

7 N& w! f9 S8 E
logical :: a,b% ]4 z- Q6 H3 J$ u" h6 [& J a = .true. 3 H, L! j |, w0 {/ `* ~- ? b = .false.1 `5 \0 U! p' v if(a)then( [# ~6 N4 Y) I . V- a9 E( K% h% B. n. z' t endif0 k$ @7 ~/ U9 M! b8 s1 c- ]$ ?( j
复数型(不常用)
complex :: a, w, G: U) K6 k& L8 [, \ a = (a1,a2)* _6 K. n& R3 X& j3 |9 ]
: A( h( j2 _" g3 M2 K; Y7 L

Fortran的数学表达

1 m# C I$ P3 W+ r

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

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

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

/ a, p" R. ~* J8 ~& Y% P8 Z

Fortran的输入输出

' C8 H7 u( g R9 e# w% Q

输入

! _7 R% T* w( h1 C5 J' x

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

方式1
integer :: a" V3 B2 @. t- M- a$ ]- \9 } wirte(*,100) a8 u6 Q0 c* x( x+ S 100 format(I4)/ q7 P T ], D! I& F
9 |6 e0 Y0 q ^; H G

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

, o3 ^0 R# g9 t# G( P% t

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

方式2,推荐
write(*,"(1X,I5)")9 Q% z+ g) E5 Q - N% G) d7 d M. M write(*,"((,i2,,,i2,))") i,j" V2 ^2 _& g5 w& k) y0 k !输出:( 1, 2) " R& g# z k- L* c& e & C/ f; N1 v- y! F/ ?% Z write(*,(5f7.2) list(4 Q( B+ j) I/ i !每行5个的形式输出list数组 0 r- T6 H# q8 y% G$ y. z ' p& Y! T0 n& ^7 O, V6 N; M write(*,(mf7.2) list(3 [% z% R) T( F K, W% ^( B3 W m>=size(list(). L: U5 w$ s2 c6 Z3 X- X: U. h !在一行内全部输出 0 o: f# K0 {( b' U6 D
方式3
character(len=10) :: fmtstring # u! H$ l% m( I. s1 A fmtstring = (I3) n) w$ X! p9 X+ Z* |- P9 g write(*,fmtstring) 3 ( D0 l' V! ^" a4 K' s6 B5 L
5 D0 N( A. h0 {4 c# r8 P

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

6 C# j1 M4 g8 L- h0 l' e

常用的格式化

4 F$ b/ @2 r" 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个字节

输出

* F% p3 s- E( Q) k* ]
write(*,*) "nice to meet you!"
write,作用是显示后面双引号中的字符串。其中括号里第一个*表示输出位置使用默认值,即屏幕,如果写成unit=6表示的也是屏幕;第二个*表示不特别设置输出格式。write每执行一次命令后,会在佛那个换到下一行准备下一次的输出其中“”是为了输出字符,若想要输出“”需要使用两层双引号双引号和单引号都可以用来输出字符,但是若想要输出双引号时,就需要连续使用两个双引号
print *, "nice to meet you too!"
4 r! @* e6 t: d) ?$ {) z

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

/ q" O ?6 E8 N; k. _1 `* P' X p ' E! i$ Y2 r& s7 T1 ~* ?8 |6 r7 w + U4 ?; D1 S( V5 s % s* X4 w: g/ S/ T$ V( K S) \# E' V) V" h4 A( u t2 d9 J8 t& C" |
回复

举报 使用道具

相关帖子

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