def plot_sed_floc_toy(pdf):
    his_name='ocean_his.nc'
    testcase='Sed_floc_toy'
    path_his_old, path_his_new=os.path.join(util_plot.get_result_path())
    path_his_old=path_his_old+'/Projects/'+testcase+'/'+his_name
    path_his_new=path_his_new+'/Projects/'+testcase+'/'+his_name

    if os.path.isfile(path_his_old) and os.path.isfile(path_his_new):
        print "------------------------------------------------------------"
        print "The history file from OLDER COAWST regression run for",testcase,"case is",path_his_old
        print "------------------------------------------------------------"
        print "The history file from NEW COAWST regression run for",testcase,"case is",path_his_new
        """ set the contour levels for zeta plots """
        colormin=-0.01
        colormax=0.01
        xmin=0.0
        xmax=30.0
        ymin=-0.01
        ymax=0.01
        """ call the function for zeta (water level) plots """
        util_plot.plot_zeta(pdf,path_his_old,path_his_new,his_name,testcase,xmin,xmax,\
                         ymin,ymax,colormin,colormax)
    else:  
        print "-------------------xx---------------------------------------"
        print "One or both HISTORY files for", testcase, "don't exist in the path"
Ejemplo n.º 2
0
def plot_sed_floc_toy(pdf):
    his_name = 'ocean_his.nc'
    testcase = 'Sed_floc_toy'
    path_his_old, path_his_new = os.path.join(util_plot.get_result_path())
    path_his_old = path_his_old + '/Projects/' + testcase + '/' + his_name
    path_his_new = path_his_new + '/Projects/' + testcase + '/' + his_name

    if os.path.isfile(path_his_old) and os.path.isfile(path_his_new):
        print "------------------------------------------------------------"
        print "The history file from OLDER COAWST regression run for", testcase, "case is", path_his_old
        print "------------------------------------------------------------"
        print "The history file from NEW COAWST regression run for", testcase, "case is", path_his_new
        """ set the contour levels for zeta plots """
        colormin = -0.01
        colormax = 0.01
        xmin = 0.0
        xmax = 30.0
        ymin = -0.01
        ymax = 0.01
        """ call the function for zeta (water level) plots """
        util_plot.plot_zeta(pdf,path_his_old,path_his_new,his_name,testcase,xmin,xmax,\
                         ymin,ymax,colormin,colormax)
    else:
        print "-------------------xx---------------------------------------"
        print "One or both HISTORY files for", testcase, "don't exist in the path"
def plot_speclight(pdf):
    his_name = 'ocean_his.nc'
    testcase = 'Spec_light'
    path_his_old, path_his_new = os.path.join(util_plot.get_result_path())
    path_his_old = path_his_old + '/Projects/' + testcase + '/' + his_name
    print "------------------------------------------------------------"
    print "The history file from older COAWST regression run for", testcase, "case is", path_his_old
    path_his_new = path_his_new + '/Projects/' + testcase + '/' + his_name
    print "------------------------------------------------------------"
    print "The history file from new COAWST regression run for", testcase, "case is", path_his_new
    print "------------------------------------------------------------"
    """ set the contour levels for zeta plots """
    colormin = -0.545
    colormax = -0.52
    xmin = -73.5
    xmax = -72.9
    ymin = -39.5
    ymax = -38.5
    """ call the function for zeta (water level) plots """
    util_plot.plot_zeta(pdf,path_his_old,path_his_new,his_name,testcase,xmin,xmax,\
                         ymin,ymax,colormin,colormax)
def plot_speclight(pdf):
    his_name='ocean_his.nc'
    testcase='Spec_light'
    path_his_old, path_his_new=os.path.join(util_plot.get_result_path())
    path_his_old=path_his_old+'/Projects/'+testcase+'/'+his_name
    print "------------------------------------------------------------"
    print "The history file from older COAWST regression run for",testcase,"case is",path_his_old
    path_his_new=path_his_new+'/Projects/'+testcase+'/'+his_name
    print "------------------------------------------------------------"
    print "The history file from new COAWST regression run for",testcase,"case is",path_his_new
    print "------------------------------------------------------------"
      
    """ set the contour levels for zeta plots """
    colormin=-0.545
    colormax=-0.52
    xmin=-73.5
    xmax=-72.9
    ymin=-39.5
    ymax=-38.5
    """ call the function for zeta (water level) plots """
    util_plot.plot_zeta(pdf,path_his_old,path_his_new,his_name,testcase,xmin,xmax,\
                         ymin,ymax,colormin,colormax)