ROMS-TOPAZ 진단변수 작게 평균되서 나오는 문제

pCO2_csurf 가 매우 작음...
~/Utility/set_diags.F90 -> DIAGS(ng)%DiaBio2d 여기에 진단변수 값을 계속 누적 더하다가 rfac으로 한번에 나눈후 write
```
.....
DIAGS(ng)%DiaBio2d(i,j,idiag)=rfac(i,j)*                &
 &                    DIAGS(ng)%DiaBio2d(i,j,idiag)
.....
```
topaz.h에서 DIAGS(ng)%DiaBio2d/3d 값을 매 타임스텝마다 넣어주는 것이 아닌 누적 더해주는 것으로 변경

또한, 특정 timestep에서 DIAGS(ng)%DiaBio2d/3d 를 리셋(=0)해주는 과정 추가 <- fennel 참조

....
   IF (((iic(ng).gt.ntsDIA(ng)).and.                                 &
 &     (MOD(iic(ng),nDIA(ng)).eq.1)).or.                            &
 &    ((iic(ng).ge.ntsDIA(ng)).and.(nDIA(ng).eq.1)).or.             &
 &    ((nrrec(ng).gt.0).and.(iic(ng).eq.ntstart(ng)))) THEN
     DO ntrc=1,NDbio2d
       DO j=Jstr,Jend
         DO i=Istr,Iend
           DIAGS(ng)%DiaBio2d(i,j,ntrc)=0.0_r8
         END DO
       END DO
     END DO
     DO ntrc=1,NDbio3d
       DO k=1,N(ng)
         DO j=Jstr,Jend
           DO i=Istr,Iend
             DIAGS(ng)%DiaBio3d(i,j,k,ntrc)=0.0_r8
           END DO
         END DO
       END DO
     END DO
   ENDIF
....
   DO ntrc = 1, NDbio3d
   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

   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,1)=DIAGS(ng)%DiaBio2d(i,j,1) + TP(ng)%co2_flux_alpha(i,j) * imask
     DIAGS(ng)%DiaBio2d(i,j,2)=DIAGS(ng)%DiaBio2d(i,j,2) + TP(ng)%co2_flux_sc_no(i,j) * imask
     DIAGS(ng)%DiaBio2d(i,j,3)=DIAGS(ng)%DiaBio2d(i,j,3) + TP(ng)%co2_flux_csurf(i,j) * imask
     DIAGS(ng)%DiaBio2d(i,j,4)=DIAGS(ng)%DiaBio2d(i,j,4) + TP(ng)%o2_flux_alpha(i,j) * imask
     DIAGS(ng)%DiaBio2d(i,j,5)=DIAGS(ng)%DiaBio2d(i,j,5) + TP(ng)%o2_flux_sc_no(i,j) * imask
     DIAGS(ng)%DiaBio2d(i,j,6)=DIAGS(ng)%DiaBio2d(i,j,6) + TP(ng)%o2_flux_csurf(i,j) * imask
     DIAGS(ng)%DiaBio2d(i,j,7)=DIAGS(ng)%DiaBio2d(i,j,7) + TP(ng)%co2_flux(i,j) * imask
     DIAGS(ng)%DiaBio2d(i,j,8)=DIAGS(ng)%DiaBio2d(i,j,8) + TP(ng)%o2_flux(i,j) * imask
     DIAGS(ng)%DiaBio2d(i,j,9)=DIAGS(ng)%DiaBio2d(i,j,9) + TP(ng)%pco2_csurf(i,j) * imask
   END DO
   END DO
.....

HCJung-jbnu / ROMS-TOPAZ

ROMS-TOPAZ 진단변수 작게 평균되서 나오는 문제 #31