Index: CONFIG/IPSLCM/IPSLCM5A/Chistorical/COMP/lmdz.card
===================================================================
--- CONFIG/IPSLCM/IPSLCM5A/Chistorical/COMP/lmdz.card	(revision 1506)
+++ CONFIG/IPSLCM/IPSLCM5A/Chistorical/COMP/lmdz.card	(revision 1507)
@@ -74,17 +74,17 @@
 
 [OutputFiles]
-List=	(histmth.nc,      ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,      Post_1M_histmth),     \
-	(histmthCOSP.nc,  ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthCOSP.nc,  Post_1M_histmthCOSP), \
-	(histmthNMC.nc,   ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthNMC.nc,   Post_1M_histmthNMC),  \
-	(dynzon.nc,       ${R_OUT_ATM_O_M}/${PREFIX}_1M_dynzon.nc,       Post_1M_dynzon),      \
-	(histday.nc,      ${R_OUT_ATM_O_D}/${PREFIX}_1D_histday.nc,      Post_1D_histday),     \
-	(histdayCOSP.nc,  ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayCOSP.nc,  Post_1D_histdayCOSP), \
-	(histdayNMC.nc,   ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayNMC.nc,   Post_1D_histdayNMC),  \
-	(histhf.nc,       ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhf.nc,       Post_HF_histhf),      \
-	(histhfCOSP.nc,   ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfCOSP.nc,   Post_HF_histhfCOSP),  \
-	(histhfNMC.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfNMC.nc,    Post_HF_histhfNMC),   \
-	(histhf3h.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3h.nc,     Post_3H_histhf3h),    \
-        (histhf3hm.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3hm.nc,    Post_3H_histhf3hm),   \
-        (paramLMDZ_phy.nc, ${R_OUT_ATM_O_M}/${PREFIX}_1M_paramLMDZ_phy.nc, Post_1M_paramLMDZ_phy)
+List=	(histmth.nc,       ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,       Post_1M_histmth),       \
+	(histmthCOSP.nc,   ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthCOSP.nc,   Post_1M_histmthCOSP),   \
+	(histmthNMC.nc,    ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthNMC.nc,    Post_1M_histmthNMC),    \
+	(dynzon.nc,        ${R_OUT_ATM_O_M}/${PREFIX}_1M_dynzon.nc,        Post_1M_dynzon),        \
+        (paramLMDZ_phy.nc, ${R_OUT_ATM_O_M}/${PREFIX}_1M_paramLMDZ_phy.nc, Post_1M_paramLMDZ_phy), \
+	(histday.nc,       ${R_OUT_ATM_O_D}/${PREFIX}_1D_histday.nc,       Post_1D_histday),       \
+	(histdayCOSP.nc,   ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayCOSP.nc,   Post_1D_histdayCOSP),   \
+	(histdayNMC.nc,    ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayNMC.nc,    Post_1D_histdayNMC),    \
+	(histhf.nc,        ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhf.nc,        Post_HF_histhf),        \
+	(histhfCOSP.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfCOSP.nc,    Post_HF_histhfCOSP),    \
+	(histhfNMC.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfNMC.nc,     Post_HF_histhfNMC),     \
+	(histhf3h.nc,      ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3h.nc,      Post_3H_histhf3h),      \
+        (histhf3hm.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3hm.nc,     Post_3H_histhf3hm)
 
 [Post_1M_histmth]
@@ -93,5 +93,5 @@
 TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, flat, ftime_con, iwp, lwp, LWdnSFC, LWdnSFCclr, LWupSFCclr, LWdn200, LWdn200clr, LWup200, LWup200clr, precip, prw, pluc, psol, q2m, rh2m, sens, sicf, slp, snow, soll0, soll, sols0, sols, SWdnSFC, SWupSFC, SWup200, SWdnTOA, SWupTOA, SWdn200, SWdnSFCclr, SWupSFCclr, SWupTOAclr, t2m, t2m_max, t2m_min, taux, taux_ter, taux_sic, taux_lic, taux_oce, tauy, tauy_ter, tauy_sic, tauy_lic, tauy_oce, lat_ter, lat_sic, lat_lic, lat_oce, nettop, pbase, ptop, topl0, topl, tops0, tops, topswad, topswai, tsol, tsol_ter, tsol_lic, tsol_oce, tsol_sic , u10m, u500, v10m, v500, w500, pourc_ter, pourc_oce, pourc_sic, pourc_lic, wind10m, od550aer, od550lt1aer, absvisaer, loadoa, loadbc, loadso4, loaddust, loadss, sconcoa, sconcbc, sconcso4, sconcdust, sconcss, reffclwtop, cldncl, swtoaas_ant, swtoacs_ant, swtoacf_ant, swsrfas_ant, swsrfcs_ant, swsrfcf_ant, swtoaas_nat, swtoacs_nat, swtoacf_nat, swsrfas_nat, swsrfcs_nat, swsrfcf_nat, solswai, swtoacf_zero, swsrfcf_zero)
 ChunckJob2D = NONE
-TimeSeriesVars3D = (upwd, lwcon, iwcon, lcc3dcon, lcc3dstra, temp, theta, ovap, ovapinit, geop, vitu, vitv, vitw, pres, rneb, rnebcon, rhum, mc, paprs, scdnc, ec550aer, concoa, concbc, concso4, concss, concdust, reffclws, reffclwc, Aga)
+TimeSeriesVars3D = (upwd, lwcon, iwcon, lcc3dcon, lcc3dstra, temp, theta, ovap, ovapinit, geop, vitu, vitv, vitw, pres, rneb, rnebcon, rhum, mc, paprs, scdnc, ec550aer, concoa, concbc, concso4, concss, concdust, reffclws, reffclwc, Aga, rlu, rsu, rld, rsd, rlucs, rsucs, rldcs, rsdcs, tnt, dtdyn, dtphy, dtlschr, tntr, tntc, tntscpbl, tnhus, dqdyn, tnhusc, dqvdf, tnhusscpbl, dqphy, evu, kz, mcd)
 ChunckJob3D = 100Y
 Seasonal=ON
@@ -131,5 +131,5 @@
 TimeSeriesVars3D = ()
 ChunckJob3D = NONE
-Seasonal=ON
+Seasonal=OFF
 
 [Post_1D_histday]
@@ -190,5 +190,5 @@
 Patches= ()
 GatherWithInternal = (lon, lat, presnivs, time_counter)
-TimeSeriesVars2D = (phis, tsol, t2m, q2m, u10m, v10m, psol, topl, bils)
+TimeSeriesVars2D = (phis, tsol, t2m, q2m, u10m, v10m, psol, topl, bils, t100, u100, v100, w100, t50, u50, v50, w50, t10, u10, v10, w10)
 ChunckJob2D = 50Y
 TimeSeriesVars3D = ()
Index: CONFIG/IPSLCM/IPSLCM5A/Chistorical/PARAM/physiq.def_L39
===================================================================
--- CONFIG/IPSLCM/IPSLCM5A/Chistorical/PARAM/physiq.def_L39	(revision 1506)
+++ CONFIG/IPSLCM/IPSLCM5A/Chistorical/PARAM/physiq.def_L39	(revision 1507)
@@ -29,8 +29,4 @@
 # Parametres fichiers de sortie
 #
-### type_run = type run par rapport aux fichiers et variables de sortie
-# - type_run = CLIM/ENSP (=1)
-# - type_run = AMIP/CFMI (=2)
-type_run=CLIM
 ### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc
 OK_journe=n
@@ -42,25 +38,86 @@
 ok_hf=n
 # Noms des fichiers
-phys_out_filenames=      histmth histday histhf  histhf3h histhf3hm
+phys_out_filenames=      histmth histday histhf histhf3h histhf3hm histstn
 # Sortir ou non les fichiers
-phys_out_filekeys=       y       y       y       y       y
+phys_out_filekeys=       y        y        y        y        y        n
 # Niveaux de sorties
