Added optional frost point trace. Turned off by default, press F to t…

…urn on.
a-urq · Nov 1, 2023 · b1ccf3a · b1ccf3a
1 parent 1b75520
commit b1ccf3a
Show file tree

Hide file tree

Showing 10 changed files with 230 additions and 145 deletions.
diff --git a/bin/com/ameliaWx/soundingViewer/Sounding.class b/bin/com/ameliaWx/soundingViewer/Sounding.class
diff --git a/bin/com/ameliaWx/soundingViewer/SoundingFrame$SoundingGraphics.class b/bin/com/ameliaWx/soundingViewer/SoundingFrame$SoundingGraphics.class
diff --git a/bin/com/ameliaWx/soundingViewer/SoundingFrame$SoundingKeyListener.class b/bin/com/ameliaWx/soundingViewer/SoundingFrame$SoundingKeyListener.class
diff --git a/bin/com/ameliaWx/soundingViewer/SoundingFrame.class b/bin/com/ameliaWx/soundingViewer/SoundingFrame.class
diff --git a/bin/com/ameliaWx/soundingViewer/res/controls.html b/bin/com/ameliaWx/soundingViewer/res/controls.html
@@ -44,6 +44,10 @@ <h2>Controls</h2>
                 <td>E</td>
                 <td>Toggles whether parcels and convective readouts consider entrainment. Entrainment is an approximation of <a href="https://arxiv.org/pdf/2301.04712.pdf">Peters et. al. 2023</a>, but not exactly what the paper describes.</td>
             </tr>
+            <tr>
+                <td>F</td>
+                <td>Toggles whether the frost point trace (dashed line, bluish-green) is shown. Off by default. Computed according to formula in <a href="https://journals.ametsoc.org/view/journals/mwre/116/11/1520-0493_1988_116_2172_tswrpa_2_0_co_2.xml?tab_body=pdf">Huffman and Norman 1988, Appendix B</a>.</td>
+            </tr>
             <tr>
                 <td>P</td>
                 <td>Shows a dialog where the user can choose the convective parcel path displayed on the Skew-T diagram.</td>

