Index: /LMDZ6/trunk/DefLists/CMIP6_ping_atmos.xml
===================================================================
--- /LMDZ6/trunk/DefLists/CMIP6_ping_atmos.xml	(revision 3111)
+++ /LMDZ6/trunk/DefLists/CMIP6_ping_atmos.xml	(revision 3112)
@@ -65,7 +65,7 @@
    <field id="CMIP6_clhcalipso"    field_ref="pclhcalipso"       /> <!-- P1 (%) cloud_area_fraction_in_atmosphere_layer : Percentage cloud cover in layer centred on 220hPa -->
    <field id="CMIP6_cli"           field_ref="iwcon"            /> <!-- P1 (kg kg-1) mass_fraction_of_cloud_ice_in_air : Includes both large-scale and convective cloud. This is calculated as the mass of cloud ice in the grid cell divided by the mass of air (including the water in all phases) in the grid cell. It includes precipitating hydrometeors ONLY if the precipitating hydrometeors affect the calculation of radiative transfer in model. -->
-   <field id="CMIP6_clic"          field_ref="dummy_XYA"        /> <!-- P2 (1.0) mass_fraction_of_convective_cloud_ice_in_air : Calculated as the mass of convective cloud ice  in the grid cell divided by the mass of air (including the water in all phases) in the grid cell.  This includes precipitating hydrometeors ONLY if the precipitating hydrometeors affect the calculation of radiative transfer in model. -->
+   <field id="CMIP6_clic"          field_ref="icc3dcon"        /> <!-- P2 (1.0) mass_fraction_of_convective_cloud_ice_in_air : Calculated as the mass of convective cloud ice  in the grid cell divided by the mass of air (including the water in all phases) in the grid cell.  This includes precipitating hydrometeors ONLY if the precipitating hydrometeors affect the calculation of radiative transfer in model. -->
    <field id="CMIP6_climodis"      field_ref="climodis"         /> <!-- P1 (%) ice_cloud_area_fraction : MODIS Ice Cloud Area Percentage -->
-   <field id="CMIP6_clis"          field_ref="dummy_XYA"        /> <!-- P2 (1.0) mass_fraction_of_stratiform_cloud_ice_in_air : Calculated as the mass of stratiform cloud ice  in the grid cell divided by the mass of air (including the water in all phases) in the grid cell.  This includes precipitating hydrometeors ONLY if the precipitating hydrometeors affect the calculation of radiative transfer in model. -->
+   <field id="CMIP6_clis"          field_ref="icc3dstra"        /> <!-- P2 (1.0) mass_fraction_of_stratiform_cloud_ice_in_air : Calculated as the mass of stratiform cloud ice  in the grid cell divided by the mass of air (including the water in all phases) in the grid cell.  This includes precipitating hydrometeors ONLY if the precipitating hydrometeors affect the calculation of radiative transfer in model. -->
    <field id="CMIP6_clisccp"       field_ref="clisccp2"         /> <!-- P1 (%) cloud_area_fraction_in_atmosphere_layer : Percentage cloud cover in optical depth categories. -->
    <field id="CMIP6_clivi"         field_ref="iwp"              /> <!-- P1 (kg m-2) atmosphere_cloud_ice_content : mass of ice water in the column divided by the area of the column (not just the area of the cloudy portion of the column). Includes precipitating frozen hydrometeors ONLY if the precipitating hydrometeor affects the calculation of radiative transfer in model. -->
@@ -168,5 +168,6 @@
    <field id="CMIP6_od550so4so"    field_ref="dummy_not_provided"         /> <!-- P1 (1.0) missing : Balkanski - LSCE -->
    <field id="CMIP6_od865dust"     field_ref="dummy_not_provided"         /> <!-- P1 (1.0) atmosphere_optical_thickness_due_to_dust_ambient_aerosol_particles : Balkanski - LSCE -->