-phys_out_filelevels=     5       5       5       5       5
-# more variables for level 5 : vitw and pres
+phys_out_filelevels=     5        5        5        5        5        5
+### Type de fichier : global (n) ou stations (y)
+phys_out_filestations =  n        n        n        n        n        y
+### Frequences des sorties
+phys_out_filetimesteps = 1.mth, 1.day, 0.25day, 0.125day, 0.125day, 1800.s
+### variables niveaux pressions 100, 50, 10 hPa a 3h-inst - F.Lott
+flag_t100       = 1, 7, 7, 5, 10, 10
+flag_u100       = 1, 7, 7, 5, 10, 10
+flag_v100       = 1, 7, 7, 5, 10, 10
+flag_w100       = 1, 7, 7, 5, 10, 10
+flag_t50        = 1, 7, 7, 5, 10, 10
+flag_u50        = 1, 7, 7, 5, 10, 10
+flag_v50        = 1, 7, 7, 5, 10, 10
+flag_w50        = 1, 7, 7, 5, 10, 10
+flag_t10        = 1, 7, 7, 5, 10, 10
+flag_u10        = 1, 7, 7, 5, 10, 10
+flag_v10        = 1, 7, 7, 5, 10, 10
+flag_w10        = 1, 7, 7, 5, 10, 10
+# more variables for level 5 histday/month: vitw and pres
 # Vertical wind
-flag_vitw         =  2, 3, 7, 6, 6
-name_vitw         =  vitw
+flag_vitw         =  2, 3, 7, 6, 6, 6
 # Air pressure
-flag_pres         =  2, 3, 7, 6, 6
-name_pres         =  pres
+flag_pres         =  2, 3, 7, 6, 6, 6
 # w500 for histday
-flag_w500         =  1, 5, 3, 10, 10
-name_w500         =  w500
+flag_w500         =  1, 5, 3, 10, 10, 10
 # t700 for histday
-flag_t700         =  1, 5, 3, 10, 10
-name_t700         =  t700
+flag_t700         =  1, 5, 3, 10, 10, 10
 # Cloud fraction
-flag_rneb         =  2, 5, 10, 5, 10
-name_rneb         =  rneb
+flag_rneb         =  2, 5, 10, 5, 10, 10
+# Aga sortie avec niveau 2
+flag_Aga          =  2, 10, 10, 10, 10, 10
+# niveaux maximum 5 mensuel pour variables CFMIP
+flag_rlu          =  5, 10, 10, 10, 10, 10
+flag_rsu          =  5, 10, 10, 10, 10, 10
+flag_rld          =  5, 10, 10, 10, 10, 10
+flag_rsd          =  5, 10, 10, 10, 10, 10
+flag_rlucs        =  5, 10, 10, 10, 10, 10
+flag_rsucs        =  5, 10, 10, 10, 10, 10
+flag_rldcs        =  5, 10, 10, 10, 10, 10
+flag_rsdcs        =  5, 10, 10, 10, 10, 10
+flag_tnt          =  5, 10, 10, 10, 10, 10
+flag_dtdyn        =  4, 10, 10, 10, 10, 10
+flag_dtphy        =  4, 10, 10, 10, 10, 10
+flag_dtlschr      =  5, 10, 10, 10, 10, 10
+flag_tntr         =  5, 10, 10, 10, 10, 10
+flag_tntc         =  5, 10, 10, 10, 10, 10
+flag_tntscpbl     =  5, 10, 10, 10, 10, 10
+flag_tnhus        =  5, 10, 10, 10, 10, 10
+flag_dqdyn        =  4, 10, 10, 10, 10, 10
+flag_tnhusc       =  5, 10, 10, 10, 10, 10
+flag_dqvdf        =  5, 10, 10, 10, 10, 10
+flag_tnhusscpbl   =  5, 10, 10, 10, 10, 10
+flag_dqphy        =  2, 10, 10, 10, 10, 10
+flag_evu          =  5, 10, 10, 10, 10, 10
+flag_kz           =  4, 10, 10, 10, 10, 10
+###
+# more variables for t tendency monthly
+###
+flag_dtvdf        =  5, 10, 10, 10, 10, 10
+flag_dteva        =  5, 10, 10, 10, 10, 10
+flag_dtcon        =  4, 10, 10, 10, 10, 10
+flag_dtajs        =  4, 10, 10, 10, 10, 10
+flag_dtlsc        =  4, 10, 10, 10, 10, 10
+flag_dtoro        =  4, 10, 10, 10, 10, 10
+flag_dtlif        =  4, 10, 10, 10, 10, 10
+flag_dthin        =  4, 10, 10, 10, 10, 10
+flag_dtec         =  4, 10, 10, 10, 10, 10
+### more variables for q tendency monthly
+flag_dqeva        =  5, 10, 10, 10, 10, 10
+flag_dqcon        =  4, 10, 10, 10, 10, 10
+flag_dqajs        =  4, 10, 10, 10, 10, 10
+flag_dqlsc        =  5, 10, 10, 10, 10, 10
+### nouvelle physique monthly
+flag_dtthe        =  4, 10, 10, 10, 10, 10
+flag_dtwak        =  4, 5, 10, 10, 10, 10
+flag_dqthe        =  4, 10, 10, 10, 10, 10
+flag_dqwak        =  4, 5, 10, 10, 10, 10
+flag_ftime_th     =  4, 10, 10, 5, 10, 10
 #
 # Parametres niveau de sorties differents fichiers
@@ -85,6 +142,4 @@
 # - lev_histmth=4 => champs tendances 3d
 lev_histmth=5
-### Aga sortie avec niveau 2
-flag_Aga       =  2, 10, 10, 10, 10
 ### parametres pour hist*NMC.nc
 #lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier
@@ -93,10 +148,4 @@
 ### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence
 ok_histNMC = y, n, n
-### ecrit_LES = frequence ecriture fichier histhf3hm.nc en jours
-ecrit_LES=0.125
-### ecrit_ins = frequence ecriture fichier histhf3h.nc en jours
-ecrit_ins=0.125
-### ecrit_hf = frequence ecriture fichier histhf.nc en jours
-ecrit_hf=0.250
 ### ecrit_day = frequence ecriture fichier histday.nc en jours
 ecrit_day=1.
Index: CONFIG/IPSLCM/IPSLCM5A/CpiControl/COMP/lmdz.card
===================================================================
--- CONFIG/IPSLCM/IPSLCM5A/CpiControl/COMP/lmdz.card	(revision 1506)
+++ CONFIG/IPSLCM/IPSLCM5A/CpiControl/COMP/lmdz.card	(revision 1507)
@@ -69,17 +69,17 @@
 
 [OutputFiles]
