source: codes/icosagcm/devel/src/icosa_init.f90 @ 992

Last change on this file since 992 was 992, checked in by rpennel, 5 years ago

devel : add reading metrics from file at init (folllowing revs @882, @887, @888)

[ !! needs to rewrite new subroutine in xios_mod using cellset strcuture ]

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