Changeset 273

Timestamp:

09/18/14 14:54:39 (10 years ago)

Author:

millour

Message:

Bug correction in the initialization of specific gas constant in Saturn physics.
EM

Location:

codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd

Files:

• inifis.F (modified) (3 diffs)
• tabfi.F (modified) (2 diffs)

codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/inifis.F

@@ -88 +88 @@
       r=pr
       rcp=r/cpp
+      ! Ehouarn debug:
+      write(*,*) "inifis: rad=",rad
+      write(*,*) "        daysec=",daysec
+      write(*,*) "        dtphys=",dtphys
+      write(*,*) "        cpp=",cpp
+      write(*,*) "        g=",g
+      write(*,*) "        r=",r
+      write(*,*) "        rcp=",rcp
 
       avocado = 6.02214179e23   ! added by RW
@@ -645 +653 @@
                STOP
            ELSE
-               write(*,*) "mugaz=",mugaz
+               write(*,*) "inifis: mugaz=",mugaz
            ENDIF
+           !Ehouarn: once mugaz has been set, r, the specific
+           ! gas constant must be computed
+           r=8.3144622/(mugaz*1.e-3)
+           write(*,*) "inifis: r=",r
+           
            cpp = -99999.
            PRINT *,'SPECIFIC HEAT CAPACITY in J K-1 kg-1 ?'
@@ -654 +667 @@
                STOP
            ELSE
-               write(*,*) "cpp=",cpp
+               write(*,*) "inifis: cpp=",cpp
            ENDIF
 !         else

codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/tabfi.F

@@ -104 +104 @@
         g=10.44
         call getin_p("g",g)
-        !mugaz=2.34 !EM: does not give cpp=11500
+        !mugaz=2.34 !EM: does not give cpp=11500)
         mugaz=2.53 ! with this value of mugaz, cpp=11500
         call getin_p("mugaz",mugaz)
+        ! recompute r=R/mugaz to be consitent
+        r=8.3144621/(mugaz*1.e-3)
         ! kappa
         rcp=0.2857143
         call getin_p("kappa",rcp)
-        cpp=(8.314511/(mugaz/1000.0))/rcp
+        cpp=(8.3144621/(mugaz/1000.0))/rcp
         call getin_p("cpp",cpp)
 !        write(*,*) "tabfi: cpp=",cpp
@@ -122 +124 @@
         ! Orbital parameters
         periastr=9.02151966094971
+        call getin_p("periastron",periastr)
         apoastr=10.054479598999
+        call getin_p("apoastron",apoastr)
         peri_day=19280
+        call getin_p("periastron_day",peri_day)
         obliquit=26.7299995422363
-        ! Other parameters some physical paréametrizations need
+        call getin_p("obliquity",obliquit)
+        ! Other parameters some physical parametrizations need
         z0=1e-2
         lmixmin=30