[MITgcm-support] Tracer Reset

[Yuan Lian]

Hi Alon,

I am not familiar with the ocean modeling. There may be some setting 
available for this.

One way to do this is to insert a timer subroutine for loading the 
initial tracer fields, e.g., somewhere
after "CALL PTRACERS_CYCLE(bi,bj,myTime,myIter,myThid)" in 
tracers_correction_step.F.
For example (see how to use the following in codes write_state.F and 
ptracers_init_varia.F):


LOGICAL  DIFFERENT_MULTIPLE
EXTERNAL DIFFERENT_MULTIPLE

_RL resetFreq
resetFreq=10.*365.*86400.

IF DIFFERENT_MULTIPLE(resetFreq,myTime,deltaTClock) THEN
     DO iTracer = 1, PTRACERS_numInUse
         IF ( PTRACERS_initialFile(iTracer) .NE. ' ' ) THEN
          CALL READ_FLD_XYZ_RL(PTRACERS_initialFile(iTracer),' ',
      &        pTracer(1-Olx,1-Oly,1,1,1,iTracer),0,myThid)
          _EXCH_XYZ_RL(pTracer(1-Olx,1-Oly,1,1,1,iTracer),myThid)
         ENDIF
        ENDDO
ENDIF

I hope someone can give you suggestion on an easier method.

Best,
Yuan


On 11/28/12 4:39 PM, Alon Stern wrote:
> Hi
>
> I am using MITGCM and I had a question:
>
> I am running an ocean model using ptracers. I am hoping to 
> reinitialize the tracers every ten years so that I can track their 
> evolution in a ten year period a number of time and gather 
> some statistics.
>
> Is there anyway to make the tracers automatically restart from an 
> initial tracer file every ten years?
>
> Thanks
> Alon
>
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support

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