Changeset 3278

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.

Location:

LMDZ6/trunk

Files:

• DefLists/config.def_LMDZ5_AGCM (modified) (1 diff)
• libf/phylmd/regr_horiz_time_climoz_m.F90 (modified) (1 diff)

LMDZ6/trunk/DefLists/config.def_LMDZ5_AGCM

@@ -94 +94 @@
 # 1: read a single ozone climatology that will be used day and night
 # 2: read two ozone climatologies, the average day and night climatology and the daylight climatology
-read_climoz=2
+read_climoz=1
 #
 # Parametres simulateur COSP (CFMIP Observational Simulator Package)

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)