[MITgcm-support] ptracer start at advanced time

[Samar Khatiwala]

Hi,

I can think of two solutions:

If you're initializing ptracers from data, as opposed to pickup files, simply set PTRACERS_Iter0 to niter0 
and read in the initial conditions by specifying PTRACERS_initialFile(1)=..., etc

On the other hand if you want to read from a pickup file, then simply hack the last few lines of ptracers_init_varia.F, 
i.e.,

C Read from a pickup file if needed
      IF ( nIter0.GT.PTRACERS_Iter0 .OR.
     &    (nIter0.EQ.PTRACERS_Iter0 .AND. pickupSuff.NE.' ')
     &   ) THEN

       CALL PTRACERS_READ_PICKUP( nIter0, myThid )
      ENDIF

Changing the IF statement above and passing PTRACERS_Iter0 instead of nIter0 to the function call 
should work. I don't offhand see what harm that would do.

Samar


> Hi Everybody,
> 
> I want to start ptracer at an advanced timestep than the model initial
> time. I tried to put "PTRACER_niter0 > niter0" for this purpose, but
> it's of no use and returned me zeros everywhere.
> 
> Is there any simple work around of this problem? I wish to use
> niter0=0 and PTRACERS_Iter0=1000, for example to start the tracers at
> 1000 time step rather than at initial time 0.
> 
> Please suggest me way to do this.
> 
> Thanks,
> 
> Suneet