この場合のフォートランのフォーマットどうすれば？

質問

今もともとあるデータをa倍してb引く値を求めそれを違うフォーマットで書き出したいのですがどういうふうに書けばいいのでしょう？

元のデータ

17346 17371 17420 17450 17409 17327 17289 17276
17237 17225 17277 17344 17365 17304 17224 17204
17223 17209 17195 17245 17274 17285 17321 17301
17302 17352 17355 17368 17408 17438 17428 17338
17300 17356 17378 17375 17350 17285 17251 17234

全て5桁の正のあたいの数値が上には40個しかのせませんでしたが8ｘ5000あります

書き出したいデータ
0.623E-001 0.108E+000 0.851E-001 0.128E+000 0.110E+000
0.125E+000 0.105E+000 0.126E+000 0.107E+000 0.114E+000
0.855E-001 0.834E-001 0.468E-001 0.326E-001 -0.105E-001
0.623E-001 0.108E+000 0.851E-001 0.128E+000 0.110E+000
0.125E+000 0.105E+000 0.126E+000 0.107E+000 0.114E+000
0.855E-001 0.834E-001 0.468E-001 0.326E-001 -0.105E-001
0.623E-001 0.108E+000 0.851E-001 0.128E+000 0.110E+000
0.125E+000 0.105E+000 0.126E+000 0.107E+000 0.114E+000

もとデータの値に定数をかけて引いた値を

空欄2つ（値が負の場合は1つ），値を0.～（小数点3桁まで），10の何乗か（10の2乗ならE+002，10のマイナス100乗ならE-100）

というフォーマットで５ｘ8000にしたいのですこのときはフォートランでどう書けばいいのでしょうか。

ソフトはMicrosoft Visual Studio 2008をつかっています。

よろしくお願いします

回答No.1

boogeyman · ベストアンサー

Visual StudioのFotranは知りませんが、一般的なＦｏｒｔｒａｎの文法では…
（符号と空白を含む）全体の桁数が１２桁、小数点以下が３桁、指数部が３桁ですので、
「Ｅ１２．３Ｅ３」（←全角で記述してます）で良いと思います。

お礼

character*30 fileout
dimension x(10000)

write(*,*)'out put file name'
read(*,*)fileout

do 100 i=1,10000
x(i)=10/i
100 continue

open(10,file=fileout,status='new')
do 110 i=1,10000
write(10,1000) x(i)
110 continue
1000 format (' ',e12.3e3)
close (10)

end

こう書いたら

0.100E+002
0.500E+001
0.300E+001
0.200E+001
0.200E+001
0.100E+001
0.100E+001
0.100E+001
0.100E+001
0.100E+001
0.000E+000
0.000E+000
0.000E+000
0.000E+000
0.000E+000
0.000E+000
0.000E+000
0.000E+000
0.000E+000

という感じで指数のマイナスが出ないのですがどうしたらいいのでしょう？

回答No.2

boogeyman

>という感じで指数のマイナスが出ないのですがどうしたらいいのでしょう？
問題は、「x(i)=10/i」ですね。

整数の定数（１０）を、整数の変数（ｉ）で割っているので、演算が整数（小数以下切り捨て）になります。
その結果、
X(1) = 10
X(2) = 5
X(3) = 3
X(4) = 2
X(5) = 1
X(6) = 1
：
X(9) = 1
X(10) = 1
X(11) = 0
X(12) = 0
　：　（以下、全てゼロ）
となります。

「ｘ（ｉ） = １０．/ｉ」、「ｘ（ｉ） = 10/ＦＬＯＡＴ（i）」、「ｘ（ｉ） = 10．/ＦＬＯＡＴ（i）」などに書き換え、実数で演算するようにして下さい。

もし、ＦＬＯＡＴ（）関数が無ければ、
do 100 i=1,10000
a=i ! (←整数ｉを実数ａに変換)
x(i)=10./a ! (←実数で演算)
100 continue
と書き換えてください．

お礼

なるほどそうなっているのですね。出来ました。ありがとうございます！！