HyunChaeJung commented 1 year ago

surface flux 관련 define 옵션

ANA_BPFLUX if analytical bottom passive tracers fluxes
- define ANA_BPFLUX => dont need tracers bottom fluxes
ANA_SPFLUX if analytical surface passive tracers fluxes
- define ANA_SPFLUX => dont need tracers surface fluxes

~/ROMS/Include/globaldefs.h

....
/* add TOPAZ hcjung 20230407 */
#if defined TOPAZ
# undef ANA_SPFLUX
# undef ANA_BPFLUX
#endif
....

HyunChaeJung commented 1 year ago

TOPAZ surface flux data

atmospheric pressure
- ROMS에서 연결 가능할 듯...
short wave fraction
- ROMS에서 연결 가능할 듯...
surface dry/wet dep : no3_wet/dry, nh4_wet/dry, lith, fed = depflux_total.mean.1860.nc
surface air-sea exchnage coef: o2_alpha, o2_sc_num, o2_csurf, co2_alpha, co2_sc_num, co2_csurf
= ocean_topaz_airsea_flux.res.nc
surface air-sea exchnage coef은 초기값을 전달하는 것이기 때문에 init 과정 접합 부분에 있음.

HyunChaeJung commented 1 year ago

ROMS 수정

TOPAZ define 시, ANA_SPFLUX, ANA_BPFLUX 삭제 -> globaldefs.h 수정
get_data.F set_data.F에서 define TOPAZ시, TOPAZ STFLUX 읽고, idTbot, idTsur set/get 삭제
- TOPAZ 시 BOTTOM FULX와 다른 SURFACE FLUX는 없음.
mod_ncparam.F90, mod_forces.F90 에 topaz stflux 등록
- mod_biology의 column_physics로 전달해야 함....
~/ROMS/Nonlinear/get_data.F : . -> ocean_var_4topaz 로 전달 : varinfo_topaz.dat에 변수 등록 -> get_data.F에서 READ
```
SUBROUTINE get_data (ng)
....
# if defined BIOLOGY || defined SEDIMENT || defined T_PASSIVE
```

!!hcjung add TOPAZ 20230412

ifdef TOPAZ

ifndef ANA_SPFLUX

! ! Soluble Fe ! CALL get_2dfld (ng, iNLM, idSolubleFe, FRCncid(idSolubleFe,ng), &

if defined PIO_LIB && defined DISTRIBUTE

 &                 FRCpioFile(idSolubleFe,ng),                       &

endif

 &                 nFfiles(ng), FRC(1,ng), update(1),              &
 &                 LBi, UBi, LBj, UBj, 2, 1,                       &

ifdef MASKING

 &                 GRID(ng) % rmask,                               &

endif

 &                 FORCES(ng) % soluble_feG)
   IF (FoundError(exit_flag, NoError, __LINE__, MyFile)) RETURN

! ! Mineral Fe ! CALL get_2dfld (ng, iNLM, idMineralFe, FRCncid(idMineralFe,ng), &

if defined PIO_LIB && defined DISTRIBUTE

 &                 FRCpioFile(idMineralFe,ng),                       &

endif

 &                 nFfiles(ng), FRC(1,ng), update(1),              &
 &                 LBi, UBi, LBj, UBj, 2, 1,                       &

ifdef MASKING

 &                 GRID(ng) % rmask,                               &

endif

 &                 FORCES(ng) % mineral_feG)
   IF (FoundError(exit_flag, NoError, __LINE__, MyFile)) RETURN

! ! Dry deposition NO3 ! CALL get_2dfld (ng, iNLM, idDryNO3, FRCncid(idDryNO3,ng), &

if defined PIO_LIB && defined DISTRIBUTE

 &                 FRCpioFile(idDryNO3,ng),                       &

endif

 &                 nFfiles(ng), FRC(1,ng), update(1),              &
 &                 LBi, UBi, LBj, UBj, 2, 1,                       &

ifdef MASKING

 &                 GRID(ng) % rmask,                               &

endif

 &                 FORCES(ng) % drydep_no3G)
   IF (FoundError(exit_flag, NoError, __LINE__, MyFile)) RETURN

! ! Wet deposition NO3 ! CALL get_2dfld (ng, iNLM, idWetNO3, FRCncid(idWetNO3,ng), &

if defined PIO_LIB && defined DISTRIBUTE

 &                 FRCpioFile(idWetNO3,ng),                       &

endif

 &                 nFfiles(ng), FRC(1,ng), update(1),              &
 &                 LBi, UBi, LBj, UBj, 2, 1,                       &

ifdef MASKING

 &                 GRID(ng) % rmask,                               &

endif

 &                 FORCES(ng) % wetdep_no3G)
   IF (FoundError(exit_flag, NoError, __LINE__, MyFile)) RETURN

