Changeset 4966 for ICOSA_LMDZ/src/phylmd

Timestamp:

Jun 7, 2024, 10:34:27 AM (6 months ago)

Author:

yann meurdesoif

Message:

Fix halo problem when computing parameters of legacy gravity wave parametrization.
YM

File:

• ICOSA_LMDZ/src/phylmd/icolmdz_param_gravity_wave.f90 (modified) (1 diff)

ICOSA_LMDZ/src/phylmd/icolmdz_param_gravity_wave.f90

@@ -401 +401 @@
     INTEGER :: i,j,ij
 
-    DO j=jj_begin,jj_end
-      DO i=ii_begin,ii_end
+    DO j=jj_begin-1,jj_end+1
+      DO i=ii_begin-1,ii_end+1
         ij=(j-1)*iim+i
         zstd(ij) = ztz(ij) - zmea(ij)*zmea(ij)
-        IF (zstd(ij)<0) zstd(ij)=0
+        IF (zstd(ij)<=0) zstd(ij)=0
         zstd(ij)=SQRT(zstd(ij))
       ENDDO