-List=	(histmth.nc,      ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,      Post_1M_histmth),     \
-	(histmthCOSP.nc,  ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthCOSP.nc,  Post_1M_histmthCOSP), \
-	(histmthNMC.nc,   ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthNMC.nc,   Post_1M_histmthNMC),  \
-	(dynzon.nc,       ${R_OUT_ATM_O_M}/${PREFIX}_1M_dynzon.nc,       Post_1M_dynzon),      \
-	(histday.nc,      ${R_OUT_ATM_O_D}/${PREFIX}_1D_histday.nc,      Post_1D_histday),     \
-	(histdayCOSP.nc,  ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayCOSP.nc,  Post_1D_histdayCOSP), \
-	(histdayNMC.nc,   ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayNMC.nc,   Post_1D_histdayNMC),  \
-	(histhf.nc,       ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhf.nc,       Post_HF_histhf),      \
-	(histhfCOSP.nc,   ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfCOSP.nc,   Post_HF_histhfCOSP),  \
-	(histhfNMC.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfNMC.nc,    Post_HF_histhfNMC),   \
-	(histhf3h.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3h.nc,     Post_3H_histhf3h),    \
-	(histhf3hm.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3hm.nc,    Post_3H_histhf3hm),   \
-        (paramLMDZ_phy.nc, ${R_OUT_ATM_O_M}/${PREFIX}_1M_paramLMDZ_phy.nc, Post_1M_paramLMDZ_phy)
+List=	(histmth.nc,       ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,       Post_1M_histmth),       \
+	(histmthCOSP.nc,   ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthCOSP.nc,   Post_1M_histmthCOSP),   \
+	(histmthNMC.nc,    ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthNMC.nc,    Post_1M_histmthNMC),    \
+	(dynzon.nc,        ${R_OUT_ATM_O_M}/${PREFIX}_1M_dynzon.nc,        Post_1M_dynzon),        \
+        (paramLMDZ_phy.nc, ${R_OUT_ATM_O_M}/${PREFIX}_1M_paramLMDZ_phy.nc, Post_1M_paramLMDZ_phy), \
+	(histday.nc,       ${R_OUT_ATM_O_D}/${PREFIX}_1D_histday.nc,       Post_1D_histday),       \
+	(histdayCOSP.nc,   ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayCOSP.nc,   Post_1D_histdayCOSP),   \
+	(histdayNMC.nc,    ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayNMC.nc,    Post_1D_histdayNMC),    \
+	(histhf.nc,        ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhf.nc,        Post_HF_histhf),        \
+	(histhfCOSP.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfCOSP.nc,    Post_HF_histhfCOSP),    \
+	(histhfNMC.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfNMC.nc,     Post_HF_histhfNMC),     \
+	(histhf3h.nc,      ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3h.nc,      Post_3H_histhf3h),      \
+        (histhf3hm.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3hm.nc,     Post_3H_histhf3hm)
 
 [Post_1M_histmth]
@@ -88,5 +88,5 @@
 TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, flat, ftime_con, iwp, lwp, LWdnSFC, LWdnSFCclr, LWupSFCclr, LWdn200, LWdn200clr, LWup200, LWup200clr, precip, prw, pluc, psol, q2m, rh2m, sens, sicf, slp, snow, soll0, soll, sols0, sols, SWdnSFC, SWupSFC, SWup200, SWdnTOA, SWupTOA, SWdn200, SWdnSFCclr, SWupSFCclr, SWupTOAclr, t2m, t2m_max, t2m_min, taux, taux_ter, taux_sic, taux_lic, taux_oce, tauy, tauy_ter, tauy_sic, tauy_lic, tauy_oce, lat_ter, lat_sic, lat_lic, lat_oce, nettop, pbase, ptop, topl0, topl, tops0, tops, topswad, topswai, tsol, tsol_ter, tsol_lic, tsol_oce, tsol_sic , u10m, u500, v10m, v500, w500, pourc_ter, pourc_oce, pourc_sic, pourc_lic, wind10m, od550aer, od550lt1aer, absvisaer, loadoa, loadbc, loadso4, loaddust, loadss, sconcoa, sconcbc, sconcso4, sconcdust, sconcss, reffclwtop, cldncl, swtoaas_ant, swtoacs_ant, swtoacf_ant, swsrfas_ant, swsrfcs_ant, swsrfcf_ant, swtoaas_nat, swtoacs_nat, swtoacf_nat, swsrfas_nat, swsrfcs_nat, swsrfcf_nat, solswai, swtoacf_zero, swsrfcf_zero)
 ChunckJob2D = NONE
-TimeSeriesVars3D = (upwd, lwcon, iwcon, lcc3dcon, lcc3dstra, temp, theta, ovap, ovapinit, geop, vitu, vitv, vitw, pres, rneb, rnebcon, rhum, mc, paprs, scdnc, ec550aer, concoa, concbc, concso4, concss, concdust, reffclws, reffclwc, Aga)
+TimeSeriesVars3D = (upwd, lwcon, iwcon, lcc3dcon, lcc3dstra, temp, theta, ovap, ovapinit, geop, vitu, vitv, vitw, pres, rneb, rnebcon, rhum, mc, paprs, scdnc, ec550aer, concoa, concbc, concso4, concss, concdust, reffclws, reffclwc, Aga, rlu, rsu, rld, rsd, rlucs, rsucs, rldcs, rsdcs, tnt, dtdyn, dtphy, dtlschr, tntr, tntc, tntscpbl, tnhus, dqdyn, tnhusc, dqvdf, tnhusscpbl, dqphy, evu, kz, mcd)
 ChunckJob3D = 100Y
 Seasonal=ON
@@ -126,11 +126,11 @@
 TimeSeriesVars3D = ()
 ChunckJob3D = NONE
-Seasonal=ON
+Seasonal=OFF
 
 [Post_1D_histday]
 Patches= ()
 GatherWithInternal = (lon, lat, presnivs, time_counter, time_counter_bnds)
-TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, flat, precip, pluc, plul, prw, psol, q2m, sens, sicf, slp, soll, sols, snow, SWdnSFC, t2m, t2m_max, t2m_min, topl, tops, tsol, t700, u10m, v10m, wind10m, wind10max, w500)
-ChunckJob2D = 200Y
+TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, flat, iwp, lwp, pbase, ptop, precip, pluc, plul, prw, psol, rh2m, rh2m_min, rh2m_max, q2m, sens, sicf, slp, soll, sols, snow, LWdnSFC, LWdnSFCclr, SWdnSFC, SWdnSFCclr, SWdnTOA, LWupSFCclr, SWupSFC, SWupSFCclr, SWupTOA, SWupTOAclr, t2m, t2m_max, t2m_min, topl, topl0, tops, tsol, t700, u10m, v10m, wind10m, wind10max, w500)
+ChunckJob2D = 50Y
 TimeSeriesVars3D = (lwcon, iwcon, mc, rhum, rneb, temp, theta, ovap, geop, vitu, vitv, vitw, pres, paprs)
 ChunckJob3D = OFF
@@ -185,5 +185,5 @@
 Patches= ()
 GatherWithInternal = (lon, lat, presnivs, time_counter)
-TimeSeriesVars2D = (phis, tsol, t2m, q2m, u10m, v10m, psol, topl, bils)
+TimeSeriesVars2D = (phis, tsol, t2m, q2m, u10m, v10m, psol, topl, bils, t100, u100, v100, w100, t50, u50, v50, w50, t10, u10, v10, w10)
 ChunckJob2D = OFF
 TimeSeriesVars3D = ()
Index: CONFIG/IPSLCM/IPSLCM5A/CpiControl/PARAM/physiq.def_L39
===================================================================
--- CONFIG/IPSLCM/IPSLCM5A/CpiControl/PARAM/physiq.def_L39	(revision 1506)
+++ CONFIG/IPSLCM/IPSLCM5A/CpiControl/PARAM/physiq.def_L39	(revision 1507)
@@ -29,8 +29,4 @@
 # Parametres fichiers de sortie
 #
-### type_run = type run par rapport aux fichiers et variables de sortie
-# - type_run = CLIM/ENSP (=1)
-# - type_run = AMIP/CFMI (=2)
-type_run=CLIM
 ### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc
 OK_journe=n
@@ -41,22 +37,79 @@
 ### OK_hf=y, ecrire sorties hautes frequence histhf.nc, =n pas de fichier histhf.nc
 ok_hf=n
-#
-# more variables for level 5 : vitw and pres
+### variables niveaux pressions 100, 50, 10 hPa a 3h-inst - F.Lott
+flag_t100       = 1, 7, 7, 5, 10, 10
+flag_u100       = 1, 7, 7, 5, 10, 10
+flag_v100       = 1, 7, 7, 5, 10, 10
+flag_w100       = 1, 7, 7, 5, 10, 10
+flag_t50        = 1, 7, 7, 5, 10, 10
+flag_u50        = 1, 7, 7, 5, 10, 10
+flag_v50        = 1, 7, 7, 5, 10, 10
+flag_w50        = 1, 7, 7, 5, 10, 10
+flag_t10        = 1, 7, 7, 5, 10, 10
+flag_u10        = 1, 7, 7, 5, 10, 10
+flag_v10        = 1, 7, 7, 5, 10, 10
+flag_w10        = 1, 7, 7, 5, 10, 10
+# more variables for level 5 histday/month: vitw and pres
 # Vertical wind
