Changeset 3231 for LMDZ6/trunk/libf

Timestamp:

Feb 21, 2018, 5:32:08 PM (7 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.

Location:

LMDZ6/trunk/libf

Files:

• phy_common/regular_lonlat_mod.F90 (modified) (1 diff)
• phylmd/regr_horiz_time_climoz_m.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phy_common/regular_lonlat_mod.F90

@@ -20 +20 @@
   REAL,SAVE,ALLOCATABLE :: lat_reg(:)      ! value of longitude cell (rad)
 
-  REAL,SAVE,ALLOCATABLE :: boundslon_reg(:,:)      ! value of boundaries cell (1=>east, 2=>west)(rad)
+  REAL,SAVE,ALLOCATABLE :: boundslon_reg(:,:)      ! value of boundaries cell (1=>west, 2=>east)(rad)
 
   REAL,SAVE,ALLOCATABLE :: boundslat_reg(:,:)      ! value of longitude cell (1=>north, 2=>south)(rad)

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)