Changeset 649 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

May 4, 2012, 5:14:30 PM (13 years ago)

Author:

jleconte

Message:

• Correction a huge bug in newstart: rcp and cpp can now be changed in start.nc files and are the same as in startfi.nc;

Even when starting from start and startfi files.

• rcp, cpp and mugaz can now be computed using gases.def in newstart

• Correction of a bug arising in gcm.F when the solar days are long (thanks Melanie V.)

File:

• trunk/LMDZ.GENERIC/libf/phystd/tabfi.F (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/phystd/tabfi.F

@@ -79 +79 @@
       INTEGER size
       CHARACTER modif*20
-
 
 c-----------------------------------------------------------------------
@@ -229 +228 @@
       write(*,*) '(31 et 32) iceradius : mean scat radius of CO2 snow'
       write(*,*) '(33 et 34) dtemisice : time scale for snow',
-     &           ' metamorphism'
+     &           '                       metamorphism'
       write(*,*) '(27)        tauvis : mean dust vis. reference ',
-     &           'opacity'
-      write(*,*) '(35)      volcapa : soil volumetric heat capacity'
+     &           '                     opacity'
+      write(*,*) '(35)     volcapa  : soil volumetric heat capacity'
       write(*,*) '(18)     obliquit : planet obliquity (deg)'
       write(*,*) '(17)     peri_day : periastron date (sols since Ls=0)'
@@ -238 +237 @@
       write(*,*) '(16)     apoastr  : max. star-planet (UA)'
       write(*,*) '(14)     year_day : length of year (in sols)'
-      write(*,*) '(5) rad : radius of the planet (m)'
-      write(*,*) '(6) omeg : planet rotation rate (rad/s)'
-      write(*,*) '(7) g : gravity (m/s2)'
-      write(*,*) '(8) mugaz : molecular mass of the atmosphere (g/mol)'
-      write(*,*) '(9) rcp : r/Cp'
-      write(*,*) '(10) daysec : length of a sol (s)'
+      write(*,*) '(5)      rad      : radius of the planet (m)'
+      write(*,*) '(6)      omeg     : planet rotation rate (rad/s)'
+      write(*,*) '(7)      g        : gravity (m/s2)'
+      write(*,*) '(8)      mugaz    : molecular mass '
+      write(*,*) '                       of the atmosphere (g/mol)'
+      write(*,*) '(9)      rcp      : r/Cp'
+      write(*,*) '(8)+(9)  calc_cpp_mugaz : r/Cp and mugaz '
+      write(*,*) '                 computed from gases.def'
+      write(*,*) '(10)     daysec   : length of a sol (s)'
       write(*,*)
  
@@ -460 +462 @@
           write(*,*) ' cpp (new value):',cpp
 
+        else if (modif(1:len_trim(modif)).eq.'calc_cpp_mugaz') then
+          write(*,*) 'current value rcp, mugaz:',rcp,mugaz
+          check_cpp_match=.false.
+          force_cpp=.false.
+          call su_gases
+          call calc_cpp_mugaz
+          write(*,*) 
+          write(*,*) ' cpp (new value):',cpp
+          write(*,*) ' mugaz (new value):',mugaz
+          r=8.314511/(mugaz/1000.0)
+          rcp=r/cpp
+          write(*,*) ' rcp (new value):',rcp
+          
         else if (modif(1:len_trim(modif)).eq.'daysec') then
           write(*,*) 'current value:',daysec