-flag_vitw         =  2, 3, 7, 6, 6
-name_vitw         =  vitw
+flag_vitw         =  2, 3, 7, 6, 6, 6
 # Air pressure
-flag_pres         =  2, 3, 7, 6, 6
-name_pres         =  pres
+flag_pres         =  2, 3, 7, 6, 6, 6
 # w500 for histday
-flag_w500         =  1, 5, 3, 10, 10
-name_w500         =  w500
+flag_w500         =  1, 5, 3, 10, 10, 10
 # t700 for histday
-flag_t700         =  1, 5, 3, 10, 10
-name_t700         =  t700
+flag_t700         =  1, 5, 3, 10, 10, 10
 # Cloud fraction
-flag_rneb         =  2, 5, 10, 5, 10
-name_rneb         =  rneb
-# Parametres niveau de sorties differents fichiers 
+flag_rneb         =  2, 5, 10, 5, 10, 10
+# Aga sortie avec niveau 2
+flag_Aga          =  2, 10, 10, 10, 10, 10
+# niveaux maximum 5 mensuel pour variables CFMIP
+flag_rlu          =  5, 10, 10, 10, 10, 10
+flag_rsu          =  5, 10, 10, 10, 10, 10
+flag_rld          =  5, 10, 10, 10, 10, 10
+flag_rsd          =  5, 10, 10, 10, 10, 10
+flag_rlucs        =  5, 10, 10, 10, 10, 10
+flag_rsucs        =  5, 10, 10, 10, 10, 10
+flag_rldcs        =  5, 10, 10, 10, 10, 10
+flag_rsdcs        =  5, 10, 10, 10, 10, 10
+flag_tnt          =  5, 10, 10, 10, 10, 10
+flag_dtdyn        =  4, 10, 10, 10, 10, 10
+flag_dtphy        =  4, 10, 10, 10, 10, 10
+flag_dtlschr      =  5, 10, 10, 10, 10, 10
+flag_tntr         =  5, 10, 10, 10, 10, 10
+flag_tntc         =  5, 10, 10, 10, 10, 10
+flag_tntscpbl     =  5, 10, 10, 10, 10, 10
+flag_tnhus        =  5, 10, 10, 10, 10, 10
+flag_dqdyn        =  4, 10, 10, 10, 10, 10
+flag_tnhusc       =  5, 10, 10, 10, 10, 10
+flag_dqvdf        =  5, 10, 10, 10, 10, 10
+flag_tnhusscpbl   =  5, 10, 10, 10, 10, 10
+flag_dqphy        =  2, 10, 10, 10, 10, 10
+flag_evu          =  5, 10, 10, 10, 10, 10
+flag_kz           =  4, 10, 10, 10, 10, 10
+###
+# more variables for t tendency monthly
+###
+flag_dtvdf        =  5, 10, 10, 10, 10, 10
+flag_dteva        =  5, 10, 10, 10, 10, 10
+flag_dtcon        =  4, 10, 10, 10, 10, 10
+flag_dtajs        =  4, 10, 10, 10, 10, 10
+flag_dtlsc        =  4, 10, 10, 10, 10, 10
+flag_dtoro        =  4, 10, 10, 10, 10, 10
+flag_dtlif        =  4, 10, 10, 10, 10, 10
+flag_dthin        =  4, 10, 10, 10, 10, 10
+flag_dtec         =  4, 10, 10, 10, 10, 10
+### more variables for q tendency monthly
+flag_dqeva        =  5, 10, 10, 10, 10, 10
+flag_dqcon        =  4, 10, 10, 10, 10, 10
+flag_dqajs        =  4, 10, 10, 10, 10, 10
+flag_dqlsc        =  5, 10, 10, 10, 10, 10
+### nouvelle physique monthly
+flag_dtthe        =  4, 10, 10, 10, 10, 10
+flag_dtwak        =  4, 5, 10, 10, 10, 10
+flag_dqthe        =  4, 10, 10, 10, 10, 10
+flag_dqwak        =  4, 5, 10, 10, 10, 10
+flag_ftime_th     =  4, 10, 10, 5, 10, 10
+#
+# Parametres niveau de sorties differents fichiers
 #
 ### lev_histhf=0-4, niveau de sortie fichier "histhf.nc" 
@@ -79,6 +132,4 @@
 # - lev_histmth=4 => champs tendances 3d
 lev_histmth=5
-### Aga sortie avec niveau 2
-flag_Aga       =  2, 10, 10, 10, 10
 ### parametres pour hist*NMC.nc
 #lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier
@@ -87,8 +138,4 @@
 ### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence
 ok_histNMC = y, n, n
-### ecrit_ins = frequence ecriture fichier histhf3h.nc en jours
-ecrit_ins=0.125
-### ecrit_hf = frequence ecriture fichier histhf.nc en jours
-ecrit_hf=0.250
 ### ecrit_day = frequence ecriture fichier histday.nc en jours
 ecrit_day=1.
Index: CONFIG/IPSLCM/IPSLCM5A/EXP00/COMP/lmdz.card
===================================================================
--- CONFIG/IPSLCM/IPSLCM5A/EXP00/COMP/lmdz.card	(revision 1506)
+++ CONFIG/IPSLCM/IPSLCM5A/EXP00/COMP/lmdz.card	(revision 1507)
@@ -67,17 +67,17 @@
 
 [OutputFiles]
-List=	(histmth.nc,      ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,      Post_1M_histmth),     \
-	(histmthCOSP.nc,  ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthCOSP.nc,  Post_1M_histmthCOSP), \
-	(histmthNMC.nc,   ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthNMC.nc,   Post_1M_histmthNMC),  \
-	(dynzon.nc,       ${R_OUT_ATM_O_M}/${PREFIX}_1M_dynzon.nc,       Post_1M_dynzon),      \
-	(histday.nc,      ${R_OUT_ATM_O_D}/${PREFIX}_1D_histday.nc,      Post_1D_histday),     \
-	(histdayCOSP.nc,  ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayCOSP.nc,  Post_1D_histdayCOSP), \
-	(histdayNMC.nc,   ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayNMC.nc,   Post_1D_histdayNMC),  \
-	(histhf.nc,       ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhf.nc,       Post_HF_histhf),      \
-	(histhfCOSP.nc,   ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfCOSP.nc,   Post_HF_histhfCOSP),  \
-	(histhfNMC.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfNMC.nc,    Post_HF_histhfNMC),   \
-	(histhf3h.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3h.nc,     Post_3H_histhf3h),    \
-	(histhf3hm.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3hm.nc,    Post_3H_histhf3hm),   \
-        (paramLMDZ_phy.nc, ${R_OUT_ATM_O_M}/${PREFIX}_1M_paramLMDZ_phy.nc, Post_1M_paramLMDZ_phy)
+List=	(histmth.nc,       ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,       Post_1M_histmth),       \
+	(histmthCOSP.nc,   ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthCOSP.nc,   Post_1M_histmthCOSP),   \
+	(histmthNMC.nc,    ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthNMC.nc,    Post_1M_histmthNMC),    \
+	(dynzon.nc,        ${R_OUT_ATM_O_M}/${PREFIX}_1M_dynzon.nc,        Post_1M_dynzon),        \
+        (paramLMDZ_phy.nc, ${R_OUT_ATM_O_M}/${PREFIX}_1M_paramLMDZ_phy.nc, Post_1M_paramLMDZ_phy), \
+	(histday.nc,       ${R_OUT_ATM_O_D}/${PREFIX}_1D_histday.nc,       Post_1D_histday),       \
+	(histdayCOSP.nc,   ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayCOSP.nc,   Post_1D_histdayCOSP),   \
+	(histdayNMC.nc,    ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayNMC.nc,    Post_1D_histdayNMC),    \
+	(histhf.nc,        ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhf.nc,        Post_HF_histhf),        \
+	(histhfCOSP.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfCOSP.nc,    Post_HF_histhfCOSP),    \
+	(histhfNMC.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfNMC.nc,     Post_HF_histhfNMC),     \
+	(histhf3h.nc,      ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3h.nc,      Post_3H_histhf3h),      \
+        (histhf3hm.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3hm.nc,     Post_3H_histhf3hm)
 
 [Post_1M_histmth]
