Index: LMDZ4/trunk/libf/phylmd/calcul_STDlev.h
===================================================================
--- LMDZ4/trunk/libf/phylmd/calcul_STDlev.h	(revision 1146)
+++ LMDZ4/trunk/libf/phylmd/calcul_STDlev.h	(revision 1352)
@@ -3,56 +3,10 @@
 c
 c
-cIM on initialise les champs en debut du jour ou du mois 
+cIM on initialise les variables 
 c
-        CALL ini_undefSTD(nlevSTD,itap,
-     $           dtime,ecrit_day,ecrit_mth,
-     $                    tnondef,tsumSTD)
-        CALL ini_undefSTD(nlevSTD,itap,
-     $           dtime,ecrit_day,ecrit_mth,
-     $                    tnondef,usumSTD)
-        CALL ini_undefSTD(nlevSTD,itap,
-     $           dtime,ecrit_day,ecrit_mth,
-     $                    tnondef,vsumSTD)
-        CALL ini_undefSTD(nlevSTD,itap,
-     $           dtime,ecrit_day,ecrit_mth,
-     $                    tnondef,wsumSTD)
-        CALL ini_undefSTD(nlevSTD,itap,
-     $           dtime,ecrit_day,ecrit_mth,
-     $                    tnondef,phisumSTD)
-        CALL ini_undefSTD(nlevSTD,itap,
-     $           dtime,ecrit_day,ecrit_mth,
-     $                    tnondef,qsumSTD)
-        CALL ini_undefSTD(nlevSTD,itap,
-     $           dtime,ecrit_day,ecrit_mth,
-     $                    tnondef,rhsumSTD)
-        CALL ini_undefSTD(nlevSTD,itap,
-     $           dtime,ecrit_day,ecrit_mth,
-     $                    tnondef,uvsumSTD)
-        CALL ini_undefSTD(nlevSTD,itap,
-     $           dtime,ecrit_day,ecrit_mth,
-     $                    tnondef,vqsumSTD)
-        CALL ini_undefSTD(nlevSTD,itap,
-     $           dtime,ecrit_day,ecrit_mth,
-     $                    tnondef,vTsumSTD)
-        CALL ini_undefSTD(nlevSTD,itap,
-     $           dtime,ecrit_day,ecrit_mth,
-     $                    tnondef,wqsumSTD)
-        CALL ini_undefSTD(nlevSTD,itap,
-     $           dtime,ecrit_day,ecrit_mth,
-     $                    tnondef,vphisumSTD)
-        CALL ini_undefSTD(nlevSTD,itap,
-     $           dtime,ecrit_day,ecrit_mth,
-     $                    tnondef,wTsumSTD)
-        CALL ini_undefSTD(nlevSTD,itap,
-     $           dtime,ecrit_day,ecrit_mth,
-     $                    tnondef,u2sumSTD)
-        CALL ini_undefSTD(nlevSTD,itap,
-     $           dtime,ecrit_day,ecrit_mth,
-     $                    tnondef,v2sumSTD)
-        CALL ini_undefSTD(nlevSTD,itap,
-     $           dtime,ecrit_day,ecrit_mth,
-     $                    tnondef,T2sumSTD)
+        CALL ini_undefSTD(itap,freq_outNMC)
 c
-cIM on interpole sur les niveaux STD de pression a chaque pas de temps de la physique
+cIM on interpole les champs sur les niveaux STD de pression 
+cIM a chaque pas de temps de la physique
 c
 c-------------------------------------------------------c
@@ -163,154 +117,16 @@
 
 c
-cIM on somme les valeurs definies a chaque pas de temps de la physique ou 
-cIM toutes les 6 heures
+cIM on somme les valeurs toutes les freq_calNMC secondes
 c
-       oknondef(1:klon,1:nlevSTD,1:nout)=.TRUE.
-       CALL undefSTD(nlevSTD,itap,tlevSTD,
-     $               dtime,ecrit_hf,
-     $               oknondef,tnondef,tsumSTD)
+       CALL undefSTD(itap,freq_calNMC)
 c
-       oknondef(1:klon,1:nlevSTD,1:nout)=.FALSE.
-       CALL undefSTD(nlevSTD,itap,ulevSTD,
-     $           dtime,ecrit_hf,
-     $               oknondef,tnondef,usumSTD)
+cIM on moyenne a la fin du mois ou du jour (toutes les freq_outNMC secondes)
 c
-       oknondef(1:klon,1:nlevSTD,1:nout)=.FALSE.
-       CALL undefSTD(nlevSTD,itap,vlevSTD,
-     $           dtime,ecrit_hf,
-     $               oknondef,tnondef,vsumSTD)
+       CALL moy_undefSTD(itap,freq_outNMC,freq_moyNMC)
 c
