Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Iceberg identifications #15

Merged
merged 3 commits into from
Jul 28, 2016
Merged

Conversation

sternalon
Copy link
Contributor

I have added a number scheme so that each iceberg gets a unique number. The number is based on a counter which counts how many icebergs have been created in a particular grid cell.

The iceberg number is called iceberg_num, and is written to the iceberg_trajectory file, to be used in post processing. It is also written and read from the iceberg.res restart file.

The iceberg grd counter is called iceberg_counter_grd. It is written and read from the calving.res restart file (since it is gridded).

The previous method used for identifying icebergs (ie: start_year, start_mass, start_lon, start_lat, start_day) has been left in place, and is still being used in several places to sort icebergs. This can be removed and replaced by the counter system at a later time.

… counts the amount of icebergs created by each cell.

This has been added to the write restart.

I tried to add it to read restart, but it is not working yet.  (And now its time for the weekend)
…erg_num and is saved in the restart file, picksup from the restart file.

The iceberg counter is now working and is saved as a gridded variable which saves to the restart file, and picksup from the restart.

The iceberg num is not yet saved in the iceberg trajectory file.
@@ -2871,9 +2887,10 @@ integer function berg_chksum(berg )
itmp(37)=berg%start_year !Changed from 29 to 37 by Alon
itmp(38)=berg%ine !Changed from 30 to 38 by Alon
itmp(39)=berg%jne !Changed from 31 to 39 by Alon
itmp(40)=berg%iceberg_num !added by Alon
Copy link
Member

Choose a reason for hiding this comment

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

itmp needs to be declared itmp(36+4) in order to avoid an out-of-bounds at this line.

@adcroft
Copy link
Member

adcroft commented Sep 26, 2015

Unfortunately I'm getting

forrtl: severe (408): fort: (2): Subscript #1 of the array SBUF has value 1 which is greater than the upper bound of 0

which has nothing to do with your commits but is due to issue #7. We just started using a patched FMS for Ben's 4d restarts to work in MOM6 and this code doesn't have the fix that allows the iceberg restarts to work. I'll ask to have this sorted out next week.

@adcroft adcroft merged commit 8c0eb53 into NOAA-GFDL:dev/master Jul 28, 2016
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