source: codes/icosagcm/trunk/src/physics.f90 @ 187

Last change on this file since 187 was 186, checked in by ymipsl, 10 years ago

Add new openMP parallelism based on distribution of domains on threads. There is no more limitation of number of threads by MPI process.

YM
File size: 2.3 KB
Line 
1MODULE physics_mod
2
3  CHARACTER(LEN=255) :: physics_type="automatic"
4!$OMP THREADPRIVATE(physics_type)
5
6CONTAINS
7
8  SUBROUTINE init_physics
9    USE mpipara
10    USE etat0_mod
11    USE icosa
12    USE physics_dcmip_mod,init_physics_dcmip=>init_physics
13!    USE physics_dry_mod
14    IMPLICIT NONE
15
16    physics_type='automatic'
17    CALL getin("physics",physics_type)
18
19    SELECT CASE(TRIM(physics_type))
20    CASE ('automatic')
21       etat0_type='jablonowsky06'
22       CALL getin("etat0",etat0_type)
23       SELECT CASE(TRIM(etat0_type))
24       CASE('held_suarez')
25       CASE DEFAULT
26          IF(is_mpi_root) PRINT*,"NO PHYSICAL PACKAGE USED"
27       END SELECT
28
29
30    CASE ('dcmip')
31       CALL init_physics_dcmip
32
33    CASE ('dry')
34!       CALL init_physics_dry
35
36    CASE DEFAULT
37       PRINT*, 'init_physics : Bad selector for variable physics <',TRIM(physics_type), &
38            '> options are <automatic>, <dcmip>, <dry>'
39       STOP
40    END SELECT
41
42  END SUBROUTINE init_physics
43
44  SUBROUTINE physics(it,jD_cur,jH_cur,f_phis, f_ps, f_theta_rhodz, f_ue, f_q)
45    USE icosa
46!    USE physics_dry_mod
47    USE physics_dcmip_mod, physics_dcmip=>physics
48    USE etat0_mod
49    USE etat0_heldsz_mod
50    IMPLICIT NONE
51    INTEGER, INTENT(IN)   :: it
52    REAL(rstd),INTENT(IN)::jD_cur,jH_cur
53    TYPE(t_field),POINTER :: f_phis(:)
54    TYPE(t_field),POINTER :: f_ps(:)
55    TYPE(t_field),POINTER :: f_theta_rhodz(:)
56    TYPE(t_field),POINTER :: f_ue(:)
57    TYPE(t_field),POINTER :: f_q(:)
58    LOGICAL:: firstcall,lastcall
59
60    SELECT CASE(TRIM(physics_type))
61    CASE ('automatic')
62
63       SELECT CASE(TRIM(etat0_type))
64       CASE('held_suarez')
65          !     CALL transfert_request(f_ps,req_i1)
66          !     CALL transfert_request(f_theta_rhodz,req_i1)
67          !     CALL transfert_request(f_ue,req_e1_vect)
68          CALL held_suarez(f_ps,f_theta_rhodz,f_ue) 
69       CASE DEFAULT
70!          PRINT*,"NO PHYSICAL PACAKAGE USED" ! FIXME MPI
71       END SELECT
72
73    CASE ('dcmip')
74       CALL physics_dcmip(it, f_phis, f_ps, f_theta_rhodz, f_ue, f_q)
75
76    CASE ('dry')
77!       CALL physics_dry(it,jD_cur,jH_cur,f_phis, f_ps, f_theta_rhodz, f_ue, f_q)
78
79    CASE DEFAULT
80       PRINT*, 'Bad selector for variable physics <',TRIM(physics_type), &
81            '> options are <automatic>, <dcmip>, <dry>'
82       STOP
83    END SELECT
84
85  END SUBROUTINE physics
86
87END MODULE physics_mod
Note: See TracBrowser for help on using the repository browser.