FENNEL pH 진단변수 출력(outputs)

[HCJung-jbnu]

• varinfo.cfs.dat : 변수 등록

  !!hcjung 20230626, add pH
  'pH'                                      ! Input/Output
  'Surface concentration of hydrogen ions'
  'nondimensional'                        ! [nondimensional]
  'pH, scalar, series'
  'ocean_time'
  'iDbio2(isfpH)'
  'r3dvar'
  1.0d0

[HCJung-jbnu]

• ~/Nonlinear/Biology/fennel_inp.h : varinfo.cfs.dat으로부터 변수 정보 read

  #ifdef CARBON
  .....
  !!!hcjung 20230626, add pH
          CASE ('Dout(isfpH)')
            IF (iDbio2(isfpH).eq.0) THEN
              IF (Master) WRITE (out,40) 'iDbio2(isfpH)'
              exit_flag=5
              RETURN
            END IF
            Npts=load_l(Nval, Cval, Ngrids, Lbio)
            i=iDbio2(isfpH)
            DO ng=1,Ngrids
              Dout(i,ng)=Lbio(ng)
            END DO
  !!!hcjung 20230626, add pH
  .....
  #endif

[HCJung-jbnu]

• /Nonlinear/Biology/fennel_mod.h : pH 변수 index 등록

  ...
  !!!hcjung 20230626, add pH
    integer  :: isfpH                       ! Surface concentration of hydrogen ions
  !!!hcjung 20230626, add pH
  ...
  # ifdef CARBON
    iCOfx=ic+1
    ipCO2=ic+2
  !!hcjung 20230627, add pH
    isfpH=ic+3
  !      ic=ic+2
    ic=ic+3
  !!hcjung 20230627, add pH
  # endif
  ...

[HCJung-jbnu]

• ~/Nonliear/Biology/fennel_var.h : pH 변수 index 지정

  ...
  # ifdef CARBON
        CASE ('iDbio2(iCOfx)')
          iDbio2(iCOfx)=varid
        CASE ('iDbio2(ipCO2)')
          iDbio2(ipCO2)=varid
  !!!hcjung 20230626, add pH
        CASE ('iDbio2(isfpH)')
          iDbio2(isfpH)=varid
  !!!hcjung 20230626, add pH
  # endif
  ...

[HCJung-jbnu]

• ~/Nonlinear/Biology/fennel.h : OCEAN%pH(i,j) -> Diabio2d(i,j,isfpH)

  ......
  !!hcjung 20230627, add pH
          DiaBio2d(i,j,isfpH)=pH(i,j)
  #  ifdef WET_DRY
          DiaBio2d(i,j,isfpH)=DiaBio2d(i,j,isfpH)*rmask_full(i,j)
  #  endif
  !!hcjung 20230627, add pH
  .....

[HCJung-jbnu]

• bio_Fennel_takYS.in : write output switch 추가

  
  ! Logical switches (TRUE/FALSE) to activate writing of time-averaged,
  ! biological processes diagnostics terms into DIAGNOSTIC output file [Ngrids].

Dout(isfpH) == T ! pH Surface concentration of hydrogen ions !!hcjung 20230627, add pH <--