Changeset 600 for trunk/LMDZ.GENERIC/libf/phystd/vdifc.F

Timestamp:

Mar 27, 2012, 5:02:10 PM (13 years ago)

Author:

jleconte

Message:

• Corrects the computation of planck function at the surface in sfluxi

so that its integral is equal to sigma Tsurf^4.

• This ensure that no flux is lost due to:

-truncation of the planck function at high/low wavenumber
-numerical error during first spectral computation of the planck function
-discrepancy between Tsurf and NTS/NTfac in sfluxi

• OLR now equal to LW net heating/cooling at equilibrium!
• As much as possible, only the value of the stephan boltzmann constant defined in racommon_h (and the

corresponding variable, sigma) should be used. Now done in physics, vdifc and turbdiff.

File:

• trunk/LMDZ.GENERIC/libf/phystd/vdifc.F (modified) (2 diffs)

trunk/LMDZ.GENERIC/libf/phystd/vdifc.F

@@ -8 +8 @@
 
       use watercommon_h, only : RLVTT, To, RCPD, mx_eau_sol
+      use radcommon_h, only : sigma
 
       implicit none
@@ -134 +135 @@
       save ivap, iice
 
-      real, parameter :: sigma=5.67e-8
       real, parameter :: karman=0.4
       real cd0, roughratio