[MITgcm-devel] updates on seaice code consolidation

[Holly Dail]

Hello -

I'm using a CVS checkout from 11/23/2010 and am using the "new" seaice code turned on with the flags set as described by Gael.  In my setup (1 degree Atlantic model from 33S to 75 N with open boundary at 33S) this seems to work fine stability-wise in both the forward and adjoint versions.  I haven't looked at the details of the sea ice extent or thickness though, so can't speak to that aspect.

Holly




On Dec 9, 2010, at Dec 9 , 6:12 PM, Gael Forget wrote:

> Hello Matt,
> 
> the code is ready to be tested. With regard to SEAICE_OPTIONS.h, the default in pkg/seaice contains 
> the currently recommended flags. At this point, the 'new code' that was merged from _if.F is inactive by default 
> though. To turn it on you want to set SEAICE_GROWTH_LEGACY to .FALSE. which would in turn set 
> SEAICE_DO_OPEN_WATER_GROWTH & SEAICE_OCN_MELT_ACT_ON_AREA to .TRUE.
> 
> With regard to the stability of the adjoint, I would be curious to know whether you see 
> problems with seaice_growth.F as it stands. I believe we made some progress along 
> those lines (at least for the thermodynamics) but we haven't gotten much feedback yet.
> 
> Cheers,
> Gael
> 
> On Dec 9, 2010, at 5:00 PM, Matthew Mazloff wrote:
> 
>> Hello,
>> 
>> Thanks for the update!
>> 
>> I'm looking forward to the new code and am happy to test it when its ready.  Is there a verification experiment or somewhere I will be able to find the recommended flags and parameters for the latest version?
>> 
>> Also, what is the status of the stability of the sea-ice adjoint -- should I be trying it with c62p?
>> 
>> Thanks,
>> -Matt
>> 
>> 
>> 
>> 
>> 
>> On Dec 9, 2010, at 1:24 PM, Patrick Heimbach wrote:
>> 
>>> Hi there,
>>> 
>>> I am summarizing a brief telecon which we had today between
>>> JPL (Dimitris, Ian, An) and MIT (Gael, Jean-Michel, myself).
>>> 
>>> Our goal is (has been for a while) to
>>> * incorporate Ian's changes into "main" code
>>> * clean up code
>>> * be able to get rid of the _if versions in pkg/seaice (and MITgcm_contrib/ )
>>> 
>>> The status as of September was:
>>> * the "legacy" code in pkg/seaice
>>> * the "_if" code in pkg/seaice
>>> * latest code from Ian in MITgcm_contrib (at my request) for the purpose of merging with code in pkg/seaice
>>> 
>>> At the time of Martin's and Dimitris' visit we had hoped to be able to
>>> run the MITgcm_contrib/ code, but attempts by several people to do so
>>> failed despite repeated iterations of fixes.
>>> Moving forward we had formulated a strategy at MIT to
>>> get a code that
>>> 1. tries to incorporate Ian's modifs,
>>> 2. do so in a controlled way (traceable via CVS),
>>> 3. cleans up many of the issues of legacy code (yneg, and many other issues)
>>> 5. a cleaned-up sequence that is more amenable to modularization, e.g. for coupling
>>> (e.g. move age calculations and similar to separate routines,
>>> deal with salinity code and salt_plume params),
>>> 5. solves non-conservation issues.
>>> 
>>> Main work on 1.-4. was taken on by Gael, and item 5. tackled by Martin.
>>> As far as I understand, starting point for Gael's merge was the "legacy"
>>> code with main merges taken from _if routines.
>>> 
>>> I guess this has gone a long way, but
>>> in order to bring it to successful completion we have the following issues:
>>> 1. extent to which code improvements in Ian's MITgcm_contrib are still missing
>>> (Ian voiced concerns that some important terms described in his thesis
>>> as adapted from McPhee are still missing)
>>> 2. make sure that major setups that we are aware of run "successfully"
>>> with the latest code, so that we can complete the clean-up
>>> 
>>> To deal with 1. Ian has agreed to document what pieces are still missing
>>> (forthcoming on mitgcm-devel).
>>> 
>>> To deal with 2. some of the major setups need to re-run (and retuned?).
>>> Especially, Matt who's probably the only one having used the _if code
>>> (I lost track which one) should re-run SOSE and see if he's happy with "latest code".
>>> Suggestion is that this be done with the next checkpoint (c62p).
>>> "full list of setups are:
>>> * Matt: SOSE
>>> * Dimitris/Hong: ECCO2
>>> * An/Gunnar: Arctic
>>> * Gael/Patrick: ECCO-GODAE codes
>>> * Holly: Atlantic setup
>>> * Martin: his setup(s)
>>> * Ian: Lab Sea setup, with special focus on adjoint behaviour
>>> 
>>> Hope I captured it all, not too confusing and not too inaccurate.
>>> Cheers
>>> -Patrick
>>> 
>>> ---
>>> Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
>>> MIT | EAPS 54-1518 | 77 Massachusetts Ave | Cambridge MA 02139 USA
>>> FON +1-617-253-5259 | FAX +1-617-253-4464 | SKYPE patrick.heimbach
>>> 
>>> 
>>> 
>>> _______________________________________________
>>> MITgcm-devel mailing list
>>> MITgcm-devel at mitgcm.org
>>> http://mitgcm.org/mailman/listinfo/mitgcm-devel
>> 
>> 
>> _______________________________________________
>> MITgcm-devel mailing list
>> MITgcm-devel at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-devel
> 
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel