163 |
#endif |
#endif |
164 |
#ifdef DYNAMIC_CHL |
#ifdef DYNAMIC_CHL |
165 |
c check Chl fields are reasonable |
c check Chl fields are reasonable |
166 |
|
#ifndef WAVEBANDS |
167 |
|
do np = 1,npmax |
168 |
|
c C:CHl minimum: chosen to be Chl:C at high light (2000uEin/m2/s) and |
169 |
|
c no temp/nutrient limitation |
170 |
|
chl2cmin(np)=chl2cmax(np)/ |
171 |
|
& (1+(chl2cmax(np)*alphachl(np)*2000. _d 0)/ |
172 |
|
& (2*pcmax(np))) |
173 |
|
chl2cmin(np)=0. _d 0 |
174 |
|
enddo |
175 |
|
#else |
176 |
do np=1,npmax |
do np=1,npmax |
177 |
chl2cmin(np)=chl2cmax(np)/ |
chl2cmin(np)=chl2cmax(np)/ |
178 |
& (1+(chl2cmax(np)* alpha_mean(np) *2000. _d 0)/ |
& (1+(chl2cmax(np)* alpha_mean(np) *2000. _d 0)/ |
179 |
& (2*pcmax(np))) |
& (2*pcmax(np))) |
180 |
enddo |
enddo |
181 |
|
#endif |
182 |
call MONOD_CHECK_CHL(myThid) |
call MONOD_CHECK_CHL(myThid) |
183 |
#endif |
#endif |
184 |
#endif |
#endif |