DIC src/sink term 진단변수로 추출

[HCJung-jbnu]

• varinfo_topaz.dat

  'srcsink_DIC'
  'Dissolved Inorganic Carbon Source Sink Term'
  'mol/m3'
  'srcsink_DIC, scalar, series'
  'ocean_time'
  'iDbio3(iSDIC)'
  'r3dvar'
  1.0d0

• bio_topaz.in

  Dout(iSDIC) == T        !Dissolved inorgaic carbon source/sink term

• topaz_inp.h

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

• topaz_mod.h

  ....
    INTEGER:: iSDIC  !"Dissolved Inorganic Carbon Source/Sink term"
  ....
    NDbio3d=topaz_diag_ntrc + 1
  ....
    iSDIC=ic+12 !"Dissolved Inorganic Carbon Source Sink Term"

• topaz_var.h

          CASE ('iDbio3(iSDIC)') !12
            iDbio3(iSDIC)=varid

• topaz.h

     DO ntrc = 1, NDbio3d - 1 !!except for DIC src/sink term
     DO k = 1, N(ng)
     DO j = Jstr,Jend
     DO i = Istr,Iend
  #ifdef MASKING
         DIAGS(ng)%DiaBio3d(i,j,N(ng)-k+1,ntrc)=DIAGS(ng)%DiaBio3d(i,j,N(ng)-k+1,ntrc) + TP(ng)%diag_trc(i,j,k,ntrc) * GRID(ng)%rmask(i,j)
  # ifdef WET_DRY
         DIAGS(ng)%DiaBio3d(i,j,N(ng)-k+1,ntrc)=DIAGS(ng)%DiaBio3d(i,j,N(ng)-k+1,ntrc) + TP(ng)%diag_trc(i,j,k,ntrc) * GRID(ng)%rmask_full(i,j)
  # endif
  #else
         DIAGS(ng)%DiaBio3d(i,j,N(ng)-k+1,ntrc)=DIAGS(ng)%DiaBio3d(i,j,N(ng)-k+1,ntrc) + TP(ng)%diag_trc(i,j,k,ntrc)
  #endif
     END DO
     END DO
     END DO
     END DO
  
     ntrc = NDbio3d !!DIC src/sink term
     DO k = 1, N(ng)
     DO j = Jstr,Jend
     DO i = Istr,Iend
         cff = TP(ng)%prog_trc(i,j,k,27) - TP(ng)%prog_trc_old(i,j,k,27)
  #ifdef MASKING
         DIAGS(ng)%DiaBio3d(i,j,N(ng)-k+1,ntrc)=DIAGS(ng)%DiaBio3d(i,j,N(ng)-k+1,ntrc) + cff * GRID(ng)%rmask(i,j)
  # ifdef WET_DRY
         DIAGS(ng)%DiaBio3d(i,j,N(ng)-k+1,ntrc)=DIAGS(ng)%DiaBio3d(i,j,N(ng)-k+1,ntrc) + cff * GRID(ng)%rmask_full(i,j)
  # endif
  #else
         DIAGS(ng)%DiaBio3d(i,j,N(ng)-k+1,ntrc)=DIAGS(ng)%DiaBio3d(i,j,N(ng)-k+1,ntrc) + cff
  #endif
     END DO
     END DO
     END DO

[HCJung-jbnu]

• ~/ROMS/Utility/def_rst.F, wrt_rst.F, get_state.F

  
  DO itrc=1,NDbio3d - 1 !!except for DIC src/sink term <-- 수정

RST(ng)%pioTrc(NT(ng)+NDbio3d-1+itrc)%dkind=PIO_FRST !!except for DIC src/sink term <-- 수정

등....