Changeset 12489 for NEMO/trunk/src/TOP/PISCES/SED
- Timestamp:
- 2020-02-28T16:55:11+01:00 (4 years ago)
- Location:
- NEMO/trunk/src/TOP/PISCES/SED
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/TOP/PISCES/SED/oce_sed.F90
r12377 r12489 18 18 USE dom_oce , ONLY : mbkt => mbkt !: vertical index of the bottom last T- ocean level 19 19 USE dom_oce , ONLY : tmask => tmask !: land/ocean mask at t-points 20 USE dom_oce , ONLY : r dt => rdt!: time step for the dynamics20 USE dom_oce , ONLY : rn_Dt => rn_Dt !: time step for the dynamics 21 21 USE dom_oce , ONLY : nyear => nyear !: Current year 22 22 USE dom_oce , ONLY : ndastp => ndastp !: time step date in year/month/day aammjj … … 50 50 USE p4zche, ONLY : sio3eq => sio3eq !: Chemical constants 51 51 USE p4zbc, ONLY : dust => dust 52 USE trc , ONLY : r 2dttrc => r2dttrc52 USE trc , ONLY : rDt_trc => rDt_trc 53 53 54 54 END MODULE oce_sed -
NEMO/trunk/src/TOP/PISCES/SED/seddta.F90
r12377 r12489 75 75 IF( kt == nitsed000 ) THEN 76 76 IF (lwp) WRITE(numsed,*) ' sed_dta : Sediment fields' 77 dtsed = r 2dttrc77 dtsed = rDt_trc 78 78 rsecday = 60.* 60. * 24. 79 79 ! conv2 = 1.0e+3 / ( 1.0e+4 * rsecday * 30. ) … … 103 103 DO_2D_11_11 104 104 ikt = mbkt(ji,jj) 105 zdep = e3t(ji,jj,ikt,Kmm) / r 2dttrc105 zdep = e3t(ji,jj,ikt,Kmm) / rDt_trc 106 106 zwsbio4(ji,jj) = MIN( 0.99 * zdep, wsbio4(ji,jj,ikt) / rday ) 107 107 zwsbio3(ji,jj) = MIN( 0.99 * zdep, wsbio3(ji,jj,ikt) / rday ) -
NEMO/trunk/src/TOP/PISCES/SED/sedini.F90
r12377 r12489 488 488 489 489 jpksedm1 = jpksed - 1 490 dtsed = r 2dttrc490 dtsed = rDt_trc 491 491 492 492 READ ( numnamsed_ref, nam_trased, IOSTAT = ios, ERR = 905) -
NEMO/trunk/src/TOP/PISCES/SED/sedrst.F90
r12377 r12489 330 330 !! In both those options, the exact duration of the experiment 331 331 !! since the beginning (cumulated duration of all previous restart runs) 332 !! is not stored in the restart and is assumed to be (nittrc000-1)*r dt.332 !! is not stored in the restart and is assumed to be (nittrc000-1)*rn_Dt. 333 333 !! This is valid is the time step has remained constant. 334 334 !! … … 381 381 ELSE 382 382 ndastp = ndate0 - 1 ! ndate0 read in the namelist in dom_nam 383 adatrj = ( REAL( nittrc000-1, wp ) * r dt ) / rday383 adatrj = ( REAL( nittrc000-1, wp ) * rn_Dt ) / rday 384 384 ! note this is wrong if time step has changed during run 385 385 ENDIF -
NEMO/trunk/src/TOP/PISCES/SED/sedstp.F90
r12377 r12489 55 55 IF(ln_sediment_offline) CALL trc_dmp_sed ( kt, Kbb, Kmm, Krhs ) 56 56 57 dtsed = r 2dttrc57 dtsed = rDt_trc 58 58 ! dtsed2 = dtsed 59 59 IF (kt /= nitsed000) THEN -
NEMO/trunk/src/TOP/PISCES/SED/sedwri.F90
r10222 r12489 94 94 DO ji = 1, jpoce 95 95 zflx(ji,jw) = ( pwcp(ji,1,jw) - pwcp_dta(ji,jw) ) & 96 & * 1.e3 / 1.e2 * dzkbot(ji) / r 2dttrc96 & * 1.e3 / 1.e2 * dzkbot(ji) / rDt_trc 97 97 ENDDO 98 98 ENDDO … … 100 100 ! Calculation of accumulation rate per dt 101 101 DO js = 1, jpsol 102 zrate = 1.0 / ( denssol * por1(jpksed) ) / r 2dttrc102 zrate = 1.0 / ( denssol * por1(jpksed) ) / rDt_trc 103 103 DO ji = 1, jpoce 104 104 zflx(ji,jpwatp1) = zflx(ji,jpwatp1) + ( tosed(ji,js) - fromsed(ji,js) ) * zrate
Note: See TracChangeset
for help on using the changeset viewer.