code: def cal_model_performance(obsl, siml):
      
      nse = MathClass.nashcoef(obsl, siml)
      r2 = MathClass.rsquare(obsl, siml)
      rmse = MathClass.rmse(obsl, siml)
      pbias = MathClass.pbias(obsl, siml)
      rsr = MathClass.rsr(obsl, siml)
      print('NSE: %.2f, R-square: %.2f, PBIAS: %.2f%%, RMSE: %.2f, RSR: %.2f' %
            (nse, r2, pbias, rmse, rsr))
docs: Calculate model performance indexes.