diff --git a/bin/com/ameliaWx/soundingViewer/unixTool/RadiosondeWrapper.class b/bin/com/ameliaWx/soundingViewer/unixTool/RadiosondeWrapper.class
diff --git a/src/com/ameliaWx/soundingViewer/Sounding.java b/src/com/ameliaWx/soundingViewer/Sounding.java
@@ -7,6 +7,7 @@ public class Sounding {
 	private double[] temperature;
 	private double[] wetbulb;
 	private double[] dewpoint;
+	private double[] frostPoint;
 	private double[] height;
 	private double[] uWind;
 	private double[] vWind;
@@ -19,6 +20,7 @@ public Sounding(double[] pressureLevels, double[] temperature, double[] dewpoint
 
 	public Sounding(double[] pressureLevels, double[] temperature, double[] dewpoint, double[] height,
 			double[] uWind, double[] vWind, double[] wWind) {
+//		System.out.println("entering constructor");
 		this.pressureLevels = pressureLevels;
 		this.temperature = temperature;
 		this.dewpoint = dewpoint;
@@ -32,6 +34,13 @@ public Sounding(double[] pressureLevels, double[] temperature, double[] dewpoint
 		for(int i = 0; i < wetbulb.length; i++) {
 			wetbulb[i] = WeatherUtils.wetBulbTemperature(temperature[i], dewpoint[i], pressureLevels[i]);
 		}
+
+		frostPoint = new double[dewpoint.length];
+
+		for(int i = 0; i < frostPoint.length; i++) {
+			frostPoint[i] = WeatherUtils.frostPointFromDewpoint(dewpoint[i]);
+//			System.out.println("assigning frost points");
+		}
 	}
 
 	public Sounding(float[] pressureLevels, float[] temperature, float[] dewpoint, float[] height,
@@ -54,6 +63,13 @@ public Sounding(float[] pressureLevels, float[] temperature, float[] dewpoint, f
 		for(int i = 0; i < wetbulb.length; i++) {
 			wetbulb[i] = WeatherUtils.wetBulbTemperature(temperature[i], dewpoint[i], pressureLevels[i]);
 		}
+
+		frostPoint = new double[dewpoint.length];
+
+		for(int i = 0; i < frostPoint.length; i++) {
+			frostPoint[i] = WeatherUtils.frostPointFromDewpoint(dewpoint[i]);
+//			System.out.println("assigning frost points");
+		}
 	}
 
 	private double[] convFloatToDouble(float[] arr) {
@@ -82,6 +98,10 @@ public double[] getDewpoint() {
 		return dewpoint;
 	}
 
+	public double[] getFrostPoint() {
+		return frostPoint;
+	}
+
 	public double[] getHeight() {
 		return height;
 	}

diff --git a/src/com/ameliaWx/soundingViewer/SoundingFrame.java b/src/com/ameliaWx/soundingViewer/SoundingFrame.java
@@ -47,6 +47,11 @@
 import com.ameliaWx.weatherUtils.WeatherUtils;
 
 public class SoundingFrame extends JFrame {
+	// IMPLEMENT WHENEVER YOU HAVE TIME
+	private void outputAsCM1Sounding() {
+		// MAKE GUI DIRECTORY SELECTOR
+	}
+
 	private static final long serialVersionUID = 396540838404275479L;
 
 	public static String dataFolder = System.getProperty("user.home") + "/Documents/SoundingViewer/data/";
@@ -68,6 +73,7 @@ public class SoundingFrame extends JFrame {
 	private ParcelPath pathType = ParcelPath.INFLOW_LAYER;
 
 	private boolean showEntrainment = false;
+	private boolean showFrostPoint = false;
 
 	private boolean parcelPathVisible = true;
 	@SuppressWarnings("unchecked")
@@ -889,7 +895,7 @@ private BufferedImage drawSoundingChart(int width, int height) {
 						g.setFont(normalFont);
 					}
 
-					drawCenteredString(timeString(time), g, (int) ((1750 * scale * (i + 1)) / 4), (int) (25 * scale));
+					drawCenteredString(timeString(time), g, (int) ((1850 * scale * (i + 1)) / 4), (int) (25 * scale));
 				}
 			}
 
@@ -913,6 +919,7 @@ private BufferedImage drawSkewT(double scale, ParcelPath pathType) {
 			double[] temperature = activeSounding.getTemperature();
 			double[] wetbulb = activeSounding.getWetbulb();
 			double[] dewpoint = activeSounding.getDewpoint();
+			double[] frostPoint = activeSounding.getFrostPoint();
 
 			double[] virtTemp = new double[dewpoint.length];
 
@@ -1095,6 +1102,16 @@ private BufferedImage drawSkewT(double scale, ParcelPath pathType) {
 					g.drawLine((int) ((x1W + skew1) * scale), (int) (y1 * scale), (int) ((x2W + skew2) * scale),
 							(int) (y2 * scale));
 
+					if(showFrostPoint) {
+						g.setStroke(thickStroke);
+						g.setColor(new Color(0, 196, 128));
+						double x1F = linScale(223.15, 323.15, 0, 800, frostPoint[i]);
+						double x2F = linScale(223.15, 323.15, 0, 800, frostPoint[i + 1]);
+
+						drawDashedLine(g, (int) ((x1F + skew1) * scale), (int) (y1 * scale), (int) ((x2F + skew2) * scale),
+								(int) (y2 * scale), 10);
+					}
+
 					g.setStroke(thickStroke);
 					g.setColor(new Color(0, 255, 0));
 					double x1D = linScale(223.15, 323.15, 0, 800, dewpoint[i]);
@@ -2729,6 +2746,10 @@ public void keyPressed(KeyEvent e) {
 				showEntrainment = !showEntrainment;
 				sg.repaint();
 				break;
+			case KeyEvent.VK_F:
+				showFrostPoint = !showFrostPoint;
+				sg.repaint();
+				break;
 			case KeyEvent.VK_P:
 				selectParcelPathType();
 				sg.repaint();

diff --git a/src/com/ameliaWx/soundingViewer/res/controls.html b/src/com/ameliaWx/soundingViewer/res/controls.html
@@ -44,6 +44,10 @@ <h2>Controls</h2>
                 <td>E</td>
                 <td>Toggles whether parcels and convective readouts consider entrainment. Entrainment is an approximation of <a href="https://arxiv.org/pdf/2301.04712.pdf">Peters et. al. 2023</a>, but not exactly what the paper describes.</td>
             </tr>
+            <tr>
+                <td>F</td>
+                <td>Toggles whether the frost point trace (dashed line, bluish-green) is shown. Off by default. Computed according to formula in <a href="https://journals.ametsoc.org/view/journals/mwre/116/11/1520-0493_1988_116_2172_tswrpa_2_0_co_2.xml?tab_body=pdf">Huffman and Norman 1988, Appendix B</a>.</td>
+            </tr>
             <tr>
                 <td>P</td>
                 <td>Shows a dialog where the user can choose the convective parcel path displayed on the Skew-T diagram.</td>