source: LMDZ5/trunk/libf/phylmd/Dust/seasalt.F

Last change on this file was 2630, checked in by fhourdin, 8 years ago

Importation du modèle d'aérosols de Boucher, Escribano et al.

File size: 1.7 KB
Line 
1c       This subroutine estimateis Sea Salt emission fluxes over
2c       Oceanic surfaces.
3c
4      SUBROUTINE seasalt(v_10m, u_10m, pct_ocean, lmt_sea_salt)
5   
6      USE dimphy
7      IMPLICIT NONE
8c
9#include "dimensions.h"
10#include "chem.h"
11#include "chem_spla.h"
12#include "YOMCST.h"
13#include "YOECUMF.h"
14c
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
23c
24c------Sea salt emission fluxes for each size bin calculated
25c------based on on parameterisation of Gong et al. (1997).
26c------Fluxes of sea salt for each size bin are given in mg/m^2/sec
27c------at wind speed of 1 m/s at 10m height (at 80% RH).
28c------Fluxes at various wind speeds (@10 m from sea
29c------surfaces are estimated using relationship: F=flux*U_10^3.14
30c
31cnhl 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
37c
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
Note: See TracBrowser for help on using the repository browser.