Changeset 1712 for LMDZ5/trunk/libf/phylmd/physiq.F

Timestamp:

Jan 18, 2013, 3:07:29 PM (12 years ago)

Author:

idelkadi

Message:

Corrections dans newmicro.F pour completer les modifications de O. Boucher sur ok_ade/ok_aie avec correction des diagnostiques CMIP5 :
Rajout d'un nouveau flag ok_cdnc (ok cloud droplet number concentration)
Dans le cas sans aérosols, nous avons flag_aerosol=0, ok_cdnc=n, ok_ade=n et ok_aie=n
Dans le cas avec aérosols, nous avons flag_aerosol=6, ok_cdnc=y, ok_ade=y et ok_aie=y

File:

• LMDZ5/trunk/libf/phylmd/physiq.F (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F

@@ -1094 +1094 @@
       ! Parameters
       LOGICAL ok_ade, ok_aie    ! Apply aerosol (in)direct effects or not
+      LOGICAL ok_cdnc          ! ok cloud droplet number concentration (O. Boucher 01-2013)
       REAL bl95_b0, bl95_b1   ! Parameter in Boucher and Lohmann (1995)
-      SAVE ok_ade, ok_aie, bl95_b0, bl95_b1
-c$OMP THREADPRIVATE(ok_ade, ok_aie, bl95_b0, bl95_b1)
+      SAVE ok_ade, ok_aie, ok_cdnc, bl95_b0, bl95_b1
+c$OMP THREADPRIVATE(ok_ade, ok_aie, ok_cdnc, bl95_b0, bl95_b1)
       LOGICAL, SAVE :: aerosol_couple ! true  : calcul des aerosols dans INCA
                                       ! false : lecture des aerosol dans un fichier
@@ -1251 +1252 @@
      .     fact_cldcon, facttemps,ok_newmicro,iflag_radia,
      .     iflag_cldcon,iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs,
-     .     ok_ade, ok_aie, aerosol_couple, 
+     .     ok_ade, ok_aie, ok_cdnc, aerosol_couple, 
      .     flag_aerosol, new_aod,
      .     bl95_b0, bl95_b1,
@@ -3145 +3146 @@
 
       if (ok_newmicro) then
-      CALL newmicro (paprs, pplay,ok_newmicro,
-     .            t_seri, cldliq, cldfra, cldtau, cldemi,
-     .            cldh, cldl, cldm, cldt, cldq,
-     .            flwp, fiwp, flwc, fiwc,
-     e            ok_aie,
-     e            mass_solu_aero, mass_solu_aero_pi,
-     e            bl95_b0, bl95_b1,
-     s            cldtaupi, re, fl, ref_liq, ref_ice)
+      CALL newmicro (ok_cdnc, bl95_b0, bl95_b1,
+     .              paprs, pplay, t_seri, cldliq, cldfra,
+     .              cldtau, cldemi, cldh, cldl, cldm, cldt, cldq,
+     e              flwp, fiwp, flwc, fiwc,
+     e              mass_solu_aero, mass_solu_aero_pi,
+     s              cldtaupi, re, fl, ref_liq, ref_ice) 
       else
       CALL nuage (paprs, pplay,
@@ -3161 +3160 @@
      e            bl95_b0, bl95_b1,
      s            cldtaupi, re, fl)
-      
       endif
 c