Skip to content

Remove unused vct_b#1067

Closed
havogt wants to merge 1 commit intomainfrom
remove_vct_b
Closed

Remove unused vct_b#1067
havogt wants to merge 1 commit intomainfrom
remove_vct_b

Conversation

@havogt
Copy link
Contributor

@havogt havogt commented Feb 18, 2026

No description provided.

@github-actions
Copy link

Mandatory Tests

Please make sure you run these tests via comment before you merge!

  • cscs-ci run default
  • cscs-ci run distributed

Optional Tests

To run benchmarks you can use:

  • cscs-ci run benchmark-bencher

To run tests and benchmarks with the DaCe backend you can use:

  • cscs-ci run dace

To run test levels ignored by the default test suite (mostly simple datatest for static fields computations) you can use:

  • cscs-ci run extra

For more detailed information please look at CI in the EXCLAIM universe.

@havogt havogt requested a review from jcanton February 18, 2026 16:20
@havogt
Copy link
Contributor Author

havogt commented Feb 18, 2026

cscs-ci run default

@havogt
Copy link
Contributor Author

havogt commented Feb 18, 2026

cscs-ci run distributed

Comment on lines 289 to 303
@@ -311,19 +300,17 @@ def _read_vct_a_and_vct_b_from_file(
file_path: Path to the vertical grid file
num_levels: number of cell levels
allocator: GT4Py field allocator
Returns: one dimensional vct_a and vct_b arrays.
Returns: one dimensional vct_a arrays.
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

maybe explain here why vct_b is ignored (even though we keep the file format including vct_b for interoperability with fortran)

Comment on lines +331 to 334
Compute vct_a.

When the thickness of lowest level grid cells is larger than 0.01:
vct_a[k] = H (2/pi arccos((k/N)^s) )^d N = num_levels s = stretch_factor in vertical configuration
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

and add here a line, something like if needed vct_b can be computed as vct_b = np.exp(-vct_a / 5000.0) (so we don't have to look it up somewhere else)

Comment on lines 519 to 520
It is also used to modify the initial vertical wind speed above surface according to a prescribed vertical profile by linearly merging the surface vertical wind with the existing vertical wind.
See init_w and adjust_w in mo_nh_init_utils.f90.
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
It is also used to modify the initial vertical wind speed above surface according to a prescribed vertical profile by linearly merging the surface vertical wind with the existing vertical wind.
See init_w and adjust_w in mo_nh_init_utils.f90.

these lines also refer to vct_b

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

(if you think they're useful, move them to the other comment with the expression for vct_b)

@havogt
Copy link
Contributor Author

havogt commented Feb 27, 2026

Conclusion was: leave it icon4py, because potentially the vct_b might be custom in the grid file (or an alternative computation might be used), but we remove it from the Fortran bindings as we pass the derived quantities.

@jcanton
Copy link
Contributor

jcanton commented Feb 27, 2026

replaced by #1082

@jcanton jcanton closed this Feb 27, 2026
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants