Skip to content

Conversation

@docguibou
Copy link
Contributor

This is the addition 7 new source terms to fully close the ice area & volume budget as requested by Clemens.

I have tried to figure out what was needed (initially we agreed on 5) but it seems if you want to understand what contributes to what completely, you need 7. It may be that some of these 7 may be inferred from other terms or that the separation of each term could be smarter, do not hesitate to come with your suggestions, but I thought by adding them all then it is easy to check whether what we think is equivalent really works or not.

The source terms are:

dci: (old) sea ice area change rate due to thermo [/day]
dci_young: young sea ice area change rate due to thermo [/day]
dci_young2old: young sea ice area transferred to 'old' due to thermo [/day]
dvi_young2old: young sea ice volume transferred to 'old' due to thermo [m/day] (Could be computed from dci_young2old ??)

dci_ridge: (old) sea ice area change rate due to dynamics (ridging) [/day]
dci_ridge_young: young sea ice area change rate due to dynamics (ridging) [/day]
dvi_ridge_young: young sea ice volume change rate due to dynamics (ridging) [m/day] (transferred to 'old')

@docguibou docguibou requested review from HR-polar and einola December 7, 2022 14:35
@docguibou
Copy link
Contributor Author

(Oh, and I forgot to mention, it obviously compiles and runs and I quickly checked that the output make sense with ncview, but haven't taken the time to make sure they all add up)

@einola einola marked this pull request as draft December 16, 2022 09:50
@einola
Copy link
Member

einola commented Dec 16, 2022

We should reconsider the names of some of the variables.

Copy link
Member

@einola einola left a comment

Choose a reason for hiding this comment

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

It looks good, except for one suspicious line delete, which I'd like you to look at.

Is it otherwise ready? It's still marked as a draft.

if ( M_h_young[i] < h_young_min*M_conc_young[i] )
{
M_conc_young[i] = M_h_young[i]/h_young_min;
young_ice_growth =M_conc_young[i] - old_conc_young ;
Copy link
Member

Choose a reason for hiding this comment

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

Why do you remove this line? It seems young_ice_growth is needed for redistributeThermoFSD.

@einola einola self-requested a review July 16, 2023 18:12
@einola
Copy link
Member

einola commented Dec 14, 2023

Ping @docguibou

@docguibou
Copy link
Contributor Author

docguibou commented Dec 14, 2023 via email

@einola
Copy link
Member

einola commented Dec 14, 2023

No worries! :D I just didn't know what the situation was.

@einola
Copy link
Member

einola commented Dec 14, 2023

And there's still the question about removing young_ice_growth from line 5000.

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.

3 participants