source: codes/icosagcm/trunk/src/etat0.f90 @ 46

Last change on this file since 46 was 37, checked in by ymipsl, 12 years ago

Implement DCMIP 3.1 etat0 testcase

YM

File size: 1.2 KB
Line 
1MODULE etat0_mod
2
3
4CONTAINS
5 
6  SUBROUTINE etat0(f_ps,f_phis,f_theta_rhodz,f_u, f_q)
7  USE icosa
8  USE etat0_jablonowsky06_mod, ONLY : etat0_jablonowsky06=>etat0
9  USE etat0_academic_mod, ONLY : etat0_academic=>etat0 
10  USE etat0_ncar_mod, ONLY : etat0_ncar=>etat0 
11  USE etat0_dcmip3_mod, ONLY : etat0_dcmip3=>etat0 
12  IMPLICIT NONE
13    TYPE(t_field),POINTER :: f_ps(:)
14    TYPE(t_field),POINTER :: f_phis(:)
15    TYPE(t_field),POINTER :: f_theta_rhodz(:)
16    TYPE(t_field),POINTER :: f_u(:)
17    TYPE(t_field),POINTER :: f_q(:)
18
19    CHARACTER(len=255) :: etat0_type
20   
21    etat0_type='jablonowsky06'
22    CALL getin("etat0",etat0_type)
23   
24    SELECT CASE (TRIM(etat0_type))
25      CASE ('jablonowsky06')
26        CALL etat0_jablonowsky06(f_ps,f_phis,f_theta_rhodz,f_u, f_q)
27      CASE ('academic')
28        CALL etat0_academic(f_ps,f_phis,f_theta_rhodz,f_u, f_q)
29      CASE ('ncar')
30        CALL etat0_ncar(f_ps,f_phis,f_theta_rhodz,f_u, f_q)
31      CASE ('dcmip3')
32        CALL etat0_dcmip3(f_ps,f_phis,f_theta_rhodz,f_u, f_q)
33      CASE DEFAULT
34        PRINT*, 'Bad selector for varaiable etat0 <',etat0_type,'> options are <jablonowsky06>, <academic>, <ncar> '
35        STOP
36    END SELECT
37
38  END SUBROUTINE etat0 
39         
40END MODULE etat0_mod
Note: See TracBrowser for help on using the repository browser.