source: trunk/src/lwnet_clark.pro

Last change on this file was 203, checked in by pinsard, 10 years ago

fix thanks to coding rules
  • Property svn:keywords set to Id URL
File size: 1.3 KB
Line 
1;+
2;
3; ===============
4; lwnet_clark.pro
5; ===============
6;
7; .. function:: lwnet_clark(sst, qair, cld, tair, ps)
8;
9; DESCRIPTION
10; ===========
11;
12; Net longwave from the Clark et al. (1974) bulk formula
13; Clark, N. E., L. Eber, R. M. Laurs, J. A. Renner, and J. F. T. Saur,
14; 1974: Heat exchange between ocean and atmosphere in the
15; eastern North Pacific for 1961-71. NOAA Tech. Rep. NMRS
16; SSRF-682, National Oceanic and Atmospheric Administration,
17; U.S. Department of Commerce, Washington, DC, 108 pp.
18;
19; See also Gill, p34.
20;
21; INPUTS
22;   sst :    Sea surface temperature (degC)
23;   qair:    specific humidity (g/kg)
24;   cld :    cloud cover
25;   tair:    surface air temperature (degC)
26;   ps  :    surface pressure (mb)
27; OUTPUT
28;   lw :    estimated longwave (W/m2)
29;
30; EVOLUTIONS
31; ==========
32;
33; $Id$
34;
35; $URL$
36;
37; - fplod 20120319
38;
39;   * minimal header
40;
41;-
42function lwnet_clark, sst,qair,cld,tair,ps
43;
44@common
45;
46eps=.97            ; emissivity of sea-water
47sig=5.67e-8        ; Stefan-Boltzman constant (W/m2/K4)
48Ts=(sst+273.15)     ; SST in K
49Ta=(tair+273.15)      ; Tair in K
50q=qair/1000.            ; specific humidity (kg/kg)
51r=q/(1-q)        ; Mixing ratio
52ea=ps*r/(r+.62197)      ; Specific humidity into vapor pressure
53;
54lw=    eps*sig*Ts^4*(.39-.05*ea^.5)*(1.-.51*cld^2) $
55   +4.*eps*sig*Ts^3*(Ts-Ta)
56;
57return, lw
58;
59end
Note: See TracBrowser for help on using the repository browser.