Major bug when activating bio-model light penetration

[cetlod]

Context

When using the simulated time varying chlorophyll by TOP biogeochemical model ( ln_qsr_bio =.true. ) to compute the light absorption in the ocean, the heat content trend due to Qsr flux is always zero everytime

Analysis

In tra_qsr_init, the level of light extinction ( nksr ) is not computed in case of BIO light penetration (np_BIO) and then set to 0 by default. The heat content trend is not computed and stay equal to its initial values

  DO jk = 1, nksr
     qsr_hc(:,:,jk) = r1_rau0_rcp * ( etot3(:,:,jk) - etot3(:,:,jk+1) )
  END DO

Recommendation

Compute the tabulated attenuation coef and the level of light extinction exactly in the same way than the RGB light penetration

     CASE( np_BIO )                   !==  BIO light penetration  ==!
         !
         IF(lwp) WRITE(numout,*) '   ==>>>   bio-model light penetration'
         IF( .NOT.lk_top )   CALL ctl_stop( 'No bio model : ln_qsr_bio = true impossible ' )
         CALL trc_oce_rgb( rkrgb )                 ! tabulated attenuation coef.
         !
         nksr = trc_oce_ext_lev( r_si2, 33._wp )   ! level of light extinction
         !
         IF(lwp) WRITE(numout,*) '        level of light extinction = ', nksr, ' ref depth = ', gdepw_1d(nksr+1), ' m'

         !
      END SELECT

[cetlod]

In 13331:

Error: Failed to load processor CommitTicketReference

No macro or processor named 'CommitTicketReference' found

[cetlod]

In 13332:

Error: Failed to load processor CommitTicketReference

No macro or processor named 'CommitTicketReference' found

[cetlod]

In 13333:

Error: Failed to load processor CommitTicketReference

No macro or processor named 'CommitTicketReference' found