source: codes/icosagcm/trunk/src/disvert.f90 @ 67

Last change on this file since 67 was 64, checked in by ymipsl, 12 years ago

bug fix for activate ncarl30 vertical discretisation

YM

File size: 1.3 KB
Line 
1MODULE disvert_mod
2  USE icosa
3  REAL(rstd), SAVE, POINTER :: ap(:)
4  REAL(rstd), SAVE, POINTER :: bp(:)
5  REAL(rstd), SAVE, POINTER :: presnivs(:)
6
7CONTAINS
8
9  SUBROUTINE init_disvert
10  USE disvert_std_mod, ONLY: ap_std=>ap, bp_std=>bp, presnivs_std=>presnivs, init_disvert_std=>init_disvert
11  USE disvert_ncar_mod, ONLY: ap_ncar=>ap, bp_ncar=>bp, presnivs_ncar=>presnivs, init_disvert_ncar=>init_disvert
12  USE disvert_ncarl30_mod, ONLY: ap_ncarl30=>ap, bp_ncarl30=>bp, presnivs_ncarl30=>presnivs, init_disvert_ncarl30=>init_disvert
13  USE icosa
14  IMPLICIT NONE
15    CHARACTER(LEN=255) :: disvert_type = 'std'
16   
17    CALL getin("disvert",disvert_type)
18   
19    SELECT CASE (TRIM(disvert_type))
20      CASE('std')
21   
22        CALL init_disvert_std
23        ap=>ap_std
24        bp=>bp_std
25        presnivs=>presnivs_std
26     
27      CASE ('ncar')
28
29        CALL init_disvert_ncar
30        ap=>ap_ncar
31        bp=>bp_ncar
32        presnivs=>presnivs_ncar
33
34      CASE ('ncarl30')
35
36        CALL init_disvert_ncarl30
37        ap=>ap_ncarl30
38        bp=>bp_ncarl30
39        presnivs=>presnivs_ncarl30
40       
41      CASE default
42        PRINT*,'Bad selector for variable disvert : <', TRIM(disvert_type),"> options are <std>, <ncar>, <ncarl30>" 
43        STOP
44       
45    END SELECT
46
47  END SUBROUTINE init_disvert 
48 
49END MODULE disvert_mod
Note: See TracBrowser for help on using the repository browser.