[MITgcm-devel] snapshots and divided.ctrl

[Matthew Mazloff]

Hello

I am having an issue where occasionally the adjoint snapshot files are  
not properly written out.  I wanted to troubleshoot this and perhaps  
just change the way the I/O for these files are performed, but I am  
having trouble figuring out where taf gets the info for snapshots.  Is  
there a place where we tell taf to write the adjoint state snapshot,  
or does it do that on its own?  I can't seem to locate in any forward  
code the call that generates the adjoint code:

C----------------------------------------------
C read snapshot
C----------------------------------------------
  9813 continue
       if (idivbeg .lt. nchklev_4) then
         open(unit=77,file='snapshot'//filen,status='old',form=
      $'unformatted',iostat=iers)
         if (iers .eq. 0) then
           read(unit=77) adapressure0,adapressure1,adaqh0,adaqh1,adarea,
...

and the later write snapshot one either.  How is this information  
provided to taf?


Matt