question about how phytoplankton bind chlorophyll in fennel

[YunPengC]

Hi dear,
I am sorry for the very basic subject, but I want to know how phytoplankton bind chlorophyll in fennel model.
I turn off the horizontal advection term just for phytoplankton in step3d_t.F and pre_step3d.F. The image below shows the result just in the first few times.
The first line is the result of turning off phytoplankton, we can see that the phytoplankton on the left has a significant reduction, the chlorophyll on the right looks like independent on phytoplankton, and the second row is the original result, the chlorophyll looks like dependent on phytoplankton.
Could anyone tell me how the relationship between phytoplankton and chlorophyll in fennel？I don't think it is just in source sink items. Anyone have better suggestions?

Thanks for any replies you have,

Yunpeng

Code: Select all

 DO j=Jstr,Jend
            DO i=Istr,Iend
              IF (itrc.eq.5) THEN            
               t(i,j,k,3,itrc)=Hz(i,j,k)*(cff1*t(i,j,k,nstp,itrc)+      &
               cff2*t(i,j,k,nnew,itrc)) 
               ELSE
              t(i,j,k,3,itrc)=Hz(i,j,k)*(cff1*t(i,j,k,nstp,itrc)+       &
     &                                   cff2*t(i,j,k,nnew,itrc))-      &
     &                        cff*pm(i,j)*pn(i,j)*                      &
     &                        (FX(i+1,j)-FX(i,j)+                       &
     &                         FE(i,j+1)-FE(i,j))
      !                   END IF              

            END DO
          END DO

Code: Select all

DO j=J_RANGE
            DO i=I_RANGE
              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.eq.5) THEN    
                Ta(i,j,k,itrc)=t(i,j,k,nnew,itrc)
                ELSE                               
               Ta(i,j,k,itrc)=t(i,j,k,nnew,itrc)-cff3
                END IF
               END DO
          END DO