@@ -86,5 +86,5 @@
 TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, flat, ftime_con, iwp, lwp, LWdnSFC, LWdnSFCclr, LWupSFCclr, LWdn200, LWdn200clr, LWup200, LWup200clr, precip, prw, pluc, psol, q2m, rh2m, sens, sicf, slp, snow, soll0, soll, sols0, sols, SWdnSFC, SWupSFC, SWup200, SWdnTOA, SWupTOA, SWdn200, SWdnSFCclr, SWupSFCclr, SWupTOAclr, t2m, t2m_max, t2m_min, taux, taux_ter, taux_sic, taux_lic, taux_oce, tauy, tauy_ter, tauy_sic, tauy_lic, tauy_oce, lat_ter, lat_sic, lat_lic, lat_oce, nettop, pbase, ptop, topl0, topl, tops0, tops, topswad, topswai, tsol, tsol_ter, tsol_lic, tsol_oce, tsol_sic , u10m, u500, v10m, v500, w500, pourc_ter, pourc_oce, pourc_sic, pourc_lic, wind10m, od550aer, od550lt1aer, absvisaer, loadoa, loadbc, loadso4, loaddust, loadss, sconcoa, sconcbc, sconcso4, sconcdust, sconcss, reffclwtop, cldncl, swtoaas_ant, swtoacs_ant, swtoacf_ant, swsrfas_ant, swsrfcs_ant, swsrfcf_ant, swtoaas_nat, swtoacs_nat, swtoacf_nat, swsrfas_nat, swsrfcs_nat, swsrfcf_nat, solswai, swtoacf_zero, swsrfcf_zero)
 ChunckJob2D = NONE
-TimeSeriesVars3D = (upwd, lwcon, iwcon, lcc3dcon, lcc3dstra, temp, theta, ovap, ovapinit, geop, vitu, vitv, vitw, pres, rneb, rnebcon, rhum, mc, paprs, scdnc, ec550aer, concoa, concbc, concso4, concss, concdust, reffclws, reffclwc, Aga)
+TimeSeriesVars3D = (upwd, lwcon, iwcon, lcc3dcon, lcc3dstra, temp, theta, ovap, ovapinit, geop, vitu, vitv, vitw, pres, rneb, rnebcon, rhum, mc, paprs, scdnc, ec550aer, concoa, concbc, concso4, concss, concdust, reffclws, reffclwc, Aga, rlu, rsu, rld, rsd, rlucs, rsucs, rldcs, rsdcs, tnt, dtdyn, dtphy, dtlschr, tntr, tntc, tntscpbl, tnhus, dqdyn, tnhusc, dqvdf, tnhusscpbl, dqphy, evu, kz, mcd)
 ChunckJob3D = OFF
 Seasonal=ON
@@ -124,5 +124,5 @@
 TimeSeriesVars3D = ()
 ChunckJob3D = NONE
-Seasonal=ON
+Seasonal=OFF
 
 [Post_1D_histday]
Index: CONFIG/IPSLCM/IPSLCM5A/EXP00/PARAM/physiq.def_L39
===================================================================
--- CONFIG/IPSLCM/IPSLCM5A/EXP00/PARAM/physiq.def_L39	(revision 1506)
+++ CONFIG/IPSLCM/IPSLCM5A/EXP00/PARAM/physiq.def_L39	(revision 1507)
@@ -29,8 +29,4 @@
 # Parametres fichiers de sortie
 #
-### type_run = type run par rapport aux fichiers et variables de sortie
-# - type_run = CLIM/ENSP (=1)
-# - type_run = AMIP/CFMI (=2)
-type_run=CLIM
 ### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc
 OK_journe=n
@@ -42,5 +38,5 @@
 ok_hf=n
 #
-# Parametres niveau de sorties differents fichiers 
+# Parametres niveau de sorties differents fichiers
 #
 ### lev_histhf=0-4, niveau de sortie fichier "histhf.nc" 
Index: CONFIG/IPSLCM/IPSLCM5A/historical/COMP/lmdz.card
===================================================================
--- CONFIG/IPSLCM/IPSLCM5A/historical/COMP/lmdz.card	(revision 1506)
+++ CONFIG/IPSLCM/IPSLCM5A/historical/COMP/lmdz.card	(revision 1507)
@@ -72,17 +72,17 @@
 
 [OutputFiles]
-List=	(histmth.nc,      ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,      Post_1M_histmth),     \
-	(histmthCOSP.nc,  ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthCOSP.nc,  Post_1M_histmthCOSP), \
-	(histmthNMC.nc,   ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthNMC.nc,   Post_1M_histmthNMC),  \
-	(dynzon.nc,       ${R_OUT_ATM_O_M}/${PREFIX}_1M_dynzon.nc,       Post_1M_dynzon),      \
-	(histday.nc,      ${R_OUT_ATM_O_D}/${PREFIX}_1D_histday.nc,      Post_1D_histday),     \
-	(histdayCOSP.nc,  ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayCOSP.nc,  Post_1D_histdayCOSP), \
-	(histdayNMC.nc,   ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayNMC.nc,   Post_1D_histdayNMC),  \
-	(histhf.nc,       ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhf.nc,       Post_HF_histhf),      \
-	(histhfCOSP.nc,   ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfCOSP.nc,   Post_HF_histhfCOSP),  \
-	(histhfNMC.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfNMC.nc,    Post_HF_histhfNMC),   \
-	(histhf3h.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3h.nc,     Post_3H_histhf3h),    \
-        (histhf3hm.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3hm.nc,    Post_3H_histhf3hm),   \
-        (paramLMDZ_phy.nc, ${R_OUT_ATM_O_M}/${PREFIX}_1M_paramLMDZ_phy.nc, Post_1M_paramLMDZ_phy)
+List=	(histmth.nc,       ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,       Post_1M_histmth),       \
+	(histmthCOSP.nc,   ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthCOSP.nc,   Post_1M_histmthCOSP),   \
+	(histmthNMC.nc,    ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthNMC.nc,    Post_1M_histmthNMC),    \
+	(dynzon.nc,        ${R_OUT_ATM_O_M}/${PREFIX}_1M_dynzon.nc,        Post_1M_dynzon),        \
+        (paramLMDZ_phy.nc, ${R_OUT_ATM_O_M}/${PREFIX}_1M_paramLMDZ_phy.nc, Post_1M_paramLMDZ_phy), \
+	(histday.nc,       ${R_OUT_ATM_O_D}/${PREFIX}_1D_histday.nc,       Post_1D_histday),       \
+	(histdayCOSP.nc,   ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayCOSP.nc,   Post_1D_histdayCOSP),   \
+	(histdayNMC.nc,    ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayNMC.nc,    Post_1D_histdayNMC),    \
+	(histhf.nc,        ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhf.nc,        Post_HF_histhf),        \
+	(histhfCOSP.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfCOSP.nc,    Post_HF_histhfCOSP),    \
+	(histhfNMC.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfNMC.nc,     Post_HF_histhfNMC),     \
+	(histhf3h.nc,      ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3h.nc,      Post_3H_histhf3h),      \
+        (histhf3hm.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3hm.nc,     Post_3H_histhf3hm)
 
 [Post_1M_histmth]
