[MITgcm-devel] TAF error

[Patrick Heimbach]

Hi Matt,

looks like you have either mixed seaice code (less likely),
or a new seaice package, with old the_main_loop.F (more likely).

Have a look in following block in latest code:
  #  ifdef SEAICE_ALLOW_DYNAMICS
  cphCADJ INIT comlev1_lsr    = COMMON,nchklev_1*2
  CADJ INIT comlev1_dynsol = COMMON,nchklev_1*MPSEUDOTIMESTEPS
  #  endif
  #  endif

(and difference to previous).
That should fix it.

Cheers
-p.

On Jun 3, 2009, at 12:56 PM, Matthew Mazloff wrote:

> Hi Patrick,
>
> I get the error:
>
> TAF ERROR   ad_input_code.f:1071033 : tape not defined
>
> from lines of the nature.
> CADJ STORE uice = comlev1_dynsol, kind=isbyte,
> CADJ &     key = ikey_dynamics + (ipseudo-1)*nchklev_1
>
> Not sure exactly what is going on but since I have seaice package  
> off in adjoint run I am planning on just commenting this store  
> out.  It shows up in
>
> grep 'CADJ &     key = ikey_dynamics' *.F
> seaice_dynsolver.F:CADJ &     key = ikey_dynamics + (ipseudo-1) 
> *nchklev_1
> seaice_dynsolver.F:CADJ &     key = ikey_dynamics + (ipseudo-1) 
> *nchklev_1
> seaice_lsr.F:CADJ &     key = ikey_dynamics + (ilcall-1)*nchklev_1
> seaice_lsr.F:CADJ &     key = ikey_dynamics + (ilcall-1)*nchklev_1
> seaice_lsr.F:CADJ &     key = ikey_dynamics + (ilcall-1)*nchklev_1
> seaice_lsr.F:CADJ &     key = ikey_dynamics + (ilcall-1)*nchklev_1
> seaice_lsr.F:CADJ &     key = ikey_dynamics + (ilcall-1)*nchklev_1
> seaice_lsr.F:CADJ &     key = ikey_dynamics + (ilcall-1)*nchklev_1
>
> I assume you have some across this issue -- is there a simple  
> remedy -- or is benign?
>
> Thanks
> Matt
>
>
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel

---
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