/[lmdze]/trunk/phylmd/ini_histrac.f
ViewVC logotype

Diff of /trunk/phylmd/ini_histrac.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/libf/phylmd/ini_histrac.f90 revision 66 by guez, Thu Sep 20 13:00:41 2012 UTC trunk/phylmd/ini_histrac.f revision 91 by guez, Wed Mar 26 17:18:58 2014 UTC
# Line 11  contains Line 11  contains
11    
12      use dimens_m, only: iim, jjm, llm      use dimens_m, only: iim, jjm, llm
13      USE calendar, only: ymds2ju      USE calendar, only: ymds2ju
14      USE histbeg_totreg_m, ONLY : histbeg_totreg      USE histbeg_totreg_m, ONLY: histbeg_totreg
15      USE histdef_m, ONLY : histdef      USE histdef_m, ONLY : histdef
16      USE histend_m, ONLY : histend      USE histend_m, ONLY : histend
17      USE histvert_m, ONLY : histvert      USE histvert_m, ONLY : histvert
18      use temps, only: annee_ref, day_ref, itau_phy      use temps, only: annee_ref, day_ref, itau_phy
19      use iniadvtrac_m, only: niadv, tnom, ttext      use iniadvtrac_m, only: tnom, ttext
20      use dimphy, only: klon      use dimphy, only: klon
21      use clesphys, only: ecrit_tra      use clesphys, only: ecrit_tra
22      use grid_change, only: gr_phy_write_2d      use grid_change, only: gr_phy_write_2d
# Line 37  contains Line 37  contains
37      REAL zx_lat(iim, jjm+1)      REAL zx_lat(iim, jjm+1)
38      INTEGER nhori, nvert      INTEGER nhori, nvert
39      REAL zsto, zout      REAL zsto, zout
40      integer it, iq, iiq      integer it, iq
41    
42      !---------------------------------------------------------      !---------------------------------------------------------
43    
# Line 45  contains Line 45  contains
45      zx_lat(:, :) = gr_phy_write_2d(rlat)      zx_lat(:, :) = gr_phy_write_2d(rlat)
46      CALL histbeg_totreg("histrac", rlon(2:iim+1), zx_lat(1, :), &      CALL histbeg_totreg("histrac", rlon(2:iim+1), zx_lat(1, :), &
47           1, iim, 1, jjm+1, itau_phy, zjulian, pdtphys, nhori, nid_tra)           1, iim, 1, jjm+1, itau_phy, zjulian, pdtphys, nhori, nid_tra)
48      CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb", llm, &      CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb", presnivs, nvert)
          presnivs, nvert)  
49    
50      zsto = pdtphys      zsto = pdtphys
51      zout = pdtphys * REAL(ecrit_tra)      zout = pdtphys * REAL(ecrit_tra)
# Line 64  contains Line 63  contains
63      DO it = 1, nq_phys      DO it = 1, nq_phys
64         ! champ 2D         ! champ 2D
65         iq=it+2         iq=it+2
66         iiq=niadv(iq)         CALL histdef(nid_tra, tnom(iq), ttext(iq), "U/kga", iim, jjm+1, &
        CALL histdef(nid_tra, tnom(iq), ttext(iiq), "U/kga", iim, jjm+1, &  
67              nhori, llm, 1, llm, nvert, "ave(X)", zsto, zout)              nhori, llm, 1, llm, nvert, "ave(X)", zsto, zout)
68         if (lessivage) THEN         if (lessivage) THEN
69            CALL histdef(nid_tra, "fl"//tnom(iq), "Flux "//ttext(iiq), &            CALL histdef(nid_tra, "fl"//tnom(iq), "Flux "//ttext(iq), &
70                 "U/m2/s", iim, jjm+1, nhori, llm, 1, llm, nvert, &                 "U/m2/s", iim, jjm+1, nhori, llm, 1, llm, nvert, &
71                 "ave(X)", zsto, zout)                 "ave(X)", zsto, zout)
72         endif         endif
73    
74         !---Ajout Olivia         !---Ajout Olivia
75         CALL histdef(nid_tra, "d_tr_th_"//tnom(iq), &         CALL histdef(nid_tra, "d_tr_th_"//tnom(iq), &
76              "tendance thermique"// ttext(iiq), "?", &              "tendance thermique"// ttext(iq), "?", &
77              iim, jjm+1, nhori, llm, 1, llm, nvert, &              iim, jjm+1, nhori, llm, 1, llm, nvert, &
78              "ave(X)", zsto, zout)              "ave(X)", zsto, zout)
79         CALL histdef(nid_tra, "d_tr_cv_"//tnom(iq), &         CALL histdef(nid_tra, "d_tr_cv_"//tnom(iq), &
80              "tendance convection"// ttext(iiq), "?", &              "tendance convection"// ttext(iq), "?", &
81              iim, jjm+1, nhori, llm, 1, llm, nvert, &              iim, jjm+1, nhori, llm, 1, llm, nvert, &
82              "ave(X)", zsto, zout)              "ave(X)", zsto, zout)
83         CALL histdef(nid_tra, "d_tr_cl_"//tnom(iq), &         CALL histdef(nid_tra, "d_tr_cl_"//tnom(iq), &
84              "tendance couche limite"// ttext(iiq), "?", &              "tendance couche limite"// ttext(iq), "?", &
85              iim, jjm+1, nhori, llm, 1, llm, nvert, &              iim, jjm+1, nhori, llm, 1, llm, nvert, &
86              "ave(X)", zsto, zout)              "ave(X)", zsto, zout)
87         !---fin Olivia         !---fin Olivia
Legend:
Removed from v.66  
changed lines
  Added in v.91
  ViewVC Help
Powered by ViewVC 1.1.21