search

orcestra-campaign / weather

Weather briefings for the ORCESTRA field campaign.
5 stars 1 forks source link

Adding the plots of surface wind convergence to the weather briefing #196

Closed dspraturi closed 1 month ago

dspraturi commented 1 month ago

This closes the issue #169 and adds the plots of forecasted convergence to the briefing.

Here is a sample plot of convergence. It looks a bit noisy (to be expected) but there is dominant convergence in the ITCZ. I have mostly repurposed the scripts that I typically use to compute convergence in nextgems data. test_sfc_conv

I also tested with generating the briefing - it works without any issues. :)

sortega87 commented 1 month ago

Thanks @dspraturi! Looks very nice. I will make some changes in the quarto template of this branch to include it as a supporting figure of the presentation. We thought it is good to include these and other figures here first.

Is it possible to also compute the vorticity in the healpix grid? That would be a great complement for the briefing.

sortega87 commented 1 month ago

@HF7weatherman I have now included the support section in the quarto template

sortega87 commented 1 month ago

I think now its ready to merge.

HF7weatherman commented 1 month ago

@dspraturi In case you have an easy way to calculate surface vorticity and would like to implement it: Please have a look at #191 👍🏻

dspraturi commented 1 month ago

@dspraturi In case you have an easy way to calculate surface vorticity and would like to implement it: Please have a look at #191 👍🏻

Yes, I have codes that compute vorticity in healpix grid. I can take a look at it tomorrow. :)

sortega87 commented 1 month ago

@HF7weatherman fixed the horizons as discussed, I think we can merge.

HF7weatherman commented 1 month ago

Looks fine, can be merged 👍🏻

HF7weatherman commented 1 month ago

@sortega87 When trying to generate the briefing, I get the following error message when the code tries to generate the new figures:

OMP: Error #15: Initializing libomp.dylib, but found libomp.dylib already initialized.
OMP: Hint This means that multiple copies of the OpenMP runtime have been linked into the program. That is dangerous, since it can degrade performance or cause incorrect results. The best thing to do is to ensure that only a single OpenMP runtime is linked into the process, e.g. by avoiding static linking of the OpenMP runtime in any library. As an unsafe, unsupported, undocumented workaround you can set the environment variable KMP_DUPLICATE_LIB_OK=TRUE to allow the program to continue to execute, but that may cause crashes or silently produce incorrect results. For more information, please see http://openmp.llvm.org/
./wbcli.sh: line 4: 17487 Abort trap: 6           python3 src/wblib/main.py "$@"
/opt/homebrew/Caskroom/miniforge/base/envs/wblib/lib/python3.12/multiprocessing/resource_tracker.py:254: UserWarning: resource_tracker: There appear to be 2 leaked semaphore objects to clean up at shutdown
  warnings.warn('resource_tracker: There appear to be %d '

Surprisingly, Python really crashes and not only throws this error. I think it has something to do with your speed-up changes, which may have not been implemented for the new convergence plots.

We should fix this tomorrow. For now, I will create a branch temporary_main from the last commit on main before the merge, which I`ll use to generate the weather briefings 👍🏻

sortega87 commented 1 month ago

Hi @HF7weatherman! Sorry for this, and thanks for letting me know. Let's take a look today, perhaps we should revert the changes regarding preloading the dask array. I tried generating all the figures in my computer and it is working without showing the warning or crashing. I wonder if it is due to differences in my environment.