- Timestamp:
- 2020-09-14T17:40:34+02:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11351_fldread_with_XIOS
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11351_fldread_with_XIOS
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEADext/AGRIF5 ^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL 8 9 # SETTE 10 ^/utils/CI/sette@13382 sette
-
- Property svn:externals
-
NEMO/branches/2019/dev_r11351_fldread_with_XIOS/src/TOP/PISCES/SED/seddta.F90
r10362 r13463 22 22 REAL(wp) :: conv2 ! [kg/m2/month]-->[g/cm2/s] ( 1 month has 30 days ) 23 23 24 !! * Substitutions 25 # include "do_loop_substitute.h90" 26 # include "domzgr_substitute.h90" 24 27 !! $Id$ 25 28 CONTAINS … … 29 32 !!--------------------------------------------------------------------------- 30 33 31 SUBROUTINE sed_dta( kt )34 SUBROUTINE sed_dta( kt, Kbb, Kmm ) 32 35 !!---------------------------------------------------------------------- 33 36 !! *** ROUTINE sed_dta *** … … 43 46 44 47 !! Arguments 45 INTEGER, INTENT(in) :: kt ! time-step 48 INTEGER, INTENT(in) :: kt ! time-step 49 INTEGER, INTENT(in) :: Kbb, Kmm ! time level indices 46 50 47 51 !! * Local declarations … … 72 76 IF( kt == nitsed000 ) THEN 73 77 IF (lwp) WRITE(numsed,*) ' sed_dta : Sediment fields' 74 dtsed = r 2dttrc78 dtsed = rDt_trc 75 79 rsecday = 60.* 60. * 24. 76 80 ! conv2 = 1.0e+3 / ( 1.0e+4 * rsecday * 30. ) … … 92 96 ! ----------------------------------------------------------- 93 97 IF (ln_sediment_offline) THEN 94 DO jj = 1, jpj 95 DO ji = 1, jpi 96 ikt = mbkt(ji,jj) 97 zwsbio4(ji,jj) = wsbio2 / rday 98 zwsbio3(ji,jj) = wsbio / rday 99 END DO 100 END DO 98 DO_2D( 1, 1, 1, 1 ) 99 ikt = mbkt(ji,jj) 100 zwsbio4(ji,jj) = wsbio2 / rday 101 zwsbio3(ji,jj) = wsbio / rday 102 END_2D 101 103 ELSE 102 DO jj = 1, jpj 103 DO ji = 1, jpi 104 ikt = mbkt(ji,jj) 105 zdep = e3t_n(ji,jj,ikt) / r2dttrc 106 zwsbio4(ji,jj) = MIN( 0.99 * zdep, wsbio4(ji,jj,ikt) / rday ) 107 zwsbio3(ji,jj) = MIN( 0.99 * zdep, wsbio3(ji,jj,ikt) / rday ) 108 END DO 109 END DO 104 DO_2D( 1, 1, 1, 1 ) 105 ikt = mbkt(ji,jj) 106 zdep = e3t(ji,jj,ikt,Kmm) / rDt_trc 107 zwsbio4(ji,jj) = MIN( 0.99 * zdep, wsbio4(ji,jj,ikt) / rday ) 108 zwsbio3(ji,jj) = MIN( 0.99 * zdep, wsbio3(ji,jj,ikt) / rday ) 109 END_2D 110 110 ENDIF 111 111 112 112 trc_data(:,:,:) = 0. 113 DO jj = 1,jpj 114 DO ji = 1, jpi 115 ikt = mbkt(ji,jj) 116 IF ( tmask(ji,jj,ikt) == 1 ) THEN 117 trc_data(ji,jj,1) = trb(ji,jj,ikt,jpsil) 118 trc_data(ji,jj,2) = trb(ji,jj,ikt,jpoxy) 119 trc_data(ji,jj,3) = trb(ji,jj,ikt,jpdic) 120 trc_data(ji,jj,4) = trb(ji,jj,ikt,jpno3) / 7.625 121 trc_data(ji,jj,5) = trb(ji,jj,ikt,jppo4) / 122. 122 trc_data(ji,jj,6) = trb(ji,jj,ikt,jptal) 123 trc_data(ji,jj,7) = trb(ji,jj,ikt,jpnh4) / 7.625 124 trc_data(ji,jj,8) = 0.0 125 trc_data(ji,jj,9) = 28.0E-3 126 trc_data(ji,jj,10) = trb(ji,jj,ikt,jpfer) 127 trc_data(ji,jj,11 ) = MIN(trb(ji,jj,ikt,jpgsi), 1E-4) * zwsbio4(ji,jj) * 1E3 128 trc_data(ji,jj,12 ) = MIN(trb(ji,jj,ikt,jppoc), 1E-4) * zwsbio3(ji,jj) * 1E3 129 trc_data(ji,jj,13 ) = MIN(trb(ji,jj,ikt,jpgoc), 1E-4) * zwsbio4(ji,jj) * 1E3 130 trc_data(ji,jj,14) = MIN(trb(ji,jj,ikt,jpcal), 1E-4) * zwsbio4(ji,jj) * 1E3 131 trc_data(ji,jj,15) = tsn(ji,jj,ikt,jp_tem) 132 trc_data(ji,jj,16) = tsn(ji,jj,ikt,jp_sal) 133 trc_data(ji,jj,17 ) = ( trb(ji,jj,ikt,jpsfe) * zwsbio3(ji,jj) + trb(ji,jj,ikt,jpbfe) & 134 & * zwsbio4(ji,jj) ) * 1E3 / ( trc_data(ji,jj,12 ) + trc_data(ji,jj,13 ) + rtrn ) 135 trc_data(ji,jj,17 ) = MIN(1E-3, trc_data(ji,jj,17 ) ) 136 ENDIF 137 ENDDO 138 ENDDO 113 DO_2D( 1, 1, 1, 1 ) 114 ikt = mbkt(ji,jj) 115 IF ( tmask(ji,jj,ikt) == 1 ) THEN 116 trc_data(ji,jj,1) = tr(ji,jj,ikt,jpsil,Kbb) 117 trc_data(ji,jj,2) = tr(ji,jj,ikt,jpoxy,Kbb) 118 trc_data(ji,jj,3) = tr(ji,jj,ikt,jpdic,Kbb) 119 trc_data(ji,jj,4) = tr(ji,jj,ikt,jpno3,Kbb) / 7.625 120 trc_data(ji,jj,5) = tr(ji,jj,ikt,jppo4,Kbb) / 122. 121 trc_data(ji,jj,6) = tr(ji,jj,ikt,jptal,Kbb) 122 trc_data(ji,jj,7) = tr(ji,jj,ikt,jpnh4,Kbb) / 7.625 123 trc_data(ji,jj,8) = 0.0 124 trc_data(ji,jj,9) = 28.0E-3 125 trc_data(ji,jj,10) = tr(ji,jj,ikt,jpfer,Kbb) 126 trc_data(ji,jj,11 ) = MIN(tr(ji,jj,ikt,jpgsi,Kbb), 1E-4) * zwsbio4(ji,jj) * 1E3 127 trc_data(ji,jj,12 ) = MIN(tr(ji,jj,ikt,jppoc,Kbb), 1E-4) * zwsbio3(ji,jj) * 1E3 128 trc_data(ji,jj,13 ) = MIN(tr(ji,jj,ikt,jpgoc,Kbb), 1E-4) * zwsbio4(ji,jj) * 1E3 129 trc_data(ji,jj,14) = MIN(tr(ji,jj,ikt,jpcal,Kbb), 1E-4) * zwsbio4(ji,jj) * 1E3 130 trc_data(ji,jj,15) = ts(ji,jj,ikt,jp_tem,Kmm) 131 trc_data(ji,jj,16) = ts(ji,jj,ikt,jp_sal,Kmm) 132 trc_data(ji,jj,17 ) = ( tr(ji,jj,ikt,jpsfe,Kbb) * zwsbio3(ji,jj) + tr(ji,jj,ikt,jpbfe,Kbb) & 133 & * zwsbio4(ji,jj) ) * 1E3 / ( trc_data(ji,jj,12 ) + trc_data(ji,jj,13 ) + rtrn ) 134 trc_data(ji,jj,17 ) = MIN(1E-3, trc_data(ji,jj,17 ) ) 135 ENDIF 136 END_2D 139 137 140 138 ! Pore water initial concentration [mol/l] in k=1 … … 167 165 CALL pack_arr ( jpoce, rainrm_dta(1:jpoce,jscal), trc_data(1:jpi,1:jpj,14), iarroce(1:jpoce) ) 168 166 rainrm_dta(1:jpoce,jscal) = rainrm_dta(1:jpoce,jscal) * 1e-4 169 ! vector temperature [ °C] and salinity167 ! vector temperature [�C] and salinity 170 168 CALL pack_arr ( jpoce, temp(1:jpoce), trc_data(1:jpi,1:jpj,15), iarroce(1:jpoce) ) 171 169 CALL pack_arr ( jpoce, salt(1:jpoce), trc_data(1:jpi,1:jpj,16), iarroce(1:jpoce) )
Note: See TracChangeset
for help on using the changeset viewer.