! ! Dry deposition NH4 ! CALL get_2dfld (ng, iNLM, idDryNH4, FRCncid(idDryNH4,ng), &

if defined PIO_LIB && defined DISTRIBUTE

 &                 FRCpioFile(idDryNH4,ng),                       &

endif

 &                 nFfiles(ng), FRC(1,ng), update(1),              &
 &                 LBi, UBi, LBj, UBj, 2, 1,                       &

ifdef MASKING

 &                 GRID(ng) % rmask,                               &

endif

 &                 FORCES(ng) % drydep_nh4G)
   IF (FoundError(exit_flag, NoError, __LINE__, MyFile)) RETURN

! ! Wet deposition NH4 ! CALL get_2dfld (ng, iNLM, idWetNH4, FRCncid(idWetNH4,ng), &

if defined PIO_LIB && defined DISTRIBUTE

 &                 FRCpioFile(idWetNH4,ng),                       &

endif

 &                 nFfiles(ng), FRC(1,ng), update(1),              &
 &                 LBi, UBi, LBj, UBj, 2, 1,                       &

ifdef MASKING

 &                 GRID(ng) % rmask,                               &

endif

 &                 FORCES(ng) % wetdep_nh4G)
   IF (FoundError(exit_flag, NoError, __LINE__, MyFile)) RETURN

endif

else

ifndef ANA_SPFLUX

...........

- ~/ROMS/Modules/mod_forcing.F : stf dry/set dep 변수 등록

..... !!hcjung add TOPAZ 20230407

if defined TOPAZ

      real(r8), pointer :: soluble_fe(:,:)
      real(r8), pointer :: soluble_feG(:,:,:)
      real(r8), pointer :: mineral_fe(:,:)
      real(r8), pointer :: mineral_feG(:,:,:)
      real(r8), pointer :: drydep_no3(:,:)
      real(r8), pointer :: drydep_no3G(:,:,:)
      real(r8), pointer :: wetdep_no3(:,:)
      real(r8), pointer :: wetdep_no3G(:,:,:)
      real(r8), pointer :: drydep_nh4(:,:)
      real(r8), pointer :: drydep_nh4G(:,:,:)
      real(r8), pointer :: wetdep_nh4(:,:)
      real(r8), pointer :: wetdep_nh4G(:,:,:)

endif

..... !!hcjung add TOPAZ

if defined TOPAZ

  allocate ( FORCES(ng) % soluble_fe(LBi:UBi,LBj:UBj) )
  Dmem(ng)=Dmem(ng)+size2d
  allocate ( FORCES(ng) % soluble_feG(LBi:UBi,LBj:UBj,2) )
  Dmem(ng)=Dmem(ng)+size2d*2
  allocate ( FORCES(ng) % mineral_fe(LBi:UBi,LBj:UBj) )
  Dmem(ng)=Dmem(ng)+size2d
  allocate ( FORCES(ng) % mineral_feG(LBi:UBi,LBj:UBj,2) )
  Dmem(ng)=Dmem(ng)+size2d*2
  allocate ( FORCES(ng) % drydep_no3(LBi:UBi,LBj:UBj) )
  Dmem(ng)=Dmem(ng)+size2d
  allocate ( FORCES(ng) % drydep_no3G(LBi:UBi,LBj:UBj,2) )
  Dmem(ng)=Dmem(ng)+size2d*2
  allocate ( FORCES(ng) % wetdep_no3(LBi:UBi,LBj:UBj) )
  Dmem(ng)=Dmem(ng)+size2d
  allocate ( FORCES(ng) % wetdep_no3G(LBi:UBi,LBj:UBj,2) )
  Dmem(ng)=Dmem(ng)+size2d*2
  allocate ( FORCES(ng) % drydep_nh4(LBi:UBi,LBj:UBj) )
  Dmem(ng)=Dmem(ng)+size2d
  allocate ( FORCES(ng) % drydep_nh4G(LBi:UBi,LBj:UBj,2) )
  Dmem(ng)=Dmem(ng)+size2d*2
  allocate ( FORCES(ng) % wetdep_nh4(LBi:UBi,LBj:UBj) )
  Dmem(ng)=Dmem(ng)+size2d
  allocate ( FORCES(ng) % wetdep_nh4G(LBi:UBi,LBj:UBj,2) )
  Dmem(ng)=Dmem(ng)+size2d*2

