Changeset 14306
- Timestamp:
- 2021-01-15T13:35:15+01:00 (3 years ago)
- 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 101 101 rcorgr = 1.E-4 ! Reactivity for refractory POC [an-1] 102 102 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] 104 104 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] 106 106 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] 108 108 xksedo2 = 1.E-6 ! Half-saturation constant for oxic remin [mol/l] 109 109 xksedno3 = 10.E-6 ! Half-saturation constant for denitrification [mol/l] 110 xksedfeo = 0. 6! Half-saturation constant for iron remin [%]110 xksedfeo = 0.007 ! Half-saturation constant for iron remin [%] 111 111 xksedso4 = 2E-3 ! Half-saturation constant for SO4 remin [mol/l] 112 112 / … … 130 130 cn_sedrst_outdir = "." ! directory to which to write output sediment restarts 131 131 / 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 158 158 ! zf1 = MIN(0.98, MAX(0., zf1 ) ) 159 159 zf1 = 0.7 160 zf2 = 0.0 4160 zf2 = 0.02 161 161 zf0 = 1.0 - zf2 - zf1 162 162 rainrm_dta(ji,jspoc) = ( zdtap(ji) + zdtag(ji) ) * 1e-4 * zf0 … … 164 164 rainrm_dta(ji,jspor) = ( zdtap(ji) + zdtag(ji) ) * 1e-4 * zf2 165 165 END DO 166 166 167 ! Sinking fluxes for Calcite in mol.m-2.s-1 ; conversion in mol.cm-2.s-1 167 168 CALL pack_arr ( jpoce, rainrm_dta(1:jpoce,jscal), trc_data(1:jpi,1:jpj,14), iarroce(1:jpoce) ) … … 177 178 rainrm_dta(1:jpoce,jsclay) = rainrm_dta(1:jpoce,jsclay) * conv2 / mol_wgt(jsclay) & 178 179 & + 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 181 183 ! rainrm_dta(1:jpoce,jsclay) = 1.0E-4 * conv2 / mol_wgt(jsclay) 182 184 -
NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/SED/sedini.F90
r14276 r14306 40 40 41 41 REAL(wp), PUBLIC :: & 42 redO2 = 1 72. , & !: Redfield coef for Oxygen42 redO2 = 133. , & !: Redfield coef for Oxygen 43 43 redNo3 = 16. , & !: Redfield coef for Nitrate 44 44 redPo4 = 1. , & !: Redfield coef for Phosphate … … 70 70 71 71 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 / 73 73 74 74 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.5 E-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 / 76 76 77 77 … … 217 217 ! -------------------------------- 218 218 IF (ln_sediment_offline) THEN 219 CALL trc_dta_ini(jptra)220 219 CALL trc_dmp_sed_ini 221 220 ENDIF -
NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/SED/sedwri.F90
r14276 r14306 103 103 ! Calculation of fluxes g/cm2/s 104 104 DO js = 1, jpsol 105 zrate = 1.0 / ( por1(jpksed) ) /r2dttrc105 zrate = 1.0 / r2dttrc 106 106 DO ji = 1, jpoce 107 107 zflx(ji,jpwat+js) = zflx(ji,jpwat+js) + ( tosed(ji,js) - fromsed(ji,js) ) * zrate … … 111 111 ! Calculation of accumulation rate per dt 112 112 DO js = 1, jpsol 113 zrate = 1.0 / ( por1(jpksed) ) /r2dttrc113 zrate = 1.0 / r2dttrc 114 114 DO ji = 1, jpoce 115 115 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 252 252 ! 253 253 ! deallocate data structure if data are not used for damping 254 #if ! defined key_sed_off 254 255 IF( .NOT.ln_trcdmp .AND. .NOT.ln_trcdmp_clo ) THEN 255 256 IF(lwp) WRITE(numout,*) 'trc_ini_state: deallocate data arrays as they are only used to initialize the run' … … 258 259 ! 259 260 ENDIF 261 #endif 260 262 ENDIF 261 263 END DO -
NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/trcnam.F90
r14276 r14306 62 62 ! 63 63 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 64 67 ! 65 68 !
Note: See TracChangeset
for help on using the changeset viewer.