## no advection of biological variables in NPZD-Powell.

Moderators: arango, robertson

Message
Author
sksdhdbfl

### no advection of biological variables in NPZD-Powell.

Hi.
I'm Yuri Oh in South Korea.

I have been conducting a 3 dimensional biological-physical coupled model using the NPZD-Powell model for a last 1 year. ROMS version is 624.
My Question is how I make an advection of biological variables (Nutient,Phyto,Zoo,Detritus) is not working??
I try to investigate the effect of nutrient supplys on ecosystem dividing into two: horizontal supply and vertical supply.
I think no advection of biological variables means vertical supply of biological biomass only.

and I find it out an advection stuff is calculated in step3d_t.F90 and pre_step3d.F90.

I edited only step3d_t.F90 for no advection as below (pre_step3d.F90 was not edited).
But model blew up with huge negative biological values. physical values were fine.

How can I make no advection for biological variables is working correctly??
Thanks

<edited step3d_t.F90 as below>
!-----------------------------------------------------------------------
!-----------------------------------------------------------------------
T_LOOP : DO itrc=1,NT(ng)
K_LOOP : DO k=1,N(ng)
.....
! of Tunits m3/s. The new tracer has units of m Tunits.
DO j=Jstr,Jend
DO i=Istr,Iend
cff=dt(ng)*pm(i,j)*pn(i,j)
cff1=cff*(FX(i+1,j)-FX(i,j))
cff2=cff*(FE(i,j+1)-FE(i,j))
cff3=cff1+cff2
if (itrc .gt. 2) then
t(i,j,k,nnew,itrc)=t(i,j,k,nnew,itrc)
else
t(i,j,k,nnew,itrc)=t(i,j,k,nnew,itrc)-cff3
endif
END DO
END DO
END DO K_LOOP
END DO T_LOOP

kate
Posts: 3804
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

### Re: no advection of biological variables in NPZD-Powell.

I suggest you print out values of T,S as well as bio tracers at every step through pre_step3d.F and step3d_t.F. You will see that time slice 3 is not say temperature but rather temperature times Hz. You might find that you also need to edit pre_step3d.F.

arango