- Timestamp:
- 2020-06-26T10:26:32+02:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement
- Files:
-
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement
- Property svn:externals
-
old new 8 8 9 9 # SETTE 10 ^/utils/CI/sette@ HEADsette10 ^/utils/CI/sette@12931 sette
-
- Property svn:externals
-
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/src/TOP/PISCES/P4Z/p4zmeso.F90
r12377 r13159 69 69 REAL(wp) :: zfact , zfood, zfoodlim, zproport, zbeta 70 70 REAL(wp) :: zmortzgoc, zfrac, zfracfe, zratio, zratio2, zfracal, zgrazcal 71 REAL(wp) :: zepsherf, zepshert, zepsherv, zgrarsig, zgraztotc, zgraztotn, zgraztotf 71 REAL(wp) :: zepsherf, zepshert, zepsherv, zepsherq 72 REAL(wp) :: zgrarsig, zgraztotc, zgraztotn, zgraztotf 72 73 REAL(wp) :: zgrarem2, zgrafer2, zgrapoc2, zprcaca, zmortz, zgrasrat, zgrasratn 73 74 REAL(wp) :: zrespz, ztortz, zgrazd, zgrazz, zgrazpof … … 156 157 zgrazing2(ji,jj,jk) = zgraztotc 157 158 158 ! Mesozooplankton efficiency 159 ! -------------------------- 159 ! Mesozooplankton efficiency. 160 ! We adopt a formulation proposed by Mitra et al. (2007) 161 ! The gross growth efficiency is controled by the most limiting nutrient. 162 ! Growth is also further decreased when the food quality is poor. This is currently 163 ! hard coded : it can be decreased by up to 50% (zepsherq) 164 ! GGE can also be decreased when food quantity is high, zepsherf (Montagnes and 165 ! Fulton, 2012) 166 ! ----------------------------------------------------------------------------------- 160 167 zgrasrat = ( zgraztotf + rtrn )/ ( zgraztotc + rtrn ) 161 168 zgrasratn = ( zgraztotn + rtrn )/ ( zgraztotc + rtrn ) … … 163 170 zbeta = MAX(0., (epsher2 - epsher2min) ) 164 171 zepsherf = epsher2min + zbeta / ( 1.0 + 0.04E6 * 12. * zfood * zbeta ) 165 zepsherv = zepsherf * zepshert 172 zepsherq = 0.5 + (1.0 - 0.5) * zepshert * ( 1.0 + 1.0 ) / ( zepshert + 1.0 ) 173 zepsherv = zepsherf * zepshert * zepsherq 166 174 167 175 zgrarem2 = zgraztotc * ( 1. - zepsherv - unass2 ) & … … 170 178 & + ferat3 * ( ( 1. - epsher2 - unass2 ) /( 1. - epsher2 ) * ztortz ) 171 179 zgrapoc2 = zgraztotc * unass2 180 172 181 173 182 ! Update the arrays TRA which contain the biological sources and sinks -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/src/TOP/PISCES/P4Z/p4zmicro.F90
r12377 r13159 67 67 REAL(wp) :: zgraze , zdenom, zdenom2 68 68 REAL(wp) :: zfact , zfood, zfoodlim, zbeta 69 REAL(wp) :: zepsherf, zepshert, zepsherv, zgrarsig, zgraztotc, zgraztotn, zgraztotf 69 REAL(wp) :: zepsherf, zepshert, zepsherv, zepsherq 70 REAL(wp) :: zgrarsig, zgraztotc, zgraztotn, zgraztotf 70 71 REAL(wp) :: zgrarem, zgrafer, zgrapoc, zprcaca, zmortz 71 72 REAL(wp) :: zrespz, ztortz, zgrasrat, zgrasratn … … 119 120 zgrazing(ji,jj,jk) = zgraztotc 120 121 121 ! Various remineralization and excretion terms 122 ! -------------------------------------------- 122 123 ! Microzooplankton efficiency. 124 ! We adopt a formulation proposed by Mitra et al. (2007) 125 ! The gross growth efficiency is controled by the most limiting nutrient. 126 ! Growth is also further decreased when the food quality is poor. This is currently 127 ! hard coded : it can be decreased by up to 50% (zepsherq) 128 ! GGE can also be decreased when food quantity is high, zepsherf (Montagnes and 129 ! Fulton, 2012) 130 ! ----------------------------------------------------------------------------- 123 131 zgrasrat = ( zgraztotf + rtrn ) / ( zgraztotc + rtrn ) 124 132 zgrasratn = ( zgraztotn + rtrn ) / ( zgraztotc + rtrn ) … … 126 134 zbeta = MAX(0., (epsher - epshermin) ) 127 135 zepsherf = epshermin + zbeta / ( 1.0 + 0.04E6 * 12. * zfood * zbeta ) 128 zepsherv = zepsherf * zepshert 136 zepsherq = 0.5 + (1.0 - 0.5) * zepshert * ( 1.0 + 1.0 ) / ( zepshert + 1.0 ) 137 zepsherv = zepsherf * zepshert * zepsherq 129 138 130 139 zgrafer = zgraztotc * MAX( 0. , ( 1. - unass ) * zgrasrat - ferat3 * zepsherv ) -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/src/TOP/PISCES/P4Z/p4zsms.F90
r12489 r13159 206 206 IF( l_trdtrc ) THEN 207 207 DO jn = jp_pcs0, jp_pcs1 208 ztrdt(:,:,:,jn) = ( tr(:,:,:,jn,Kbb) - ztrdt(:,:,:,jn) ) * rfact 2r208 ztrdt(:,:,:,jn) = ( tr(:,:,:,jn,Kbb) - ztrdt(:,:,:,jn) ) * rfactr 209 209 CALL trd_trc( tr(:,:,:,jn,Krhs), jn, jptra_sms, kt, Kmm ) ! save trends 210 210 END DO -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/src/TOP/PISCES/SED/sedchem.F90
r12377 r13159 577 577 saltprac(:) = salt(:) * 35.0 / 35.16504 578 578 ELSE 579 saltprac(:) = temp(:)579 saltprac(:) = salt(:) 580 580 ENDIF 581 581 -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/src/TOP/PISCES/SED/sedinorg.F90
r10225 r13159 89 89 zsolcpcl = zsolcpcl + solcp(ji,jk,jsclay) * dz(jk) 90 90 END DO 91 zsolcpsi = MAX( zsolcpsi, rtrn ) 91 92 zsieq(ji) = sieqs(ji) * MAX(0.25, 1.0 - (0.045 * zsolcpcl / zsolcpsi )**0.58 ) 92 93 zsieq(ji) = MAX( rtrn, sieqs(ji) ) -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/src/TOP/PISCES/SED/sedrst.F90
r12489 r13159 80 80 IF(lwp) WRITE(numsed,*) & 81 81 ' open sed restart.output NetCDF file: ',TRIM(clpath)//clname 82 CALL iom_open( TRIM(clpath)//TRIM(clname), numrsw, ldwrt = .TRUE., kdlev = jpksed )82 CALL iom_open( TRIM(clpath)//TRIM(clname), numrsw, ldwrt = .TRUE., kdlev = jpksed, cdcomp = 'SED' ) 83 83 lrst_sed = .TRUE. 84 84 ENDIF -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/src/TOP/trcbc.F90
r12489 r13159 151 151 IF(trcdta_bdy(jn,ib)%cn_obc == 'frs' .AND. nn_trcdmp_bdy(ib) /= 0 ) & 152 152 & CALL ctl_stop( 'trc_bc_ini: Use FRS OR relaxation' ) 153 IF( .NOT.( 0 < nn_trcdmp_bdy(ib) .AND. nn_trcdmp_bdy(ib) <= 2 ) ) &153 IF( .NOT.( 0 <= nn_trcdmp_bdy(ib) .AND. nn_trcdmp_bdy(ib) <= 2 ) ) & 154 154 & CALL ctl_stop( 'trc_bc_ini: Not a valid option for nn_trcdmp_bdy. Allowed: 0,1,2.' ) 155 155 END DO -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/src/TOP/trcstp.F90
r12489 r13159 142 142 ! 143 143 ! Define logical parameter ton control dirunal cycle in TOP 144 l_trcdm2dc = ln_dm2dc .OR. ( ln_cpl .AND. ncpl_qsr_freq /= 1 ) 145 l_trcdm2dc = l_trcdm2dc .AND. .NOT. l_offline 144 l_trcdm2dc = ln_dm2dc .OR. ( ln_cpl .AND. ncpl_qsr_freq /= 1 .AND. ncpl_qsr_freq /= 0 ) 145 l_trcdm2dc = l_trcdm2dc .AND. .NOT. l_offline 146 ! 146 147 IF( l_trcdm2dc .AND. lwp ) CALL ctl_warn( 'Coupling with passive tracers and used of diurnal cycle.', & 147 148 & 'Computation of a daily mean shortwave for some biogeochemical models ' )
Note: See TracChangeset
for help on using the changeset viewer.