source: codes/icosagcm/trunk/src/icosa_init.f90 @ 713

Last change on this file since 713 was 710, checked in by ymipsl, 6 years ago

Remove conflict with OpenMP for LMDZ physic initialization.

At this stage, LDZ physic initialization must be called outside a OpenMP parallel section.

YM

File size: 2.4 KB
Line 
1MODULE icosa_init_mod
2
3CONTAINS
4
5  SUBROUTINE icosa_init
6  USE icosa
7  USE timeloop_gcm_mod
8  USE disvert_mod
9  USE etat0_mod
10  USE wind_mod
11  USE mpipara
12  USE omp_para
13  USE vertical_interp_mod
14  USE trace
15  USE output_field_mod
16  USE xios_mod
17  USE write_field_mod
18  USE physics_mod
19  USE tracer_mod
20  USE restart_mod
21  USE etat0_mod
22  USE diagflux_mod
23  USE profiling_mod
24  IMPLICIT NONE
25 
26    CALL init_profiling
27    CALL init_mpipara
28    CALL trace_off
29    CALL xios_init
30    CALL init_earth_const 
31    CALL init_grid_param(is_mpi_master)
32    CALL init_omp_para(is_mpi_master)
33    CALL compute_metric
34    CALL compute_domain
35    CALL init_transfert
36    CALL init_writefield
37    CALL init_trace
38
39  !$OMP PARALLEL 
40    CALL switch_omp_no_distrib_level
41    CALL compute_geometry
42    CALL check_total_area
43 
44    CALL init_disvert 
45    CALL init_vertical_interp
46    CALL init_tracer
47   
48    CALL WriteField("Ai",geom%Ai)
49 
50    IF (is_mpi_root) CALL write_apbp
51    CALL init_time
52    CALL init_restart
53    CALL init_etat0
54
55    CALL output_field_init
56    CALL init_timeloop
57
58  !$OMP END PARALLEL
59
60    CALL init_physics
61
62  !$OMP PARALLEL
63 
64    CALL init_diagflux
65    CALL zero_du_phys
66    CALL timeloop
67    CALL switch_omp_no_distrib_level
68  !$OMP END PARALLEL
69
70    CALL output_field_finalize
71    CALL close_files
72    CALL close_time_counter
73    CALL finalize_mpipara
74
75  CONTAINS
76 
77    SUBROUTINE check_total_area
78    IMPLICIT NONE
79    TYPE(t_field),POINTER,SAVE :: sum_ne(:)
80    TYPE(t_field),POINTER,SAVE :: sum_ne_glo(:)
81    REAL(rstd),POINTER :: pt_sum_ne(:)
82    INTEGER :: ind,i,j,k,n
83    REAL(rstd) :: vect(3,6)
84    REAL(rstd) :: centr(3),dist   
85    REAL(rstd) :: tot_sum=0
86
87      CALL allocate_field(sum_ne,field_T,type_real)
88   
89    !$OMP BARRIER
90    !$OMP MASTER
91      DO ind=1,ndomain
92       
93        pt_sum_ne=sum_ne(ind)
94       
95        CALL swap_dimensions(ind)
96        CALL swap_geometry(ind)
97       
98        DO j=jj_begin,jj_end
99          DO i=ii_begin,ii_end
100            n=(j-1)*iim+i
101            pt_sum_ne(n)=0
102            DO k=1,6
103              pt_sum_ne(n)=pt_sum_ne(n)+ne(n,k)
104            ENDDO
105            IF (domain(ind)%own(i,j)) tot_sum=tot_sum+Ai(n)
106          ENDDO
107        ENDDO
108      ENDDO
109   
110      IF (is_mpi_root) PRINT *," Diff surf",1-tot_sum/(4*Pi*radius*radius)
111    !$OMP END MASTER
112   
113    END SUBROUTINE check_total_area
114 
115  END SUBROUTINE icosa_init
116
117END MODULE icosa_init_mod
Note: See TracBrowser for help on using the repository browser.