Changeset 517 for LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F

Timestamp:

Apr 16, 2004, 5:43:38 PM (20 years ago)

Author:

lmdzadmin

Message:

Inclusion des modifications de O. Boucher et de J. Quaas pour le calcul des
premiers effets directs et indirects dus aux aerosols
LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F (modified) (12 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F

@@ -1 @@
-c
+C
 c $Header$
 c
@@ -143 +143 @@
       LOGICAL ok_region ! sortir le fichier regional
       PARAMETER (ok_region=.FALSE.)
+c
+c
+      LOGICAL ok_polder ! sortir échantillonné de manière POLDER
+      save ok_polder
 c======================================================================
 c
@@ -877 +881 @@
       REAL zx_lon(iim,jjmp1), zx_lat(iim,jjmp1)
 c
-      INTEGER nid_day, nid_mth, nid_ins, nid_nmc
-      SAVE nid_day, nid_mth, nid_ins, nid_nmc
+      INTEGER nid_day, nid_mth, nid_ins, nid_nmc, nid_pol
+      SAVE nid_day, nid_mth, nid_ins, nid_nmc, nid_pol
 c
       INTEGER nhori, nvert
@@ -928 +932 @@
       REAL zu10m(klon), zv10m(klon)           !vents a 10m moyennes s/1 maille
       CHARACTER*40 t2mincels, t2maxcels       !t2m min., t2m max
+cjq   Aerosol effects (Johannes Quaas, 27/11/2003)
+      REAL sulfate(klon, klev) ! SO4 aerosol concentration [ug/m3]
+      REAL sulfate_pi(klon, klev) ! SO4 aerosol concentration [ug/m3] (pre-industrial value)
+      SAVE sulfate_pi
+
+      REAL cldtaupi(klon,klev)  ! Cloud optical thickness for pre-industrial (pi) aerosols
+
+      REAL re(klon, klev)       ! Cloud droplet effective radius
+      REAL fl(klon, klev)  ! denominator of re
+
+      REAL re_top(klon), fl_top(klon) ! CDR at top of liquid water clouds
+
+      ! Aerosol optical properties
+      REAL tau_ae(klon,klev,2), piz_ae(klon,klev,2)
+      REAL cg_ae(klon,klev,2)
+
+      REAL topswad(klon), solswad(klon) ! Aerosol direct effect.
+      ! ok_ade=T -ADE=topswad-topsw
+
+      REAL topswai(klon), solswai(klon) ! Aerosol indirect effect.
+      ! ok_aie=T ->
+      !        ok_ade=T -AIE=topswai-topswad
+      !        ok_ade=F -AIE=topswai-topsw
+
+      ! For POLDER swath
+      INTEGER pyr, pmo, pday    ! Year, month and day
+      INTEGER poldermask(klon)  ! POLDER swath mask (0 or 1)
+
+      REAL aerindex(klon)       ! POLDER aerosol index
+     
+      ! Parameters
+      LOGICAL ok_ade, ok_aie    ! Apply aerosol (in)direct effects or not
+      REAL bl95_b0, bl95_b1   ! Parameter in Boucher and Lohmann (1995)
+cjq-end
 c
 c Declaration des constantes et des fonctions thermodynamiques
@@ -991 +1029 @@
          call conf_phys(ocean, ok_veget, ok_journe, ok_mensuel,
      .                  ok_instan, fact_cldcon, facttemps,ok_newmicro,
-     .                  iflag_cldcon,ratqsbas,ratqshaut, if_ebil)
+     .                  iflag_cldcon,ratqsbas,ratqshaut, if_ebil,
+     .                  ok_ade, ok_aie, 
+     .                  bl95_b0, bl95_b1)
 cIM  .                  , RI0)
 
@@ -1145 +1185 @@
 
 #undef histmthNMC
-#define histmthNMC
+cccccccc#define histmthNMC
 #ifdef histmthNMC
 #include "ini_histmthNMC.h"
@@ -1151 +1191 @@
 
 #include "ini_histins.h"
+#include "ini_histpol.h"
 
 #ifdef histREGDYN
@@ -1912 +1953 @@
      .            t_seri, convliq, convfra, dtau_c, dem_c,
      .            cldh_c, cldl_c, cldm_c, cldt_c, cldq_c,
-     .            flwp_c, fiwp_c, flwc_c, fiwc_c)
+     .            flwp_c, fiwp_c, flwc_c, fiwc_c,
+     e            ok_aie,
+     e            sulfate, sulfate_pi,
+     e            bl95_b0, bl95_b1,
+     s            cldtaupi, re, fl)
 c
 cIM calcul tau. emi nuages startiformes
@@ -1918 +1963 @@
      .            t_seri, cldliq, cldfra, dtau_s, dem_s,
      .            cldh_s, cldl_s, cldm_s, cldt_s, cldq_s,
-     .            flwp_s, fiwp_s, flwc_s, fiwc_s)
+     .            flwp_s, fiwp_s, flwc_s, fiwc_s,
+     e            ok_aie,
+     e            sulfate, sulfate_pi,
+     e            bl95_b0, bl95_b1,
+     s            cldtaupi, re, fl)
 c
       cldtot(:,:)=min(max(cldfra(:,:),rnebcon(:,:)),1.)
@@ -2185 +2234 @@
       ENDDO
       ENDDO
-c
+cjq - introduce the aerosol direct and first indirect radiative forcings
+cjq - Johannes Quaas, 27/11/2003 (quaas@lmd.jussieu.fr)
+      IF (ok_ade.OR.ok_aie) THEN
+         ! Get sulfate aerosol distribution
+         CALL readsulfate(rjourvrai, debut, sulfate)
+         CALL readsulfate_preind(rjourvrai, debut, sulfate_pi)
+
+         ! Calculate aerosol optical properties (Olivier Boucher)
+         CALL aeropt(pplay, paprs, t_seri, sulfate, rhcl,
+     .        tau_ae, piz_ae, cg_ae, aerindex)
+      ENDIF
+
+c     
 c Calculer les parametres optiques des nuages et quelques
 c parametres pour diagnostiques:
@@ -2193 +2254 @@
      .            t_seri, cldliq, cldfra, cldtau, cldemi,
      .            cldh, cldl, cldm, cldt, cldq,
-     .            flwp, fiwp, flwc, fiwc)
+     .            flwp, fiwp, flwc, fiwc,
+     e            ok_aie,
+     e            sulfate, sulfate_pi,
+     e            bl95_b0, bl95_b1,
+     s            cldtaupi, re, fl)
       else
       CALL nuage (paprs, pplay,
      .            t_seri, cldliq, cldfra, cldtau, cldemi,
-     .            cldh, cldl, cldm, cldt, cldq)
+     .            cldh, cldl, cldm, cldt, cldq,
+     e            ok_aie,
+     e            sulfate, sulfate_pi,
+     e            bl95_b0, bl95_b1,
+     s            cldtaupi, re, fl)
+      
       endif
 c
@@ -2233 +2303 @@
      s             topsw0,toplw0,solsw0,sollw0,
      s             lwdn0, lwdn, lwup0, lwup, 
-     s             swdn0, swdn, swup0, swup     )
+     s             swdn0, swdn, swup0, swup,
+     e             ok_ade, ok_aie, ! new for aerosol radiative effects
+     e             tau_ae, piz_ae, cg_ae, ! ="=
+     s             topswad, solswad, ! ="=
+     e             cldtaupi, ! ="=
+     s             topswai, solswai) ! ="=
       itaprad = 0
       ENDIF