340 |
|
|
341 |
#else |
#else |
342 |
|
|
343 |
z_dvm_regr = 398. _d 0 - 0.56 _d 0*min(300. _d 0, |
z_dvm_regr = 398. _d 0 |
344 |
& max(-10. _d 0,(o2_upper - o2_lower))) - |
& - 0.56 _d 0*min(300. _d 0,max(-10. _d 0,(o2_upper - o2_lower))) |
345 |
& 115. _d 0*min(0.85 _d 0,max(-1.80 _d 0,log10(chl(i,j,1,bi,bj)))) |
& - 115. _d 0*min(0.85 _d 0,max(-1.80 _d 0, |
346 |
& + 0.36 _d 0*min(500. _d 0,max(epsln,mld(i,j))) - |
& log10(max(chl(i,j,1,bi,bj),chl_min)))) |
347 |
& 2.4 _d 0*min(20. _d 0,max(-3. _d 0, (temp_upper - temp_lower))) |
& + 0.36 _d 0*min(500. _d 0,max(epsln,mld(i,j))) |
348 |
|
& - 2.40 _d 0*min(20. _d 0,max(-3. _d 0,(temp_upper-temp_lower))) |
349 |
|
|
350 |
! Limit the depth of migration in polar winter. |
c ! Limit the depth of migration in polar winter. |
351 |
! Use irr_mem since this is averaged over multiple days, dampening the diurnal cycle. |
c ! Use irr_mem since this is averaged over multiple days, dampening the diurnal cycle. |
352 |
! Tapers Z_DVM to the minimum when surface irradince is below a given threshold (here 10 W/m2). |
c ! Tapers Z_DVM to the minimum when surface irradince is below a given threshold (here 10 W/m2). |
353 |
|
|
354 |
if ( irr_mem(i,j,1,bi,bj) .lt. 10. ) then |
if ( irr_mem(i,j,1,bi,bj) .lt. 10. ) then |
355 |
z_dvm_regr = 150. _d 0 + (z_dvm_regr - 150. _d 0) * |
z_dvm_regr = 150. _d 0 + (z_dvm_regr - 150. _d 0) * |