-   <field id="CMIP6_parasolRefl_sea" field_ref="parasolRefl" parasolRefl*fract_oce </field> <!-- P1 (1.0) toa_bidirectional_reflectance : Simulated reflectance from PARASOL as seen at the top of the atmosphere for 5 solar zenith angles. Valid only over ocean and for one viewing direction (viewing zenith angle of 30 degrees and relative azimuth angle 320 degrees). -->
+   <field id="CMIP6_parasolRefl_sea" field_ref="parasol_refl" parasol_refl*fract_oce </field> <!-- P1 (1.0) toa_bidirectional_reflectance : Simulated reflectance from PARASOL as seen at the top of the atmosphere for 5 solar zenith angles. Valid only over ocean and for one viewing direction (viewing zenith angle of 30 degrees and relative azimuth angle 320 degrees). -->
+   <field id="HomeCMIP6_parasolCRefl" field_ref="parasol_crefl" /> <!-- P4 (1.0) toa_bidirectional_reflectance (integral) : Simulated reflectance from PARASOL as seen at the top of the atmosphere integrated over all 5 solar zenith angles. -->
    <field id="CMIP6_pctisccp"      field_ref="ctpisccp"         /> <!-- P1 (Pa) air_pressure_at_cloud_top : ISCCP Mean Cloud Top Pressure. Time-means are weighted by the ISCCP Total Cloud Fraction {:cltisccp} - see  http://cfmip.metoffice.com/COSP.html -->
    <field id="CMIP6_pfull"         field_ref="pres"             /> <!-- P1 (Pa) air_pressure : Air pressure on model levels -->
@@ -230,6 +231,6 @@
    <field id="CMIP6_rsdscsafbnd"   field_ref="dummy_not_provided"         /> <!-- P1 (W m-2) surface_downwelling_shortwave_flux_in_air_assuming_clean_clear_sky : Calculated in the absence of aerosols and clouds, following Ghan (2013, ACP). This requires a double-call in the radiation code with precisely the same meteorology. -->
    <field id="CMIP6_rsdscsbnd"     field_ref="dummy_not_provided"         /> <!-- P1 (W m-2) surface_downwelling_shortwave_flux_in_air_assuming_clear_sky : Calculated with aerosols but without clouds. This is a standard clear-sky calculation -->
-   <field id="CMIP6_rsdscsdiff"    field_ref="dummy_XY"         /> <!-- P1 (W m-2) surface_diffuse_downwelling_shortwave_flux_in_air_assuming_clear_sky : unset -->
-   <field id="CMIP6_rsdsdiff"      field_ref="dummy_XY"         /> <!-- P1 (W m-2) surface_diffuse_downwelling_shortwave_flux_in_air : unset -->
+   <field id="CMIP6_rsdscsdiff"    field_ref="rsdscsdiff"         /> <!-- P1 (W m-2) surface_diffuse_downwelling_shortwave_flux_in_air_assuming_clear_sky : unset -->
+   <field id="CMIP6_rsdsdiff"      field_ref="rsdsdiff"         /> <!-- P1 (W m-2) surface_diffuse_downwelling_shortwave_flux_in_air : unset -->
    <field id="CMIP6_rsdt"          field_ref="SWdnTOA"          /> <!-- P1 (W m-2) toa_incoming_shortwave_flux : Shortwave radiation incident at the top of the atmosphere -->
    <field id="CMIP6_rss"           field_ref="sols"             /> <!-- P1 (W m-2) surface_net_downward_shortwave_flux : Net downward shortwave radiation at the surface -->
Index: /LMDZ6/trunk/DefLists/field_def_lmdz.xml
===================================================================
--- /LMDZ6/trunk/DefLists/field_def_lmdz.xml	(revision 3111)
+++ /LMDZ6/trunk/DefLists/field_def_lmdz.xml	(revision 3112)
@@ -568,4 +568,6 @@
         <field id="lcc3dcon"    long_name="Convective cloud liquid fraction"    unit="1" />
         <field id="lcc3dstra"    long_name="Stratiform cloud liquid fraction"    unit="1" />
