Changeset 5910

Timestamp:

Dec 4, 2025, 3:14:38 PM (2 hours ago)

Author:

evignon

Message:

commentaire de lmdz_cloud_optics_prop suite au dernier atelier nuages-rayonnement

File:

• LMDZ6/trunk/libf/phylmd/lmdz_cloud_optics_prop.f90 (modified) (8 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_cloud_optics_prop.f90

@@ -18 +18 @@
     IF (ANY(first)) THEN
       IF (novlp==1) THEN
-        WRITE (*, *) 'Hypothese de recouvrement: MAXIMUM_ &
+        WRITE (*, *) 'Hypothesis on cloud overlap: MAXIMUM_ &
                 &                                             &
                 &                                          RANDOM'
         first_first = .FALSE.
       ELSEIF (novlp==2) THEN
-        WRITE (*, *) 'Hypothese de recouvrement: MAXIMUM'
+        WRITE (*, *) 'Hypothesis on cloud overlap: MAXIMUM'
         first_first = .FALSE.
       ELSEIF (novlp==3) THEN
-        WRITE (*, *) 'Hypothese de recouvrement: RANDOM'
+        WRITE (*, *) 'Hypothesis on cloud overlap: RANDOM'
         first_first = .FALSE.
       ENDIF
@@ -75 +75 @@
   !          M. Coulon-Decorzens: replayisation of the routine + cleaning
   !                               and commentaries
+  !          "atelier optics of clouds": careful reading and comments 
+  !                                   things to address -> see ?aoc
   !
-  ! Aim: compute condensate optical properties,
-  !      cloud optical depth and emissivity
+  ! Aim: compute condensates' optical properties
+  !      (cloud optical depth,) and emissivity ?aoc
   ! ======================================================================
   
   ! List of arguments
   !------------------
+
+  ! ?aoc all the pi arguments should be optional
 
   ! input:
@@ -92 +96 @@
   REAL, INTENT(IN) :: rnebcon(klon,klev) ! convection cloud fraction [-]
   REAL, INTENT(IN) :: ccwcon(klon,klev)  ! condensed water from deep convection [kg/kg]
-  ! jq for the aerosol indirect effect
-  ! jq introduced by Johannes Quaas (quaas@lmd.jussieu.fr), 27/11/2003
   REAL, INTENT(IN) :: mass_solu_aero(klon, klev)    ! total mass concentration for all soluble aerosols [ug m-3]
   REAL, INTENT(IN) :: mass_solu_aero_pi(klon, klev) ! - (pre-industrial value)
   REAL, INTENT(IN)  :: dNovrN(klon)         ! enhancement factor for cdnc
-  REAL, INTENT(OUT) :: distcltop(klon,klev) ! distance from large scale cloud top [m]
-  REAL, INTENT(OUT) :: temp_cltop(klon,klev)!temperature at large scale cloud top [K]
+  REAL, INTENT(OUT) :: distcltop(klon,klev) ! distance from large scale cloud top [m] ?aoc
+  REAL, INTENT(OUT) :: temp_cltop(klon,klev)! temperature at large scale cloud top [K] ?aoc
 
   LOGICAL, INTENT(IN) :: ptconv(klon, klev) ! flag for grid points affected by deep convection
 
   ! inout:
-  REAL, INTENT(INOUT) :: pclc(klon, klev) ! cloud fraction for radiation [-]
+  REAL, INTENT(INOUT) :: pclc(klon, klev) ! cloud fraction for radiation [-] ?aoc
+  REAL, INTENT(INOUT) :: icefrac_optics(klon, klev)! ice fraction in clouds seen by radiation [-]
+
 
   ! out:
+  ! only diagnostics:
   REAL, INTENT(OUT) :: pct(klon)      ! 2D total cloud cover [-]
   REAL, INTENT(OUT) :: pcl(klon)      ! 2D low cloud cover [-]
@@ -119 +124 @@
                                       ! introduced to avoid problems in the averaging of the output
                                       ! water clouds within a grid cell
-
+ 
+  ! diagnostics or properties if one uses oldrad
   REAL, INTENT(OUT) :: pcltau(klon, klev) ! cloud optical depth [m]
   REAL, INTENT(OUT) :: pclemi(klon, klev) ! cloud emissivity [-]
@@ -125 +131 @@
                                              ! values of optical thickness that does not account
                                              ! for aerosol effects on cloud droplet radius [m]
-
-  REAL, INTENT(OUT) :: reliq(klon, klev)   ! liquid droplet effective radius [m] 
-  REAL, INTENT(OUT) :: reice(klon, klev)   ! ice effective radius [m]
-  REAL, INTENT(OUT) :: reliq_pi(klon, klev)! liquid droplet effective radius [m], pre-industrial
-  REAL, INTENT(OUT) :: reice_pi(klon, klev)! ice effective radius [m], pre-industrial
   REAL, INTENT(OUT) :: scdnc(klon, klev)   ! cloud droplet number concentration, mean over the whole mesh [m-3]
   REAL, INTENT(OUT) :: cldncl(klon)        ! cloud droplet number concentration at top of cloud [m-3]
@@ -142 +143 @@
   REAL, INTENT(OUT) :: icc3dcon(klon, klev)! cloud fraction for liquid part only, convective clouds [-]
   REAL, INTENT(OUT) :: icc3dstra(klon, klev)! cloud fraction for ice part only, stratiform clouds [-] 
-  REAL, INTENT(INOUT) :: icefrac_optics(klon, klev)! ice fraction in clouds seen by radiation [-]
-
+
+  ! actual optical variables used as inputs for the radiative scheme
+  REAL, INTENT(OUT) :: reliq(klon, klev)   ! liquid droplet effective radius [m] 
+  REAL, INTENT(OUT) :: reice(klon, klev)   ! ice effective radius [m]
+  REAL, INTENT(OUT) :: reliq_pi(klon, klev)! liquid droplet effective radius [m], pre-industrial
+  REAL, INTENT(OUT) :: reice_pi(klon, klev)! ice effective radius [m], pre-industrial
+  
+  
+  
   ! Local variables
   !----------------
@@ -158 +166 @@
   REAL rel, tc, rei, iwc, dei, deimin, deimax
   REAL k_ice
-
-  ! jq for the aerosol indirect effect
-  ! jq introduced by Johannes Quaas (quaas@lmd.jussieu.fr), 27/11/2003
   REAL cdnc(klon, klev) ! cloud droplet number concentration [m-3]
   REAL cdnc_pi(klon, klev) ! cloud droplet number concentration [m-3] (pi value)
@@ -222 +227 @@
     DO k = 1, klev
 
-
-!!$      IF (ok_new_lscp) THEN
-!!$          CALL icefrac_lscp(klon,temp(:,k),iflag_ice_thermo,distcltop(:,k),temp_cltop(:,k),icefrac_optics(:,k),dzfice(:,k))
-!!$      ELSE
-!!$          CALL icefrac_lsc(klon,temp(:,k),pplay(:,k)/paprs(:,1),icefrac_optics(:,k))
-!!$      ENDIF
 
       DO i = 1, klon