text: #!/usr/bin/env python
  import sys 
  import numpy as np
  from netCDF4 import Dataset
  
  with Dataset(sys.argv[1]) as nc1, Dataset(sys.argv[2]) as nc2:
    # Check if the list of variables are the same
    if nc1.variables.keys()!=nc2.variables.keys():
      print("Variables are different")
      sys.exit(2)
  
    for varname in nc1.variables.keys():
      if varname == "time_iso":
        continue
      # First check if each variable has the same dimension
      if np.shape(nc1[varname][:])!=np.shape(nc2[varname][:]):
        print(varname,"dimension is different")
        sys.exit(2)
      # If dimension is the same, compare data
      else:
        diff = nc2[varname][:]-nc1[varname][:]
  
        if (np.abs(diff)).max() != 0:
          print(varname,"is different")
          sys.exit(2)
  