實數乘法

Question

我有一個包含一些實數乘法事件的代碼，並且無論輸入變量如何，我總能得到該特定子例程的答案零。

所有變量都是real*4類型聲明並具有像例如， a = 5，b = 3.5 & c = 0.0005。

我使用柏拉圖IDE

`

Subroutine Load_Modulation_100() 

Real*4 :: Te  
Real*4 :: Ee  
Real*4 :: Et  
Real*4 :: QLHV  
Real*4 :: mfg  
Real*4 :: Qfeul 
Real*4 :: mf  
Real*4 :: Pe  
Real*4 :: Pt 

Te  = 100 
Ee  = 28.2 
mf  = 0.0005 
QLHV = 47141 
Et  = 66.8  
Cpfg = 4.5 

      Qfeul = (mf * QLHV) 
      Pe = (Qfeul * (Ee/100)) 
      Pt = (Qfeul * (Et/100)) 

Write (*,*) Pe, 'kW' 
Write (*,*) Pt, 'kW' 

End Subroutine Load_Modulation_100` 

Answer 1

可能並非如此，但整數截斷可能是一個問題...... 我會更加明確，可能使用一些調試。

Subroutine Load_Modulation_100() 
IMPLICIT NONE 
Real*4 :: Te  
Real*4 :: Ee  
Real*4, PARAMETER :: Et = 66.8E0 
Real*4 :: QLHV  
Real*4 :: mfg  
Real*4 :: Qfeul 
Real*4 :: mf  
Real*4 :: Pe  
Real*4 :: Pt 
LOGICAL*4 :: IO = .TRUE. 

! Te  = 100 ! This is a float, it just does;t jump off the page as one 
Te  = 100.0E0 ! 
Ee  = 28.2 
mf  = 0.0005 
QLHV = 47141 ! This is a float, it just does't jump off the page as one 
!Moved to parameter Et  = 66.8  
Cpfg = 4.5 

    IF(IO) WRITE(*,*)' (Ee/100)=', (Ee/100) 
     Qfeul = (mf * QLHV) 
!  Pe = (Qfeul * (Ee/100)) 
     Pe = (Qfeul * (Ee/FLOAT(100))) 
!  Pt = (Qfeul * (Et/100)) 
     Pt = (Qfeul * (Et/100.0E0)) 

Write (*,*) Pe, 'kW' 
Write (*,*) Pt, 'kW' 

RETURN 
End Subroutine Load_Modulation_100!`