[5246] | 1 | ! This subroutine estimateis Sea Salt emission fluxes over |
---|
| 2 | ! Oceanic surfaces. |
---|
| 3 | ! |
---|
| 4 | SUBROUTINE seasalt(v_10m, u_10m, pct_ocean, lmt_sea_salt) |
---|
[2630] | 5 | |
---|
[5246] | 6 | USE dimphy |
---|
[5271] | 7 | USE dimensions_mod, ONLY: iim, jjm, llm, ndm |
---|
[5274] | 8 | USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO & |
---|
| 9 | , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA & |
---|
| 10 | , R_ecc, R_peri, R_incl & |
---|
| 11 | , RA, RG, R1SA & |
---|
| 12 | , RSIGMA & |
---|
| 13 | , R, RMD, RMV, RD, RV, RCPD & |
---|
| 14 | , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12 & |
---|
| 15 | , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w & |
---|
| 16 | , RCW, RCS & |
---|
| 17 | , RLVTT, RLSTT, RLMLT, RTT, RATM & |
---|
| 18 | , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS & |
---|
| 19 | , RALPD, RBETD, RGAMD |
---|
[5271] | 20 | IMPLICIT NONE |
---|
[5246] | 21 | ! |
---|
[5271] | 22 | |
---|
[5246] | 23 | INCLUDE "chem.h" |
---|
| 24 | INCLUDE "chem_spla.h" |
---|
[5274] | 25 | |
---|
[5246] | 26 | INCLUDE "YOECUMF.h" |
---|
| 27 | ! |
---|
| 28 | INTEGER :: i, bin !local variables |
---|
| 29 | REAL :: pct_ocean(klon) !hfraction of Ocean in each grid |
---|
| 30 | REAL :: v_10m(klon), u_10m(klon) !V&H components of wind @10 m |
---|
| 31 | REAL :: w_speed_10m(klon) !wind speed at 10m from surface |
---|
| 32 | REAL :: lmt_sea_salt(klon,ss_bins)!sea salt emission flux - mg/m2/s |
---|
| 33 | REAL :: sea_salt_flux(ss_bins) !sea salt emission flux per unit wind speed |
---|
[2630] | 34 | |
---|
[5246] | 35 | REAL :: wind, ocean |
---|
| 36 | ! |
---|
| 37 | !------Sea salt emission fluxes for each size bin calculated |
---|
| 38 | !------based on on parameterisation of Gong et al. (1997). |
---|
| 39 | !------Fluxes of sea salt for each size bin are given in mg/m^2/sec |
---|
| 40 | !------at wind speed of 1 m/s at 10m height (at 80% RH). |
---|
| 41 | !------Fluxes at various wind speeds (@10 m from sea |
---|
| 42 | !------surfaces are estimated using relationship: F=flux*U_10^3.14 |
---|
| 43 | ! |
---|
| 44 | !nhl for size bin of 0.03-0.5 and 0.5-20 |
---|
| 45 | DATA sea_salt_flux/4.5E-09,8.7E-7/ |
---|
| 46 | |
---|
| 47 | DO i=1, klon |
---|
| 48 | w_speed_10m(i)= (v_10m(i)**2.0+u_10m(i)**2.0)**0.5 |
---|
| 49 | ENDDO |
---|
| 50 | ! |
---|
| 51 | DO bin=1,ss_bins |
---|
| 52 | wind=0.0 |
---|
| 53 | ocean=0.0 |
---|
| 54 | DO i=1, klon |
---|
| 55 | lmt_sea_salt(i,bin)=sea_salt_flux(bin)*(w_speed_10m(i)**3.41) & |
---|
| 56 | *pct_ocean(i)*1.e-4*1.e-3 !g/cm2/s |
---|
| 57 | wind=wind+w_speed_10m(i) |
---|
| 58 | ocean=ocean+pct_ocean(i) |
---|
| 59 | ENDDO |
---|
| 60 | ! print *,'Sea Salt flux = ',sea_salt_flux(bin) |
---|
| 61 | ENDDO |
---|
| 62 | ! print *,'SUM OF WIND = ',wind |
---|
| 63 | ! print *,'SUM OF OCEAN SURFACE = ',ocean |
---|
| 64 | RETURN |
---|
| 65 | END SUBROUTINE seasalt |
---|