Changeset 149 for codes/icosagcm/trunk/src/etat0.f90

Timestamp:

04/03/13 12:05:12 (11 years ago)

Author:

sdubey

Message:

Added few new routines to read NC files and compute diagnostics to r145.

Few routines of dry physics including radiation module, surface process and convective adjustment in new routine phyparam.f90. dynetat to read start files for dynamics.

check_conserve routine to compute conservation of quatities like mass, energy etc.etat0_heldsz.f90 for held-suarez test case initial conditions.

new Key time_style=lmd or dcmip to use day_step, ndays like in LMDZ

File:

• codes/icosagcm/trunk/src/etat0.f90 (modified) (5 diffs)

codes/icosagcm/trunk/src/etat0.f90

@@ -1 +1 @@
 MODULE etat0_mod
+    CHARACTER(len=255),SAVE :: etat0_type
 
 CONTAINS
@@ -12 +13 @@
     USE etat0_dcmip4_mod, ONLY : etat0_dcmip4=>etat0  
     USE etat0_dcmip5_mod, ONLY : etat0_dcmip5=>etat0  
+    USE etat0_heldsz_mod, ONLY : etat0_heldsz=>etat0  
+    USE dynetat0_gcm_mod, ONLY : dynetat0_start=>etat0 
+    USE dynetat0_hz_mod,  ONLY : dynetat0_hz=>etat0 
+
     IMPLICIT NONE
     TYPE(t_field),POINTER :: f_ps(:)
@@ -19 +24 @@
     TYPE(t_field),POINTER :: f_q(:)
     
-    CHARACTER(len=255) :: etat0_type
     etat0_type='jablonowsky06'
     CALL getin("etat0",etat0_type)
@@ -28 +32 @@
     CASE ('academic')
        CALL etat0_academic(f_ps,f_phis,f_theta_rhodz,f_u, f_q)
+    CASE ('heldsz')
+        print*,"heldsz test case"
+       CALL etat0_heldsz(f_ps,f_phis,f_theta_rhodz,f_u, f_q)
     CASE ('dcmip1')
        CALL etat0_dcmip1(f_ps,f_phis,f_theta_rhodz,f_u, f_q)
@@ -38 +45 @@
      CASE ('dcmip5')
        CALL etat0_dcmip5(f_ps,f_phis,f_theta_rhodz,f_u, f_q)
+     CASE ('readnf_start') 
+          print*,"readnf_start used"    
+       CALL dynetat0_start(f_ps,f_phis,f_theta_rhodz,f_u,f_q) 
+        CASE ('readnf_hz') 
+          print*,"readnf_hz used"
+       CALL dynetat0_hz(f_ps,f_phis,f_theta_rhodz,f_u,f_q) 
    CASE DEFAULT
        PRINT*, 'Bad selector for variable etat0 <',etat0_type, &