@@ -91,5 +91,5 @@
 TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, flat, ftime_con, iwp, lwp, LWdnSFC, LWdnSFCclr, LWupSFCclr, LWdn200, LWdn200clr, LWup200, LWup200clr, precip, prw, pluc, psol, q2m, rh2m, sens, sicf, slp, snow, soll0, soll, sols0, sols, SWdnSFC, SWupSFC, SWup200, SWdnTOA, SWupTOA, SWdn200, SWdnSFCclr, SWupSFCclr, SWupTOAclr, t2m, t2m_max, t2m_min, taux, taux_ter, taux_sic, taux_lic, taux_oce, tauy, tauy_ter, tauy_sic, tauy_lic, tauy_oce, lat_ter, lat_sic, lat_lic, lat_oce, nettop, pbase, ptop, topl0, topl, tops0, tops, topswad, topswai, tsol, tsol_ter, tsol_lic, tsol_oce, tsol_sic , u10m, u500, v10m, v500, w500, pourc_ter, pourc_oce, pourc_sic, pourc_lic, wind10m, od550aer, od550lt1aer, absvisaer, loadoa, loadbc, loadso4, loaddust, loadss, sconcoa, sconcbc, sconcso4, sconcdust, sconcss, reffclwtop, cldncl, swtoaas_ant, swtoacs_ant, swtoacf_ant, swsrfas_ant, swsrfcs_ant, swsrfcf_ant, swtoaas_nat, swtoacs_nat, swtoacf_nat, swsrfas_nat, swsrfcs_nat, swsrfcf_nat, solswai, swtoacf_zero, swsrfcf_zero)
 ChunckJob2D = NONE
-TimeSeriesVars3D = (upwd, lwcon, iwcon, lcc3dcon, lcc3dstra, temp, theta, ovap, ovapinit, geop, vitu, vitv, vitw, pres, rneb, rnebcon, rhum, mc, paprs, scdnc, ec550aer, concoa, concbc, concso4, concss, concdust, reffclws, reffclwc, Aga)
+TimeSeriesVars3D = (upwd, lwcon, iwcon, lcc3dcon, lcc3dstra, temp, theta, ovap, ovapinit, geop, vitu, vitv, vitw, pres, rneb, rnebcon, rhum, mc, paprs, scdnc, ec550aer, concoa, concbc, concso4, concss, concdust, reffclws, reffclwc, Aga, rlu, rsu, rld, rsd, rlucs, rsucs, rldcs, rsdcs, tnt, dtdyn, dtphy, dtlschr, tntr, tntc, tntscpbl, tnhus, dqdyn, tnhusc, dqvdf, tnhusscpbl, dqphy, evu, kz, mcd)
 ChunckJob3D = 100Y
 Seasonal=ON
@@ -129,5 +129,5 @@
 TimeSeriesVars3D = ()
 ChunckJob3D = NONE
-Seasonal=ON
+Seasonal=OFF
 
 [Post_1D_histday]
@@ -188,5 +188,5 @@
 Patches= ()
 GatherWithInternal = (lon, lat, presnivs, time_counter)
-TimeSeriesVars2D = (phis, tsol, t2m, q2m, u10m, v10m, psol, topl, bils)
+TimeSeriesVars2D = (phis, tsol, t2m, q2m, u10m, v10m, psol, topl, bils, t100, u100, v100, w100, t50, u50, v50, w50, t10, u10, v10, w10)
 ChunckJob2D = 50Y
 TimeSeriesVars3D = ()
Index: CONFIG/IPSLCM/IPSLCM5A/historical/PARAM/physiq.def_L39
===================================================================
--- CONFIG/IPSLCM/IPSLCM5A/historical/PARAM/physiq.def_L39	(revision 1506)
+++ CONFIG/IPSLCM/IPSLCM5A/historical/PARAM/physiq.def_L39	(revision 1507)
@@ -29,8 +29,4 @@
 # Parametres fichiers de sortie
 #
-### type_run = type run par rapport aux fichiers et variables de sortie
-# - type_run = CLIM/ENSP (=1)
-# - type_run = AMIP/CFMI (=2)
-type_run=CLIM
 ### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc
 OK_journe=n
@@ -42,25 +38,86 @@
 ok_hf=n
 # Noms des fichiers
-phys_out_filenames=      histmth histday histhf  histhf3h histhf3hm
+phys_out_filenames=      histmth histday histhf histhf3h histhf3hm histstn
 # Sortir ou non les fichiers
-phys_out_filekeys=       y       y       y       y       y
+phys_out_filekeys=       y        y        y        y        y        n
 # Niveaux de sorties
-phys_out_filelevels=     5       5       5       5       5
-# more variables for level 5 : vitw and pres
+phys_out_filelevels=     5        5        5        5        5        5
+### Type de fichier : global (n) ou stations (y)
+phys_out_filestations =  n        n        n        n        n        y
+### Frequences des sorties
+phys_out_filetimesteps = 1.mth, 1.day, 0.25day, 0.125day, 0.125day, 1800.s
+### variables niveaux pressions 100, 50, 10 hPa a 3h-inst - F.Lott
+flag_t100       = 1, 7, 7, 5, 10, 10
+flag_u100       = 1, 7, 7, 5, 10, 10
+flag_v100       = 1, 7, 7, 5, 10, 10
+flag_w100       = 1, 7, 7, 5, 10, 10
+flag_t50        = 1, 7, 7, 5, 10, 10
+flag_u50        = 1, 7, 7, 5, 10, 10
+flag_v50        = 1, 7, 7, 5, 10, 10
+flag_w50        = 1, 7, 7, 5, 10, 10
+flag_t10        = 1, 7, 7, 5, 10, 10
+flag_u10        = 1, 7, 7, 5, 10, 10
+flag_v10        = 1, 7, 7, 5, 10, 10
+flag_w10        = 1, 7, 7, 5, 10, 10
+# more variables for level 5 histday/month: vitw and pres
 # Vertical wind
-flag_vitw         =  2, 3, 7, 6, 6
-name_vitw         =  vitw
+flag_vitw         =  2, 3, 7, 6, 6, 6
 # Air pressure
-flag_pres         =  2, 3, 7, 6, 6
-name_pres         =  pres
+flag_pres         =  2, 3, 7, 6, 6, 6
 # w500 for histday
-flag_w500         =  1, 5, 3, 10, 10
-name_w500         =  w500
+flag_w500         =  1, 5, 3, 10, 10, 10
 # t700 for histday
-flag_t700         =  1, 5, 3, 10, 10
-name_t700         =  t700
+flag_t700         =  1, 5, 3, 10, 10, 10
 # Cloud fraction
