New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 14306 – NEMO

Changeset 14306


Ignore:
Timestamp:
2021-01-15T13:35:15+01:00 (3 years ago)
Author:
aumont
Message:

small bug fixes and improvements in the diagenetic module

Location:
NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/cfgs/SHARED/namelist_sediment_ref

    r14276 r14306  
    101101   rcorgr   =  1.E-4   ! Reactivity for refractory POC [an-1] 
    102102   rcnh4    =  1E7     ! Reactivity for O2/NH4 [l.mol-1.an-1] 
    103    rch2s    =  1.6E5   ! Reactivity for O2/H2S [l.mol-1.an-1] 
     103   rch2s    =  2E8     ! Reactivity for O2/H2S [l.mol-1.an-1] 
    104104   rcfe2    =  1E8     ! Reactivity for O2/Fe2+ [l.mol-1.an-1] 
    105    rcfeh2s  =  8E3     ! Reactivity for FEOH/H2S [l.mol-1.an-1] 
     105   rcfeh2s  =  1E4     ! Reactivity for FEOH/H2S [l.mol-1.an-1] 
    106106   rcfes    =  1E6     ! Reactivity for FE2+/H2S [l.mol-1.an-1] 
    107    rcfeso   =  3E5     ! Reactivity for FES/O2 [l.mol-1.an-1] 
     107   rcfeso   =  2E7     ! Reactivity for FES/O2 [l.mol-1.an-1] 
    108108   xksedo2  =  1.E-6   ! Half-saturation constant for oxic remin [mol/l] 
    109109   xksedno3 =  10.E-6  ! Half-saturation constant for denitrification [mol/l]  
    110    xksedfeo =  0.   ! Half-saturation constant for iron remin [%] 
     110   xksedfeo =  0.007   ! Half-saturation constant for iron remin [%] 
    111111   xksedso4 =  2E-3    ! Half-saturation constant for SO4 remin [mol/l] 
    112112/ 
     
    130130   cn_sedrst_outdir = "."          !  directory to which to write output sediment restarts 
    131131/ 
     132!----------------------------------------------------------------------- 
     133&nam_output     !   parameters for outputing the sediment module 
     134!----------------------------------------------------------------------- 
     135   ldefsedpis_avg = .true.            !  write averaged output variables 
     136   cn_sedwri_out  = "output_sed.nc"   !  name of the input restart file name of the sediment module 
     137   nrpfsedpis_avg = 0                 ! Frequency of the averaged outputs 
     138   nwrtsedpis_avg = 24                ! Frequency of the averaged outputs 
     139   ntssedpis_avg  =  1                ! ??? 
     140/ 
  • NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/SED/seddta.F90

    r14276 r14306  
    158158!        zf1    = MIN(0.98, MAX(0., zf1 ) ) 
    159159         zf1    = 0.7 
    160          zf2    = 0.04 
     160         zf2    = 0.02 
    161161         zf0    = 1.0 - zf2 - zf1 
    162162         rainrm_dta(ji,jspoc) =   ( zdtap(ji) +  zdtag(ji) ) * 1e-4 * zf0 
     
    164164         rainrm_dta(ji,jspor) =   ( zdtap(ji) +  zdtag(ji) ) * 1e-4 * zf2 
    165165      END DO 
     166 
    166167      !  Sinking fluxes for Calcite in mol.m-2.s-1 ; conversion in mol.cm-2.s-1 
    167168      CALL pack_arr ( jpoce,  rainrm_dta(1:jpoce,jscal), trc_data(1:jpi,1:jpj,14), iarroce(1:jpoce) ) 
     
    177178      rainrm_dta(1:jpoce,jsclay) = rainrm_dta(1:jpoce,jsclay) * conv2 / mol_wgt(jsclay)   & 
    178179      &                            + wacc(1:jpoce) * por1(2) * denssol / mol_wgt(jsclay) / ( rsecday * 365.0 ) 
    179       rainrm_dta(1:jpoce,jsfeo)  = rainrm_dta(1:jpoce,jsclay) * mol_wgt(jsclay) / mol_wgt(jsfeo) * 0.035 * 0.5 
    180       rainrm_dta(1:jpoce,jsclay) = rainrm_dta(1:jpoce,jsclay) * 0.965 
     180      rainrm_dta(1:jpoce,jsfeo)  = rainrm_dta(1:jpoce,jsclay) * mol_wgt(jsclay) / mol_wgt(jsfeo) * 0.035 * 0.5 * 0.333 
     181      rainrm_dta(1:jpoce,jsclay) = rainrm_dta(1:jpoce,jsclay) * (1.0 - 0.035 * 0.5 * 0.333 ) 
     182 
    181183!    rainrm_dta(1:jpoce,jsclay) = 1.0E-4 * conv2 / mol_wgt(jsclay) 
    182184 
  • NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/SED/sedini.F90

    r14276 r14306  
    4040 
    4141   REAL(wp), PUBLIC    ::  & 
    42       redO2    =  172.  ,  &  !: Redfield coef for Oxygen 
     42      redO2    =  133.  ,  &  !: Redfield coef for Oxygen 
    4343      redNo3   =   16.  ,  &  !: Redfield coef for Nitrate 
    4444      redPo4   =    1.  ,  &  !: Redfield coef for Phosphate 
     
    7070 
    7171   REAL(wp), DIMENSION(jpwat), PUBLIC  :: diff1 
    72    DATA diff1/4.59E-6, 1.104E-5, 4.81E-6 , 9.78E-6, 3.58E-6, 4.01E-6, 9.8E-6, 9.73E-6, 5.0E-6, 3.31E-6 / 
     72   DATA diff1/4.59E-6, 1.104E-5, 4.81E-6 , 9.78E-6, 3.58E-6, 4.81E-6, 9.8E-6, 9.73E-6, 5.0E-6, 3.31E-6 / 
    7373 
    7474   REAL(wp), DIMENSION(jpwat), PUBLIC  :: diff2 
    75    DATA diff2/1.74E-7, 4.47E-7, 2.51E-7, 3.89E-7, 1.77E-7, 2.5E-7, 3.89E-7, 3.06E-7, 2.5E-7, 1.5E-7 / 
     75   DATA diff2/1.74E-7, 4.47E-7, 2.51E-7, 3.89E-7, 1.77E-7, 2.51E-7, 3.89E-7, 3.06E-7, 2.5E-7, 1.5E-7 / 
    7676 
    7777 
     
    217217      ! -------------------------------- 
    218218      IF (ln_sediment_offline) THEN 
    219          CALL trc_dta_ini(jptra) 
    220219         CALL trc_dmp_sed_ini 
    221220      ENDIF 
  • NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/SED/sedwri.F90

    r14276 r14306  
    103103      ! Calculation of fluxes g/cm2/s 
    104104      DO js = 1, jpsol 
    105          zrate =  1.0 / ( por1(jpksed) ) / r2dttrc 
     105         zrate =  1.0 / r2dttrc 
    106106         DO ji = 1, jpoce 
    107107            zflx(ji,jpwat+js) = zflx(ji,jpwat+js) + ( tosed(ji,js) - fromsed(ji,js) ) * zrate 
     
    111111      ! Calculation of accumulation rate per dt 
    112112      DO js = 1, jpsol 
    113          zrate =  1.0 / ( por1(jpksed) ) / r2dttrc 
     113         zrate =  1.0 / r2dttrc 
    114114         DO ji = 1, jpoce 
    115115            zflx(ji,jptrased+1) = zflx(ji,jptrased+1) + ( tosed(ji,js) - fromsed(ji,js) ) * zrate 
  • NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/trcini.F90

    r12537 r14306  
    252252                  ! 
    253253                  ! deallocate data structure if data are not used for damping 
     254#if ! defined key_sed_off 
    254255                  IF( .NOT.ln_trcdmp .AND. .NOT.ln_trcdmp_clo ) THEN 
    255256                     IF(lwp) WRITE(numout,*) 'trc_ini_state: deallocate data arrays as they are only used to initialize the run' 
     
    258259                     ! 
    259260                  ENDIF 
     261#endif 
    260262               ENDIF 
    261263            END DO 
  • NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/trcnam.F90

    r14276 r14306  
    6262      ! 
    6363      IF( ln_trcdmp .OR. ln_trcdmp_clo )   ln_trcdta = .TRUE.    ! damping : need to have clim data 
     64#if defined key_sed_off 
     65      ln_trcdta = .TRUE. 
     66#endif 
    6467      ! 
    6568      ! 
Note: See TracChangeset for help on using the changeset viewer.