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 件のコメント:
コメントを投稿