Changeset 6810 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM
- Timestamp:
- 2016-07-20T14:54:06+02:00 (8 years ago)
- Location:
- branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/IDTRA/trcsms_idtra.F90
r6509 r6810 133 133 DO ji = 2,jpim1 134 134 135 !IF (trn(ji,jj,jk,jn) > 0.0) THEN135 IF (trn(ji,jj,jk,jn) > 0.0) THEN 136 136 WTEMP = trn(ji,jj,jk,jn) * (1. - FDEC ) 137 137 tra(ji,jj,jk,jn) = (tra(ji,jj,jk,jn) - WTEMP/rdt ) * & 138 138 tmask(ji,jj,jk) 139 ! ENDIF 139 ELSE 140 trn(ji,jj,jk,jn) = 0.0 141 tra(ji,jj,jk,jn) = 0.0 142 ENDIF 140 143 ENDDO 141 144 ENDDO -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/sms_medusa.F90
r6715 r6810 260 260 !! 261 261 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: dust !: dust parameter 1 262 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zirondep !! Fe deposition 262 263 263 264 !!---------------------------------------------------------------------- … … 488 489 !* 2D fields of miscellaneous parameters 489 490 ALLOCATE( ocal_ccd(jpi,jpj) , dust(jpi,jpj) , & 491 & zirondep(jpi,jpj) , & 490 492 & riv_n(jpi,jpj) , & 491 493 & riv_si(jpi,jpj) , riv_c(jpi,jpj) , & -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90
r6806 r6810 2509 2509 !! aeolian iron deposition 2510 2510 if (jk.eq.1) then 2511 !! dust is in g Fe / m2 / month2512 !! ffetop is in mmol/ m3 / day2513 ffetop = (((dust(ji,jj) * 1.e3 * xfe_sol) / xfe_mass) / fthk) / 30.2511 !! zirondep is in mmol-Fe / m2 / day 2512 !! ffetop is in mmol-dissolved-Fe / m3 / day 2513 ffetop = zirondep(ji,jj) * xfe_sol / fthk 2514 2514 else 2515 2515 ffetop = 0.0 -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcsed_medusa.F90
r6744 r6810 98 98 99 99 !! AXY (10/02/09) 100 REAL(wp), DIMENSION(jpi,jpj) :: zsidep !! Si deposition101 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zirondep !! Fe deposition100 !! REAL(wp), DIMENSION(jpi,jpj) :: zsidep !! Si deposition 101 !! REAL(wp), DIMENSION(jpi,jpj,jpk) :: zirondep !! Fe deposition 102 102 REAL(wp) :: rfact2 103 103 … … 106 106 !! JPALM - 26-11-2015 -add iom_use for diagnostic 107 107 REAL(wp), POINTER, DIMENSION(:,: ) :: zw2d 108 !! JPALM -- 20-07-2016 -- 109 REAL(wp), PARAMETER :: Fe_dust_mratio = 0.035 !! Fe:dust mass ratio = 0.035 108 110 !!--------------------------------------------------------------------- 109 111 !! … … 144 146 !! AXY (10/02/09) 145 147 !!IF( (jnt == 1) .and. (bdustfer) ) CALL trc_sed_medusa_sbc( kt ) 148 146 149 !! JPALM -- 31-03-2016 -- rewrite trc_sed_medusa_sbc. 147 150 !! IF (kt == nittrc000 ) CALL trc_sed_medusa_sbc 151 152 !! JPALM -- 20-07-2016 -- change the dust dep conversion. 153 !! now all dust dep forcings and coupling are in SI units : 154 !! kg-dust/m2/s instead of g-Fe/m2/month. 155 !! and we convert them here into mmol-Fe/m2/day 156 !! what will be changed into mmol-Fe/m3/d in trcbio_med 157 !! 148 158 IF( bdustfer ) THEN 149 159 IF( kt == nittrc000 .OR. ( kt /= nittrc000 .AND. ntimes_dust > 1 ) ) THEN … … 157 167 ENDIF 158 168 !! 159 160 !! 161 zirondep(:,:,:) = 0.e0 !! Initialisation of deposition variables 162 zsidep (:,:) = 0.e0 163 !! 164 !! Iron and Si deposition at the surface 169 170 zirondep(:,:) = 0.e0 !! Initialisation of deposition variables 171 zirondep(:,:) = dust(:,:) * Fe_dust_mratio / xfe_mass * 1.e6 * 86400 !! mmol-Fe/m2/d 172 173 !! 174 !! JPALM -- 20-07-2016 -- previous Fe and Si deposition. 175 !! not used --> commented. 176 !! but may want to use it later on. 177 !! 178 !! Iron and Si deposition at the surface 165 179 !! ------------------------------------- 166 180 !! 167 DO jj = 1, jpj168 DO ji = 1, jpi169 zirondep(ji,jj,1) = (dustsolub * dust(ji,jj) / (55.85 * rmtss) + 3.e-10 / ryyss) &170 & * rfact2 / fse3t(ji,jj,1)171 zsidep (ji,jj) = 8.8 * 0.075 * dust(ji,jj) * rfact2 / &172 & (fse3t(ji,jj,1) * 28.1 * rmtss)173 END DO174 END DO181 !! DO jj = 1, jpj 182 !! DO ji = 1, jpi 183 !! zirondep(ji,jj,1) = (dustsolub * dust(ji,jj) / (55.85 * rmtss) + 3.e-10 / ryyss) & 184 !! & * rfact2 / fse3t(ji,jj,1) 185 !! zsidep (ji,jj) = 8.8 * 0.075 * dust(ji,jj) * rfact2 / & 186 !! & (fse3t(ji,jj,1) * 28.1 * rmtss) 187 !! END DO 188 !! END DO 175 189 176 190 ! sedimentation of detrital nitrogen : upstream scheme
Note: See TracChangeset
for help on using the changeset viewer.