Changeset 4966


Ignore:
Timestamp:
Jun 7, 2024, 10:34:27 AM (7 months ago)
Author:
yann meurdesoif
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • ICOSA_LMDZ/src/phylmd/icolmdz_param_gravity_wave.f90

    r4858 r4966  
    401401    INTEGER :: i,j,ij
    402402
    403     DO j=jj_begin,jj_end
    404       DO i=ii_begin,ii_end
     403    DO j=jj_begin-1,jj_end+1
     404      DO i=ii_begin-1,ii_end+1
    405405        ij=(j-1)*iim+i
    406406        zstd(ij) = ztz(ij) - zmea(ij)*zmea(ij)
    407         IF (zstd(ij)<0) zstd(ij)=0
     407        IF (zstd(ij)<=0) zstd(ij)=0
    408408        zstd(ij)=SQRT(zstd(ij))
    409409      ENDDO
Note: See TracChangeset for help on using the changeset viewer.