[MITgcm-cvs] MITgcm/pkg/seaice CVS Commit

[Gael Forget]

Update of /u/gcmpack/MITgcm/pkg/seaice
In directory faulks.csail.mit.edu:/u/u2/gforget/MITgcm/pkg/seaice

Modified Files:
	seaice_growth.F 
Log Message:

Merging seaice_growth codes -- part 3 -- the merging, at last.

- re-arranging the blocks of code in a more logical order. That is
c =================PART 0: constants and initializations=============
c ===========PART 1: treat pathological cases (post advdiff)===========
c ================PART 2: determine heat fluxes/reservoirs===========
c =========PART 3: determine effective thicknesses increments========
c ==========PART 4: determine ice cover fraction increments=========-
c =============PART 5: determine ice salinity increments=============
c =====LEGACY PART 5.5: treat pathological cases, then do flooding ======
c ===============PART 6: determine ice age increments================
c ==============PART 7: determine ocean model forcing================

This is only cosmetic change, by virtue of the previous revisions.
It does not affect fwd results, which I wont update. On Faulks, with g77:
Y Y Y Y>16<16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  . pass  lab_sea
Y Y Y Y>16<16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  . pass  lab_sea.fd
Y Y Y Y>16<16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  . pass  lab_sea.hb87
Y Y Y Y>16<16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  . pass  lab_sea.salt_plume
Y Y Y Y -- 16 16 16 16 16 16 16 16 22 22 22 22 22 22 22 22 16 16 16>16< .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  . pass  offline_exf_seaice
Y Y Y Y -- 16 16 16 16 16 16 16 16 22 22 22 22 22 22 22 22 16 16 16 16 16 16 16>16< .  .  .  .  .  .  .  .  .  .  .  . pass  offline_exf_seaice.seaicetd
Y Y Y Y>16<16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  . pass  seaice_obcs
In ad tests, the revision did not affect the foward run, but it did the gradients
Y Y Y Y 16> 4<FAIL  lab_sea
Y Y Y Y 16> 4<FAIL  lab_sea.evp
Y Y Y Y 16>16<pass  lab_sea.noseaice
Y Y Y Y 16> 4<FAIL  lab_sea.noseaicedyn
Y Y Y Y 16>16<pass  offline_exf_seaice
The gradient check were precise to 4 digits. They remain so, and improved slightly. 
The gradient change due to the revision is at the noise level. There is no additional 
recomputations. So we are fine. I will update those results.