Changeset 1506 for trunk/LMDZ.VENUS/libf/phyvenus/physiq.F

Timestamp:

Dec 24, 2015, 2:10:58 PM (9 years ago)

Author:

emillour

Message:

Venus GCM:
Some minor modifications to compile/run with gfortran:

• removed reference to intrinsic "ieee_arithmetic" module (only available with gfortran v5 and higher), and replaced calls to intrinsic "ieee_is_nan" with equivalent tests in physiq.F
• added a LAPACK preprocessing flag around calls to dgesv (Lapack subroutine) in new_photochemistry_venus.F90 and adapted MESU and OCCIGEN arch files (GNOME and ADA arch files already include this flag).

EM

File:

• trunk/LMDZ.VENUS/libf/phyvenus/physiq.F (modified) (2 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/physiq.F

@@ -69 +69 @@
       USE compo_hedin83_mod2
       use moyzon_mod, only: tmoy
-      use ieee_arithmetic
+!      use ieee_arithmetic
       use temps_mod, only: itau_phy,day_ref,annee_ref
       use logic_mod, only: iflag_trac
@@ -962 +962 @@
 c--------------------
 c   Ce test est necessaire pour eviter Xliq=NaN   
-        IF (ieee_is_nan(d_tr_sed(i,k,1)).OR.
-     &  ieee_is_nan(d_tr_sed(i,k,2))) THEN
+!        IF (ieee_is_nan(d_tr_sed(i,k,1)).OR.
+!     &  ieee_is_nan(d_tr_sed(i,k,2))) THEN
+        IF ((d_tr_sed(i,k,1).ne.d_tr_sed(i,k,1)).OR.
+     &      (d_tr_sed(i,k,2).ne.d_tr_sed(i,k,2))) THEN
         PRINT*,'sedim NaN PROBLEM'
         PRINT*,'d_tr_sed Nan?',d_tr_sed(i,k,:),'Temp',t_seri(i,k)