Changeset 1680 for LMDZ5/trunk/libf/dyn3dpar/filtreg_p.F

Timestamp:

Nov 15, 2012, 2:43:45 PM (12 years ago)

Author:

Ehouarn Millour

Message:

• bug correction in filtreg_p.F on array boundaries (was a problem only in cases where BLAS was not used).
• modified "makelmdz_fcm" script: added '-mem' option to trigger using "dyn3dmem" dynamics instead of (default) "dyn3dpar".

EM

File:

• LMDZ5/trunk/libf/dyn3dpar/filtreg_p.F (modified) (6 diffs)

LMDZ5/trunk/libf/dyn3dpar/filtreg_p.F

@@ -214 +214 @@
      &                    champ_fft(1,j-jdfil+1,1), iip1*nlat)
 #else
-                     champ_fft(:,j-jdfil+1,:)
+                     champ_fft(:iim,j-jdfil+1,:)
      &                    =matmul(matrinvn(:,:,j),champ_loc(:iim,j,:))
 #endif
@@ -227 +227 @@
      &                    champ_fft(1,j-jdfil+1,1), iip1*nlat)
 #else
-                     champ_fft(:,j-jdfil+1,:)
+                     champ_fft(:iim,j-jdfil+1,:)
      &                    =matmul(matriceun(:,:,j),champ_loc(:iim,j,:))
 #endif
@@ -240 +240 @@
      &                    champ_fft(1,j-jdfil+1,1), iip1*nlat)
 #else
-                     champ_fft(:,j-jdfil+1,:)
+                     champ_fft(:iim,j-jdfil+1,:)
      &                    =matmul(matricevn(:,:,j),champ_loc(:iim,j,:))
 #endif
@@ -257 +257 @@
      &                    champ_fft(1,j-jdfil+1,1), iip1*nlat)
 #else
-                     champ_fft(:,j-jdfil+1,:)
+                     champ_fft(:iim,j-jdfil+1,:)
      &                    =matmul(matrinvs(:,:,j-jfiltsu+1),
      &                            champ_loc(:iim,j,:))
@@ -272 +272 @@
      &                    champ_fft(1,j-jdfil+1,1), iip1*nlat)
 #else
-                     champ_fft(:,j-jdfil+1,:)
+                     champ_fft(:iim,j-jdfil+1,:)
      &                    =matmul(matriceus(:,:,j-jfiltsu+1),
      &                            champ_loc(:iim,j,:))
@@ -287 +287 @@
      &                    champ_fft(1,j-jdfil+1,1), iip1*nlat)
 #else
-                     champ_fft(:,j-jdfil+1,:)
+                     champ_fft(:iim,j-jdfil+1,:)
      &                    =matmul(matricevs(:,:,j-jfiltsv+1),
      &                            champ_loc(:iim,j,:))