Changeset 2215


Ignore:
Timestamp:
09/06/21 10:53:04 (3 months ago)
Author:
jderouillat
Message:

Manage null elements computing relative error in XIOS_COUPLING test suite

File:
1 edited

Legend:

Unmodified
Added
Removed
  • XIOS/dev/dev_ym/XIOS_COUPLING/xios_test_suite/TEST_SUITE/step2.py

    r2214 r2215  
    9494                        validated = 1 
    9595                        for var in res.variables: 
    96                             if (not (var.startswith('lon_'))) and (not (var.startswith('lat_'))) and (not (var.startswith('time_'))): 
     96                            if (not (var.startswith('lon_'))) and (not (var.startswith('lat_'))) and (not (var.startswith('time_'))) and (not (var.startswith('atm__'))): 
    9797                                ref_interp = ref.variables[var] 
    9898                                ref_array = ref_interp[:] 
    9999                                res_interp = res.variables[var] 
    100100                                res_array = res_interp[:] 
    101                                 diff = (ref_array-res_array)/ref_array 
    102                                 if ( np.max(np.abs(diff)) >  2*10**-3 ): 
    103                                     print( var,  ", max relative error : ", np.max(diff) ) 
     101                                diff = np.zeros_like( ref_array ) 
     102                                np.divide(ref_array-res_array,ref_array,diff,where=(ref_array[:]>10**-15)) 
     103                                if ( np.max(np.abs(diff)) >  1*10**-3 ): 
     104                                    validated = -1 
     105                                diff = np.zeros_like( ref_array ) 
     106                                np.divide(ref_array-res_array,res_array,diff,where=(ref_array[:]>10**-15)) 
     107                                if ( np.max(np.abs(diff)) >  1*10**-3 ): 
    104108                                    validated = -1 
    105109                        report.write(folder_name+" "+folder_name+"@"+config_name+" "+folder_name+"@"+config_name+"@"+checkfile+" "+str(validated)+"\n") 
Note: See TracChangeset for help on using the changeset viewer.