Changeset 3278 for LMDZ6/trunk/libf

Timestamp:

Mar 16, 2018, 4:03:21 PM (6 years ago)

Author:

dcugnet

Message:

Fix of the rev. 3231 for ozone horizontal interpolation.
Fields pre-r3231 were one point shifted in longitude in the easter part of the globe due to an index offset.
Difference in small, thanks to the quasi-zonal symetry of the 3D ozone fields.

File:

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

LMDZ6/trunk/libf/phylmd/regr_horiz_time_climoz_m.F90

@@ -190 +190 @@
     !--- Ensure first input longitudes interval contains first output point boundslon_reg(1,west)
     dx1=locate(v1,boundslon_reg(1,west))-1
-    v1=CSHIFT(v1,SHIFT=dx1,DIM=1); v1(nlon_in-dx1+1:)=v1(nlon_in-dx1+1:)+2.*pi
+    v1=CSHIFT(v1,SHIFT=dx1,DIM=1); v1(nlon_in-dx1+2:)=v1(nlon_in-dx1+2:)+2.*pi
 
     !--- Extend input longitudes vector until last interval contains boundslon_reg(nlon_ou,east)