Hacktoberfest 2026: những issue maintainer đã đánh dấu cho tháng Mười, đang mở và phù hợp người mới. Xem issue Hacktoberfest

smb_input confusing logic

Đang mở
  • Có gì không ổn?
#76 0 bình luận 0 reaction 0 người được giao Xem trên GitHub

Chưa có ai nhận issue này.

Đánh giá

Độ khó
5/5
Thời gian dự kiến
Hơn một tuần
Mức phù hợp với người mới
25/100
Loại issue
Tái cấu trúc
Độ rõ ràng
Cần làm rõ
Mức độ hoạt động
Đình trệ
Công nghệ
fortran
Lĩnh vực
backend

Hướng nghiên cứu

Theo dõi smb_input qua trình điều khiển Glissade, bắt đầu từ mô tả issue và các entry point được tham chiếu là glad_main.F90, glad_timestep.F90 và glide_io.F90. So sánh đầu vào acab/smb standalone với chuyển đổi CESM coupled, sau đó xác định hành vi hoặc tên tùy chọn rõ ràng hơn và đầu ra log hữu ích. Hoàn tất khi người dùng có thể biết trường nào cung cấp SMB và các trường hợp standalone và coupled sử dụng nguồn dự kiến.

Do mô hình lập chỉ mục viết ra từ nội dung của issue.

Mô tả

As of cism_main_2.02.003 the logic surrounding the smb_input field has some confusing and unintended side effects. There is a cism_wrapper issue describing the problems HERE. For ease of reference, here is the text from Bill's email (feel free to add more details in the comments @whlipscomb !)

I think the problem may be with the config setting smb_input. This setting tells CISM (in the Glissade driver) whether the input SMB has units of mm/yr water equivalent (the units of model%climate%smb) or m/s ice (the units of model%climate%acab). If smb_input = 1, Glissade will compute acab by copying from smb and changing the units. If smb_input = 0, Glissade will ignore the smb field and use the existing value of acab.

In standalone CISM, this means we should have smb_input = 0 when reading acab from the input file, and smb_input = 1 when reading smb from the input file.

But in coupled CESM, it's not so simple. That's because we have Glad code in between. While the incoming qsmb from the coupled has units of mm/yr w.e., Glad converts to acab units when it's filling the arrays params%instances%acab and model%climate%acab. So when we get to Glissade, we should not compute model%climate%acab from model%climate%smb, but rather should go with the existing model%climate%acab, i.e. we want smb_input = 0.

[If you want to see the details, take a look at glad_main.F90, starting around l. 800. This is where Glad converts units. The qsmb field comes from the coupler with units of mm/yr w.e. (or equivalently, kg/m2/s) and is copied into params%instances(i)%acab. Then at l. 822, the units are converted to m/yr ice just before the call to glad_i_tstep_gcm. Subroutine glad_i_tstep_gcm is in glad_timestep.F90. At l. 200 there is a call to glide_set_acab. This is where the acab value in the params%instances derived type is copied into model%climate%acab. Subroutine glide_set_acab is in glide_io.F90. It converts the units from m/yr ice to m/s ice, which are the units CISM now uses.]

I think either the logic around these options or the names of the options (or both) should be cleaned up so it is more obvious if you are about to change your settings such that you are not getting SMB where you expect to get it from. This is a high-impact change, and it should be clear to the user what they are doing. Perhaps some log output would be good to have for this as well.

Ngôn ngữ chính
Fortran
Star
13
Fork
17
Merge trung bình
3 ngày 15 giờ
Pull request đã merge (30 ngày)
1

Chuẩn bị môi trường

Dự án này không cung cấp dev container, Dockerfile hay hướng dẫn đóng góp, nên bạn cần tự thiết lập môi trường: hãy bắt đầu từ README và xem hướng dẫn đóng góp lần đầu của chúng tôi để biết các bước chung.

Bắt đầu từ đâu

  1. Đọc hết issue, rồi đọc hướng dẫn đóng góp của dự án.
  2. Bình luận trên issue rằng bạn sẽ nhận — tránh hai người làm cùng một việc.
  3. Fork repository và làm thay đổi trên một nhánh.
  4. Mở pull request có tham chiếu số hiệu của issue.

Issue khác của ESCOMP/CISM

Tất cả issue của ESCOMP/CISM

Issue tương tự

Thêm issue về Backend & API Design

Nhận issue mới trong hộp thư của bạn

Bản tóm tắt ngắn những issue GitHub phù hợp với người mới.