ROMS-FENNEL river runoff 분석

[HCJung-jbnu]

단위

• ROMS-FENNEL 변수 단위: mMol/m3 <-> RUNOFF 단위: mMol/m3
• TOPAZ 변수 단위 Mol/kg <-> river runoff 단위 Mol/m3
• ROMS-FENNEL을 분석해 보면 생지화학모델 내부 단위와 같은 단위의 RUNOFF를 넣어줌.
• 그러므로 TOPAZ도 RIVER RUNOFF 만들 때 내부 변수 단위인 Mol/kg으로 환산해서 넣어주는 것이 맞음.

~/ROMS/Module/mod_source.F

!
!  Allocate structure.
!
      IF (ng.eq.1) allocate ( SOURCES(Ngrids) )

[HCJung-jbnu]

ROMS-FENNEL RIVER RUNOFF 데이터 READ

• 두 파일에서 source forcing data read & time interp.

• get_data.F

  !
  !  Tracer Sources/Sinks.
  !
    DO i=1,NT(ng)
      IF (LtracerSrc(i,ng)) THEN
        CALL get_ngfld (ng, iNLM, idRtrc(i), SSF(ng)%ncid,            &
  #  if defined PIO_LIB && defined DISTRIBUTE
   &                    SSF(ng)%pioFile,                              &
  #  endif
   &                    1, SSF(ng), recordless, update(1),            &
   &                    1, Nsrc(ng), N(ng), 2, 1, Nsrc(ng), N(ng),    &
   &                    SOURCES(ng) % TsrcG(:,:,:,i))
        IF (FoundError(exit_flag, NoError, __LINE__, MyFile)) RETURN
      END IF
    END DO

• set_data.F

  !
  !  Tracer Sources/Sinks.
  !
      DO itrc=1,NT(ng)
        IF (LtracerSrc(itrc,ng)) THEN
          CALL set_ngfld (ng, iNLM, idRtrc(itrc), 1, Nsrc(ng), N(ng), &
   &                      1, Nsrc(ng), N(ng),                         &
   &                      SOURCES(ng) % TsrcG(:,:,:,itrc),            &
   &                      SOURCES(ng) % Tsrc(:,:,itrc),               &
   &                      update)
          IF (FoundError(exit_flag, NoError, __LINE__, MyFile)) RETURN
        END IF
      END DO

[HCJung-jbnu]

• ~/Nonlinear/pre_step3d.F : advection과 함께 river runoff(source) 적용

  .....
  !
  !  Apply tracers point sources to the horizontal advection terms,
  !  if any.
  !
  !    Dsrc(is) = 0,  flow across grid cell u-face (positive or negative)
  !    Dsrc(is) = 1,  flow across grid cell v-face (positive or negative)
  !
        IF (LuvSrc(ng)) THEN
          DO is=1,Nsrc(ng)
            Isrc=SOURCES(ng)%Isrc(is)
            Jsrc=SOURCES(ng)%Jsrc(is)
            IF (((Istr.le.Isrc).and.(Isrc.le.Iend+1)).and.            &
   &            ((Jstr.le.Jsrc).and.(Jsrc.le.Jend+1))) THEN
              IF (INT(SOURCES(ng)%Dsrc(is)).eq.0) THEN
                IF (LtracerSrc(itrc,ng)) THEN
                  FX(Isrc,Jsrc)=Huon(Isrc,Jsrc,k)*                    &
   &                            SOURCES(ng)%Tsrc(is,k,itrc)
                ELSE
                  FX(Isrc,Jsrc)=0.0_r8
                END IF
              ELSE IF (INT(SOURCES(ng)%Dsrc(is)).eq.1) THEN
                IF (LtracerSrc(itrc,ng)) THEN
                  FE(Isrc,Jsrc)=Hvom(Isrc,Jsrc,k)*                    &
   &                            SOURCES(ng)%Tsrc(is,k,itrc)
                ELSE
                   FE(Isrc,Jsrc)=0.0_r8
                END IF
              END IF
            END IF
          END DO
        END IF
  .....