Ignore:
Timestamp:
Feb 2, 2016, 10:35:57 AM (9 years ago)
Author:
ymipsl
Message:

Add global index array in order to be able of correctly restart when changing data distribution between 2 jobs (not same number of mpi process or openmp threads)

YM

File:
1 edited

Legend:

Unmodified
Added
Removed

  • dynamico_lmdz/aquaplanet/LMDZ5/libf/dynlonlat_phylonlat/phylmd/iniphysiq.F90

    r3838 r3901  
    8787  REAL,ALLOCATABLE,SAVE :: boundslonfi(:,:)
    8888  REAL,ALLOCATABLE,SAVE :: boundslatfi(:,:)
    89 !$OMP THREADPRIVATE (latfi,lonfi,cufi,cvfi,airefi,boundslonfi,boundslatfi)
     89  INTEGER,ALLOCATABLE,SAVE :: ind_cell_glo_fi(:)
     90!$OMP THREADPRIVATE (latfi,lonfi,cufi,cvfi,airefi,boundslonfi,boundslatfi,ind_cell_glo_fi)
    9091 
    9192
    … …  
    211212  ALLOCATE(boundslonfi(klon_omp,4))
    212213  ALLOCATE(boundslatfi(klon_omp,4))
     214  ALLOCATE(ind_cell_glo_fi(klon_omp))
    213215
    214216
    … …  
    221223  boundslonfi(1:klon_omp,:) = boundslonfi_glo(offset+klon_omp_begin:offset+klon_omp_end,:)
    222224  boundslatfi(1:klon_omp,:) = boundslatfi_glo(offset+klon_omp_begin:offset+klon_omp_end,:)
    223 
     225  ind_cell_glo_fi(1:klon_omp)=(/ (i,i=offset+klon_omp_begin,offset+klon_omp_end) /)
     226 
    224227  ! copy over global grid longitudes and latitudes
    225   CALL init_geometry(lonfi, latfi, boundslonfi, boundslatfi, airefi, cufi, cvfi)
     228  CALL init_geometry(lonfi, latfi, boundslonfi, boundslatfi, airefi, ind_cell_glo_fi, cufi, cvfi)
    226229
    227230 
Note: See TracChangeset for help on using the changeset viewer.