[MITgcm-devel] lab_sea.hb87 restart problem

[Martin Losch]

OK, I will give it a try, but not right now.
M.
On 4 Nov 2007, at 21:19, Jean-Michel Campin wrote:

> 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
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel