Changeset 1764 for LMDZ5/trunk/libf

Timestamp:

Jun 10, 2013, 3:40:50 PM (11 years ago)

Author:

Laurent Fairhead

Message:

Inclusion d'une routine qui lit des champs d'aérosols stratosphériques
mensuels, prescrit des propriétés optiques et modifie le rayonnement en
conséquence. Pour le moment, seule l'interaction avec le rayonnement ondes
courtes est pris en compte. Les fichiers d'input doivent être au format des
fichiers de sortie. Contrôlé par la variable logique: flag_aerosol_strat
(false par défaut dans DefLists?/config.def)

15. Boucher

---

A new routine has been added to the code that reads in monthly stratospheric
aerosols, prescribes optical properties and modifies radiation accordingly.
Presently, only the interaction with short wave radiation is taken into account.
Input files must be formatted as are the aerosol output fields. Control is by
the logical flag: flag_aerosol_strat (which is false by default and included
DefLists?/config.def)

15. Boucher

Location:

LMDZ5/trunk/libf/phylmd

Files:

• aero_mod.F90 (modified) (5 diffs)
• conf_phys_m.F90 (modified) (7 diffs)
• etat0_netcdf.F90 (modified) (2 diffs)
• phys_output_mod.F90 (modified) (7 diffs)
• phys_output_write.h (modified) (1 diff)
• physiq.F (modified) (6 diffs)
• radlwsw_m.F90 (modified) (5 diffs)
• readaerosolstrato.F90 (added)
• sw_aeroAR4.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/aero_mod.F90

@@ -5 +5 @@
 
   ! Total number of aerosols
-  INTEGER, PARAMETER :: naero_tot = 10 
+!  INTEGER, PARAMETER :: naero_tot = 10 
+!--STRAT AER
+  INTEGER, PARAMETER :: naero_tot = 11
 
   ! Identification number used in aeropt_2bands and aeropt_5wv
@@ -19 +21 @@
   INTEGER, PARAMETER :: id_AIBCM    = 9
   INTEGER, PARAMETER :: id_AIPOMM   = 10
+!--STRAT AER
+  INTEGER, PARAMETER :: id_strat   = 11
+
 
   ! Total number of aerosols actually used in LMDZ 
@@ -31 +36 @@
   ! 9 =  AIBCM
   !10 =  AIPOMM
-  INTEGER, PARAMETER :: naero_spc = 10
+!--STRAT AER
+  !11 = aerosols stratos
+!  INTEGER, PARAMETER :: naero_spc = 10
+  INTEGER, PARAMETER :: naero_spc = 11
 
   ! Corresponding names for the aerosols
@@ -44 +52 @@
        "CIDUSTM", &
        "AIBCM  ", &
-       "AIPOMM " /)
+!       "AIPOMM " /)
+       "AIPOMM ", &
+       "STRAT  " /)
 
 
@@ -65 +75 @@
   INTEGER, parameter :: nbands = 2
 
-
 END MODULE aero_mod

LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

@@ -18 +18 @@
                        iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs, &
                        ok_ade, ok_aie, ok_cdnc, aerosol_couple, &
-                       flag_aerosol, new_aod, &
+                       flag_aerosol, flag_aerosol_strat, new_aod, &
                        bl95_b0, bl95_b1,&
                        read_climoz, &
@@ -60 +60 @@
 ! ok_ade, ok_aie: apply or not aerosol direct and indirect effects
 ! ok_cdnc, ok cloud droplet number concentration
+! flag_aerosol_strat : flag pour les aerosols stratos
 ! bl95_b*: parameters in the formula to link CDNC to aerosol mass conc 
 !
@@ -72 +73 @@
   LOGICAL              :: ok_ade, ok_aie, ok_cdnc, aerosol_couple
   INTEGER              :: flag_aerosol
+  LOGICAL              :: flag_aerosol_strat
   LOGICAL              :: new_aod
   REAL                 :: bl95_b0, bl95_b1
@@ -87 +89 @@
   LOGICAL,SAVE        :: ok_ade_omp, ok_aie_omp, ok_cdnc_omp, aerosol_couple_omp
   INTEGER, SAVE       :: flag_aerosol_omp
+  INTEGER, SAVE       :: flag_aerosol_strat_omp
   LOGICAL, SAVE       :: new_aod_omp
   REAL,SAVE           :: bl95_b0_omp, bl95_b1_omp
@@ -307 +310 @@
   flag_aerosol_omp = 0
   CALL getin('flag_aerosol',flag_aerosol_omp)
+!
+!Config Key  = flag_aerosol_strat
+!Config Desc = use stratospheric aerosols T/F
+!Config Def  = false
+!Config Help = Used in physiq.F
+!
+!
+  flag_aerosol_strat_omp = .false.
+  CALL getin('flag_aerosol_strat',flag_aerosol_strat_omp)
 
 ! Temporary variable for testing purpose!!
@@ -1717 +1729 @@
     aerosol_couple = aerosol_couple_omp
     flag_aerosol=flag_aerosol_omp
+    flag_aerosol_strat=flag_aerosol_strat_omp
     new_aod=new_aod_omp
     aer_type = aer_type_omp
@@ -1901 +1914 @@
   write(lunout,*)' aerosol_couple = ', aerosol_couple
   write(lunout,*)' flag_aerosol = ', flag_aerosol
+  write(lunout,*)' flag_aerosol_strat = ', flag_aerosol_strat
   write(lunout,*)' new_aod = ', new_aod
   write(lunout,*)' aer_type = ',aer_type

LMDZ5/trunk/libf/phylmd/etat0_netcdf.F90

@@ -101 +101 @@
   LOGICAL :: ok_LES, ok_ade, ok_aie, ok_cdnc, aerosol_couple, new_aod, callstats
   INTEGER :: iflag_radia, flag_aerosol
+  LOGICAL :: flag_aerosol_strat
   REAL    :: bl95_b0, bl95_b1, fact_cldcon, facttemps, ratqsbas, ratqshaut
   REAL    :: tau_ratqs
@@ -137 +138 @@
                    iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs,            &
                    ok_ade, ok_aie, ok_cdnc, aerosol_couple,             &
-                   flag_aerosol, new_aod,                               &
+                   flag_aerosol, flag_aerosol_strat, new_aod,           &
                    bl95_b0, bl95_b1,                                    &
                    read_climoz,                                         &

LMDZ5/trunk/libf/phylmd/phys_output_mod.F90

@@ -410 +410 @@
   type(ctrl_out),save :: o_solswai      = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'solswai')
 
-  type(ctrl_out),save,dimension(10) :: o_tausumaero  = (/ ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_ASBCM'), &
+!  type(ctrl_out),save,dimension(10) :: o_tausumaero  = (/ ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_ASBCM'), &
+  type(ctrl_out),save,dimension(11) :: o_tausumaero  = & 
+    (/ ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_ASBCM'), &
        ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_ASPOMM'), &
        ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_ASSO4M'), &
@@ -419 +421 @@
        ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_CIDUSTM'), &
        ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_AIBCM'), &
-       ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_AIPOMM') /)
+       ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_AIPOMM'), &
+       ctrl_out((/ 2, 2, 10, 10, 10, 10 /),'OD550_STRAT') /)
 
   type(ctrl_out),save :: o_od550aer     = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'od550aer')
@@ -579 +582 @@
   type(ctrl_out),save :: o_e_th         = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'e_th')
   type(ctrl_out),save :: o_w_th         = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'w_th')
+  type(ctrl_out),save :: o_lambda_th    = ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'lambda_th')
   type(ctrl_out),save :: o_ftime_th     = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'ftime_th')
   type(ctrl_out),save :: o_q_th         = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'q_th')
@@ -663 +667 @@
        ok_hf,ok_instan,ok_LES,ok_ade,ok_aie, read_climoz, &
        phys_out_filestations, &
-       new_aod, aerosol_couple)   
+       new_aod, aerosol_couple, flag_aerosol_strat)   
 
     USE iophy 
@@ -693 +697 @@
     integer                               :: nbteta, nlevSTD, radpas
     logical                               :: ok_mensuel, ok_journe, ok_hf, ok_instan
-    logical                               :: ok_LES,ok_ade,ok_aie
+    logical                               :: ok_LES,ok_ade,ok_aie,flag_aerosol_strat
     logical                               :: new_aod, aerosol_couple
     integer, intent(in)::  read_climoz ! read ozone climatology
@@ -1165 +1169 @@
                 CALL histdef2d(iff,clef_stations(iff), &
                      o_loaddust%flag,o_loaddust%name,"Column Load of Dust ","kg/m2")
-
+!--STRAT AER
+             ENDIF
+             IF (ok_ade.OR.ok_aie.OR.flag_aerosol_strat) THEN
                 DO naero = 1, naero_spc
                    CALL histdef2d(iff,clef_stations(iff), &
@@ -1650 +1656 @@
              CALL histdef3d(iff,clef_stations(iff),o_e_th%flag,o_e_th%name,"Thermal plume entrainment","K/s")
              CALL histdef3d(iff,clef_stations(iff),o_w_th%flag,o_w_th%name,"Thermal plume vertical velocity","m/s")
+             CALL histdef3d(iff,clef_stations(iff), &
+                  o_lambda_th%flag,o_lambda_th%name,"Thermal plume vertical velocity","m/s")
              CALL histdef2d(iff,clef_stations(iff), &
                   o_ftime_th%flag,o_ftime_th%name,"Fraction of time Shallow convection occurs"," ")

LMDZ5/trunk/libf/phylmd/phys_output_write.h

@@ -1305 +1305 @@
           ENDIF
           
+c--STRAT AER
+          endif
+          IF (ok_ade.OR.ok_aie.OR.flag_aerosol_strat) THEN
           DO naero = 1, naero_spc
             IF (o_tausumaero(naero)%flag(iff)<=lev_files(iff)) THEN

LMDZ5/trunk/libf/phylmd/physiq.F

@@ -1120 +1120 @@
       LOGICAL, SAVE :: new_aod
 c$OMP THREADPRIVATE(new_aod) 
-   
+c
+c--STRAT AEROSOL
+      LOGICAL, SAVE :: flag_aerosol_strat
+c$OMP THREADPRIVATE(flag_aerosol_strat)
+cc-fin STRAT AEROSOL
 c
 c Declaration des constantes et des fonctions thermodynamiques
@@ -1270 +1274 @@
      .     iflag_cldcon,iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs,
      .     ok_ade, ok_aie, ok_cdnc, aerosol_couple, 
-     .     flag_aerosol, new_aod,
+     .     flag_aerosol, flag_aerosol_strat, new_aod,
      .     bl95_b0, bl95_b1,
 c     nv flags pour la convection et les poches froides
@@ -1595 +1599 @@
      &                       ok_hf,ok_instan,ok_LES,ok_ade,ok_aie, 
      &                       read_climoz, phys_out_filestations,
-     &                       new_aod, aerosol_couple
-     &                        )
+     &                       new_aod, aerosol_couple,
+     &                       flag_aerosol_strat )
 c$OMP END MASTER
 c$OMP BARRIER
@@ -2972 +2976 @@
          cg_aero(:,:,:,:)  = 0.
       ENDIF
+c
+c--STRAT AEROSOL
+c--updates tausum_aero,tau_aero,piz_aero,cg_aero
+      IF (flag_aerosol_strat) THEN
+         PRINT *,'appel a readaerosolstrat', mth_cur
+         CALL readaerosolstrato(debut)
+      ENDIF
+c--fin STRAT AEROSOL
 
 cIM calcul nuages par le simulateur ISCCP
@@ -3359 +3371 @@
      e        t_seri,q_seri,wo,
      e        cldfrarad, cldemirad, cldtaurad,
-     e        ok_ade, ok_aie, flag_aerosol,
+     e        ok_ade.OR.flag_aerosol_strat, ok_aie, flag_aerosol,
+     e        flag_aerosol_strat,
      e        tau_aero, piz_aero, cg_aero,
      e        cldtaupirad,new_aod,
@@ -3401 +3414 @@
      e        t_seri,q_seri,wo,
      e        cldfra, cldemi, cldtau,
-     e        ok_ade, ok_aie, flag_aerosol,
+     e        ok_ade.OR.flag_aerosol_strat, ok_aie, flag_aerosol,
+     e        flag_aerosol_strat,
      e        tau_aero, piz_aero, cg_aero,
      e        cldtaupi,new_aod,

LMDZ5/trunk/libf/phylmd/radlwsw_m.F90

@@ -11 +11 @@
    cldfra, cldemi, cldtaupd,&
    ok_ade, ok_aie, flag_aerosol,&
+   flag_aerosol_strat,&
    tau_aero, piz_aero, cg_aero,&
    cldtaupi, new_aod, &
@@ -57 +58 @@
   ! ok_aie---input-L- apply the Aerosol Indirect Effect or not?
   ! flag_aerosol-input-I- aerosol flag from 0 to 6
+  ! flag_aerosol_strat-input-I- use stratospheric aerosols flag (T/F)
   ! tau_ae, piz_ae, cg_ae-input-R- aerosol optical properties (calculated in aeropt.F)
   ! cldtaupi-input-R- epaisseur optique des nuages dans le visible
@@ -121 +123 @@
   LOGICAL, INTENT(in)  :: ok_ade, ok_aie                                 ! switches whether to use aerosol direct (indirect) effects or not
   INTEGER, INTENT(in)  :: flag_aerosol                                   ! takes value 0 (no aerosol) or 1 to 6 (aerosols)
+  LOGICAL, INTENT(in)  :: flag_aerosol_strat                             ! use stratospheric aerosols
   REAL,    INTENT(in)  :: cldfra(KLON,KLEV), cldemi(KLON,KLEV), cldtaupd(KLON,KLEV)
   REAL,    INTENT(in)  :: tau_aero(KLON,KLEV,9,2)                        ! aerosol optical properties (see aeropt.F)
@@ -360 +363 @@
                ztopswadaero,zsolswadaero,&
                ztopswaiaero,zsolswaiaero,& 
-               ok_ade, ok_aie, flag_aerosol) 
+               ok_ade, ok_aie) 
           
        ELSE ! new_aod=T         
@@ -379 +382 @@
                zsolsw_aero,zsolsw0_aero,&
                ztopswcf_aero,zsolswcf_aero, & 
-               ok_ade, ok_aie, flag_aerosol) 
+               ok_ade, ok_aie, flag_aerosol,flag_aerosol_strat) 
        ENDIF
 

LMDZ5/trunk/libf/phylmd/sw_aeroAR4.F90

@@ -18 +18 @@
      PSOLSWAERO,PSOLSW0AERO,&
      PTOPSWCFAERO,PSOLSWCFAERO,&
-     ok_ade, ok_aie, flag_aerosol )
+     ok_ade, ok_aie, flag_aerosol, flag_aerosol_strat )
 
   USE dimphy
@@ -138 +138 @@
 
   LOGICAL ok_ade, ok_aie    ! use aerosol forcings or not?
+  LOGICAL flag_aerosol_strat ! use stratospehric aerosols
   INTEGER flag_aerosol      ! global flag for aerosol 0 (no aerosol) or 1-5 (aerosols)
   REAL(KIND=8) tauaero(kdlon,kflev,9,2)  ! aerosol optical properties
@@ -307 +308 @@
      ENDIF ! swaero_diag .or. .not. AEROSOLFEEDBACK_ACTIVE
 
-     IF (flag_aerosol .GT. 0 ) THEN
+     IF (flag_aerosol .GT. 0 .OR. flag_aerosol_strat) THEN
 
      IF (ok_ade.and.swaero_diag .or. .not. ok_ade) THEN
@@ -498 +499 @@
      ENDIF ! ok_aie      
 
-     ENDIF !--if flag_aerosol GT 0
+     ENDIF !--if flag_aerosol GT 0 OR flag_aerosol_strat
 
      itapsw = 0
@@ -504 +505 @@
   itapsw = itapsw + 1
 
-  IF  ( AEROSOLFEEDBACK_ACTIVE .AND. flag_aerosol .GT. 0 ) THEN
+  IF  ( AEROSOLFEEDBACK_ACTIVE .AND. (flag_aerosol .GT. 0 .OR. flag_aerosol_strat) ) THEN
   IF ( ok_ade .and. ok_aie  ) THEN
     ZFSUP(:,:) =    ZFSUP_AERO(:,:,4)