Changeset 2668

Timestamp:

Apr 21, 2022, 5:34:49 PM (3 years ago)

Author:

aslmd

Message:

Fixed floating overflow that could occur in setspi.F90

File:

• trunk/LMDZ.GENERIC/libf/phystd/setspi.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/setspi.F90

@@ -158 +158 @@
          bpa = (b+a)/2.0D0
          bma = (b-a)/2.0D0
+         ! if (nw .eq. 25) then !LT debug
+         !    print*, "a = ",a 
+         !    print*, "b= ",b 
+         !    print*,"bpa = ",bpa 
+         !    print*, "bma = ",bma
+         ! endif
          do nt=NTstart,NTstop
             T   = dble(NT)/NTfac
             ans = 0.0D0
-
             do mm=1,12
                y    = bma*x(mm)+bpa
-               ans  = ans + w(mm)*c1/(y**5*(exp(c2/(y*T))-1.0D0))
+               !to avoid floating overflow when T is low and optical wavelength
+               if ((c2/(y*T)) .lt. 700.0D0) then 
+                  ans  = ans + w(mm)*c1/(y**5*(exp(c2/(y*T))-1.0D0))
+               else 
+                  ans = ans +0.0D0
+               endif
             end do
-
             planckir(NW,nt-NTstart+1) = ans*bma/(PI*DWNI(NW))
          end do