ROMS-TOPAZ surface tracer flux 변수들 출력

[HCJung-jbnu]

• ALK, LITH, NH4, NO3, FED Dry/Wet deposition output으로 뽑기

[HCJung-jbnu]

• varinfo_topaz.dat

  
  'flux_FED'
  'FED surface flux'
  'mol/m2/s'
  'flux_FED, scalar, series'
  'ocean_time'
  'iDbio2(ifFED)'
  'r2dvar'
  1.0d0

'flux_LITH' 'LITH surface flux' 'g/m2/s' 'flux_LITH, scalar, series' 'ocean_time' 'iDbio2(ifLTH)' 'r2dvar' 1.0d0

'flux_NO3' 'NO3 surface flux' 'mol/m2/s' 'flux_NO3, scalar, series' 'ocean_time' 'iDbio2(ifNO3)' 'r2dvar' 1.0d0

'flux_NH4' 'NH4 surface flux' 'mol/m2/s' 'flux_NH4, scalar, series' 'ocean_time' 'iDbio2(ifNH4)' 'r2dvar' 1.0d0

- bio_topaz.in

Dout(ifFED) == T !FED surface flux Dout(ifLTH) == T !LITH surface flux Dout(ifNO3) == T !NO3 surface flux Dout(ifNH4) == T !NH4 surface flux

- topaz_inp.h

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

- topaz_mod.h

...... INTEGER:: ifFED !"FED surface flux" "mol/m2/sec" INTEGER:: ifLTH !"LITH surface flux" "g/m2/sec" INTEGER:: ifNO3 !"NO3 surface flux" "mol/m2/sec" INTEGER:: ifNH4 !"NH4 surface flux" "mol/m2/sec" ...... NDbio2d=13 ...... ifFED=ic+10 !"FED surface flux" "mol/m2/sec" ifLTH=ic+11!"LITH surface flux" "g/m2/sec" ifNO3=ic+12!"NO3 surface flux" "mol/m2/sec" ifNH4=ic+13!"NH4 surface flux" "mol/m2/sec" ......

- topaz_var.h

        CASE ('iDbio2(ifFED)') !10
          iDbio2(ifFED)=varid
        CASE ('iDbio2(ifLTH)') !11
          iDbio2(ifLTH)=varid
        CASE ('iDbio2(ifNO3)') !12
          iDbio2(ifNO3)=varid
        CASE ('iDbio2(ifNH4)') !13
          iDbio2(ifNH4)=varid

- topaz.h

       DIAGS(ng)%DiaBio2d(i,j,10)=DIAGS(ng)%DiaBio2d(i,j,10) + TP(ng)%flux_fed(i,j) * imask
       DIAGS(ng)%DiaBio2d(i,j,11)=DIAGS(ng)%DiaBio2d(i,j,11) + TP(ng)%flux_lith(i,j) * imask
       DIAGS(ng)%DiaBio2d(i,j,12)=DIAGS(ng)%DiaBio2d(i,j,12) + (TP(ng)%flux_no3_wet(i,j) + TP(ng)%flux_no3_dry(i,j) ) * imask
       DIAGS(ng)%DiaBio2d(i,j,13)=DIAGS(ng)%DiaBio2d(i,j,13) + (TP(ng)%flux_nh4_wet(i,j) + TP(ng)%flux_nh4_dry(i,j) ) * imask

[HCJung-jbnu]

[HCJung-jbnu]