-flag_rneb         =  2, 5, 10, 5, 10
-name_rneb         =  rneb
+flag_rneb         =  2, 5, 10, 5, 10, 10
+# Aga sortie avec niveau 2
+flag_Aga          =  2, 10, 10, 10, 10, 10
+# niveaux maximum 5 mensuel pour variables CFMIP
+flag_rlu          =  5, 10, 10, 10, 10, 10
+flag_rsu          =  5, 10, 10, 10, 10, 10
+flag_rld          =  5, 10, 10, 10, 10, 10
+flag_rsd          =  5, 10, 10, 10, 10, 10
+flag_rlucs        =  5, 10, 10, 10, 10, 10
+flag_rsucs        =  5, 10, 10, 10, 10, 10
+flag_rldcs        =  5, 10, 10, 10, 10, 10
+flag_rsdcs        =  5, 10, 10, 10, 10, 10
+flag_tnt          =  5, 10, 10, 10, 10, 10
+flag_dtdyn        =  4, 10, 10, 10, 10, 10
+flag_dtphy        =  4, 10, 10, 10, 10, 10
+flag_dtlschr      =  5, 10, 10, 10, 10, 10
+flag_tntr         =  5, 10, 10, 10, 10, 10
+flag_tntc         =  5, 10, 10, 10, 10, 10
+flag_tntscpbl     =  5, 10, 10, 10, 10, 10
+flag_tnhus        =  5, 10, 10, 10, 10, 10
+flag_dqdyn        =  4, 10, 10, 10, 10, 10
+flag_tnhusc       =  5, 10, 10, 10, 10, 10
+flag_dqvdf        =  5, 10, 10, 10, 10, 10
+flag_tnhusscpbl   =  5, 10, 10, 10, 10, 10
+flag_dqphy        =  2, 10, 10, 10, 10, 10
+flag_evu          =  5, 10, 10, 10, 10, 10
+flag_kz           =  4, 10, 10, 10, 10, 10
+###
+# more variables for t tendency monthly
+###
+flag_dtvdf        =  5, 10, 10, 10, 10, 10
+flag_dteva        =  5, 10, 10, 10, 10, 10
+flag_dtcon        =  4, 10, 10, 10, 10, 10
+flag_dtajs        =  4, 10, 10, 10, 10, 10
+flag_dtlsc        =  4, 10, 10, 10, 10, 10
+flag_dtoro        =  4, 10, 10, 10, 10, 10
+flag_dtlif        =  4, 10, 10, 10, 10, 10
+flag_dthin        =  4, 10, 10, 10, 10, 10
+flag_dtec         =  4, 10, 10, 10, 10, 10
+### more variables for q tendency monthly
+flag_dqeva        =  5, 10, 10, 10, 10, 10
+flag_dqcon        =  4, 10, 10, 10, 10, 10
+flag_dqajs        =  4, 10, 10, 10, 10, 10
+flag_dqlsc        =  5, 10, 10, 10, 10, 10
+### nouvelle physique monthly
+flag_dtthe        =  4, 10, 10, 10, 10, 10
+flag_dtwak        =  4, 5, 10, 10, 10, 10
+flag_dqthe        =  4, 10, 10, 10, 10, 10
+flag_dqwak        =  4, 5, 10, 10, 10, 10
+flag_ftime_th     =  4, 10, 10, 5, 10, 10
 #
 # Parametres niveau de sorties differents fichiers
@@ -85,6 +142,4 @@
 # - lev_histmth=4 => champs tendances 3d
 lev_histmth=5
-### Aga sortie avec niveau 2
-flag_Aga       =  2, 10, 10, 10, 10
 ### parametres pour hist*NMC.nc
 #lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier
@@ -93,10 +148,4 @@
 ### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence
 ok_histNMC = y, n, n
-### ecrit_LES = frequence ecriture fichier histhf3hm.nc en jours
-ecrit_LES=0.125
-### ecrit_ins = frequence ecriture fichier histhf3h.nc en jours
-ecrit_ins=0.125
-### ecrit_hf = frequence ecriture fichier histhf.nc en jours
-ecrit_hf=0.250
 ### ecrit_day = frequence ecriture fichier histday.nc en jours
 ecrit_day=1.
Index: CONFIG/IPSLCM/IPSLCM5A/piControl/COMP/lmdz.card
===================================================================
--- CONFIG/IPSLCM/IPSLCM5A/piControl/COMP/lmdz.card	(revision 1506)
+++ CONFIG/IPSLCM/IPSLCM5A/piControl/COMP/lmdz.card	(revision 1507)
@@ -66,17 +66,17 @@
 
 [OutputFiles]
-List=	(histmth.nc,      ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,      Post_1M_histmth),     \
-	(histmthCOSP.nc,  ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthCOSP.nc,  Post_1M_histmthCOSP), \
-	(histmthNMC.nc,   ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthNMC.nc,   Post_1M_histmthNMC),  \
-	(dynzon.nc,       ${R_OUT_ATM_O_M}/${PREFIX}_1M_dynzon.nc,       Post_1M_dynzon),      \
-	(histday.nc,      ${R_OUT_ATM_O_D}/${PREFIX}_1D_histday.nc,      Post_1D_histday),     \
-	(histdayCOSP.nc,  ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayCOSP.nc,  Post_1D_histdayCOSP), \
-	(histdayNMC.nc,   ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayNMC.nc,   Post_1D_histdayNMC),  \
-	(histhf.nc,       ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhf.nc,       Post_HF_histhf),      \
-	(histhfCOSP.nc,   ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfCOSP.nc,   Post_HF_histhfCOSP),  \
-	(histhfNMC.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfNMC.nc,    Post_HF_histhfNMC),   \
-	(histhf3h.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3h.nc,     Post_3H_histhf3h),    \
-	(histhf3hm.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3hm.nc,    Post_3H_histhf3hm),   \
-        (paramLMDZ_phy.nc, ${R_OUT_ATM_O_M}/${PREFIX}_1M_paramLMDZ_phy.nc, Post_1M_paramLMDZ_phy)
+List=	(histmth.nc,       ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,       Post_1M_histmth),       \
+	(histmthCOSP.nc,   ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthCOSP.nc,   Post_1M_histmthCOSP),   \
+	(histmthNMC.nc,    ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthNMC.nc,    Post_1M_histmthNMC),    \
+	(dynzon.nc,        ${R_OUT_ATM_O_M}/${PREFIX}_1M_dynzon.nc,        Post_1M_dynzon),        \
+        (paramLMDZ_phy.nc, ${R_OUT_ATM_O_M}/${PREFIX}_1M_paramLMDZ_phy.nc, Post_1M_paramLMDZ_phy), \
+	(histday.nc,       ${R_OUT_ATM_O_D}/${PREFIX}_1D_histday.nc,       Post_1D_histday),       \
+	(histdayCOSP.nc,   ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayCOSP.nc,   Post_1D_histdayCOSP),   \
+	(histdayNMC.nc,    ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayNMC.nc,    Post_1D_histdayNMC),    \
+	(histhf.nc,        ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhf.nc,        Post_HF_histhf),        \
+	(histhfCOSP.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfCOSP.nc,    Post_HF_histhfCOSP),    \
+	(histhfNMC.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfNMC.nc,     Post_HF_histhfNMC),     \
+	(histhf3h.nc,      ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3h.nc,      Post_3H_histhf3h),      \
+        (histhf3hm.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3hm.nc,     Post_3H_histhf3hm)
 
 [Post_1M_histmth]
@@ -85,5 +85,5 @@
 TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, flat, ftime_con, iwp, lwp, LWdnSFC, LWdnSFCclr, LWupSFCclr, LWdn200, LWdn200clr, LWup200, LWup200clr, precip, prw, pluc, psol, q2m, rh2m, sens, sicf, slp, snow, soll0, soll, sols0, sols, SWdnSFC, SWupSFC, SWup200, SWdnTOA, SWupTOA, SWdn200, SWdnSFCclr, SWupSFCclr, SWupTOAclr, t2m, t2m_max, t2m_min, taux, taux_ter, taux_sic, taux_lic, taux_oce, tauy, tauy_ter, tauy_sic, tauy_lic, tauy_oce, lat_ter, lat_sic, lat_lic, lat_oce, nettop, pbase, ptop, topl0, topl, tops0, tops, topswad, topswai, tsol, tsol_ter, tsol_lic, tsol_oce, tsol_sic , u10m, u500, v10m, v500, w500, pourc_ter, pourc_oce, pourc_sic, pourc_lic, wind10m, od550aer, od550lt1aer, absvisaer, loadoa, loadbc, loadso4, loaddust, loadss, sconcoa, sconcbc, sconcso4, sconcdust, sconcss, reffclwtop, cldncl, swtoaas_ant, swtoacs_ant, swtoacf_ant, swsrfas_ant, swsrfcs_ant, swsrfcf_ant, swtoaas_nat, swtoacs_nat, swtoacf_nat, swsrfas_nat, swsrfcs_nat, swsrfcf_nat, solswai, swtoacf_zero, swsrfcf_zero)
 ChunckJob2D = NONE
