| 1 | |
|---|
| 2 | |
|---|
| 3 | |
|---|
| 4 | |
|---|
| 5 | |
|---|
| 6 | |
|---|
| 7 | |
|---|
| 8 | |
|---|
| 9 | |
|---|
| 10 | |
|---|
| 11 | |
|---|
| 12 | SUBROUTINE init_inca_geometry( & |
|---|
| 13 | longitude_lmdz, latitude_lmdz, & |
|---|
| 14 | boundslon_lmdz, boundslat_lmdz, & |
|---|
| 15 | dx_lmdz, dy_lmdz, cell_area_lmdz, & |
|---|
| 16 | ind_cell_glo_lmdz) |
|---|
| 17 | |
|---|
| 18 | USE INCA_DIM, ONLY : PLON |
|---|
| 19 | USE MOD_GEOMETRY_INCA |
|---|
| 20 | USE MOD_GRID_INCA, ONLY: nvertex |
|---|
| 21 | USE CONST_MOD, ONLY: PI |
|---|
| 22 | USE PRINT_INCA |
|---|
| 23 | IMPLICIT NONE |
|---|
| 24 | |
|---|
| 25 | REAL,INTENT(IN) :: longitude_lmdz(PLON) |
|---|
| 26 | REAL,INTENT(IN) :: latitude_lmdz(PLON) |
|---|
| 27 | REAL,INTENT(IN) :: boundslon_lmdz(PLON,nvertex) |
|---|
| 28 | REAL,INTENT(IN) :: boundslat_lmdz(PLON,nvertex) |
|---|
| 29 | REAL,INTENT(IN) :: cell_area_lmdz(PLON) |
|---|
| 30 | INTEGER,OPTIONAL,INTENT(IN) :: ind_cell_glo_lmdz(PLON) |
|---|
| 31 | REAL,OPTIONAL,INTENT(IN) :: dx_lmdz(PLON) |
|---|
| 32 | REAL,OPTIONAL,INTENT(IN) :: dy_lmdz(PLON) |
|---|
| 33 | |
|---|
| 34 | ALLOCATE(longitude(PLON)) |
|---|
| 35 | ALLOCATE(latitude(PLON)) |
|---|
| 36 | ALLOCATE(longitude_deg(PLON)) |
|---|
| 37 | ALLOCATE(latitude_deg(PLON)) |
|---|
| 38 | ALLOCATE(boundslon(PLON,nvertex)) |
|---|
| 39 | ALLOCATE(boundslat(PLON,nvertex)) |
|---|
| 40 | ALLOCATE(cell_area(PLON)) |
|---|
| 41 | IF (PRESENT(ind_cell_glo_lmdz)) ALLOCATE(ind_cell_glo(PLON)) |
|---|
| 42 | IF (PRESENT(dx_lmdz)) ALLOCATE(dx(PLON)) |
|---|
| 43 | IF (PRESENT(dy_lmdz))ALLOCATE(dy(PLON)) |
|---|
| 44 | |
|---|
| 45 | |
|---|
| 46 | |
|---|
| 47 | longitude(:) = longitude_lmdz(:) |
|---|
| 48 | latitude(:) = latitude_lmdz(:) |
|---|
| 49 | longitude_deg(:) = longitude(:)*180./PI |
|---|
| 50 | latitude_deg(:) = latitude(:)*180./PI |
|---|
| 51 | boundslon(:,:) = boundslon_lmdz(:,:) |
|---|
| 52 | boundslat(:,:) = boundslat_lmdz(:,:) |
|---|
| 53 | cell_area(:) = cell_area_lmdz(:) |
|---|
| 54 | IF (PRESENT(ind_cell_glo_lmdz)) ind_cell_glo(:) = ind_cell_glo_lmdz(:) |
|---|
| 55 | IF (PRESENT(dx_lmdz)) dx(:) = dx_lmdz(:) |
|---|
| 56 | IF (PRESENT(dy_lmdz)) dy(:) = dy_lmdz(:) |
|---|
| 57 | |
|---|
| 58 | |
|---|
| 59 | ! write(lunout,*) ' dans init_inca_geometry longitude = ' |
|---|
| 60 | ! write(lunout,*) longitude(1:PLON) |
|---|
| 61 | ! write(lunout,*) ' dans init_inca_geometry latitude = ' |
|---|
| 62 | ! write(lunout,*) latitude(1:PLON) |
|---|
| 63 | ! write(lunout,*) 'PLON = ', PLON |
|---|
| 64 | |
|---|
| 65 | |
|---|
| 66 | |
|---|
| 67 | |
|---|
| 68 | END SUBROUTINE init_inca_geometry |
|---|
