Changeset 1123 for trunk/LMDZ.GENERIC/libf

Timestamp:

Dec 6, 2013, 4:48:32 PM (11 years ago)

Author:

jleconte

Message:

correct a copy/paste bug from last commit and gfortran bug in watercommon

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• sugas_corrk.F90 (modified) (1 diff)
• watercommon_h.F90 (modified) (1 diff)

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

@@ -103 +103 @@
       do igas=1,ngas
          read(111,*) gastype(igas)
-         print*,'Gas ',igas,' is ',gastype(igas+      real kappa_IR, kappa_VI, IR_VI_wnlimit
-)
+         print*,'Gas ',igas,' is ',gastype(igas)
       enddo
 

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

@@ -214 +214 @@
          else if (T.lt.Tmin) then
             print*, "careful, T<Tmin in psat water"
-            psat = Pref_solid_liquid*Exp(r3sublimation*(Tmin-T_h2O_ice_liq)/(Tmin-Trefsublimation)) ! min psat  
+         !   psat = Pref_solid_liquid*Exp(r3sublimation*(Tmin-T_h2O_ice_liq)/(Tmin-Trefsublimation)) ! min psat  
+         ! Ehouarn: gfortran says: Error: Result of EXP underflows its kind,
+         !          so set psat to the smallest possible value instead
+            psat=tiny(psat)
          else                 
             psat = Pref_solid_liquid*Exp(r3sublimation*(T-T_h2O_ice_liq)/(T-Trefsublimation)) ! solid / vapour