Changeset 526 for trunk/LMDZ.GENERIC/libf/phystd/aeropacity.F90

Timestamp:

Feb 13, 2012, 6:30:26 PM (13 years ago)

Author:

jleconte

Message:

13/02/2012 == JL + AS

• All outputs are now in netCDF format. Even in 1D (No more G1D)
• Clean up of the call to callcorrk when CLFvarying=true
• Corrects a bug in writediagspecIR/VI. Output are now in W/m2/cm-1 as a function of the wavenumber in cm-1
• Enable writediagspecIR/V to work in the CLFvarying=true case (output now done in Physiq after writediagfi)
• Add a simple treatment for the supersaturation of CO2 (see forget et al 2012)
• corrects a small bug when no clouds are present in aeropacity

File:

• trunk/LMDZ.GENERIC/libf/phystd/aeropacity.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/aeropacity.F90

@@ -235 +235 @@
                           (  0.75 * QREFvis3d(ig,l,iaer) /        &
                           ( rho_ice * reffrad(ig,l,iaer) )  ) *   &
-                          ( pq(ig,l,i_h2oice) + 1.E-9 ) *         &
-                          ( pplev(ig,l) - pplev(ig,l+1) ) / g /   &
-                            CLF1
+                            pq(ig,l,i_h2oice) *                   & !JL I dropped the +1e-9 here to have the same
+                          ( pplev(ig,l) - pplev(ig,l+1) ) / g /   & !   opacity in the clearsky=true and the 
+                            CLF1                                    !   clear=false/pq=0 case
 
                      aerosol(ig,l,iaer) = -log(1 - CLF2 + CLF2*exp(-aerosol0))