Changeset 243 for codes/icosagcm/branches

Timestamp:

07/21/14 12:02:38 (10 years ago)

Author:

ymipsl

Message:

calling iniphysics directly from dynamico

YM

Location:

codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src

Files:

• earth_const.f90 (modified) (2 diffs)
• physics_lmdz_generic.f90 (modified) (3 diffs)

codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src/earth_const.f90

@@ -13 +13 @@
   REAL(rstd),SAVE :: scale_height=8000. ! atmospheric scale height (m)
   REAL(rstd),SAVE :: scale_factor=1.
+  REAL(rstd),SAVE :: gas_constant = 8.3144621 
+  REAL(rstd),SAVE :: mu                 ! molar mass of the atmosphere
 
   LOGICAL, SAVE :: boussinesq
@@ -32 +34 @@
     CALL getin("scale_height",scale_height)
     
+    mu=kappa/cpp
     boussinesq=.FALSE.
     CALL getin("boussinesq",boussinesq)  

codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src/physics_lmdz_generic.f90

@@ -42 +42 @@
   USE mpi_mod
   USE mpipara
+  USE disvert_mod
+  
   IMPLICIT NONE
   INTEGER  :: distrib(0:mpi_size-1)
@@ -49 +51 @@
   REAL(rstd),ALLOCATABLE :: lonfi(:)
   REAL(rstd),ALLOCATABLE :: airefi(:)
-   
+
     start_day=0
     day_length=86400
@@ -108 +110 @@
     
     CALL init_gcm_lmdz(nbp_phys,mpi_size,distrib,latfi,lonfi,airefi)
+    CALL iniphysiq(llm,day_length,start_day,itau_physics*dt,     &
+                   latfi,lonfi,airefi,radius,g, gas_constant/mu, cpp,     &
+                   ap, bp)
+    
+    
     
   END SUBROUTINE init_physics