source: codes/icosagcm/trunk/src/icosa_sw.f90 @ 12

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

dynamico tree creation

YM
File size: 1.1 KB
Line 
1PROGRAM ICOSA_SW
2  USE spherical_geom_mod
3  USE metric
4  USE domain_mod
5  USE field_mod
6  USE write_field
7  USE geometry
8  USE dimensions
9  USE timeloop_sw_mod
10  USE transfert_mod
11  USE dissip_mod
12  USE disvert_mod
13  USE transfert_mod
14  IMPLICIT NONE
15 
16
17  TYPE(t_field),POINTER :: sum_ne(:)
18  REAL(rstd),POINTER :: pt_sum_ne(:)
19 
20  INTEGER :: ind,i,j,k,n
21  REAL(rstd) :: tot_sum=0
22 
23  CALL compute_metric
24  CALL compute_domain
25  CALL compute_geometry
26  CALL init_transfert
27 
28  CALL allocate_field(sum_ne,field_T,type_real)
29 
30
31  DO ind=1,ndomain
32     
33    pt_sum_ne=sum_ne(ind)
34     
35    CALL swap_dimensions(ind)
36    CALL swap_geometry(ind)
37   
38    DO j=jj_begin,jj_end
39      DO i=ii_begin,ii_end
40        n=(j-1)*iim+i
41        pt_sum_ne(n)=0
42        DO k=1,6
43          pt_sum_ne(n)=pt_sum_ne(n)+ne(n,k)
44        ENDDO
45        IF (domain(ind)%own(i,j)) tot_sum=tot_sum+Ai(n)
46      ENDDO
47    ENDDO
48  ENDDO
49
50  PRINT *," Diff surf",1-tot_sum/(4*Pi*radius*radius)
51
52 
53!  CALL WriteField("Ai",geom%Ai)
54!  CALL WriteField("sum_ne",sum_ne)
55
56 CALL timeloop
57 
58END PROGRAM ICOSA_SW 
Note: See TracBrowser for help on using the repository browser.