endif

...... !!hcjung add TOPAZ 20230407

ifdef TOPAZ

        FORCES(ng) % soluble_fe(i,j) = IniVal
        FORCES(ng) % soluble_feG(i,j,1) = IniVal
        FORCES(ng) % soluble_feG(i,j,2) = IniVal
        FORCES(ng) % mineral_fe(i,j) = IniVal
        FORCES(ng) % mineral_feG(i,j,1) = IniVal
        FORCES(ng) % mineral_feG(i,j,2) = IniVal
        FORCES(ng) % drydep_no3(i,j) = IniVal
        FORCES(ng) % drydep_no3G(i,j,1) = IniVal
        FORCES(ng) % drydep_no3G(i,j,2) = IniVal
        FORCES(ng) % wetdep_no3(i,j) = IniVal
        FORCES(ng) % wetdep_no3G(i,j,1) = IniVal
        FORCES(ng) % wetdep_no3G(i,j,2) = IniVal
        FORCES(ng) % drydep_nh4(i,j) = IniVal
        FORCES(ng) % drydep_nh4G(i,j,1) = IniVal
        FORCES(ng) % drydep_nh4G(i,j,2) = IniVal
        FORCES(ng) % wetdep_nh4(i,j) = IniVal
        FORCES(ng) % wetdep_nh4G(i,j,1) = IniVal
        FORCES(ng) % wetdep_nh4G(i,j,2) = IniVal

endif

.......

- ~/ROMS/Modules/mod_ncparam.F : varinfo_topaz.dat -> topaz_var.h -> 변수 indices 등록

.... !!hcjung add topaz 20230407

ifdef TOPAZ

  integer  :: idSolubleFe   ! dust deposition to iron
  integer  :: idMineralFe   ! dust deposition to lith
  integer  :: idDryNO3      ! dry dep no3
  integer  :: idWetNO3      ! wet dep no3
  integer  :: idDryNH4      ! dry dep nh4
  integer  :: idWetNH4      ! wet dep nh4

endif

..... !hcjung add topaz 20230407

ifdef TOPAZ

      CASE ('idSolubleFe')
          idSolubleFe=varid
      CASE ('idMineralFe')
          idMineralFe=varid
      CASE ('idDryNO3')
          idDryNO3=varid
      CASE ('idWetNO3')
          idWetNO3=varid
      CASE ('idDryNH4')
          idDryNH4=varid
      CASE ('idWetNH4')
          idWetNH4=varid

endif

.....

- ~/ROMS/Nonlinear/set_data.F

SUBROUTINE set_data_tile ..... !!hcjung add TOPAZ 20230407

if defined TOPAZ

! ! Soluble Fe ! CALL set_2dfld_tile (ng, tile, iNLM, idSolubleFe, & & LBi, UBi, LBj, UBj, & & FORCES(ng)%soluble_feG, & & FORCES(ng)%soluble_fe, & & update) IF (FoundError(exit_flag, NoError, LINE, MyFile)) RETURN

! ! Mineral Fe ! CALL set_2dfld_tile (ng, tile, iNLM, idMineralFe, & & LBi, UBi, LBj, UBj, & & FORCES(ng)%mineral_feG, & & FORCES(ng)%mineral_fe, & & update) IF (FoundError(exit_flag, NoError, LINE, MyFile)) RETURN .....

HCJung-jbnu / ROMS-TOPAZ

ROMS-TOPAZ surface flux 처방 #8

surface flux 관련 define 옵션

define ANA_BPFLUX => dont need tracers bottom fluxes

define ANA_SPFLUX => dont need tracers surface fluxes

TOPAZ surface flux data

ROMS 수정

ifdef TOPAZ

ifndef ANA_SPFLUX

if defined PIO_LIB && defined DISTRIBUTE

endif

ifdef MASKING

endif

if defined PIO_LIB && defined DISTRIBUTE

endif

ifdef MASKING

endif

if defined PIO_LIB && defined DISTRIBUTE

endif

ifdef MASKING

endif

if defined PIO_LIB && defined DISTRIBUTE

endif

ifdef MASKING

endif

if defined PIO_LIB && defined DISTRIBUTE

endif

ifdef MASKING

endif

if defined PIO_LIB && defined DISTRIBUTE

endif

ifdef MASKING

endif

endif

else

ifndef ANA_SPFLUX

if defined TOPAZ

endif

if defined TOPAZ

endif

ifdef TOPAZ

endif

ifdef TOPAZ

endif

ifdef TOPAZ

endif

if defined TOPAZ