ROMS-TOPAZ pCO2csurf, CO2_flux, O2_flux 2차원 진단변수 output 추가

[HCJung-jbnu]

• ~/Nonlinear/Biology/topaz_inp.h

          CASE ('Dout(iCOfx)')
            IF (iDbio2(iCOfx).eq.0) THEN
              IF (Master) WRITE (out,40) 'iDbio2(iCOfx)'
              exit_flag=5
              RETURN
            END IF
            Npts=load_l(Nval, Cval, Ngrids, Lbio)
            i=iDbio2(iCOfx)
            DO ng=1,Ngrids
              Dout(i,ng)=Lbio(ng)
            END DO
          CASE ('Dout(iO2fx)')
            IF (iDbio2(iO2fx).eq.0) THEN
              IF (Master) WRITE (out,40) 'iDbio2(iO2fx)'
              exit_flag=5
              RETURN
            END IF
            Npts=load_l(Nval, Cval, Ngrids, Lbio)
            i=iDbio2(iO2fx)
            DO ng=1,Ngrids
              Dout(i,ng)=Lbio(ng)
            END DO
          CASE ('Dout(ipCO2)')
            IF (iDbio2(ipCO2).eq.0) THEN
              IF (Master) WRITE (out,40) 'iDbio2(ipCO2)'
              exit_flag=5
              RETURN
            END IF
            Npts=load_l(Nval, Cval, Ngrids, Lbio)
            i=iDbio2(ipCO2)
            DO ng=1,Ngrids
              Dout(i,ng)=Lbio(ng)
            END DO

• ~/Nonlinear/Biology/topaz_mod.h

    INTEGER:: iCOfx  !"CO2 air-sea flux" "mol/m2/sec"
    INTEGER:: iO2fx  !"O2 air-sea flux" "mol/m2/sec"
    INTEGER:: ipCO2  !"Ocean surface pCO2" "ppmv"
  ...........
    iCOfx=ic+7 !"CO2 air-sea flux" "mol/m2/sec"
    iO2fx=ic+8 !"O2 air-sea flux" "mol/m2/sec"
    ipCO2=ic+9 !"Ocean surface pCO2" "ppmv"
  ...........
    NDbio2d=9
  .........

• ~/Nonlinear/Biology/topaz_var.h

          CASE ('iDbio2(iCOfx)') !7
            iDbio2(iCOfx)=varid
          CASE ('iDbio2(iO2fx)') !8
            iDbio2(iO2fx)=varid
          CASE ('iDbio2(ipCO2)') !9
            iDbio2(ipCO2)=varid

• ~/Nonlinear/Biology/topaz.h

     DO j = Jstr,Jend
     DO i = Istr,Iend
  #ifdef MASKING
       imask=GRID(ng)%rmask(i,j)
  # ifdef WET_DRY
       imask=GRID(ng)%rmask_full(i,j)
  # endif
  #else
       imask=1 !nomasking
  #endif
  ................
       DIAGS(ng)%DiaBio2d(i,j,7)=TP(ng)%co2_flux(i,j) * imask
       DIAGS(ng)%DiaBio2d(i,j,8)=TP(ng)%o2_flux(i,j) * imask
       DIAGS(ng)%DiaBio2d(i,j,9)=TP(ng)%pco2_csurf(i,j) * imask
     END DO
     END DO

• bio_topaz.in

  ! Logical switches (TRUE/FALSE) to activate writing of time-averaged,
  ! biological processes diagnostics terms into DIAGNOSTIC output file [Ngrids].
  ..................
  Dout(iCOfx) == T        !CO2 air-sea flux mol/m2/sec
  Dout(iO2fx) == T        !O2 air-sea flux mol/m2/sec
  Dout(ipCO2) == T        !Ocean surface pCO2 ppmv

• varinfo_topaz.in

  
  'CO2_airsea'
  'CO2 air-sea flux'
  'mol/m2/sec'
  'CO2_airsea, scalar, series'
  'ocean_time'
  'iDbio2(iCOfx)'
  'r2dvar'
  1.0d0

'O2_airsea' 'O2 air-sea flux' 'mol/m2/sec' 'O2_airsea, scalar, series' 'ocean_time' 'iDbio2(iO2fx)' 'r2dvar' 1.0d0

'pCO2csurf' 'Ocean surface pCO2' 'ppmv' 'pCO2csurf, scalar, series' 'ocean_time' 'iDbio2(ipCO2)' 'r2dvar' 1.0d0

- ~/Nonlinear/Utility/get_state.F : iCOfx, iO2fx, ipCO2 는 초기장 READ시 제외

....... !!!iDbio2d DO itrc=1, 6 !!NDbio2d, only iCOAL, iCOSC, iCOCS, iO2AL, iO2SC, iO2CS foundit=find_string(var_name, n_var, & & TRIM(Vname(1,iDbio2(itrc))), varid) .......

[HCJung-jbnu]

• ~/Nonlinear/Biology/TOPAZ/generic_TOPAZ_model.F90

    if ( co2flux ) then !!turn on CO2 air-sea exchange
      do j = jsc, jec
      do i = isc, iec
        !co2
        cair(i,j)=TP(ng)%co2_flux_alpha(i,j)*TP(ng)%pCO2air*TP(ng)%apr(i,j)*9.7561e-006
        flux_co2(i,j)=(pvel(i,j)*sqrt(660/(TP(ng)%co2_flux_sc_no(i,j)+epsln))* &
                                  (TP(ng)%co2_flux_csurf(i,j)-cair(i,j)))*INTP(ng)%grid_tmask(i,j,1)
      end do
      end do
      flux_co2=-flux_co2
      call g_tracer_set_values(tracer_list,'dic','stf',flux_co2,isd,jsd)
  
      TP(ng)%co2_flux(:,:)=flux_co2(:,:)   <----
  
    end if
  
    if ( o2flux ) then !!turn on O2 air-sea exchange
      do j = jsc, jec
      do i = isc, iec
        pvel(i,j)=9.36e-007*(TP(ng)%wndm(i,j)**2)*(1.0-TP(ng)%fr_i(i,j))
        !o2
        cair(i,j)=TP(ng)%o2_flux_alpha(i,j)*TP(ng)%pO2air*TP(ng)%apr(i,j)*9.7561e-006
        flux_o2(i,j)=(pvel(i,j)*sqrt(660/(TP(ng)%o2_flux_sc_no(i,j)+epsln))* &
                                 (TP(ng)%o2_flux_csurf(i,j)-cair(i,j)))*INTP(ng)%grid_tmask(i,j,1)
      end do
      end do
      flux_o2=-flux_o2
      call g_tracer_set_values(tracer_list,'o2','stf',flux_o2,isd,jsd)
  
      TP(ng)%o2_flux(:,:)=flux_o2(:,:) <----
  
    end if

• ~/Nonlinear/Biology/TOPAZ/generic_TOPAZ.F90 : pCO2csurf 값 가져오기

  subroutine generic_TOPAZ_update_from_source(tracer_list,Temp,Salt,rho_dzt,dzt,hblt_depth,&
     ilb,jlb,tau,dt,grid_dat,model_time,nbands,max_wavelength_band,sw_pen_band,opacity_band, rad)
  ....
  !!hcjung 20230630, get pco2_csurf
  TP(tp_ng)%pco2_csurf(:,:)=topaz%pco2_csurf(:,:)