Changeset 207

Timestamp:

07/10/14 11:14:42 (10 years ago)

Author:

ymipsl

Message:

Add new vertical discretisation by reading ap bp in an input file (ehouarn subroutines)

YM

Location:

codes/icosagcm/trunk/src

Files:

• disvert.f90 (modified) (2 diffs)
• disvert_apbp.f90 (added)
• earth_const.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/disvert.f90

@@ -23 +23 @@
   SUBROUTINE init_disvert
   USE disvert_std_mod, ONLY: ap_std=>ap, bp_std=>bp, presnivs_std=>presnivs, init_disvert_std=>init_disvert
+  USE disvert_apbp_mod, ONLY: ap_apbp=>ap, bp_apbp=>bp, presnivs_apbp=>presnivs, init_disvert_apbp=>init_disvert
   USE disvert_ncar_mod, ONLY: ap_ncar=>ap, bp_ncar=>bp, presnivs_ncar=>presnivs, init_disvert_ncar=>init_disvert
   USE disvert_ncarl30_mod, ONLY: ap_ncarl30=>ap, bp_ncarl30=>bp, presnivs_ncarl30=>presnivs, init_disvert_ncarl30=>init_disvert
@@ -64 +65 @@
         presnivs=>presnivs_std
       
+      CASE('read_apbp')
+    
+        CALL init_disvert_apbp
+        ap=>ap_apbp
+        bp=>bp_apbp
+        presnivs=>presnivs_apbp
+        
       CASE ('ncar')
 

codes/icosagcm/trunk/src/earth_const.f90

@@ -11 +11 @@
   REAL(rstd),SAVE :: preff=101325.
   REAL(rstd),SAVE :: pa=50000.
+  REAL(rstd),SAVE :: scaleheight=8. ! atmospheric scale height (km)
   REAL(rstd),SAVE :: scale_factor=1.
 
@@ -29 +30 @@
     CALL getin("cpp",cpp)  
     CALL getin("preff",preff)  
+    CALL getin("scaleheight",scaleheight)
     
     boussinesq=.FALSE.