Changeset 2773 for trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zrem.F90
- Timestamp:
- 2011-06-01T12:49:52+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zrem.F90
r2715 r2773 61 61 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 62 62 USE wrk_nemo, ONLY: ztempbac => wrk_2d_1 63 USE wrk_nemo, ONLY: zdepbac => wrk_3d_2 , z fesatur => wrk_3d_2 , zolimi => wrk_3d_463 USE wrk_nemo, ONLY: zdepbac => wrk_3d_2 , zolimi => wrk_3d_3 64 64 ! 65 65 INTEGER, INTENT(in) :: kt ! ocean time step … … 67 67 INTEGER :: ji, jj, jk 68 68 REAL(wp) :: zremip, zremik , zlam1b 69 REAL(wp) :: zkeq , zfeequi, zsiremin 69 REAL(wp) :: zkeq , zfeequi, zsiremin, zfesatur 70 70 REAL(wp) :: zsatur, zsatur2, znusil 71 71 REAL(wp) :: zbactfer, zorem, zorem2, zofer … … 78 78 !!--------------------------------------------------------------------- 79 79 80 IF( wrk_in_use(2, 1) .OR. wrk_in_use(3, 2,3 ,4) ) THEN80 IF( wrk_in_use(2, 1) .OR. wrk_in_use(3, 2,3) ) THEN 81 81 CALL ctl_stop('p4z_rem: requested workspace arrays unavailable') ; RETURN 82 82 ENDIF … … 84 84 ! Initialisation of temprary arrys 85 85 zdepbac (:,:,:) = 0._wp 86 zfesatur(:,:,:) = 0._wp87 86 zolimi (:,:,:) = 0._wp 88 87 ztempbac(:,:) = 0._wp … … 276 275 zsatur = ( sio3eq(ji,jj,jk) - trn(ji,jj,jk,jpsil) ) / ( sio3eq(ji,jj,jk) + rtrn ) 277 276 zsatur = MAX( rtrn, zsatur ) 278 zsatur2 = zsatur * ( 1. + t n(ji,jj,jk) / 400.)**4279 znusil = 0.225 * ( 1. + t n(ji,jj,jk) / 15.) * zsatur + 0.775 * zsatur2**9277 zsatur2 = zsatur * ( 1. + tsn(ji,jj,jk,jp_tem) / 400.)**4 278 znusil = 0.225 * ( 1. + tsn(ji,jj,jk,jp_tem) / 15.) * zsatur + 0.775 * zsatur2**9 280 279 zsiremin = xsirem * zstep * znusil 281 280 zosil = zsiremin * trn(ji,jj,jk,jpdsi) … … 294 293 ENDIF 295 294 296 zfesatur (:,:,:)= 0.6e-9295 zfesatur = 0.6e-9 297 296 !CDIR NOVERRCHK 298 297 DO jk = 1, jpkm1 … … 323 322 ! ------------------------------------------------------------ 324 323 zkeq = fekeq(ji,jj,jk) 325 zfeequi = ( -( 1. + zfesatur (ji,jj,jk)* zkeq - zkeq * trn(ji,jj,jk,jpfer) ) &326 & + SQRT( ( 1. + zfesatur (ji,jj,jk)* zkeq - zkeq * trn(ji,jj,jk,jpfer) )**2 &324 zfeequi = ( -( 1. + zfesatur * zkeq - zkeq * trn(ji,jj,jk,jpfer) ) & 325 & + SQRT( ( 1. + zfesatur * zkeq - zkeq * trn(ji,jj,jk,jpfer) )**2 & 327 326 & + 4. * trn(ji,jj,jk,jpfer) * zkeq) ) / ( 2. * zkeq ) 328 327 … … 395 394 ! 396 395 IF( wrk_not_released(2, 1) .OR. & 397 wrk_not_released(3, 2,3 ,4) ) CALL ctl_stop('p4z_rem: failed to release workspace arrays')396 wrk_not_released(3, 2,3) ) CALL ctl_stop('p4z_rem: failed to release workspace arrays') 398 397 ! 399 398 END SUBROUTINE p4z_rem
Note: See TracChangeset
for help on using the changeset viewer.