Index: LMDZ6/trunk/DefLists/field_def_lmdz.xml
===================================================================
--- LMDZ6/trunk/DefLists/field_def_lmdz.xml	(revision 3109)
+++ LMDZ6/trunk/DefLists/field_def_lmdz.xml	(revision 3110)
@@ -24,4 +24,5 @@
         <field id="contfracATM"     long_name="% sfce ter+lic"  unit="-"    />
         <field id="contfracOR"      long_name="% sfce terre OR" unit="-" />
+        <field id="sza"     long_name="Solar zenithal angle "   unit="deg" />
         <field id="flat"    long_name="Latent heat flux"        unit="W/m2" />
         <field id="slp"     long_name="Sea Level Pressure"      unit="Pa" />
Index: LMDZ6/trunk/DefLists/file_def_histLES_lmdz.xml
===================================================================
--- LMDZ6/trunk/DefLists/file_def_histLES_lmdz.xml	(revision 3109)
+++ LMDZ6/trunk/DefLists/file_def_histLES_lmdz.xml	(revision 3110)
@@ -18,4 +18,5 @@
                 <field field_ref="flat" level="5" />
                 <field field_ref="slp" level="10" />
+                <field field_ref="sza" level="10" />
                 <field field_ref="ptstar" level="10" />
                 <field field_ref="pt0" level="10" />
Index: LMDZ6/trunk/DefLists/file_def_histday_lmdz.xml
===================================================================
--- LMDZ6/trunk/DefLists/file_def_histday_lmdz.xml	(revision 3109)
+++ LMDZ6/trunk/DefLists/file_def_histday_lmdz.xml	(revision 3110)
@@ -18,4 +18,5 @@
                 <field field_ref="flat" level="1" />
                 <field field_ref="slp" level="1" />
+                <field field_ref="sza" level="1" />
                 <field field_ref="ptstar" level="1" />
                 <field field_ref="pt0" level="1" />
Index: LMDZ6/trunk/DefLists/file_def_histhf_lmdz.xml
===================================================================
--- LMDZ6/trunk/DefLists/file_def_histhf_lmdz.xml	(revision 3109)
+++ LMDZ6/trunk/DefLists/file_def_histhf_lmdz.xml	(revision 3110)
@@ -18,4 +18,5 @@
                 <field field_ref="flat" level="5" />
                 <field field_ref="slp" level="10" />
+                <field field_ref="sza" level="10" />
                 <field field_ref="ptstar" level="10" />
                 <field field_ref="pt0" level="10" />
Index: LMDZ6/trunk/DefLists/file_def_histins_lmdz.xml
===================================================================
--- LMDZ6/trunk/DefLists/file_def_histins_lmdz.xml	(revision 3109)
+++ LMDZ6/trunk/DefLists/file_def_histins_lmdz.xml	(revision 3110)
@@ -18,4 +18,5 @@
                 <field field_ref="flat" level="10" />
                 <field field_ref="slp" level="1" />
+                <field field_ref="sza" level="1" />
                 <field field_ref="ptstar" level="1" />
                 <field field_ref="pt0" level="1" />
Index: LMDZ6/trunk/DefLists/file_def_histmth_lmdz.xml
===================================================================
--- LMDZ6/trunk/DefLists/file_def_histmth_lmdz.xml	(revision 3109)
+++ LMDZ6/trunk/DefLists/file_def_histmth_lmdz.xml	(revision 3110)
@@ -47,4 +47,5 @@
                 <field field_ref="flat"         level="5" />
                 <field field_ref="slp"          level="1" />
+                <field field_ref="sza"          level="1" />
                 <field field_ref="ptstar"       level="1" />
                 <field field_ref="pt0"          level="1" />
Index: LMDZ6/trunk/DefLists/file_def_histstn_lmdz.xml
===================================================================
--- LMDZ6/trunk/DefLists/file_def_histstn_lmdz.xml	(revision 3109)
+++ LMDZ6/trunk/DefLists/file_def_histstn_lmdz.xml	(revision 3110)
@@ -18,4 +18,5 @@
                 <field field_ref="flat" level="10" />
                 <field field_ref="slp" level="10" />
