@@ -797,18 +797,22 @@ def calc_vol_flux(name):
797
797
# Actual cross-sectional area (sq metres)
798
798
# Calculated by interpolating across my bathymetry map
799
799
# Wellington_cross_3 accounts for bearing (6687758 m2 x sin(313.7))
800
+ # Wellington_CTD_1 accounts for ADCP track not following line of longitude
801
+
800
802
actual_area = dict (
801
- wellington_cross_1 = 4629513 ,
802
- wellington_cross_2 = 6378842 ,
803
- wellington_cross_3 = 4835029 ,
804
- wellington_CTD_1 = 7562074 ,
805
- wellington_CTD_2 = 4564666 ,
806
- wellington_CTD_3 = 5188899 ,
807
- barrow_CTD = 12452423 )
803
+ wellington_cross_1 = 4629e3 ,
804
+ wellington_cross_2 = 6378e3 ,
805
+ wellington_cross_3 = 4835e3 ,
806
+ wellington_CTD_1 = 6972e3 ,
807
+ wellington_CTD_2 = 4564e3 ,
808
+ wellington_CTD_3 = 5188e3 ,
809
+ barrow_CTD = 12452e3 )
808
810
809
811
area_percent = 100 * area / actual_area [name ]
810
812
print (' Percent area covered by ADCP: {0:2.0f}\n \n ' .format (area_percent ))
811
813
812
814
813
815
# if __name__ == '__main__':
814
- # froude_plots()
816
+ # plot_ctd_as_section(np.r_[102:107], 'wellington_CTD_1')
817
+ # plot_ctd_as_section(np.r_[107:112], 'wellington_CTD_2')
818
+ # plot_ctd_as_section(np.r_[112:117], 'wellington_CTD_3')
0 commit comments