fortran90

0熱度

1回答

相同的操作數在fortran90語言中產生與REAL類型單精度不同的結果，並且它似乎是編譯器選項給出不同的結果。 real , dimension(n_pft) :: sapwood_ratio real , dimension(n_pft) :: qsw real , dimension(n_pft) :: SLA !these 4 variables are also real type

1熱度

1回答

在Fortran中使用遞歸函數時出現隨機大浮點錯誤和NaN

我正在嘗試使用b樣條和遞歸函數創建箔片部分。問題在於它有時會返回一些非常大的浮點數或NaN甚至無窮大的值。出於某種原因，哪些值在重新編譯後似乎會改變，有時甚至比其他值更大。它通常返回的值是正確的。有沒有人有任何建議可能是什麼問題？任何幫助深表感謝。我知道我的代碼很長，但我認爲一切都是必要的。首先是程序，然後是具有子程序和函數的模塊。正確的答案應該是： xb = 0 0.0147 0.0529

0熱度

1回答

用ifort編譯Fortran 77代碼：libg2c缺失

我正在編譯一些fortran 77代碼，編譯器爲ifort。該代碼需要鏈接到libg2c，但該庫缺失。我在哪裏可以找到它？我在遠程計算機羣集上工作，有沒有辦法只在本地安裝這個庫？

0熱度

1回答

如何聲明來自文件的未知行數

我必須使用可以更改由於前面的操作（來自python腳本）的行數的文件，並且該行數將用於聲明其他變量，喜歡這個。 integer NumberOfLine !This is not work,Of course. real(8) F(FixedDimension,NumberOfLine) integer, parameter :: NewDimension = ANumber*NumberOf

0熱度

1回答

輸出文件包含額外的零

我有一個小問題，這個簡單的Fortran程序 PROGRAM Calc_U0 Implicit None real, dimension(:,:), allocatable :: x,y integer :: n=16589,m=9,nm=6,I,J,r open (unit = 7, file = "U0.txt") open (unit=99, file='slice06.txt'

0熱度

1回答

設置隨機種子的問題

我正在嘗試編寫一個Montecarlo算法來模擬羣體中的代理之間的交互。該算法需要在每次迭代時調用兩個隨機數（例如，10^9次迭代）。這裏我的問題是，每次我改變種子（以獲得不同的實現）時，RNG給我相同的輸出（相同的蒙特卡洛事件）。我嘗試了不同的方式來實施它，但無濟於事。我完全是Fortran的新手，從MATLAB複製這段代碼。我在執行此代碼的方式上做錯了什麼？下面是我的嘗試： progra

0熱度

1回答

如何模塊

我有很多的模塊中聲明的變量，如 module test use othermod, only: n integer, dimension(n) :: var0 real, dimension(n) :: var1 real, dimension(n) :: var2 ..... real, dimension(n) :: var1000 end module test 中創

0熱度

2回答

Fortran 90中的DATA語句的語法錯誤

我必須計算一些復積分，爲此我從我的導師舊程序中編寫Fortran 77.但是我幾乎沒有遇到任何問題。大多數與DATA語句的語法錯誤有關。這是代碼的函數計算實積分的一部分： FUNCTION CAUSSA(F,A,B,EPS) IMPLICIT DOUBLE PRECISION (A-H,O-Z) external f REAL :: W(12),X(12)

1熱度

3回答

在long fortran代碼中確定if-endif和do-enddo結構

我想修改一個非常大的Fortran遺留代碼（fortran77），但由於代碼的作者並不在意寫評論，也沒有使用適當的縮進等等，我在修改代碼時遇到了很大的困難。我需要的是某種方式（一個包或一個命令行工具），它使我能夠識別給定結構的末尾，如if-then或do-enddo，如果我將它展示在起始位置。例如： if(x .eq. 0)then if (y .eq. 0)then prin

0熱度

1回答

如果數組的等級/大小可用，是否有任何理由使用假定大小的數組？

如果我知道傳遞給函數或子例程的數組的等級和/或大小，是否有理由使用假定形狀或假定大小的數組？例如，如果我可以代替 function f(a,m,n) real,dimension(*),intent(inout) :: a ! ... end function 與 function f(a,m,n) real,dimension(m,n),intent(inou