Changeset 10788 for NEMO/trunk/src/TOP/PISCES
- Timestamp:
- 2019-03-21T12:15:14+01:00 (5 years ago)
- Location:
- NEMO/trunk/src/TOP/PISCES
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/TOP/PISCES/P4Z/p4zsbc.F90
r10522 r10788 27 27 LOGICAL , PUBLIC :: ln_ironsed !: boolean for Fe input from sediments 28 28 LOGICAL , PUBLIC :: ln_hydrofe !: boolean for Fe input from hydrothermal vents 29 LOGICAL , PUBLIC :: ln_ironice !: boolean for Fe input from sea ice30 29 REAL(wp), PUBLIC :: sedfeinput !: Coastal release of Iron 31 30 REAL(wp), PUBLIC :: dustsolub !: Solubility of the dust … … 109 108 IF( kt == nit000 .OR. ( kt /= nit000 .AND. ntimes_dust > 1 ) ) THEN 110 109 CALL fld_read( kt, 1, sf_dust ) 111 IF( nn_ice_tr == -1 .AND. .NOT.ln_ironice ) THEN ; dust(:,:) = MAX( rtrn, sf_dust(1)%fnow(:,:,1) ) 112 ELSE ; dust(:,:) = MAX( rtrn, sf_dust(1)%fnow(:,:,1) * ( 1.-fr_i(:,:) ) ) 113 ENDIF 110 dust(:,:) = MAX( rtrn, sf_dust(1)%fnow(:,:,1) ) * ( 1.0 - fr_i(ji,jj) ) 114 111 ENDIF 115 112 ENDIF … … 265 262 END IF 266 263 267 IF( nn_ice_tr >= 0 .AND. ln_ironice ) THEN268 IF(lwp) THEN269 WRITE(numout,*) ' ==>>> ln_ironice incompatible with nn_ice_tr = ', nn_ice_tr270 WRITE(numout,*) ' Specify your sea ice iron concentration in nampisice instead '271 WRITE(numout,*) ' ln_ironice is forced to .FALSE. '272 ENDIF273 ln_ironice = .FALSE.274 ENDIF275 276 264 IF( ln_dust .OR. ln_river .OR. ln_ndepo ) THEN ; ll_sbc = .TRUE. 277 265 ELSE ; ll_sbc = .FALSE. -
NEMO/trunk/src/TOP/PISCES/P4Z/p4zsed.F90
r10425 r10788 53 53 INTEGER :: ji, jj, jk, ikt 54 54 REAL(wp) :: zrivalk, zrivsil, zrivno3 55 REAL(wp) :: zwflux, zfminus, zfplus 56 REAL(wp) :: zlim, zfact, zfactcal 55 REAL(wp) :: zwflux, zlim, zfact, zfactcal 57 56 REAL(wp) :: zo2, zno3, zflx, zpdenit, z1pdenit, zolimit 58 57 REAL(wp) :: zsiloss, zcaloss, zws3, zws4, zwsc, zdep … … 104 103 zdep = rfact2 / e3t_n(ji,jj,1) 105 104 zwflux = fmmflx(ji,jj) / 1000._wp 106 zfminus = MIN( 0._wp, -zwflux ) * trb(ji,jj,1,jpfer) * zdep 107 zfplus = MAX( 0._wp, -zwflux ) * icefeinput * zdep 108 zironice(ji,jj) = zfplus + zfminus 105 zironice(ji,jj) = MAX( -0.99 * trb(ji,jj,1,jpfer), -zwflux * icefeinput * zdep ) 109 106 END DO 110 107 END DO -
NEMO/trunk/src/TOP/PISCES/sms_pisces.F90
r10425 r10788 67 67 !!* Mass conservation 68 68 LOGICAL :: ln_check_mass !: Flag to check mass conservation 69 LOGICAL , PUBLIC :: ln_ironice !: boolean for Fe input from sea ice 69 70 70 71 !!* Biological fluxes for light : variables shared by pisces & lobster -
NEMO/trunk/src/TOP/PISCES/trcice_pisces.F90
r10213 r10788 266 266 zpisc(jpgon,4) = 1.05e-8_wp 267 267 zpisc(jpgop,4) = 1.05e-8_wp 268 268 ! 269 ! ln_ironice and tracers in seaice are redundant. Thus, if tracers in ice 270 ! is activated, ln_ironice should be set to false 271 ! ------------------------------------------------------------------------ 272 IF( nn_ice_tr >= 0 .AND. ln_ironice ) THEN 273 IF(lwp) THEN 274 WRITE(numout,*) ' ==>>> ln_ironice incompatible with nn_ice_tr = ', nn_ice_tr 275 WRITE(numout,*) ' Specify your sea ice iron concentration in nampisice instead ' 276 WRITE(numout,*) ' ln_ironice is forced to .FALSE. ' 277 ENDIF 278 ln_ironice = .FALSE. 279 ENDIF 280 ! 269 281 DO jn = jp_pcs0, jp_pcs1 270 282 IF( cn_trc_o(jn) == 'GL ' ) trc_o(:,:,jn) = zpisc(jn,1) ! Global case … … 310 322 trc_i(:,:,jn) = trc_ice_prescr(jn) 311 323 ENDIF 312 313 324 !-- Baltic 314 325 IF( cn_cfg == "orca" .OR. cn_cfg == "ORCA" ) THEN … … 318 329 trc_i(:,:,jn) = zratio(jn,2) * trc_o(:,:,jn) 319 330 END WHERE 320 ELSE ! prescribed tracer concentration in ice 321 WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND. & 322 54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp ) 323 trc_i(:,:,jn) = trc_ice_prescr(jn) 324 END WHERE 325 ENDIF ! trc_ice_ratio 331 ENDIF 326 332 ENDIF 327 333 !
Note: See TracChangeset
for help on using the changeset viewer.