[MITgcm-devel] lab_sea.hb87 restart problem

[Jean-Michel Campin]

Hi Martin,

On Sun, Nov 04, 2007 at 09:14:48AM +0100, Martin Losch wrote:
> Hi Jean-Michel,
> 
> no, I did not check that, because the Hibler+Bryan stress formulation  
> should not require any restarts ... I can't think of a reason why it  
> should:
> for useHB87stressCoupling=.true., the stress on the ocean is computed  
> differently at the end of seaice_dynamics (in the case of the  
> verification experiment it's based on the internal stress divergence  
> computed earlier in seaice_evp.F). So if the restart does not work,  
> it should be related to the evp solver rather than the hb87 stress ...
> 
> Strange. Where do I find you restart test script, so I can start  
> playing with this, too.
> Martin
> 

I updated the little script that I use:
MITgcm_contrib/jmc_script/tst1+1
 
and the sequence of command I typed to check the restart was:
% cd verification
% ./testreport -dd -t lab_sea
% cd lab_sea/tr_run.hb87
% mkdir temp
% mv pickup* temp
% cp data data.tst
% rm data

edit data.tst:
(useSingleCpuIO is commented because of limitation of the tst1+1 script)
% diff data data.tst
53c53
<  useSingleCpuIO=.TRUE.,
---
> #useSingleCpuIO=.TRUE.,
86a87,88
>  nIter0=10,
>  nTimeSteps=4,
88,89c90,91
<  startTime=0.0,
<  endTime=36000.,
---
> #startTime=0.0,
> #endTime=36000.,

% cp ../../../../MITgcm_contrib/jmc_script/tst1+1 .
% ./tst1+1 1
% ./tst1+1 2
And then you should see some differences in cg2d_ini output.
All the results are put in separated dir, so you can make some
plots and see where the differences are located.

the last step is to make a diff of pickup files:
% ./tst1+1 3

Jean-Michel