Skip to content
This repository has been archived by the owner on Nov 1, 2023. It is now read-only.

Support trajectory profile DSG data type #6

Open
noaaroland opened this issue Dec 8, 2020 · 1 comment
Open

Support trajectory profile DSG data type #6

noaaroland opened this issue Dec 8, 2020 · 1 comment

Comments

@noaaroland
Copy link
Collaborator

I have enhanced LAS 9 to include products to make a trajectory plot and a profile plot from a trajectoryprofile data sets. When these products are selected they produce the symbol scripts in the attached zip, but the current products that run to make each of those types of plots fail with errors since the DSG type is not as expected.

We need to either modify the script sequence that gets invoked or create new ones for these two types of products.

Completing this sequence of steps will give us very basic support for this data type where in we can create a trajectory plot (which would if we use the default PyFerret behavior, show the average value of the profile between the selected depths) and a ZT profile ribbon plot of the selected profiles.

Here are the steps for this to be included in an LAS release:

  1. Release PyFerret 7.6.3 (or higher) with native support for trajectory profiles. Already done if you don't include the -c conda-forge on the conda update pyferret nothing happens to pyferret. Should be conda update -c conda-forge pyferret.
    2.Using the symbol set tjp_trajectory.jnl (see attached tjp.zip) (and any additional inputs needed) create a trajectory plot (probably using PLOT/ALONG=xy) This runs to completion but the plot doesn't show much because the range is too small.
  2. Using the symbol set tjp_profile.jnl (see attached tjp.zip) create a ribbon plot in ZT showing the profiles in the selected XYZT region. This appears to work just fine.

tjp.zip

The Trajectory plot:
b12976f9cf262c9f21f07b84f3f95cba_plot_image

The profile plot:
c16ba57ca446dda89c513174a7a2c0db_plot_image

@noaaroland
Copy link
Collaborator Author

Here's the current state:

It runs on whatever the default pyferret build is on dunkel:

[rhs@dunkel scripts]$ pyferret
        NOAA/PMEL TMAP
        PyFerret v7.63 (optimized)
        Linux 2.6.32-754.33.1.el6.x86_64 - 09/28/20
         8-Dec-20 10:57

yes? use c16ba57ca446dda89c513174a7a2c0db_netcdf.nc
yes? define variable allgood = temperature[z=@ngd]
yes? define symbol ngood_data = `allgood[e=@sum]`
 !-> define symbol ngood_data = 8
yes?

However that same command fails on the build that currently ships with conda:

(pyferret) schweitzer@GREEN:/mnt/c/Users/schweitzer/Documents/tjp$ pyferret
        NOAA/PMEL TMAP
        PyFerret v7.63 (optimized)
        Linux 4.15.0-1096-azure - 10/13/20
         8-Dec-20 11:00

yes? use c16ba57ca446dda89c513174a7a2c0db_netcdf.nc
yes? define variable allgood = temperature[z=@ngd]
yes? define symbol ngood_data = `allgood[e=@sum]`
 **ERROR: invalid command: grave accent doesnt evaluate to scalar
          allgood[e=@sum]
yes?

@noaaroland noaaroland transferred this issue from NOAA-PMEL/las9-old Jan 25, 2022
Sign up for free to subscribe to this conversation on GitHub. Already have an account? Sign in.
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant