Changeset 5193 for branches/2015/dev_r5171_CNRS_LIM3_seaicebgc/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsed.F90
- Timestamp:
- 2015-04-03T10:48:51+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
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.