Changeset 3231 for LMDZ6/trunk/libf/phylmd/regr_horiz_time_climoz_m.F90

Timestamp:

Feb 21, 2018, 5:32:08 PM (6 years ago)

Author:

dcugnet

Message:

Fix in regr_horiz_time_climoz_m: <= instead of < in an index search loop to avoid
a possible error in the 3D ozone interpolation when input and output longitudes
numbers are equal.
Few typos in the comments corrected.

File:

• LMDZ6/trunk/libf/phylmd/regr_horiz_time_climoz_m.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/regr_horiz_time_climoz_m.F90

@@ -192 +192 @@
     v1=CSHIFT(v1,SHIFT=dx1,DIM=1); v1(nlon_in-dx1+1:)=v1(nlon_in-dx1+1:)+2.*pi
 
-    !--- Extend input longitudes vector until last interval contains boundslon_reg(nlat_ou,east)
-    dx2=0; DO WHILE(v1(1+dx2)+2.*pi<boundslon_reg(nlon_ou,east)); dx2=dx2+1; END DO
+    !--- Extend input longitudes vector until last interval contains boundslon_reg(nlon_ou,east)
+    dx2=0; DO WHILE(v1(1+dx2)+2.*pi<=boundslon_reg(nlon_ou,east)); dx2=dx2+1; END DO
 
     !--- Final edges longitudes vector (with margin and end point)