- Timestamp:
- 2018-01-15T17:49:43+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90
r7522 r9238 118 118 CALL fld_read( kt, 1, sf_dust ) 119 119 IF( nn_ice_tr == -1 .AND. .NOT. ln_ironice ) THEN 120 dust(:,:) = sf_dust(1)%fnow(:,:,1)120 dust(:,:) = MAX( rtrn, sf_dust(1)%fnow(:,:,1) ) 121 121 ELSE 122 dust(:,:) = sf_dust(1)%fnow(:,:,1) * ( 1.0 - fr_i(:,:) )122 dust(:,:) = MAX( rtrn, sf_dust(1)%fnow(:,:,1) ) * ( 1.0 - fr_i(:,:) ) 123 123 ENDIF 124 124 ENDIF … … 158 158 ! Compute N deposition at nit000 or only if there is more than 1 time record in N deposition file 159 159 IF( ln_ndepo ) THEN 160 ! from kg m-2 s-1 to molC l-1 s-1 160 161 IF( kt == nit000 .OR. ( kt /= nit000 .AND. ntimes_ndep > 1 ) ) THEN 161 zcoef = rno3 * 14E6 * ryyss162 zcoef = 14. * rno3 162 163 CALL fld_read( kt, 1, sf_ndepo ) 163 nitdep(:,:) = sf_ndepo(1)%fnow(:,:,1) / zcoef / fse3t(:,:,1)164 nitdep(:,:) = MAX( rtrn, sf_ndepo(1)%fnow(:,:,1) ) / zcoef / fse3t(:,:,1) 164 165 ENDIF 165 166 IF( lk_vvl ) THEN 166 zcoef = rno3 * 14E6 * ryyss167 nitdep(:,:) = sf_ndepo(1)%fnow(:,:,1) / zcoef / fse3t(:,:,1)167 zcoef = 14. * rno3 168 nitdep(:,:) = MAX ( rtrn, sf_ndepo(1)%fnow(:,:,1) ) / zcoef / fse3t(:,:,1) 168 169 ENDIF 169 170 ENDIF … … 229 230 WRITE(numout,*) ' Specify your sea ice iron concentration in nampisice instead ' 230 231 WRITE(numout,*) ' ln_ironice is forced to .FALSE. ' 231 ln_ironice = .FALSE.232 ENDIF232 ENDIF 233 ln_ironice = .FALSE. 233 234 ENDIF 234 235 … … 283 284 ! Get total input dust ; need to compute total atmospheric supply of Si in a year 284 285 CALL iom_open ( TRIM( sn_dust%clname ) , numdust ) 285 CALL iom_gettime( numdust, zsteps, kntime=ntimes_dust) ! get number of record in file 286 ALLOCATE( zdust(jpi,jpj,ntimes_dust) ) 287 DO jm = 1, ntimes_dust 288 CALL iom_get( numdust, jpdom_data, TRIM( sn_dust%clvar ), zdust(:,:,jm), jm ) 289 END DO 290 CALL iom_close( numdust ) 291 ztimes_dust = 1._wp / FLOAT( ntimes_dust ) 292 sumdepsi = 0.e0 293 DO jm = 1, ntimes_dust 294 sumdepsi = sumdepsi + glob_sum( zdust(:,:,jm) * e1e2t(:,:) * tmask(:,:,1) * ztimes_dust ) 295 ENDDO 296 sumdepsi = sumdepsi / ( nyear_len(1) * rday ) * 12. * 8.8 * 0.075 * mfrac / 28.1 297 DEALLOCATE( zdust) 298 ENDIF 299 ELSE 300 sumdepsi = 0._wp 286 CALL iom_gettime( numdust, zsteps, kntime=ntimes_dust) ! get number of record in file` 287 ENDIF 301 288 END IF 302 289 … … 395 382 CALL iom_open ( TRIM( sn_ndepo%clname ), numdepo ) 396 383 CALL iom_gettime( numdepo, zsteps, kntime=ntimes_ndep) 397 ALLOCATE( zndepo(jpi,jpj,ntimes_ndep) ) 398 DO jm = 1, ntimes_ndep 399 CALL iom_get( numdepo, jpdom_data, TRIM( sn_ndepo%clvar ), zndepo(:,:,jm), jm ) 400 END DO 401 CALL iom_close( numdepo ) 402 ztimes_ndep = 1._wp / FLOAT( ntimes_ndep ) 403 nitdepinput = 0._wp 404 DO jm = 1, ntimes_ndep 405 nitdepinput = nitdepinput + glob_sum( zndepo(:,:,jm) * e1e2t(:,:) * tmask(:,:,1) * ztimes_ndep ) 406 ENDDO 407 nitdepinput = nitdepinput / rno3 / 14E6 408 DEALLOCATE( zndepo) 409 ENDIF 410 ELSE 411 nitdepinput = 0._wp 384 ENDIF 412 385 ENDIF 413 386 … … 496 469 WRITE(numout,*) ' DIC Supply : ', rivdicinput*1E3*12./1E12,'TgC/yr' 497 470 WRITE(numout,*) 498 WRITE(numout,*) ' Total input of elements from atmospheric supply'499 WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~'500 WRITE(numout,*) ' N Supply : ', nitdepinput*rno3*1E3/1E12*14.,' TgN/yr'501 WRITE(numout,*)502 471 ENDIF 503 472 !
Note: See TracChangeset
for help on using the changeset viewer.