-       oknondef(1:klon,1:nlevSTD,1:nout)=.FALSE.
-       CALL undefSTD(nlevSTD,itap,wlevSTD,
-     $           dtime,ecrit_hf,
-     $               oknondef,tnondef,wsumSTD)
-c
-       oknondef(1:klon,1:nlevSTD,1:nout)=.FALSE.
-       CALL undefSTD(nlevSTD,itap,philevSTD,
-     $           dtime,ecrit_hf,
-     $               oknondef,tnondef,phisumSTD)
-c
-       oknondef(1:klon,1:nlevSTD,1:nout)=.FALSE.
-       CALL undefSTD(nlevSTD,itap,qlevSTD,
-     $           dtime,ecrit_hf,
-     $               oknondef,tnondef,qsumSTD)
-c
-       oknondef(1:klon,1:nlevSTD,1:nout)=.FALSE.
-       CALL undefSTD(nlevSTD,itap,rhlevSTD,
-     $           dtime,ecrit_hf,
-     $               oknondef,tnondef,rhsumSTD)
-c
-       oknondef(1:klon,1:nlevSTD,1:nout)=.FALSE.
-       CALL undefSTD(nlevSTD,itap,uvSTD,
-     $           dtime,ecrit_hf,
-     $               oknondef,tnondef,uvsumSTD)
-c
-       oknondef(1:klon,1:nlevSTD,1:nout)=.FALSE.
-       CALL undefSTD(nlevSTD,itap,vqSTD,
-     $           dtime,ecrit_hf,
-     $               oknondef,tnondef,vqsumSTD)
-c
-       oknondef(1:klon,1:nlevSTD,1:nout)=.FALSE.
-       CALL undefSTD(nlevSTD,itap,vTSTD,
-     $           dtime,ecrit_hf,
-     $               oknondef,tnondef,vTsumSTD)
-c
-       oknondef(1:klon,1:nlevSTD,1:nout)=.FALSE.
-       CALL undefSTD(nlevSTD,itap,wqSTD,
-     $           dtime,ecrit_hf,
-     $               oknondef,tnondef,wqsumSTD)
-c
-       oknondef(1:klon,1:nlevSTD,1:nout)=.FALSE.
-       CALL undefSTD(nlevSTD,itap,vphiSTD,
-     $           dtime,ecrit_hf,
-     $               oknondef,tnondef,vphisumSTD)
-c
-       oknondef(1:klon,1:nlevSTD,1:nout)=.FALSE.
-       CALL undefSTD(nlevSTD,itap,wTSTD,
-     $           dtime,ecrit_hf,
-     $               oknondef,tnondef,wTsumSTD)
-c
-       oknondef(1:klon,1:nlevSTD,1:nout)=.FALSE.
-       CALL undefSTD(nlevSTD,itap,u2STD,
-     $           dtime,ecrit_hf,
-     $               oknondef,tnondef,u2sumSTD)
-c
-       oknondef(1:klon,1:nlevSTD,1:nout)=.FALSE.
-       CALL undefSTD(nlevSTD,itap,v2STD,
-     $           dtime,ecrit_hf,
-     $               oknondef,tnondef,v2sumSTD)
-c
-       oknondef(1:klon,1:nlevSTD,1:nout)=.FALSE.
-       CALL undefSTD(nlevSTD,itap,T2STD,
-     $           dtime,ecrit_hf,
-     $               oknondef,tnondef,T2sumSTD)
-c
-cIM on moyenne a la fin du jour ou du mois 
-c
-       CALL moy_undefSTD(nlevSTD,itap,
-     $      dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
-     $                   tnondef,tsumSTD)
-c
-       CALL moy_undefSTD(nlevSTD,itap,
-     $      dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
-     $                   tnondef,usumSTD)
-c
-       CALL moy_undefSTD(nlevSTD,itap,
-     $      dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
-     $                   tnondef,vsumSTD)
-c
-       CALL moy_undefSTD(nlevSTD,itap,
-     $      dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
-     $                   tnondef,wsumSTD)
-c
-       CALL moy_undefSTD(nlevSTD,itap,
-     $      dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
-     $                   tnondef,phisumSTD)
-c
-       CALL moy_undefSTD(nlevSTD,itap,
-     $      dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
-     $                   tnondef,qsumSTD)
-c
-       CALL moy_undefSTD(nlevSTD,itap,
-     $      dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
-     $                   tnondef,rhsumSTD)
-c
-       CALL moy_undefSTD(nlevSTD,itap,
-     $      dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
-     $                   tnondef,uvsumSTD)
-c
-       CALL moy_undefSTD(nlevSTD,itap,
-     $      dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
-     $                   tnondef,vqsumSTD)
-c
-       CALL moy_undefSTD(nlevSTD,itap,
-     $      dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
-     $                   tnondef,vTsumSTD)
-c
-       CALL moy_undefSTD(nlevSTD,itap,
-     $      dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
-     $                   tnondef,wqsumSTD)
-c
-       CALL moy_undefSTD(nlevSTD,itap,
-     $      dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
-     $                   tnondef,vphisumSTD)
-c
-       CALL moy_undefSTD(nlevSTD,itap,
-     $      dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
-     $                   tnondef,wTsumSTD)
-c
-       CALL moy_undefSTD(nlevSTD,itap,
-     $      dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
-     $                   tnondef,u2sumSTD)
-c
-       CALL moy_undefSTD(nlevSTD,itap,
-     $      dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
-     $                   tnondef,v2sumSTD)
-c
-       CALL moy_undefSTD(nlevSTD,itap,
-     $      dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
-     $                   tnondef,T2sumSTD)
-c
-cIM interpolation a chaque pas de temps du SWup(clr) et SWdn(clr) a 200 hPa
+       CALL plevel(klon,klev,.true.,pplay,50000.,
+     &              zphi/RG,geo500)
+
+cIM on interpole a chaque pas de temps le SWup(clr) et SWdn(clr) a 200 hPa
 c
       CALL plevel(klon,klevp1,.true.,paprs,20000.,
