Ignore:
Timestamp:
Jan 8, 2007, 3:36:39 PM (17 years ago)
Author:
sexton
Message:

report updates and fixes to check_list

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/abstract_2d_finite_volumes/util.py

    r4131 r4145  
    833833                   
    834834    import sys
    835     if not sys.version.startswith('2.4'):
     835    #if not sys.version.startswith('2.4'):
    836836        # Backwards compatibility
    837         from sets import Set as set
    838        
    839    
    840     if sys.platform == 'win32':
    841         for i,j in enumerate(quantity):
    842             quantity[i] = quantity[i].lower()
    843         p = list(set(quantity).difference(set(all_quantity)))
    844         if len(p) <> 0:
    845             msg = 'Quantities %s do not exist - please try again' %p
    846             raise Exception, msg
     837    #   from sets import Set as set
     838
     839    from sets import Set as set
     840           
     841    for i,j in enumerate(quantity):
     842        quantity[i] = quantity[i].lower()
     843    p = list(set(quantity).difference(set(all_quantity)))
     844    if len(p) <> 0:
     845        msg = 'Quantities %s do not exist - please try again' %p
     846        raise Exception, msg
    847847       
    848848    return
     
    928928    max_momentums = []
    929929    max_speeds = []
    930     c = 0
     930    max_depths = []
    931931    model_time_plot3d = zeros((n0,m), Float)
    932932    stages_plot3d = zeros((n0,m), Float)
     
    944944            max_stage = 0
    945945            max_momentum = 0
    946             max_speed = 0   
     946            max_speed = 0
     947            max_depth = 0
    947948            gaugeloc = locations[k]
    948949            thisfile = file_loc[j]+sep+'gauges_time_series'+'_'+gaugeloc+'.csv'
     
    983984                    if m > max_momentum: max_momentum = m
    984985                    if vel > max_speed: max_speed = vel
     986                    if z > 0 and depth > max_depth: max_depth = depth
    985987                   
    986988                   
     
    990992            min_stages.append(min_stage)
    991993            max_momentums.append(max_momentum)
    992             max_speeds.append(max_speed)     
     994            max_speeds.append(max_speed)
     995            max_depths.append(max_depth)
    993996            #### finished generating quantities for each swwfile #####
    994997       
     
    10261029        profilefig = 'solution_xprofile'
    10271030        savefig('profilefig')
    1028                
    1029     #stage_axis = axis([time_min/60.0, time_max/60.0, min(min_stages), max(max_stages)*1.1])
    1030     #stage_axis = axis([time_min/60.0, time_max/60.0, -3.0, 3.0])   
    1031     #vel_axis = axis([time_min/60.0, time_max/60.0, min(max_speeds), max(max_speeds)*1.1])
    1032     #mom_axis = axis([time_min/60.0, time_max/60.0, min(max_momentums), max(max_momentums)*1.1]) 
     1031
     1032    depth_axis = axis([time_min/60.0, time_max/60.0, 0, max(max_depths)*1.1])
     1033    stage_axis = axis([time_min/60.0, time_max/60.0, min(min_stages), max(max_stages)*1.1])
     1034    stage_axis = axis([time_min/60.0, time_max/60.0, -3.0, 3.0])   
     1035    vel_axis = axis([time_min/60.0, time_max/60.0, min(max_speeds), max(max_speeds)*1.1])
     1036    mom_axis = axis([time_min/60.0, time_max/60.0, min(max_momentums), max(max_momentums)*1.1]) 
    10331037   
    10341038    cstr = ['g', 'r', 'b', 'c', 'm', 'y', 'k']
     
    10371041    elev_output = []
    10381042    if len(label_id) > 1: graphname_report = []
     1043    pp = 1
     1044    div = 11.
    10391045    for k in gauge_index:
    10401046        g = gauges[k]
    10411047        count1 = 0
    10421048        if report == True and len(label_id) > 1:
    1043             s = '\\begin{figure}[hbt] \n \\centering \n \\begin{tabular}{cc} \n'
     1049            s = '\\begin{figure}[ht] \n \\centering \n \\begin{tabular}{cc} \n'
    10441050            fid.write(s)
    10451051        if len(label_id) > 1: graphname_report = []
     
    10631069                    plot(model_time[0:n[j]-1,k,j], depths[0:n[j]-1,k,j], '-', c = cstr[j])
    10641070                    units = 'm'
     1071                    axis(depth_axis)
    10651072                if which_quantity == 'stage':
    1066                     if elevations[0:n[j]-1,k,j] < 0:
     1073                    if elevations[0,k,j] < 0:
    10671074                        plot(model_time[0:n[j]-1,k,j], stages[0:n[j]-1,k,j], '-', c = cstr[j])
     1075                        axis(stage_axis)
    10681076                    else:
    10691077                        plot(model_time[0:n[j]-1,k,j], depths[0:n[j]-1,k,j], '-', c = cstr[j])
    1070                     #axis(stage_axis)
     1078                        axis(depth_axis)                 
    10711079                    units = 'm'
    10721080                if which_quantity == 'momentum':
    10731081                    plot(model_time[0:n[j]-1,k,j], momenta[0:n[j]-1,k,j], '-', c = cstr[j])
    1074                     #axis(mom_axis)
     1082                    axis(mom_axis)
    10751083                    units = 'm^2 / sec'
    10761084                if which_quantity == 'xmomentum':
    10771085                    plot(model_time[0:n[j]-1,k,j], xmom[0:n[j]-1,k,j], '-', c = cstr[j])
    1078                     #axis(mom_axis)
     1086                    axis(mom_axis)
    10791087                    units = 'm^2 / sec'
    10801088                if which_quantity == 'ymomentum':
    10811089                    plot(model_time[0:n[j]-1,k,j], ymom[0:n[j]-1,k,j], '-', c = cstr[j])
    1082                     #axis(mom_axis)
     1090                    axis(mom_axis)
    10831091                    units = 'm^2 / sec'
    10841092                if which_quantity == 'speed':
    10851093                    plot(model_time[0:n[j]-1,k,j], speed[0:n[j]-1,k,j], '-', c = cstr[j])
    1086                     #axis(vel_axis)
     1094                    axis(vel_axis)
    10871095                    units = 'm / sec'
    10881096                if which_quantity == 'bearing':
     
    11701178                s = '\end{tabular} \n \\caption{%s} \n \label{fig:%s} \n \end{figure} \n \n' %(caption, label)
    11711179                fid.write(s)
    1172                 c += 1
    1173                 if c % 6 == 0: fid.write('\\clearpage \n')
     1180                cc += 1
     1181                if cc % 6 == 0: fid.write('\\clearpage \n')
    11741182                savefig(graphname_latex)               
    11751183               
     
    11861194                    else:
    11871195                        word_quantity += plot_quantity[i]
    1188                 print 'hello', elevations[0,k,j]
    11891196                where1 = 0
    11901197                count1 += 1
     
    12131220            s = '\end{tabular} \n \\caption{%s} \n \label{fig:%s} \n \end{figure} \n \n' %(caption, label)
    12141221            fid.write(s)
    1215             c += 1
    1216             if c % 6 == 0: fid.write('\\clearpage \n')         
     1222            if float((k+1)/div - pp) == 0.:
     1223                fid.write('\\clearpage \n')
     1224                pp += 1
    12171225           
    12181226            #### finished generating figures ###
Note: See TracChangeset for help on using the changeset viewer.