+        <field id="icc3dcon"    long_name="Mass Fraction of Convective Cloud Ice"    unit="kg/kg" />
+        <field id="icc3dstra"    long_name="Mass Fraction of Stratiform Cloud Ice "    unit="kg/kg" />
         <field id="clwcon"    long_name="Convective Cloud Liquid water content"    unit="kg/kg" />
         <field id="Ma"    long_name="undilute adiab updraft"    unit="kg/m2/s" />
@@ -804,5 +806,7 @@
 
     <field_group id="fields_solbnd_3D" grid_ref="grid_glo_spectband" >
-      <field id="solbnd"    long_name="Top-of-Atmosphere Solar Insolation for each band"    unit="W m-2" />
+      <field id="solbnd"    long_name="Top-of-Atmosphere Solar Insolation for each band" unit="W m-2" />
+      <field id="rsdsdiff"  long_name="Surface Diffuse Downwelling Shortwave Radiation"  unit="W m-2" />
+      <field id="rsdscsdiff" long_name="Surface Diffuse Downwelling Clear Sky Shortwave Radiation" unit="W m-2" />
     </field_group>
 
Index: /LMDZ6/trunk/DefLists/file_def_histday_lmdz.xml
===================================================================
--- /LMDZ6/trunk/DefLists/file_def_histday_lmdz.xml	(revision 3111)
+++ /LMDZ6/trunk/DefLists/file_def_histday_lmdz.xml	(revision 3112)
@@ -512,4 +512,6 @@
                 <field field_ref="lcc3dcon" level="6" />
                 <field field_ref="lcc3dstra" level="6" />
+                <field field_ref="icc3dcon" level="6" />
+                <field field_ref="icc3dstra" level="6" />
                 <field field_ref="clwcon" level="10" />
                 <field field_ref="Ma" level="10" />
@@ -622,4 +624,12 @@
                 <field field_ref="rldcs4co2" level="10" />
             </field_group>
+
+            <!-- VARS 3D -->
+            <field_group operation="average" freq_op="1ts" grid_ref="grid_glo_spectband">
+                <field field_ref="solbnd" level="4" />
+                <field field_ref="rsdsdiff" level="4" />
+                <field field_ref="rsdscsdiff" level="4" />
+            </field_group>
+
         </file>
     </file_group>
Index: /LMDZ6/trunk/DefLists/file_def_histhf_lmdz.xml
===================================================================
--- /LMDZ6/trunk/DefLists/file_def_histhf_lmdz.xml	(revision 3111)
+++ /LMDZ6/trunk/DefLists/file_def_histhf_lmdz.xml	(revision 3112)
@@ -531,6 +531,8 @@
                 <field field_ref="reffclwc" level="10" />
                 <field field_ref="lcc3d" level="10" />
-                <field field_ref="lcc3dcon" level="10" />
-                <field field_ref="lcc3dstra" level="10" />
+                <field field_ref="lcc3dcon" level="5" />
+                <field field_ref="lcc3dstra" level="5" />
+                <field field_ref="icc3dcon" level="5" />
+                <field field_ref="icc3dstra" level="5" />
                 <field field_ref="clwcon" level="10" />
                 <field field_ref="Ma" level="10" />
@@ -643,4 +645,12 @@
                 <field field_ref="rldcs4co2" level="10" />
             </field_group>
+
+            <!-- VARS 3D -->
+            <field_group operation="average" freq_op="1ts" grid_ref="grid_glo_spectband">
+                <field field_ref="solbnd" level="4" />
+                <field field_ref="rsdsdiff" level="4" />
+                <field field_ref="rsdscsdiff" level="4" />
+            </field_group>
+
         </file>
     </file_group>
Index: /LMDZ6/trunk/DefLists/file_def_histmth_lmdz.xml
===================================================================
--- /LMDZ6/trunk/DefLists/file_def_histmth_lmdz.xml	(revision 3111)
+++ /LMDZ6/trunk/DefLists/file_def_histmth_lmdz.xml	(revision 3112)
@@ -554,4 +554,6 @@
                 <field field_ref="lcc3dcon" level="2" />
                 <field field_ref="lcc3dstra" level="2" />
