SUBROUTINE init_inca_geometry(         &
       longitude_lmdz, latitude_lmdz,    &
       boundslon_lmdz, boundslat_lmdz,   &
       dx_lmdz, dy_lmdz, cell_area_lmdz, &
       ind_cell_glo_lmdz) 

    USE INCA_DIM, ONLY : PLON
    USE MOD_GEOMETRY_INCA
    USE MOD_GRID_INCA, ONLY: nvertex
    USE CONST_MOD, ONLY: PI 
    USE PRINT_INCA
  IMPLICIT NONE

    REAL,INTENT(IN) :: longitude_lmdz(PLON)
    REAL,INTENT(IN) :: latitude_lmdz(PLON)
    REAL,INTENT(IN) :: boundslon_lmdz(PLON,nvertex)
    REAL,INTENT(IN) :: boundslat_lmdz(PLON,nvertex)
    REAL,INTENT(IN) :: cell_area_lmdz(PLON)
    INTEGER,OPTIONAL,INTENT(IN) :: ind_cell_glo_lmdz(PLON)
    REAL,OPTIONAL,INTENT(IN) :: dx_lmdz(PLON)
    REAL,OPTIONAL,INTENT(IN) :: dy_lmdz(PLON)

    ALLOCATE(longitude(PLON))
    ALLOCATE(latitude(PLON))
    ALLOCATE(longitude_deg(PLON))
    ALLOCATE(latitude_deg(PLON))
    ALLOCATE(boundslon(PLON,nvertex))
    ALLOCATE(boundslat(PLON,nvertex))
    ALLOCATE(cell_area(PLON))
    IF (PRESENT(ind_cell_glo_lmdz)) ALLOCATE(ind_cell_glo(PLON))
    IF (PRESENT(dx_lmdz)) ALLOCATE(dx(PLON))
    IF (PRESENT(dy_lmdz))ALLOCATE(dy(PLON))



    longitude(:) = longitude_lmdz(:)
    latitude(:) = latitude_lmdz(:)
    longitude_deg(:) = longitude(:)*180./PI
    latitude_deg(:) = latitude(:)*180./PI
    boundslon(:,:) = boundslon_lmdz(:,:)
    boundslat(:,:) = boundslat_lmdz(:,:)
    cell_area(:) = cell_area_lmdz(:)
    IF (PRESENT(ind_cell_glo_lmdz)) ind_cell_glo(:) = ind_cell_glo_lmdz(:)
    IF (PRESENT(dx_lmdz)) dx(:) = dx_lmdz(:)
    IF (PRESENT(dy_lmdz)) dy(:) = dy_lmdz(:)


!    write(lunout,*) ' dans init_inca_geometry longitude = '
!    write(lunout,*) longitude(1:PLON) 
!    write(lunout,*) ' dans init_inca_geometry latitude = '
!    write(lunout,*) latitude(1:PLON) 
!    write(lunout,*) 'PLON = ', PLON




  END SUBROUTINE init_inca_geometry