+                <field field_ref="sza" level="10" />
                 <field field_ref="ptstar" level="10" />
                 <field field_ref="pt0" level="10" />
Index: LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90	(revision 3109)
+++ LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90	(revision 3110)
@@ -51,4 +51,6 @@
 
 !!! 2D
+  TYPE(ctrl_out), SAVE :: o_sza = ctrl_out((/ 1, 1, 10, 10, 5, 10, 11, 11, 11, 11/), &
+    'sza', 'Solar zenithal angle', 'degrees', (/ ('', i=1, 10) /))
 
 ! Marine
Index: LMDZ6/trunk/libf/phylmd/phys_output_var_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/phys_output_var_mod.F90	(revision 3109)
+++ LMDZ6/trunk/libf/phylmd/phys_output_var_mod.F90	(revision 3110)
@@ -14,4 +14,6 @@
   REAL, SAVE, ALLOCATABLE :: snow_o(:), zfra_o(:)
   !$OMP THREADPRIVATE(snow_o, zfra_o)
+  REAL, SAVE, ALLOCATABLE :: sza_o(:) ! solar zenithal angle
+  !$OMP THREADPRIVATE(sza_o)
   INTEGER, SAVE, ALLOCATABLE ::  itau_con(:)       ! Nombre de pas ou rflag <= 1
   !$OMP THREADPRIVATE(itau_con)
@@ -144,4 +146,5 @@
 
     allocate(snow_o(klon), zfra_o(klon))
+    allocate(sza_o(klon) )
     allocate(itau_con(klon))
     allocate(sens_prec_liq_o(klon,2))
@@ -198,4 +201,5 @@
 
     deallocate(snow_o,zfra_o,itau_con)
+    deallocate(sza_o)
     deallocate (bils_ec,bils_ech,bils_tke,bils_diss,bils_kinetic,bils_enthalp,bils_latent)
     deallocate (d_qw_col, d_ql_col, d_qs_col, d_qt_col, d_ek_col, d_h_dair_col &
Index: LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90	(revision 3109)
+++ LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90	(revision 3110)
@@ -178,4 +178,5 @@
          o_lat_prec_liq_oce, o_lat_prec_liq_sic, &
          o_lat_prec_sol_oce, o_lat_prec_sol_sic, &
+         o_sza, &
 ! Marine
          o_map_prop_hc, o_map_prop_hist, o_map_emis_hc, o_map_iwp_hc, &
@@ -319,5 +320,5 @@
          itau_con, nfiles, clef_files, nid_files, dryaod_diag, &
          zustr_gwd_hines, zvstr_gwd_hines,zustr_gwd_rando, zvstr_gwd_rando, &
-         zustr_gwd_front, zvstr_gwd_front,     &
+         zustr_gwd_front, zvstr_gwd_front, sza_o,    &
          sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o, &
          cloudth_sth,cloudth_senv,cloudth_sigmath,cloudth_sigmaenv, &
@@ -614,4 +615,5 @@
      ENDIF
 
+       CALL histwrite_phy(o_sza, sza_o)
        CALL histwrite_phy(o_flat, zxfluxlat)
        CALL histwrite_phy(o_ptstar, ptstar)
Index: LMDZ6/trunk/libf/phylmd/physiq_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/physiq_mod.F90	(revision 3109)
+++ LMDZ6/trunk/libf/phylmd/physiq_mod.F90	(revision 3110)
@@ -1122,4 +1122,8 @@
     logical, parameter :: mass_fixer=.false.
     real qql1(klon),qql2(klon),corrqql
+
+    REAL pi
+
+    pi = 4. * ATAN(1.)
 
     ! Ehouarn: set value of jjmp1 since it is no longer a "fixed parameter"
@@ -2190,4 +2194,5 @@
        END SELECT
     ENDIF
+    sza_o = ACOS (rmu0) *180./pi
 
     IF (mydebug) THEN