+                <field field_ref="icc3dcon" level="2" />
+                <field field_ref="icc3dstra" level="2" />
                 <field field_ref="clwcon" level="4" />
                 <field field_ref="Ma" level="4" />
@@ -668,4 +670,6 @@
             <field_group operation="average" freq_op="1ts" grid_ref="grid_glo_spectband">
                 <field field_ref="solbnd" level="4" />
+                <field field_ref="rsdsdiff" level="4" />
+                <field field_ref="rsdscsdiff" level="4" />
             </field_group>
 
Index: /LMDZ6/trunk/libf/phylmd/newmicro.F90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/newmicro.F90	(revision 3111)
+++ /LMDZ6/trunk/libf/phylmd/newmicro.F90	(revision 3112)
@@ -10,5 +10,5 @@
   USE dimphy
   USE phys_local_var_mod, ONLY: scdnc, cldncl, reffclwtop, lcc, reffclws, &
-    reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra
+    reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra
   USE phys_state_var_mod, ONLY: rnebcon, clwcon
   USE icefrac_lsc_mod ! computes ice fraction (JBM 3/14)
@@ -586,4 +586,11 @@
         lcc3dstra(i, k) = lcc3dstra(i, k) - lcc3dcon(i, k) ! eau liquide stratiforme
         lcc3dstra(i, k) = max(lcc3dstra(i,k), 0.0)
