ROMS-FENNEL carbon tendency -> diag output 추가

[HCJung-jbnu]

• varinfo.cfs.in

  
  !!hcjung 20230706, add TIC tendencies
  'TIC_uptake'
  'TIC uptake during phytoplankton growth'
  'millimole_carbon meter-3 day-1'               ! [millimole/m3/day]
  'TIC_uptake, scalar, series'
  'ocean_time'
  'iDbio3(iTICu)'
  'r3dvar'
  1.0d0

'TIC_zoo' 'TIC zooplankton basal metabolism' 'millimole_carbon meter-3 day-1' ! [millimole/m3/day] 'TIC_zoo, scalar, series' 'ocean_time' 'iDbio3(iTICz)' 'r3dvar' 1.0d0

'TIC_re' 'TIC remineralization' 'millimole_carbon meter-3 day-1' ! [millimole/m3/day] 'TIC_re, scalar, series' 'ocean_time' 'iDbio3(iTICr)' 'r3dvar' 1.0d0

[HCJung-jbnu]

• ~/Nonlinear/Biology/fennel_inp.h

  !!hcjung 20230706, add TIC tendencies
          CASE ('Dout(iTICu)')
            IF (iDbio3(iTICu).eq.0) THEN
              IF (Master) WRITE (out,40) 'iDbio3(iTICu)'
              exit_flag=5
              RETURN
            END IF
            Npts=load_l(Nval, Cval, Ngrids, Lbio)
            i=iDbio3(iTICu)
            DO ng=1,Ngrids
              Dout(i,ng)=Lbio(ng)
            END DO
          CASE ('Dout(iTICz)')
            IF (iDbio3(iTICz).eq.0) THEN
              IF (Master) WRITE (out,40) 'iDbio3(iTICz)'
              exit_flag=5
              RETURN
            END IF
            Npts=load_l(Nval, Cval, Ngrids, Lbio)
            i=iDbio3(iTICz)
            DO ng=1,Ngrids
              Dout(i,ng)=Lbio(ng)
            END DO
          CASE ('Dout(iTICr)')
            IF (iDbio3(iTICr).eq.0) THEN
              IF (Master) WRITE (out,40) 'iDbio3(iTICr)'
              exit_flag=5
              RETURN
            END IF
            Npts=load_l(Nval, Cval, Ngrids, Lbio)
            i=iDbio3(iTICr)
            DO ng=1,Ngrids
              Dout(i,ng)=Lbio(ng)
            END DO

[HCJung-jbnu]

• ~/Nonlinear/Biology/fennel_mod.h

  #ifdef CARBON
  !!hcjung 20230706, add TIC tendencies
    integer  :: iTICu = 4                   ! TIC uptake during phytoplankton growth
    integer  :: iTICz = 5                   ! TIC zooplankton basal metabolism
    integer  :: iTICr = 6                   ! TIC remineralization
  #endif
  .....
  # if defined CARBON
  !!hcjung 20230706, add TIC tendencies
    NDbio3d=6
  # else
    NDbio3d=3
  # endif

[HCJung-jbnu]

• ~/Nonlinear/Biology/fennel_var.h

  # ifdef CARBON
  !!hcjung 20230706, add TIC tendencies
        CASE ('iDbio3(iTICu)')
          iDbio3(iTICu)=varid
        CASE ('iDbio3(iTICz)')
          iDbio3(iTICz)=varid
        CASE ('iDbio3(iTICr)')
          iDbio3(iTICr)=varid
  # endif

[HCJung-jbnu]

• ~/Nonlinear/Biology/fennel.h

  #ifdef DIAGNOSTICS_BIO
  !!hcjung 20230706, add TIC tendencies
              DiaBio3d(i,j,k,iTICu)=DiaBio3d(i,j,k,iTICu)-            &
  # ifdef WET_DRY
   &                                rmask_full(i,j)*                  &
  # endif
   &                                cff1*fiter
  #endif
  .....
  #ifdef DIAGNOSTICS_BIO
  !!hcjung 20230706, add TIC tendencies
              DiaBio3d(i,j,k,iTICz)=DiaBio3d(i,j,k,iTICz)+            &
  # ifdef WET_DRY
   &                                rmask_full(i,j)*                  &
  # endif
   &                                (ZooCN(ng)*(N_Flux_Zmetabo+N_Flux_Zexcret))*fiter
  #endif
  .....
  #ifdef DIAGNOSTICS_BIO
  !!hcjung 20230706, add TIC tendencies
              DiaBio3d(i,j,k,iTICr)=DiaBio3d(i,j,k,iTICr)+            &
  # ifdef WET_DRY
   &                                rmask_full(i,j)*                  &
  # endif
   &                                (C_Flux_RemineS+C_Flux_RemineL)*fiter
  #endif