Changeset 366 for codes/icosagcm/trunk/src/earth_const.f90

Timestamp:

10/30/15 15:41:06 (9 years ago)

Author:

dubos

Message:

Progress towards NH

File:

• codes/icosagcm/trunk/src/earth_const.f90 (modified) (3 diffs)

codes/icosagcm/trunk/src/earth_const.f90

@@ -16 +16 @@
   REAL(rstd),SAVE :: mu                 ! molar mass of the atmosphere
 
-  LOGICAL, SAVE :: boussinesq
+  LOGICAL, SAVE :: boussinesq, hydrostatic
 
 CONTAINS
@@ -34 +34 @@
     CALL getin("scale_height",scale_height)
     
-    mu=kappa/cpp
     boussinesq=.FALSE.
     CALL getin("boussinesq",boussinesq)  
+    hydrostatic=.TRUE.
+    CALL getin("hydrostatic",hydrostatic)  
+    IF(boussinesq .AND. .NOT. hydrostatic) THEN
+       PRINT *, 'boussinesq=.TRUE. and hydrostatic=.FALSE. : Non-hydrostatic boussinesq equations are not supported'
+       STOP
+    END IF
     
+    mu=kappa/cpp
     radius=radius/scale_factor
     omega=omega*scale_factor
@@ -45 +51 @@
   
 END MODULE earth_const
-