Changeset 3097 for LMDZ6/trunk/libf


Ignore:
Timestamp:
Dec 2, 2017, 6:19:58 PM (7 years ago)
Author:
oboucher
Message:

minor change to diag to avoid copy of arrays

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/StratAer/interp_sulf_input.F90

    r2752 r3097  
    3232  INTEGER, SAVE :: n_lev   ! number of levels in the input data
    3333  INTEGER n_mth   ! number of months in the input data
    34   REAL OCS_tmp, SO2_tmp
    3534  INTEGER, SAVE :: mth_pre
    3635!$OMP THREADPRIVATE(mth_pre)
     
    238237    DO k=1, klev
    239238      !
    240       OCS_tmp=tr_seri(i,k,id_OCS_strat)
    241       SO2_tmp=tr_seri(i,k,id_SO2_strat)
    242       !--OCS and SO2 prescribed below p_bound
     239      !--OCS and SO2 prescribed back to their clim values below p_bound
    243240      IF (paprs(i,k).GT.p_bound) THEN
     241        budg_3D_backgr_ocs(i,k)=OCS_clim(i,k)-tr_seri(i,k,id_OCS_strat)
     242        budg_3D_backgr_so2(i,k)=SO2_clim(i,k)-tr_seri(i,k,id_SO2_strat)
    244243        tr_seri(i,k,id_OCS_strat)=OCS_clim(i,k)
    245244        tr_seri(i,k,id_SO2_strat)=SO2_clim(i,k)
    246245      ENDIF
    247       budg_3D_backgr_ocs(i,k)=tr_seri(i,k,id_OCS_strat)-OCS_tmp
    248       budg_3D_backgr_so2(i,k)=tr_seri(i,k,id_SO2_strat)-SO2_tmp
    249246    ENDDO
    250247  ENDDO
Note: See TracChangeset for help on using the changeset viewer.