2016年9月29日木曜日

time step for tracers in MY_TRC

PISCES/trcini_pisces.F90

rfact = rdttrc(1)

HH this should equal rdt (time step for dynamics) unless different tracer time step is given.

rfact2 = rfact / float (nrdttrc)

HH this would still be rfact unless nrdttrc (time step for biology) is different from tracers.

xstep = rfact2 / rday

HH rday = 86400 sec/day (defined in OPA_SRC/DOM/phycst.F90), so xstep gives the time step in days.. PISCES uses xstep when the rate constant is given in units of d-1. I will avoid using xstep and instead convert the rate constant to s-1 during the initialization as in FABM.

0 件のコメント:

コメントを投稿