Changeset 535

Timestamp:

Feb 17, 2012, 11:34:42 PM (13 years ago)

Author:

rwordsworth

Message:

fixed kappa bug

File:

• trunk/LMDZ.GENERIC/libf/dyn3d/newstart.F (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/dyn3d/newstart.F

@@ -176 +176 @@
 !      CHARACTER*80 :: visu_file
 
-      cpp    = 744.499 ! for Mars, instead of 1004.70885 (Earth)
-      preff  = 610.    ! for Mars, instead of 101325. (Earth)
-      pa     = 20      ! for Mars, instead of 500 (Earth)
+!      cpp    = 744.499 ! for Mars, instead of 1004.70885 (Earth)
+!      preff  = 610.    ! for Mars, instead of 101325. (Earth)
+!      pa     = 20      ! for Mars, instead of 500 (Earth)
+      cpp    = 0.
+      preff  = 0.
+      pa     = 0. ! to ensure disaster rather than minor error if we don`t
+                  ! make deliberate choice of these values elsewhere.
 
 c=======================================================================
@@ -343 +347 @@
       daysec = p_daysec
 
-!      print*,'daysec=',p_daysec
-!      stop
-
-!      write(*,*) 'aire',aire
+      kappa = 8.314*1000./(p_mugaz*p_cpp) ! added by RDW
 
 
@@ -1726 +1727 @@
       end if
 
-
 c=======================================================================
 c=======================================================================