[MITgcm-support] problem running optim.x

[Charlotte Breitkreuz]

Hello everybody,

I recently started to work with the adjoint of the MITgcm and having a 
problem using optim.x.
The problem comes up in a configuration similar to 
./verification/global_ocean.cs32x15, but also in the unchanged(!) 
./verification/tutorial_global_ocean_optim.

Compiling and running mitgcmuv_ad as described in the manual in section 
3.18 works fine. Compiling optim.x also works fine.
Then when I try to use optim.x after running mitgcmuv_ad I get the 
following message:

----------------------------------------------------------------------------------
   ==================================================
   Large Scale Optimization with off-line capability.
   ==================================================

                  Version 2.1.0

   OPTIM_NUMBMOD: Control options have been read.
   OPTIM_NUMBMOD: Minimization options have been read.
  pathei-lsopt in optim_readdata

   OPTIM_READDATA: Reading control vector
              for optimization cycle:            0

  dfile= ecco_ctrl
  yctrlid= MIT_CE_000
  nopt=            0
  fname= ecco_ctrl_MIT_CE_000.opt0000
  opened file ecco_ctrl_MIT_CE_000.opt0000
At line 1717 of file optim_readdata.f (unit = 20, file = 
'ecco_ctrl_MIT_CE_000.opt0000')
Fortran runtime error: End of file

-----------------------------------------------------------------------------------

The lines in optim_readdata.f before the error look like this:

-----------------------------------------------------------------------
       write(fname(1:128),'(4a,i4.4)')
      &     dfile,'_',yctrlid(1:10),'.opt', nopt

       open( funit, file   = fname,
      &     status = 'old',
      &     form   = 'unformatted',
      &     access = 'sequential'   )
       print*, 'opened file ', fname

c--   Read the header.
       read( funit ) nvartype
       read( funit ) nvarlength <-- 1717
-----------------------------------------------------------------

So the file 'ecco_ctrl_MIT_CE_000.opt0000' is opened and the first value 
nvartype is actually read, but then the file does not contain a second one.
Possibly, someone else can reproduce this error just compiling and 
running mitgcmuv_ad in verification/tutorial_global_oce_optim/code_ad 
and input_ad, respectively, and then looking at the file 
'ecco_ctrl_MIT_CE_000.opt0000'.


I guess that this is not a problem of optim.x but more a problem of how 
the ecco_ctrl file is written by mitgcmuv_ad. This is done in 
ctrl_pack.F, but everything looks fine there..

I work with checkpoint 65v, use the adjoint_options file which is 
located in the ./tutorial_global_ocean_optim/code_ad folder. My compiler 
for the mitgcmuv_ad is  mpif77.

If someone else does not encounter this problem, could it be a problem 
of my compiler or my machine?

Any help is appreciated! Let me know if you need any more information.

Thanks a lot!

Regards,
Charlotte


-- 
*Charlotte Breitkreuz*, M. Sc. Technomathematik
PhD Student
Research Group 'Geosystem Modeling'

*MARUM* - Center for Marine Environmental Sciences, University of Bremen
*FB5* - Department of Geosciences, University of Bremen
Klagenfurter Straße 2
28359 Bremen, Germany
Room 5420

Tel: +49(0)421 - 218 65448

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20160726/df109f2e/attachment.htm>