Changeset 1725 for trunk/LMDZ.TITAN/libf

Timestamp:

Jul 21, 2017, 5:30:44 PM (7 years ago)

Author:

jvatant

Message:

Optimization of the optci/cv routines

• The repeated calls to huge matrices gasi/v increased dramatically the execution time because of memory access
• Added a tmpk variable
• Save ~ 50% time on the RT, ~30% on the whole code on the tested simulations

JVO

Location:

trunk/LMDZ.TITAN/libf/phytitan

Files:

• optci.F90 (modified) (2 diffs)
• optcv.F90 (modified) (2 diffs)

trunk/LMDZ.TITAN/libf/phytitan/optci.F90

@@ -74 +74 @@
 
   real*8  KCOEF(4)
+   
+  ! temporary variable to reduce memory access time to gasi
+  real*8 tmpk(2,2)
   
   ! temporary variables for multiple aerosol calculation
@@ -216 +219 @@
            ! Now compute TAUGAS
 
-           KCOEF(1) = GASI(MT(K),MP(K),1,NW,NG)
-           KCOEF(2) = GASI(MT(K),MP(K)+1,1,NW,NG)
-           KCOEF(3) = GASI(MT(K)+1,MP(K)+1,1,NW,NG)
-           KCOEF(4) = GASI(MT(K)+1,MP(K),1,NW,NG)
+           ! JVO 2017 : added tmpk because the repeated calls to gasi/v increased dramatically
+           ! the execution time of optci/v -> ~ factor 2 on the whole radiative
+           ! transfer on the tested simulations !
+
+           tmpk = GASI(MT(K):MT(K)+1,MP(K):MP(K)+1,1,NW,NG)
+
+           KCOEF(1) = tmpk(1,1) ! KCOEF(1) = GASI(MT(K),MP(K),1,NW,NG)
+           KCOEF(2) = tmpk(1,2) ! KCOEF(2) = GASI(MT(K),MP(K)+1,1,NW,NG)
+           KCOEF(3) = tmpk(2,2) ! KCOEF(3) = GASI(MT(K)+1,MP(K)+1,1,NW,NG)
+           KCOEF(4) = tmpk(2,1) ! KCOEF(4) = GASI(MT(K)+1,MP(K),1,NW,NG)
+
 
            ! Interpolate the gaseous k-coefficients to the requested T,P values

trunk/LMDZ.TITAN/libf/phytitan/optcv.F90

@@ -83 +83 @@
 
   real*8  KCOEF(4)
+  
+  ! temporary variable to reduce memory access time to gasv
+  real*8 tmpk(2,2)
 
   ! temporary variables for multiple aerosol calculation
@@ -233 +236 @@
            ! Now compute TAUGAS
 
-           KCOEF(1) = GASV(MT(K),MP(K),1,NW,NG)
-           KCOEF(2) = GASV(MT(K),MP(K)+1,1,NW,NG)
-           KCOEF(3) = GASV(MT(K)+1,MP(K)+1,1,NW,NG)
-           KCOEF(4) = GASV(MT(K)+1,MP(K),1,NW,NG)
+           ! JVO 2017 : added tmpk because the repeated calls to gasi/v increased dramatically
+           ! the execution time of optci/v -> ~ factor 2 on the whole radiative
+           ! transfer on the tested simulations !
+
+           tmpk = GASV(MT(K):MT(K)+1,MP(K):MP(K)+1,1,NW,NG)
+              
+           KCOEF(1) = tmpk(1,1) ! KCOEF(1) = GASV(MT(K),MP(K),1,NW,NG)
+           KCOEF(2) = tmpk(1,2) ! KCOEF(2) = GASV(MT(K),MP(K)+1,1,NW,NG)
+           KCOEF(3) = tmpk(2,2) ! KCOEF(3) = GASV(MT(K)+1,MP(K)+1,1,NW,NG)
+           KCOEF(4) = tmpk(2,1) ! KCOEF(4) = GASV(MT(K)+1,MP(K),1,NW,NG)
 
            ! Interpolate the gaseous k-coefficients to the requested T,P values