c
      dt_cum = dtvr*day_step

      zan = annee_ref
      dayref = day_ref
      CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)
      tau0 = itau_dyn
c
       pi = 4.0 * ATAN(1.0)
       degres = 180./pi
       rlong = rlonu * degres
       rlatg = rlatu * degres
c
      CALL histbeg("paramLMDZ_dyn.nc", 
     .                 iip1,rlong, jjp1,rlatg,
     .                 1,1,1,1,
     .                 tau0, zjulian, dt_cum,
     .                 thoriid, nid_ctesGCM)
c
         CALL histdef(nid_ctesGCM, "prt_level", 
     .        "Niveau impression debuggage dynamique",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "dayref", 
     .        "Jour de l etat initial ( = 350  si 20 Decembre par ex.)",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "anneeref", 
     .        "Annee de l etat initial",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "anneelim", 
     .        "Annee du fichier limitxxxx.nc  si  ok_limitvrai =y",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "raz_date", 
     .   "Remise a zero (raz) date init.: 0 pas de raz;1=date gcm.def",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "nday", 
     .   "Nombre de jours d integration",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "day_step", 
     .   "nombre de pas par jour pour dt = 1 min",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "iperiod", 
     .   "periode pour le pas Matsuno (en pas de temps)",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "iapp_tracvl", 
     .   "frequence du groupement des flux (en pas de temps)",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "iconser", 
     .  "periode de sortie des variables de controle (en pas de temps)",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "iecri", 
     .  "periode d ecriture du fichier histoire (en jour)",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "periodav", 
     .  "periode de stockage fichier histmoy (en jour)",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "idissip", 
     .  "periode de la dissipation (en pas) ... a completer",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "lstardis", 
     .  "choix de l operateur de dissipation: 1= star,0=non-star ??",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "nitergdiv", 
     .  "nombre d iterations de l operateur de dissipation gradiv",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "nitergrot", 
     .  "nombre d iterations de l operateur de dissipation nxgradrot",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "niterh", 
     .  "nombre d iterations de l operateur de dissipation divgrad",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "tetagdiv", 
     ."temps dissipation des + petites long. d ondes pour u,v (gradiv)",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "tetagrot", 
     ."temps diss. des + petites long. d ondes pour u,v (nxgradrot)",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "tetatemp", 
     ."temps diss. des + petites long. d ondes pour h (divgrad)",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "coefdis", 
     ."coefficient pour gamdissip",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "purmats", 
     ."Choix schema integration temporel: 1=Matsuno,0=Matsuno-leapfrog",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "ok_guide", 
     ."Guidage: 1=true ,0=false",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "true_calendar", 
     ."Choix du calendrier: 1=gregorien ,0=calen. a 360 j",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "guide_calend", 
     ."Guidage calendrier gregorien: 1=oui ,0=non",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "iflag_phys", 
     ."Permet de faire tourner le modele sans physique: 1=avec ,0=sans",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "iphysiq", 
     ."Periode de la physique en pas de temps de la dynamique",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "ecritphy", 
     ."frequence  de l ecriture du fichier histphy en jours",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "cycle_diurne", 
     ."Avoir un cycle diurne: 1=true, 0=false",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "soil_model", 
     ."Choix du modele de sol (Thermique ?): 1=true, 0=false",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "new_oliq", 
     ."Nouvelle eau liquide: 1=true, 0=false",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "ok_orodr", 
     ."Orography drag : 1=true, 0=false",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "ok_orolf", 
     ."Orography lift : 1=true, 0=false",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "ok_limitvrai", 
     ."Force la lecture de la bonne annee: 1=true, 0=false",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "nbapp_rad", 
     ."Nombre  d appels des routines de rayonnements par jour",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "iflag_con", 
     ."Flag  pour la convection: 1=LMD, 2=Tiedtke, 3=KE, 4=KE vect.",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "clon", 
     ."longitude en degres du centre du zoom",
     .                "deg",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "clat", 
     ."latitude en degres du centre du zoom",
     .                "deg",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "grossismx", 
     ."facteur de grossissement du zoom, selon la longitude",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "grossismy", 
     ."facteur de grossissement du zoom, selon la latitude",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "fxyhypb", 
     ."Fonction f(y) hyperbolique  si true=1, sinusoidale si false=0",
     .                "-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "dzoomx", 
     ."extension en longitude de la zone du zoom (fraction zone totale)"
     .                ,"-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "dzoomy", 
     ."extension en latitude de la zone du zoom (fraction zone totale)"
     .                ,"-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "taux", 
     ."raideur du zoom en  X"
     .                ,"-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "tauy", 
     ."raideur du zoom en  Y"
     .                ,"-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "ysinus", 
     ."ysinus=1: Ftion f(y) avec y=Sin(latit.)/ ysinus=0: y = latit"
     .                ,"-",iip1,jjp1,thoriid, 1,1,1, -99, 32, 
     .                "once", dt_cum,dt_cum)
c
         CALL histdef(nid_ctesGCM, "ip_ebil_dyn", 
     ."PRINTlevel for energy conservation diag.; 0/1= pas de print,
     . 2= print","-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
     .                "once", dt_cum,dt_cum)
c
c=================================================================
c
         CALL histend(nid_ctesGCM)
c
c=================================================================