+        !
+        icc3dcon(i, k) = rnebcon(i, k)*(1-phase3d(i, k))*clwcon(i, k) ! glace
+                                                                      ! convective
+        icc3dstra(i, k) = pclc(i, k)*pqlwp(i, k)*(1-phase3d(i, k))
+        icc3dstra(i, k) = icc3dstra(i, k) - icc3dcon(i, k) ! glace stratiforme
+        icc3dstra(i, k) = max( icc3dstra(i, k), 0.0) 
+        !
         ! Compute cloud droplet radius as above in meter
         radius = 1.1*((pqlwp(i,k)*pplay(i,k)/(rd*t(i,k)))/(4./3*rpi*1000.* &
@@ -626,4 +633,6 @@
         IF (lcc3dcon(i,k)<=0.0) lcc3dcon(i, k) = 0.0
         IF (lcc3dstra(i,k)<=0.0) lcc3dstra(i, k) = 0.0
+        IF (icc3dcon(i,k)<=0.0) icc3dcon(i, k) = 0.0
+        IF (icc3dstra(i,k)<=0.0) icc3dstra(i, k) = 0.0 
       END DO
       IF (reffclwtop(i)<=0.0) reffclwtop(i) = 0.0
Index: /LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90	(revision 3111)
+++ /LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90	(revision 3112)
@@ -140,4 +140,6 @@
       REAL, SAVE, ALLOCATABLE :: ZFSDN0(:,:),  ZFSUP0(:,:)      ! diag
       !$OMP THREADPRIVATE(ZFSDN0,ZFSUP0)
+      REAL, SAVE, ALLOCATABLE :: PSFSWDIR(:,:), PSFSWDIF(:,:)
+      !$OMP THREADPRIVATE(PSFSWDIR, PSFSWDIF)
 !
       REAL, SAVE, ALLOCATABLE :: stratomask(:,:) 
@@ -169,4 +171,8 @@
       REAL, SAVE, ALLOCATABLE :: lcc3dstra(:,:)
       !$OMP THREADPRIVATE(lcc3dstra)
+      REAL, SAVE, ALLOCATABLE :: icc3dstra(:,:)
+      !$OMP THREADPRIVATE(icc3dstra)
+      REAL, SAVE, ALLOCATABLE :: icc3dcon(:,:)
+      !$OMP THREADPRIVATE(icc3dcon)
       REAL, SAVE, ALLOCATABLE :: od443aer(:) 
       !$OMP THREADPRIVATE(od443aer) 
@@ -516,5 +522,6 @@
 USE phys_state_var_mod
 
-IMPLICIT NONE
+      IMPLICIT NONE
+      include "clesphys.h"
       ALLOCATE(t_seri(klon,klev),q_seri(klon,klev),ql_seri(klon,klev),qs_seri(klon,klev))
       ALLOCATE(u_seri(klon,klev),v_seri(klon,klev))
@@ -564,4 +571,5 @@
       ALLOCATE(ZLWFT0_i(klon,klev+1),ZSWFT0_i(klon,klev+1),ZFLDN0(klon,klev+1))
       ALLOCATE(ZFLUP0(klon,klev+1),ZFSDN0(klon,klev+1),ZFSUP0(klon,klev+1))
+      ALLOCATE(PSFSWDIR(klon,NSW), PSFSWDIF(klon,NSW))
 !
       ALLOCATE(topswad_aero(klon), solswad_aero(klon))
@@ -593,4 +601,6 @@
       ALLOCATE(lcc3dcon(klon, klev))
       ALLOCATE(lcc3dstra(klon, klev))
+      ALLOCATE(icc3dcon(klon, klev))
+      ALLOCATE(icc3dstra(klon, klev))
       ALLOCATE(od443aer(klon))
       ALLOCATE(od550aer(klon))
@@ -845,4 +855,7 @@
       DEALLOCATE(d_u_lif,d_v_lif)
       DEALLOCATE(d_ts, d_tr)
+      DEALLOCATE(ZLWFT0_i,ZSWFT0_i,ZFLDN0)
+      DEALLOCATE(ZFLUP0,ZFSDN0,ZFSUP0)
+      DEALLOCATE(PSFSWDIR, PSFSWDIF)
       DEALLOCATE(topswad_aero,solswad_aero)
       DEALLOCATE(topswai_aero,solswai_aero)
@@ -868,4 +881,6 @@
       DEALLOCATE(lcc3dcon)
       DEALLOCATE(lcc3dstra)
+      DEALLOCATE(icc3dcon)
+      DEALLOCATE(icc3dstra)
       DEALLOCATE(od443aer)
       DEALLOCATE(od550aer)
Index: /LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90	(revision 3111)
+++ /LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90	(revision 3112)
@@ -1447,6 +1447,16 @@
   TYPE(ctrl_out), SAVE :: o_lcc3dstra = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lcc3dstra', 'Stratiform cloud liquid fraction', '1', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_icc3dcon = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'icc3dcon', 'Mass Fraction of Convective Cloud Ice', &
+    'kg kg-1', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_icc3dstra = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'icc3dstra', 'Mass Fraction of Stratiform Cloud Ice', &
+    'kg kg-1', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_solbnd = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'solbnd', 'Top-of-Atmosphere Solar Insolation for each band', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_rsdsdiff = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'rsdsdiff', 'Surface Diffuse Downwelling Shortwave Radiation', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_rsdscsdiff = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'rsdscsdiff', 'Surface Diffuse Downwelling Clear Sky Shortwave Radiation', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_stratomask = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'stratomask', 'Stratospheric fraction', '1', (/ ('', i=1, 10) /))
Index: /LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90	(revision 3111)
+++ /LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90	(revision 3112)
@@ -120,7 +120,8 @@
          o_swtoacf_zero, o_swsrfcf_zero, &
          o_topswai, o_solswai, o_toplwai, o_sollwai, o_scdnc, &
-         o_cldncl, o_reffclws, o_reffclwc, o_solbnd, o_stratomask,&
+         o_cldncl, o_reffclws, o_reffclwc, o_solbnd, &
+         o_rsdsdiff, o_rsdscsdiff, o_stratomask,&
          o_cldnvi, o_lcc, o_lcc3d, o_lcc3dcon, &
-         o_lcc3dstra, o_reffclwtop, o_ec550aer, &
+         o_lcc3dstra, o_icc3dcon, o_icc3dstra, o_reffclwtop, o_ec550aer, &
          o_lwcon, o_iwcon, o_temp, o_theta, &
          o_ovapinit, o_ovap, o_oliq, o_ocond, o_geop, &
@@ -277,6 +278,7 @@
          toplwad_aero, toplwad0_aero, sollwad_aero, &
          sollwad0_aero, toplwai_aero, sollwai_aero, &
-         scdnc, cldncl, reffclws, reffclwc, cldnvi, stratomask,&
-         lcc, lcc3d, lcc3dcon, lcc3dstra, reffclwtop, &
+         scdnc, cldncl, reffclws, reffclwc, cldnvi, &
+         PSFSWDIR, PSFSWDIF, stratomask,&
+         lcc, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra, reffclwtop, &
          ec550aer, flwc, fiwc, t_seri, theta, q_seri, &
          ql_seri, qs_seri, tr_seri, &
@@ -1406,4 +1408,6 @@
           CALL histwrite_phy(o_lcc3dcon, lcc3dcon)
           CALL histwrite_phy(o_lcc3dstra, lcc3dstra)
+          CALL histwrite_phy(o_icc3dcon, icc3dcon)
+          CALL histwrite_phy(o_icc3dstra, icc3dstra)
           CALL histwrite_phy(o_reffclwtop, reffclwtop)
        ENDIF
@@ -1432,4 +1436,5 @@
 #ifdef CPP_RRTM
       IF (iflag_rrtm.EQ.1) THEN
+
        IF (vars_defined) THEN
         DO ISW=1, NSW
@@ -1438,4 +1443,13 @@
         CALL histwrite_phy(o_solbnd, zx_tmp_fi3dsp)
        ENDIF
+
+       IF (vars_defined) THEN
+        CALL histwrite_phy(o_rsdscsdiff, PSFSWDIR)
+       ENDIF
+
+       IF (vars_defined) THEN
+        CALL histwrite_phy(o_rsdsdiff, PSFSWDIF)
+       ENDIF
+
       ENDIF
 #endif
Index: /LMDZ6/trunk/libf/phylmd/physiq_mod.F90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/physiq_mod.F90	(revision 3111)
+++ /LMDZ6/trunk/libf/phylmd/physiq_mod.F90	(revision 3112)
@@ -80,4 +80,5 @@
        ZLWFT0_i,ZSWFT0_i,ZFLDN0,  &
        ZFLUP0,ZFSDN0,ZFSUP0,      &
+       PSFSWDIR, PSFSWDIF, &
        !
        topswad_aero,solswad_aero,   &
@@ -3792,5 +3793,6 @@
                                 !-end
                ZLWFT0_i, ZFLDN0, ZFLUP0, &
-               ZSWFT0_i, ZFSDN0, ZFSUP0)
+               ZSWFT0_i, ZFSDN0, ZFSUP0, &
+               PSFSWDIR, PSFSWDIF)
 
           !lwoff=y, betalwoff=1. : offset LW CRE for radiation code and other
@@ -3876,5 +3878,6 @@
                                 !-end
                      ZLWFT0_i, ZFLDN0, ZFLUP0, &
-                     ZSWFT0_i, ZFSDN0, ZFSUP0)
+                     ZSWFT0_i, ZFSDN0, ZFSUP0, &
+                     PSFSWDIR, PSFSWDIF)
           endif !ok_4xCO2atm
        ENDIF ! aerosol_couple
Index: /LMDZ6/trunk/libf/phylmd/radlwsw_m.F90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/radlwsw_m.F90	(revision 3111)
+++ /LMDZ6/trunk/libf/phylmd/radlwsw_m.F90	(revision 3112)
@@ -42,5 +42,6 @@
 !-end
    ZLWFT0_i, ZFLDN0, ZFLUP0,&
-   ZSWFT0_i, ZFSDN0, ZFSUP0)
+   ZSWFT0_i, ZFSDN0, ZFSUP0,&
+   PSFSWDIR, PSFSWDIF)