| 1 | SUBROUTINE nightingale(u, v, u_10m, v_10m, paprs, pplay, & |
|---|
| 2 | cdragh, cdragm, t, q, ftsol, tsol, & |
|---|
| 3 | pctsrf, lmt_dmsconc, lmt_dms) |
|---|
| 4 | ! |
|---|
| 5 | USE dimphy |
|---|
| 6 | USE indice_sol_mod |
|---|
| 7 | USE yomcst_mod_h |
|---|
| 8 | IMPLICIT NONE |
|---|
| 9 | ! |
|---|
| 10 | |
|---|
| 11 | |
|---|
| 12 | ! |
|---|
| 13 | REAL :: u(klon,klev), v(klon,klev) |
|---|
| 14 | REAL :: u_10m(klon), v_10m(klon) |
|---|
| 15 | REAL :: ftsol(klon,nbsrf) |
|---|
| 16 | REAL :: tsol(klon) |
|---|
| 17 | REAL :: paprs(klon,klev+1), pplay(klon,klev) |
|---|
| 18 | REAL :: t(klon,klev) |
|---|
| 19 | REAL :: q(klon,klev) |
|---|
| 20 | REAL :: cdragh(klon), cdragm(klon) |
|---|
| 21 | REAL :: pctsrf(klon,nbsrf) |
|---|
| 22 | REAL :: lmt_dmsconc(klon) ! concentration oceanique DMS |
|---|
| 23 | REAL :: lmt_dms(klon) ! flux de DMS |
|---|
| 24 | ! |
|---|
| 25 | REAL :: ustar(klon), obklen(klon) |
|---|
| 26 | REAL :: u10(klon), u10n(klon) |
|---|
| 27 | REAL :: tvelocity, schmidt_corr |
|---|
| 28 | REAL :: t1, t2, t3, t4, viscosity_kin, diffusivity, schmidt |
|---|
| 29 | INTEGER :: i |
|---|
| 30 | ! |
|---|
| 31 | CALL bl_for_dms(u, v, paprs, pplay, cdragh, cdragm, & |
|---|
| 32 | t, q, tsol, ustar, obklen) |
|---|
| 33 | ! |
|---|
| 34 | DO i=1,klon |
|---|
| 35 | u10(i)=SQRT(u_10m(i)**2+v_10m(i)**2) |
|---|
| 36 | ENDDO |
|---|
| 37 | ! |
|---|
| 38 | CALL neutral(u10, ustar, obklen, u10n) |
|---|
| 39 | ! |
|---|
| 40 | DO i=1,klon |
|---|
| 41 | ! |
|---|
| 42 | ! tvelocity - transfer velocity, also known as kw (cm/s) |
|---|
| 43 | ! schmidt_corr - Schmidt number correction factor (dimensionless) |
|---|
| 44 | ! Reference: Nightingale, P.D., G. Malin, C. S. Law, J. J. Watson, P.S. Liss |
|---|
| 45 | ! M. I. Liddicoat, J. Boutin, R.C. Upstill-Goddard. 'In situ evaluation |
|---|
| 46 | ! of air-sea gas exchange parameterizations using conservative and |
|---|
| 47 | ! volatile tracers.' Glob. Biogeochem. Cycles, 14:373-387, 2000. |
|---|
| 48 | ! compute transfer velocity using u10neutral |
|---|
| 49 | ! |
|---|
| 50 | tvelocity = 0.222*u10n(i)*u10n(i) + 0.333*u10n(i) |
|---|
| 51 | ! |
|---|
| 52 | ! above expression gives tvelocity in cm/hr. convert to cm/s. 1hr =3600 sec |
|---|
| 53 | |
|---|
| 54 | tvelocity = tvelocity / 3600. |
|---|
| 55 | |
|---|
| 56 | ! compute the correction factor, which for Nightingale parameterization is |
|---|
| 57 | ! based on how different the schmidt number is from 600. |
|---|
| 58 | ! correction factor based on temperature in Kelvin. good |
|---|
| 59 | ! only for t<=30 deg C. for temperatures above that, set correction factor |
|---|
| 60 | ! equal to value at 30 deg C. |
|---|
| 61 | |
|---|
| 62 | IF (ftsol(i,is_oce) .LE. 303.15) THEN |
|---|
| 63 | t1 = ftsol(i,is_oce) |
|---|
| 64 | ELSE |
|---|
| 65 | t1 = 303.15 |
|---|
| 66 | ENDIF |
|---|
| 67 | |
|---|
| 68 | t2 = t1 * t1 |
|---|
| 69 | t3 = t2 * t1 |
|---|
| 70 | t4 = t3 * t1 |
|---|
| 71 | viscosity_kin = 3.0363e-9*t4 - 3.655198e-6*t3 + 1.65333e-3*t2 & |
|---|
| 72 | - 3.332083e-1*t1 + 25.26819 |
|---|
| 73 | diffusivity = 0.01922 * exp(-2177.1/t1) |
|---|
| 74 | schmidt = viscosity_kin / diffusivity |
|---|
| 75 | schmidt_corr = (schmidt/600.)**(-.5) |
|---|
| 76 | ! |
|---|
| 77 | lmt_dms(i) = tvelocity * pctsrf(i,is_oce) & |
|---|
| 78 | * lmt_dmsconc(i)/1.0e12 * schmidt_corr * RNAVO |
|---|
| 79 | ! |
|---|
| 80 | IF (lmt_dmsconc(i).LE.1.e-20) lmt_dms(i)=0.0 |
|---|
| 81 | ! |
|---|
| 82 | ENDDO |
|---|
| 83 | ! |
|---|
| 84 | END SUBROUTINE nightingale |
|---|