-TimeSeriesVars3D = (upwd, lwcon, iwcon, lcc3dcon, lcc3dstra, temp, theta, ovap, ovapinit, geop, vitu, vitv, vitw, pres, rneb, rnebcon, rhum, mc, paprs, scdnc, ec550aer, concoa, concbc, concso4, concss, concdust, reffclws, reffclwc, Aga)
+TimeSeriesVars3D = (upwd, lwcon, iwcon, lcc3dcon, lcc3dstra, temp, theta, ovap, ovapinit, geop, vitu, vitv, vitw, pres, rneb, rnebcon, rhum, mc, paprs, scdnc, ec550aer, concoa, concbc, concso4, concss, concdust, reffclws, reffclwc, Aga, rlu, rsu, rld, rsd, rlucs, rsucs, rldcs, rsdcs, tnt, dtdyn, dtphy, dtlschr, tntr, tntc, tntscpbl, tnhus, dqdyn, tnhusc, dqvdf, tnhusscpbl, dqphy, evu, kz, mcd)
 ChunckJob3D = 100Y
 Seasonal=ON
@@ -123,12 +123,12 @@
 TimeSeriesVars3D = ()
 ChunckJob3D = NONE
-Seasonal=ON
+Seasonal=OFF
 
 [Post_1D_histday]
 Patches= ()
 GatherWithInternal = (lon, lat, presnivs, time_counter, time_counter_bnds)
-TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, flat, precip, pluc, plul, prw, psol, q2m, sens, sicf, slp, soll, sols, snow, SWdnSFC, t2m, t2m_max, t2m_min, topl, tops, tsol, t700, u10m, v10m, wind10m, wind10max, w500)
-ChunckJob2D = 200Y
-TimeSeriesVars3D = (lwcon, iwcon, mc, rhum,  rneb, temp, theta, ovap, geop, vitu, vitv, vitw, pres, paprs)
+TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, flat, iwp, lwp, pbase, ptop, precip, pluc, plul, prw, psol, rh2m, rh2m_min, rh2m_max, q2m, sens, sicf, slp, soll, sols, snow, LWdnSFC, LWdnSFCclr, SWdnSFC, SWdnSFCclr, SWdnTOA, LWupSFCclr, SWupSFC, SWupSFCclr, SWupTOA, SWupTOAclr, t2m, t2m_max, t2m_min, topl, topl0, tops, tsol, t700, u10m, v10m, wind10m, wind10max, w500)
+ChunckJob2D = 50Y
+TimeSeriesVars3D = (lwcon, iwcon, mc, rhum, rneb, temp, theta, ovap, geop, vitu, vitv, vitw, pres, paprs)
 ChunckJob3D = OFF
 Seasonal=OFF
Index: CONFIG/IPSLCM/IPSLCM5A/piControl/PARAM/physiq.def_L39
===================================================================
--- CONFIG/IPSLCM/IPSLCM5A/piControl/PARAM/physiq.def_L39	(revision 1506)
+++ CONFIG/IPSLCM/IPSLCM5A/piControl/PARAM/physiq.def_L39	(revision 1507)
@@ -29,8 +29,4 @@
 # Parametres fichiers de sortie
 #
-### type_run = type run par rapport aux fichiers et variables de sortie
-# - type_run = CLIM/ENSP (=1)
-# - type_run = AMIP/CFMI (=2)
-type_run=CLIM
 ### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc
 OK_journe=n
@@ -44,19 +40,64 @@
 # more variables for level 5 : vitw and pres
 # Vertical wind
-flag_vitw         =  2, 3, 7, 6, 6
-name_vitw         =  vitw
+flag_vitw         =  2, 3, 7, 6, 6, 6
 # Air pressure
-flag_pres         =  2, 3, 7, 6, 6
-name_pres         =  pres
+flag_pres         =  2, 3, 7, 6, 6, 6
 # w500 for histday
-flag_w500         =  1, 5, 3, 10, 10
-name_w500         =  w500
+flag_w500         =  1, 5, 3, 10, 10, 10
 # t700 for histday
-flag_t700         =  1, 5, 3, 10, 10
-name_t700         =  t700
+flag_t700         =  1, 5, 3, 10, 10, 10
 # Cloud fraction
-flag_rneb         =  2, 5, 10, 5, 10
-name_rneb         =  rneb
-# Parametres niveau de sorties differents fichiers 
+flag_rneb         =  2, 5, 10, 5, 10, 10
+# Aga sortie avec niveau 2
+flag_Aga          =  2, 10, 10, 10, 10, 10
+# niveaux maximum 5 mensuel pour variables CFMIP
+flag_rlu          =  5, 10, 10, 10, 10, 10
+flag_rsu          =  5, 10, 10, 10, 10, 10
+flag_rld          =  5, 10, 10, 10, 10, 10
+flag_rsd          =  5, 10, 10, 10, 10, 10
+flag_rlucs        =  5, 10, 10, 10, 10, 10
+flag_rsucs        =  5, 10, 10, 10, 10, 10
+flag_rldcs        =  5, 10, 10, 10, 10, 10
+flag_rsdcs        =  5, 10, 10, 10, 10, 10
+flag_tnt          =  5, 10, 10, 10, 10, 10
+flag_dtdyn        =  4, 10, 10, 10, 10, 10
+flag_dtphy        =  4, 10, 10, 10, 10, 10
+flag_dtlschr      =  5, 10, 10, 10, 10, 10
+flag_tntr         =  5, 10, 10, 10, 10, 10
+flag_tntc         =  5, 10, 10, 10, 10, 10
+flag_tntscpbl     =  5, 10, 10, 10, 10, 10
+flag_tnhus        =  5, 10, 10, 10, 10, 10
+flag_dqdyn        =  4, 10, 10, 10, 10, 10
+flag_tnhusc       =  5, 10, 10, 10, 10, 10
+flag_dqvdf        =  5, 10, 10, 10, 10, 10
+flag_tnhusscpbl   =  5, 10, 10, 10, 10, 10
+flag_dqphy        =  2, 10, 10, 10, 10, 10
+flag_evu          =  5, 10, 10, 10, 10, 10
+flag_kz           =  4, 10, 10, 10, 10, 10
+###
+# more variables for t tendency monthly
+###
+flag_dtvdf        =  5, 10, 10, 10, 10, 10
+flag_dteva        =  5, 10, 10, 10, 10, 10
+flag_dtcon        =  4, 10, 10, 10, 10, 10
+flag_dtajs        =  4, 10, 10, 10, 10, 10
+flag_dtlsc        =  4, 10, 10, 10, 10, 10
+flag_dtoro        =  4, 10, 10, 10, 10, 10
+flag_dtlif        =  4, 10, 10, 10, 10, 10
+flag_dthin        =  4, 10, 10, 10, 10, 10
+flag_dtec         =  4, 10, 10, 10, 10, 10
+### more variables for q tendency monthly
+flag_dqeva        =  5, 10, 10, 10, 10, 10
+flag_dqcon        =  4, 10, 10, 10, 10, 10
+flag_dqajs        =  4, 10, 10, 10, 10, 10
+flag_dqlsc        =  5, 10, 10, 10, 10, 10
+### nouvelle physique monthly
+flag_dtthe        =  4, 10, 10, 10, 10, 10
+flag_dtwak        =  4, 5, 10, 10, 10, 10
+flag_dqthe        =  4, 10, 10, 10, 10, 10
+flag_dqwak        =  4, 5, 10, 10, 10, 10
+flag_ftime_th     =  4, 10, 10, 5, 10, 10
+#
+# Parametres niveau de sorties differents fichiers
 #
 ### lev_histhf=0-4, niveau de sortie fichier "histhf.nc" 
@@ -79,6 +120,4 @@
 # - lev_histmth=4 => champs tendances 3d
 lev_histmth=5
-### Aga sortie avec niveau 2
-flag_Aga       =  2, 10, 10, 10, 10
 ### parametres pour hist*NMC.nc
 #lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier
@@ -87,8 +126,4 @@
 ### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence
 ok_histNMC = y, n, n
-### ecrit_ins = frequence ecriture fichier histhf3h.nc en jours
-ecrit_ins=0.125
-### ecrit_hf = frequence ecriture fichier histhf.nc en jours
-ecrit_hf=0.250
 ### ecrit_day = frequence ecriture fichier histday.nc en jours
 ecrit_day=1.
