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

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

[复制链接]
5 [ a0 O8 z2 o( g$ V" l

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

9 Z% j$ V" ?" f: D2 [. i* z

书面格式

' \9 R4 g) x- |3 x {# }8 y

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

. f( V, j, p' N7 s8 T5 u$ h* k+ s

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

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

数据格式

6 h7 e* i# o+ m# l( G

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

) h7 o; y0 F M) @

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

) n6 l. i, S/ E
real(kind=4) :: distance
/ J7 U8 f. O' T) r

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

0 H; J6 I: f6 M: U* G& X
program test2 o' M& [+ `5 ]( f+ e1 i" f implicit none2 h2 H; F# @2 x5 A) b real(KIND=4) :: a8 } `6 R, Q4 D% s6 H5 I% u real(KIND=8) :: b% W8 Y; R% M0 `; A# ^ a=0.123456789123456 , `& `1 n2 G" W b=0.123456789123456$ z4 E8 u' ?2 J: b/ J# Q5 i write(*,*) a,b=,a,b4 @- V* R( c& L' ~9 ~* r end 0 U/ I+ \( `- ]8 w5 n !!!输出结果: a,b = 0.123456791 0.12345679104328156 ) u/ T* `, X! {5 V' A
( E. N# l7 m5 J8 v! D1 l

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

$ B0 O6 @- Z5 h
int(7.8) = 7 !忽略小数部分1 r$ N8 @% a* G. b& X# e1 y nint(7.8) = 8 !四舍五入后的结果 B- R+ L1 F. G ceiling(7.8) = 8 !大于等于的整数) U G0 u, ~! w/ [/ U6 m floor(7.8) = 7 !小于等于的整数
字符型
charcter(len=16) :: atmosphere 5 X; {& a, k a5 h% @6 s% z* s charcter(len=20) :: ocean. h/ U4 N/ K4 L" y2 G + z8 b$ e0 f, G1 M. ^0 x- _ atmosphere = 1234567890* V! S! @8 y+ p( o( I/ Q4 f" z4 B ocean = 123456789123456 !字符赋值注意使用引号 / U* E: E; W) H/ o) O/ f . H @$ P8 c/ }
逻辑型

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

3 M, G) w6 ^2 q3 d
logical :: a,b4 U+ T. D% x7 o5 \9 X3 c a = .true. , H% R* \0 P" J% {# g8 h/ T b = .false. ' @6 y% G9 M0 z v+ B* }6 Z if(a)then # W4 h$ P0 G1 z) o" F3 O8 x) }) s% M : q' I& |4 u9 o1 y& I endif & x' ]! M7 i5 h; D5 `
复数型(不常用)
complex :: a* U# x) X2 T$ E3 I& D2 Y a = (a1,a2) 5 `. D5 N/ ^6 {* z% L X
* P- S' F' h. { K/ h% ~% ^; q( T

Fortran的数学表达

! D$ q3 ~8 W) ?: Y& p. K, l# h; {! J8 |

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

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

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

$ K: V: T7 e* R0 ~- z- I( r

Fortran的输入输出

P. }7 \5 F6 @, p2 C. [

输入

- Q- k& f# F% R2 [" @: t6 v

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

方式1
integer :: a 6 t9 s5 m* @% [ wirte(*,100) a3 o+ }# L5 a7 Q. F 100 format(I4) 8 w. ^' p( B1 }' n# _3 y/ T" i K
0 L o! b% S$ @1 j4 J/ [! Q

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

+ p s( s: M8 G/ t% b

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

方式2,推荐
write(*,"(1X,I5)") . \* j: m! o$ X- ^+ o$ X- p/ [6 L+ k; a9 @ t3 T. A$ K write(*,"((,i2,,,i2,))") i,j $ x0 d, b0 i# w% J0 Y, @ !输出:( 1, 2)1 ]* m6 Y1 l! L! Z, L8 R9 g ' Q' y8 w: h% s9 h% l& \: y7 }2 d* } write(*,(5f7.2) list( 1 H* x8 B/ ^, V# c! ~3 t !每行5个的形式输出list数组 3 I+ J# n- u' M" F, Z! C* { + s: a H, j: N7 S; v7 i write(*,(mf7.2) list( 8 A9 V$ a% T; o0 g) e- z) v m>=size(list()0 U) z: i2 m+ V1 Q: O !在一行内全部输出" n, @- l* y$ y, D
方式3
character(len=10) :: fmtstring& M# B+ L/ U+ U g5 c0 } fmtstring = (I3) ) Y2 r$ N2 _) p write(*,fmtstring) 3$ k1 Z& D- f; P& D5 u2 X1 [! T
) }; S, r c, o F- N7 J2 T2 e

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

1 ~7 ~ K" E7 Y' T( D8 l

常用的格式化

9 e" k: ~7 u1 P4 F4 Y& ~2 G# g+ ?

一般来说,最常用的几种格式为“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个字节

输出

6 p( H8 h, {% i+ Y) d* ?$ P
write(*,*) "nice to meet you!"
write,作用是显示后面双引号中的字符串。其中括号里第一个*表示输出位置使用默认值,即屏幕,如果写成unit=6表示的也是屏幕;第二个*表示不特别设置输出格式。write每执行一次命令后,会在佛那个换到下一行准备下一次的输出其中“”是为了输出字符,若想要输出“”需要使用两层双引号双引号和单引号都可以用来输出字符,但是若想要输出双引号时,就需要连续使用两个双引号
print *, "nice to meet you too!"
, t! \. M% n( {6 |9 D7 r9 L2 f B

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

! c+ q, V, d) {! l! a$ d0 [ + X) [3 Q3 n. d0 X( J% _2 v) T: R# P " }' O9 Z; Y j/ D7 ? 6 T5 F7 D2 W! @, S$ o 7 {; T# |$ m, p5 c N
回复

举报 使用道具

相关帖子

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