| 1 | c This subroutine estimateis Sea Salt emission fluxes over |
|---|
| 2 | c Oceanic surfaces. |
|---|
| 3 | c |
|---|
| 4 | SUBROUTINE seasalt(v_10m, u_10m, pct_ocean, lmt_sea_salt) |
|---|
| 5 | |
|---|
| 6 | USE dimphy |
|---|
| 7 | IMPLICIT NONE |
|---|
| 8 | c |
|---|
| 9 | INCLUDE "dimensions.h" |
|---|
| 10 | INCLUDE "chem.h" |
|---|
| 11 | INCLUDE "chem_spla.h" |
|---|
| 12 | INCLUDE "YOMCST.h" |
|---|
| 13 | INCLUDE "YOECUMF.h" |
|---|
| 14 | c |
|---|
| 15 | INTEGER i, bin !local variables |
|---|
| 16 | REAL pct_ocean(klon) !hfraction of Ocean in each grid |
|---|
| 17 | REAL v_10m(klon), u_10m(klon) !V&H components of wind @10 m |
|---|
| 18 | REAL w_speed_10m(klon) !wind speed at 10m from surface |
|---|
| 19 | REAL lmt_sea_salt(klon,ss_bins)!sea salt emission flux - mg/m2/s |
|---|
| 20 | REAL sea_salt_flux(ss_bins) !sea salt emission flux per unit wind speed |
|---|
| 21 | |
|---|
| 22 | REAL wind, ocean |
|---|
| 23 | c |
|---|
| 24 | c------Sea salt emission fluxes for each size bin calculated |
|---|
| 25 | c------based on on parameterisation of Gong et al. (1997). |
|---|
| 26 | c------Fluxes of sea salt for each size bin are given in mg/m^2/sec |
|---|
| 27 | c------at wind speed of 1 m/s at 10m height (at 80% RH). |
|---|
| 28 | c------Fluxes at various wind speeds (@10 m from sea |
|---|
| 29 | c------surfaces are estimated using relationship: F=flux*U_10^3.14 |
|---|
| 30 | c |
|---|
| 31 | cnhl for size bin of 0.03-0.5 and 0.5-20 |
|---|
| 32 | DATA sea_salt_flux/4.5E-09,8.7E-7/ |
|---|
| 33 | |
|---|
| 34 | DO i=1, klon |
|---|
| 35 | w_speed_10m(i)= (v_10m(i)**2.0+u_10m(i)**2.0)**0.5 |
|---|
| 36 | ENDDO |
|---|
| 37 | c |
|---|
| 38 | DO bin=1,ss_bins |
|---|
| 39 | wind=0.0 |
|---|
| 40 | ocean=0.0 |
|---|
| 41 | DO i=1, klon |
|---|
| 42 | lmt_sea_salt(i,bin)=sea_salt_flux(bin)*(w_speed_10m(i)**3.41) |
|---|
| 43 | . *pct_ocean(i)*1.e-4*1.e-3 !g/cm2/s |
|---|
| 44 | wind=wind+w_speed_10m(i) |
|---|
| 45 | ocean=ocean+pct_ocean(i) |
|---|
| 46 | ENDDO |
|---|
| 47 | ! print *,'Sea Salt flux = ',sea_salt_flux(bin) |
|---|
| 48 | ENDDO |
|---|
| 49 | ! print *,'SUM OF WIND = ',wind |
|---|
| 50 | ! print *,'SUM OF OCEAN SURFACE = ',ocean |
|---|
| 51 | RETURN |
|---|
| 52 | END |
|---|
