Fix time derivative

demo/app/convert-time-data-to-json.f90

@@ -0,0 +1,36 @@
+! Copyright (c) 2022-2024, The Regents of the University of California and Sourcery Institute
+! Terms of use are as specified in LICENSE.txt
+program convert_time_data_to_json
+  !! Convert icar time-step output file to JSON format
+  use julienne_m, only : string_t, file_t, command_line_t
+  use time_data_m, only : time_data_t, icar_output_file_t
+  implicit none
+
+  character(len=*), parameter :: usage_info = &
+    new_line('') // new_line('') // &
+    'Usage:  ./build/run-fpm.sh run convert-time-data-to-json --  --input-file <string> --output-file <string>' // new_line('') // &
+    'where angular brackets (<>) denote user-provided input.' // &
+    new_line('')
+
+  associate(time_data => time_data_t(icar_output_file_t(file_t(file_path("--input-file")))))
+    associate(json_file => time_data%to_json())
+      call json_file%write_lines(file_path("--output-file"))
+    end associate
+  end associate
+
+  print *,new_line('') // "______ convert-time-data-to-json done _______"
+
+contains
+
+  function file_path(flag) result(string)
+    character(len=*), intent(in) :: flag
+    type(string_t) string
+    character(len=:), allocatable :: flag_value
+    type(command_line_t) command_line
+
+    flag_value = command_line%flag_value(flag)
+    if (len(flag_value)==0) error stop usage_info
+    string = string_t(flag_value)
+  end function
+
+end program

demo/app/plot-normalized-histograms.gnu

@@ -2,10 +2,10 @@ set logscale y 10
 set key font ",20"
 set title "Input-Tensor Histograms Mapped to the Interval [0,1]" font "Arial,20"
 
-plot for [col=2:8] 'training_inputs_stats.plt' using 1:col with linespoints title columnheader
+plot for [col=2:8] 'inputs_stats.plt' using 1:col with linespoints title columnheader
 pause mouse any "Plot 1 of 2. Press any key to continue.\n"
 
 set title "Output-Tensor Histograms Mapped to the Interval [0,1]" font "Arial,20"
 
-plot for [col=2:6] 'training_outputs_stats.plt' using 1:col with linespoints title columnheader
+plot for [col=2:6] 'outputs_stats.plt' using 1:col with linespoints title columnheader
 pause mouse any "Plot 2 of 2. Press any key to exit.\n"

demo/app/plot-raw-histograms.gnu

@@ -3,7 +3,7 @@ set term png size 2600, 1200
 
 set output "input-tensor-histograms.png"
 set multiplot layout 3,3
-do for [name in "potential-temperature pressure qc qr qs qv temperature"] {
+do for [name in "potential_temperature pressure qc qr qs qv temperature"] {
   filename = name . ".plt"
   set title sprintf("Histogram for %s",name)
   plot filename using 1:2 with linespoints title columnheader
@@ -14,7 +14,7 @@ unset multiplot
 
 set output "output-tensor-histograms.png"
 set multiplot layout 3,2
-do for [name in "dptdt dqcdt dqrdt dqsdt dqvdt"] {
+do for [name in "dpotential_temperature_dt dqc_dt dqr_dt dqs_dt dqv_dt"] {
   filename = name . ".plt"
   set title sprintf("Histogram for %s",name)
   plot filename using 1:2 with linespoints title columnheader