Changeset 11822 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zpoc.F90
- Timestamp:
- 2019-10-29T11:41:36+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zpoc.F90
r10975 r11822 103 103 ! particle scheme. Does not work with Kriest parameterization. 104 104 ! ----------------------------------------------------------------------- 105 ztremint(:,:,:) = zremigoc(:,:,:) 105 106 DO jk = 2, jpkm1 106 107 DO jj = 1, jpj … … 167 168 & + zsizek1 ) ) * zpoc + ( prodgoc(ji,jj,jk-1) / tgfunc(ji,jj,jk-1) * ( 1. & 168 169 & - exp( -reminp(jn) * zsizek1 ) ) * exp( -reminp(jn) * zsizek ) + prodgoc(ji,jj,jk) & 169 & / tgfunc(ji,jj,jk) * ( 1. - exp( -reminp(jn) * zsizek ) ) ) * rday / rfact2 / reminp(jn) 170 & / tgfunc(ji,jj,jk) * ( 1. - exp( -reminp(jn) * zsizek ) ) ) * rday / rfact2 / reminp(jn) * alphan(jn) 170 171 alphat = alphat + alphag(ji,jj,jk,jn) 171 172 remint = remint + alphag(ji,jj,jk,jn) * reminp(jn) … … 286 287 ! layer, this spectrum is supposed to be uniform. 287 288 ! --------------------------------------------------------------------- 289 ztremint(:,:,:) = zremipoc(:,:,:) 288 290 DO jk = 1, jpkm1 289 291 DO jj = 1, jpj … … 498 500 REWIND( numnatp_ref ) ! Namelist nampisrem in reference namelist : Pisces remineralization 499 501 READ ( numnatp_ref, nampispoc, IOSTAT = ios, ERR = 901) 500 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampispoc in reference namelist' , lwp)502 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampispoc in reference namelist' ) 501 503 REWIND( numnatp_cfg ) ! Namelist nampisrem in configuration namelist : Pisces remineralization 502 504 READ ( numnatp_cfg, nampispoc, IOSTAT = ios, ERR = 902 ) 503 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'nampispoc in configuration namelist' , lwp)505 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'nampispoc in configuration namelist' ) 504 506 IF(lwm) WRITE( numonp, nampispoc ) 505 507
Note: See TracChangeset
for help on using the changeset viewer.