- Timestamp:
- 2015-04-03T10:48:51+02:00 (9 years ago)
- Location:
- branches/2015/dev_r5171_CNRS_LIM3_seaicebgc/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5171_CNRS_LIM3_seaicebgc/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90
r4996 r5193 216 216 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampissbc in configuration namelist', lwp ) 217 217 IF(lwm) WRITE ( numonp, nampissbc ) 218 219 IF ( ( nn_ice_tr >= 0 ) .AND. ln_ironice ) THEN 220 WRITE(numout,*) ' ln_ironice incompatible with nn_ice_tr = ', nn_ice_tr 221 WRITE(numout,*) ' Specify your sea ice iron concentration in nampisice instead ' 222 WRITE(numout,*) ' ln_ironice is forced to .FALSE. ' 223 ln_ironice = .FALSE. 224 ENDIF 218 225 219 226 IF(lwp) THEN -
branches/2015/dev_r5171_CNRS_LIM3_seaicebgc/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsed.F90
r5184 r5193 73 73 CHARACTER (len=25) :: charout 74 74 REAL(wp), POINTER, DIMENSION(:,: ) :: zpdep, zsidep, zwork1, zwork2, zwork3, zwork4 75 REAL(wp), POINTER, DIMENSION(:,: ) :: zdenit2d, z bureff75 REAL(wp), POINTER, DIMENSION(:,: ) :: zdenit2d, zironice, zbureff 76 76 REAL(wp), POINTER, DIMENSION(:,: ) :: zwsbio3, zwsbio4, zwscal 77 77 REAL(wp), POINTER, DIMENSION(:,:,:) :: znitrpot, zirondep, zsoufer … … 97 97 zwork3 (:,:) = 0.e0 98 98 zwork4 (:,:) = 0.e0 99 100 ! Iron input/uptake due to sea ice : Crude parameterization based on Lancelot et al. 101 ! ---------------------------------------------------- 102 IF( ln_ironice ) THEN 103 ! 104 CALL wrk_alloc( jpi, jpj, zironice ) 105 ! 106 DO jj = 1, jpj 107 DO ji = 1, jpi 108 zdep = rfact2 / fse3t(ji,jj,1) 109 zwflux = fmmflx(ji,jj) / 1000._wp 110 zfminus = MIN( 0._wp, -zwflux ) * trn(ji,jj,1,jpfer) * zdep 111 zfplus = MAX( 0._wp, -zwflux ) * icefeinput * zdep 112 zironice(ji,jj) = zfplus + zfminus 113 END DO 114 END DO 115 ! 116 trn(:,:,1,jpfer) = trn(:,:,1,jpfer) + zironice(:,:) 117 ! 118 IF( lk_iomput .AND. jnt == nrdttrc .AND. iom_use( "Ironice" ) ) & 119 & CALL iom_put( "Ironice", zironice(:,:) * 1.e+3 * rfact2r * fse3t(:,:,1) * tmask(:,:,1) ) ! iron flux from ice 120 ! 121 CALL wrk_dealloc( jpi, jpj, zironice ) 122 ! 123 ENDIF 124 99 125 100 126 ! Add the external input of nutrients from dust deposition
Note: See TracChangeset
for help on using the changeset viewer.