- Timestamp:
- 2016-11-28T17:04:10+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90
r5836 r7351 25 25 PUBLIC p4z_sbc_init 26 26 27 !! * Shared module variables28 27 LOGICAL , PUBLIC :: ln_dust !: boolean for dust input from the atmosphere 29 28 LOGICAL , PUBLIC :: ln_solub !: boolean for variable solubility of atmospheric iron … … 45 44 LOGICAL , PUBLIC :: ll_sbc 46 45 47 !! * Module variables48 46 LOGICAL :: ll_solub 49 47 … … 80 78 REAL(wp), PUBLIC :: rivdininput, rivdipinput, rivdsiinput 81 79 82 83 80 !! * Substitutions 84 # include "domzgr_substitute.h90"85 81 # include "vectopt_loop_substitute.h90" 86 82 !!---------------------------------------------------------------------- … … 89 85 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 90 86 !!---------------------------------------------------------------------- 91 92 87 CONTAINS 93 88 … … 160 155 IF( ln_ndepo ) THEN 161 156 IF( kt == nit000 .OR. ( kt /= nit000 .AND. ntimes_ndep > 1 ) ) THEN 162 CALL fld_read( kt, 1, sf_ndepo ) 163 DO jj = 1, jpj 164 DO ji = 1, jpi 165 nitdep(ji,jj) = sf_ndepo(1)%fnow(ji,jj,1) / rno3 / ( 14E6 * ryyss * fse3t(ji,jj,1) + rtrn ) 166 END DO 167 END DO 157 zcoef = rno3 * 14E6 * ryyss 158 CALL fld_read( kt, 1, sf_ndepo ) 159 nitdep(:,:) = sf_ndepo(1)%fnow(:,:,1) / zcoef / e3t_n(:,:,1) 160 ENDIF 161 IF( .NOT.ln_linssh ) THEN 162 zcoef = rno3 * 14E6 * ryyss 163 nitdep(:,:) = sf_ndepo(1)%fnow(:,:,1) / zcoef / e3t_n(:,:,1) 168 164 ENDIF 169 165 ENDIF … … 267 263 IF( lk_offline ) THEN 268 264 nk_rnf(:,:) = 1 269 h_rnf (:,:) = fsdept(:,:,1)265 h_rnf (:,:) = gdept_n(:,:,1) 270 266 ENDIF 271 267 … … 456 452 DO jj = 1, jpj 457 453 DO ji = 1, jpi 458 zexpide = MIN( 8.,( fsdept(ji,jj,jk) / 500. )**(-1.5) )454 zexpide = MIN( 8.,( gdept_n(ji,jj,jk) / 500. )**(-1.5) ) 459 455 zdenitide = -0.9543 + 0.7662 * LOG( zexpide ) - 0.235 * LOG( zexpide )**2 460 456 zcmask(ji,jj,jk) = zcmask(ji,jj,jk) * MIN( 1., EXP( zdenitide ) / 0.5 ) … … 466 462 ironsed(:,:,jpk) = 0._wp 467 463 DO jk = 1, jpkm1 468 ironsed(:,:,jk) = sedfeinput * zcmask(:,:,jk) / ( fse3t(:,:,jk) * rday )464 ironsed(:,:,jk) = sedfeinput * zcmask(:,:,jk) / ( e3t_0(:,:,jk) * rday ) 469 465 END DO 470 466 DEALLOCATE( zcmask) … … 484 480 CALL iom_close( numhydro ) 485 481 ! 486 hydrofe(:,:,:) = ( hydrofe(:,:,:) * hratio ) / ( cvol(:,:,:) * ryyss + rtrn ) / 1000._wp 482 DO jk = 1, jpk 483 hydrofe(:,:,jk) = ( hydrofe(:,:,jk) * hratio ) / ( e1e2t(:,:) * e3t_0(:,:,jk) * ryyss + rtrn ) / 1000._wp 484 ENDDO 487 485 ! 488 486 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.