Ignore:
Timestamp:
Jan 5, 2017, 4:40:31 PM (8 years ago)
Author:
oboucher
Message:

Adding threadprivate for mth_pre variable
This should fix some reproducibility effects for runs longer than 1 months....

Location:
LMDZ5/trunk/libf/phylmd/rrtm
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phylmd/rrtm/readaerosolstrato1_rrtm.F90

    r2694 r2744  
    3535    REAL, POINTER:: lev(:)
    3636    INTEGER k, band, wave, i
    37     INTEGER, SAVE :: mth_pre
     37    INTEGER, SAVE :: mth_pre=1
     38!$OMP THREADPRIVATE(mth_pre)
    3839
    3940    REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: tau_aer_strat
    … …  
    4445    REAL, ALLOCATABLE:: tauaerstrat_mois(:, :, :)
    4546    REAL, ALLOCATABLE:: tauaerstrat_mois_glo(:, :)
    46 
    47     REAL, ALLOCATABLE:: sum_tau_aer_strat(:)
    4847
    4948! For NetCDF:
    … …  
    7776
    7877    IF (.not.ALLOCATED(tau_aer_strat)) ALLOCATE(tau_aer_strat(klon,klev))
    79     IF (.not.ALLOCATED(sum_tau_aer_strat)) ALLOCATE(sum_tau_aer_strat(klon))
    8078
    8179!--we only read monthly strat aerosol data

  • LMDZ5/trunk/libf/phylmd/rrtm/readaerosolstrato2_rrtm.F90

    r2725 r2744  
    4141    REAL, POINTER:: wav(:)
    4242    INTEGER i,k,wave,band
    43     INTEGER, SAVE :: mth_pre
     43    INTEGER, SAVE :: mth_pre=1
     44!$OMP THREADPRIVATE(mth_pre)
    4445
    4546    REAL, ALLOCATABLE, DIMENSION(:,:,:), SAVE :: tau_aer_strat
Note: See TracChangeset for help on using the changeset viewer.