Last change
on this file since 2422 was
2351,
checked in by Ehouarn Millour, 9 years ago
|
More on physics/dynamics separation and cleanup:
- Set things up so that all physics-related initializations are done via iniphysiq.
- Created a "geometry_mod.F90" module in phy_common to store information on the loacl grid (i.e. replaces comgeomphy) and moreover give these variables more obvious names (e.g.: rlond => longitude, rlatd => latitude, airephy => cell_area).
- removed obsolete comgeomphy.h and comgeomphy.F90
EM
|
-
Property copyright set to
Name of program: LMDZ Creation date: 1984 Version: LMDZ5 License: CeCILL version 2 Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539 See the license file in the root directory
-
Property svn:eol-style set to
native
-
Property svn:keywords set to
Author Date Id Revision
|
File size:
1.1 KB
|
Line | |
---|
1 | ! |
---|
2 | ! $Id: mod_interface_dyn_phys.F90 2351 2015-08-25 15:14:59Z emillour $ |
---|
3 | ! |
---|
4 | MODULE mod_interface_dyn_phys |
---|
5 | INTEGER,SAVE,dimension(:),allocatable :: index_i |
---|
6 | INTEGER,SAVE,dimension(:),allocatable :: index_j |
---|
7 | |
---|
8 | |
---|
9 | CONTAINS |
---|
10 | |
---|
11 | #ifdef CPP_PARA |
---|
12 | ! Interface with parallel physics, |
---|
13 | SUBROUTINE Init_interface_dyn_phys |
---|
14 | USE mod_phys_lmdz_mpi_data |
---|
15 | IMPLICIT NONE |
---|
16 | include 'dimensions.h' |
---|
17 | |
---|
18 | INTEGER :: i,j,k |
---|
19 | |
---|
20 | ALLOCATE(index_i(klon_mpi)) |
---|
21 | ALLOCATE(index_j(klon_mpi)) |
---|
22 | |
---|
23 | k=1 |
---|
24 | IF (is_north_pole) THEN |
---|
25 | index_i(k)=1 |
---|
26 | index_j(k)=1 |
---|
27 | k=2 |
---|
28 | ELSE |
---|
29 | DO i=ii_begin,iim |
---|
30 | index_i(k)=i |
---|
31 | index_j(k)=jj_begin |
---|
32 | k=k+1 |
---|
33 | ENDDO |
---|
34 | ENDIF |
---|
35 | |
---|
36 | DO j=jj_begin+1,jj_end-1 |
---|
37 | DO i=1,iim |
---|
38 | index_i(k)=i |
---|
39 | index_j(k)=j |
---|
40 | k=k+1 |
---|
41 | ENDDO |
---|
42 | ENDDO |
---|
43 | |
---|
44 | IF (is_south_pole) THEN |
---|
45 | index_i(k)=1 |
---|
46 | index_j(k)=jj_end |
---|
47 | ELSE |
---|
48 | DO i=1,ii_end |
---|
49 | index_i(k)=i |
---|
50 | index_j(k)=jj_end |
---|
51 | k=k+1 |
---|
52 | ENDDO |
---|
53 | ENDIF |
---|
54 | |
---|
55 | END SUBROUTINE Init_interface_dyn_phys |
---|
56 | #else |
---|
57 | SUBROUTINE Init_interface_dyn_phys |
---|
58 | ! dummy routine for seq case |
---|
59 | END SUBROUTINE Init_interface_dyn_phys |
---|
60 | #endif |
---|
61 | ! of #ifdef CPP_PARA |
---|
62 | END MODULE mod_interface_dyn_phys |
---|
Note: See
TracBrowser
for help on using the repository browser.