Changeset 3097 for LMDZ6

Timestamp:

Dec 2, 2017, 6:19:58 PM (7 years ago)

Author:

oboucher

Message:

minor change to diag to avoid copy of arrays

File:

• LMDZ6/trunk/libf/phylmd/StratAer/interp_sulf_input.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/StratAer/interp_sulf_input.F90

@@ -32 +32 @@
   INTEGER, SAVE :: n_lev   ! number of levels in the input data
   INTEGER n_mth   ! number of months in the input data
-  REAL OCS_tmp, SO2_tmp
   INTEGER, SAVE :: mth_pre
 !$OMP THREADPRIVATE(mth_pre)
@@ -238 +237 @@
     DO k=1, klev
       !
-      OCS_tmp=tr_seri(i,k,id_OCS_strat)
-      SO2_tmp=tr_seri(i,k,id_SO2_strat)
-      !--OCS and SO2 prescribed below p_bound
+      !--OCS and SO2 prescribed back to their clim values below p_bound
       IF (paprs(i,k).GT.p_bound) THEN
+        budg_3D_backgr_ocs(i,k)=OCS_clim(i,k)-tr_seri(i,k,id_OCS_strat)
+        budg_3D_backgr_so2(i,k)=SO2_clim(i,k)-tr_seri(i,k,id_SO2_strat)
         tr_seri(i,k,id_OCS_strat)=OCS_clim(i,k)
         tr_seri(i,k,id_SO2_strat)=SO2_clim(i,k)
       ENDIF
-      budg_3D_backgr_ocs(i,k)=tr_seri(i,k,id_OCS_strat)-OCS_tmp
-      budg_3D_backgr_so2(i,k)=tr_seri(i,k,id_SO2_strat)-SO2_tmp
     ENDDO
   ENDDO