Changeset 2215 for XIOS

Timestamp:

09/06/21 10:53:04 (3 years ago)

Author:

jderouillat

Message:

Manage null elements computing relative error in XIOS_COUPLING test suite

File:

• XIOS/dev/dev_ym/XIOS_COUPLING/xios_test_suite/TEST_SUITE/step2.py (modified) (1 diff)

XIOS/dev/dev_ym/XIOS_COUPLING/xios_test_suite/TEST_SUITE/step2.py

@@ -94 +94 @@
                         validated = 1
                         for var in res.variables:
-                            if (not (var.startswith('lon_'))) and (not (var.startswith('lat_'))) and (not (var.startswith('time_'))):
+                            if (not (var.startswith('lon_'))) and (not (var.startswith('lat_'))) and (not (var.startswith('time_'))) and (not (var.startswith('atm__'))):
                                 ref_interp = ref.variables[var]
                                 ref_array = ref_interp[:]
                                 res_interp = res.variables[var]
                                 res_array = res_interp[:]
-                                diff = (ref_array-res_array)/ref_array
-                                if ( np.max(np.abs(diff)) >  2*10**-3 ):
-                                    print( var,  ", max relative error : ", np.max(diff) )
+                                diff = np.zeros_like( ref_array )
+                                np.divide(ref_array-res_array,ref_array,diff,where=(ref_array[:]>10**-15))
+                                if ( np.max(np.abs(diff)) >  1*10**-3 ):
+                                    validated = -1
+                                diff = np.zeros_like( ref_array )
+                                np.divide(ref_array-res_array,res_array,diff,where=(ref_array[:]>10**-15))
+                                if ( np.max(np.abs(diff)) >  1*10**-3 ):
                                     validated = -1
                         report.write(folder_name+" "+folder_name+"@"+config_name+" "+folder_name+"@"+config_name+"@"+checkfile+" "+str(validated)+"\n")