Changeset 989 for trunk/LMDZ.GENERIC/libf

Timestamp:

Jun 14, 2013, 9:17:44 AM (12 years ago)

Author:

emillour

Message:

Generic GCM:

• Some minor changes so that gcm compiles with gfortran:
  • Added option to compile "long lines" (>132 characters) in makegcm_gfortran
  • Removed use of isnan() in physiq.F90 (it is not a standard function)
  • Avoid possible underflow of psat in watercommon_h.F90
  • Adapted the checks on the *IR and *VI band files to be more strict

EM

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• physiq.F90 (modified) (1 diff)
• setspi.F90 (modified) (3 diffs)
• setspv.F90 (modified) (2 diffs)
• watercommon_h.F90 (modified) (1 diff)

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

@@ -1103 +1103 @@
                   enddo
                   dEtot=cpp*SUM(massarea(:,:)*(dtlscale(:,:)))/totarea
-                  if(isnan(dEtot)) then
-                     print*,'Nan in largescale, abort'
-                     STOP
-                  endif
+!                 if(isnan(dEtot)) then ! NB: isnan() is not a standard function...
+!                    print*,'Nan in largescale, abort'
+!                     STOP
+!                 endif
                   print*,'In largescale atmospheric energy change =',dEtot,' W m-2'
 

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

@@ -33 +33 @@
       logical file_ok
       integer nw, nt, m, mm, file_entries
-      real*8 a, b, ans, y, bpa, bma, T
+      real*8 a, b, ans, y, bpa, bma, T, dummy
 
       character(len=30)  :: temp1
@@ -48 +48 @@
 
       !! used to count lines
-      integer :: nb=-1  !because first line is not an actual value
+      integer :: nb=0
       integer :: ierr=0
 
@@ -94 +94 @@
       ! check that the file contains the right number of bands 
       open(131,file=file_path,form='formatted')
+      read(131,*,iostat=ierr) file_entries
       do while (ierr==0)
-        read(131,*,iostat=ierr) file_entries
+        read(131,*,iostat=ierr) dummy
+!        write(*,*) 'setspi: file_entries:',dummy,'ierr=',ierr
         if (ierr==0) nb=nb+1
       enddo
       close(131)
+
       write(*,*) 'setspi: L_NSPECTI = ',L_NSPECTI, 'in the model '
       write(*,*) '        there are   ',nb, 'entries in ',TRIM(file_path)

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

@@ -44 +44 @@
 
       real*8 STELLAR(L_NSPECTV)
-      real*8 sum
+      real*8 sum, dummy
 
       !! used to count lines
-      integer :: nb=-1  !because first line is not an actual value
+      integer :: nb=0 
       integer :: ierr=0
 
@@ -72 +72 @@
       ! check that the file contains the right number of bands 
       open(131,file=file_path,form='formatted')
+      read(131,*,iostat=ierr) file_entries
       do while (ierr==0)
-        read(131,*,iostat=ierr) file_entries
+        read(131,*,iostat=ierr) dummy
         if (ierr==0) nb=nb+1
       enddo
       close(131)
+
       write(*,*) 'setspv: L_NSPECTV = ',L_NSPECTV, 'in the model '
       write(*,*) '        there are   ',nb, 'entries in ',TRIM(file_path)

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

@@ -61 +61 @@
          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