source: LMDZ6/trunk/libf/phylmd/ecrad/practical/compare_output_scalar.py @ 5322

#!/usr/bin/env python3

def warn(*args, **kwargs):
    pass
    
import os, warnings
warnings.warn = warn

from ecradplot import plot as eplt

def main(input_srcfile, reference_output_srcfile, output_srcfiles, dstdir):
    """
    Plot input files
    """
    
    import os
    if not os.path.isdir(dstdir):
        os.makedirs(dstdir)

    import seaborn as sns
    name_string  = os.path.splitext(os.path.basename(input_srcfile))[0]
    outputs_string = "_".join([os.path.splitext(os.path.basename(f))[0] for f in output_srcfiles])
    reference_string  = os.path.splitext(os.path.basename(reference_output_srcfile))[0]
            
    styles = [{'lw':2, 'color':'k', 'ls':'-', 'zorder':10},
              {'lw':3, 'color':sns.color_palette()[0], 'ls':'-'},
              {'lw':3, 'color':sns.color_palette()[2], 'ls':'-'},
              {'lw':2, 'color':sns.color_palette()[3], 'ls':'--'},
              {'lw':2, 'color':sns.color_palette()[5], 'ls':'-.'},
              {'lw':4, 'color':sns.color_palette()[6], 'ls':'-'},
              {'lw':4, 'color':sns.color_palette()[9], 'ls':'-'}]
        
    reference_name_string = os.path.splitext(os.path.basename(reference_output_srcfile))[0]
    dstfile = f"{dstdir}/{name_string}_{outputs_string}_vs_{reference_name_string}_scalar.png"
    print(f"Plotting integrated, surface and TOA outputs to {dstfile}")
    eplt.compare_output_scalar(input_srcfile, output_srcfiles, reference_output_srcfile, styles, dstfile=dstfile)

    
if __name__ == "__main__":
    import argparse
    parser = argparse.ArgumentParser(description="Plot radiative fluxes and heating rates from ecRAD output file.")
    parser.add_argument("input",    help="ecRAD input file")
    parser.add_argument("reference", help="ecRAD output file to use as a reference")
    parser.add_argument("outputs",  help="ecRAD output files", nargs='+')
    parser.add_argument("--dstdir", help="Destination directory for plots", default="./")
    args = parser.parse_args()
        
    main(args.input,  args.reference, args.outputs, args.dstdir)