/[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 34 by guez, Wed Jun 2 11:01:12 2010 UTC trunk/phylmd/ini_histrac.f revision 92 by guez, Wed Mar 26 18:16:05 2014 UTC
# Line 4  module ini_histrac_m Line 4  module ini_histrac_m
4    
5  contains  contains
6    
7    
8    subroutine ini_histrac(nid_tra, pdtphys, nq_phys, lessivage)    subroutine ini_histrac(nid_tra, pdtphys, nq_phys, lessivage)
9    
10      ! From phylmd/ini_histrac.h, version 1.10 2006/02/21 08:08:30      ! From phylmd/ini_histrac.h, version 1.10 2006/02/21 08:08:30
11    
12        use clesphys, only: ecrit_tra
13      use dimens_m, only: iim, jjm, llm      use dimens_m, only: iim, jjm, llm
14      USE calendar, only: ymds2ju      use disvert_m, only: presnivs
     use histcom, only: histbeg_totreg, histvert, histend, histdef  
     use temps, only: annee_ref, day_ref, itau_phy  
     use iniadvtrac_m, only: niadv, tnom, ttext  
15      use dimphy, only: klon      use dimphy, only: klon
     use clesphys, only: ecrit_tra  
16      use grid_change, only: gr_phy_write_2d      use grid_change, only: gr_phy_write_2d
17        USE histbeg_totreg_m, ONLY: histbeg_totreg
18        USE histdef_m, ONLY : histdef
19        USE histend_m, ONLY : histend
20        USE histvert_m, ONLY : histvert
21        use iniadvtrac_m, only: tnom, ttext
22      use phyetat0_m, only: rlon, rlat      use phyetat0_m, only: rlon, rlat
23      use comvert, only: presnivs      use temps, only: annee_ref, day_ref, itau_phy
24        USE ymds2ju_m, only: ymds2ju
25    
26      INTEGER, intent(out):: nid_tra      INTEGER, intent(out):: nid_tra
27      real, intent(in):: pdtphys  ! pas d'integration pour la physique (s)      real, intent(in):: pdtphys  ! pas d'integration pour la physique (s)
# Line 33  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 41  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 60  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.34  
changed lines
  Added in v.92
  ViewVC Help
Powered by ViewVC 1.1.21