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

Diff of /trunk/Sources/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/Sources/phylmd/ini_histrac.f revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC
# Line 4  module ini_histrac_m Line 4  module ini_histrac_m
4    
5  contains  contains
6    
   
7    subroutine ini_histrac(nid_tra, pdtphys, nq_phys, lessivage)    subroutine ini_histrac(nid_tra, pdtphys, nq_phys, lessivage)
8    
9      ! 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
10    
11        use clesphys, only: ecrit_tra
12      use dimens_m, only: iim, jjm, llm      use dimens_m, only: iim, jjm, llm
13      USE calendar, only: ymds2ju      use disvert_m, only: presnivs
14      USE histbeg_totreg_m, ONLY : histbeg_totreg      use dimphy, only: klon
15        use dynetat0_m, only: day_ref, annee_ref
16        use grid_change, only: gr_phy_write_2d
17        USE histbeg_totreg_m, ONLY: histbeg_totreg
18      USE histdef_m, ONLY : histdef      USE histdef_m, ONLY : histdef
19      USE histend_m, ONLY : histend      USE histend_m, ONLY : histend
20      USE histvert_m, ONLY : histvert      USE histvert_m, ONLY : histvert
21      use temps, only: annee_ref, day_ref, itau_phy      use iniadvtrac_m, only: tname, ttext
     use iniadvtrac_m, only: niadv, tnom, ttext  
     use dimphy, only: klon  
     use clesphys, only: ecrit_tra  
     use grid_change, only: gr_phy_write_2d  
22      use phyetat0_m, only: rlon, rlat      use phyetat0_m, only: rlon, rlat
23      use disvert_m, only: presnivs      use temps, only: 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 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, tname(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"//tname(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_"//tname(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_"//tname(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_"//tname(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.134

  ViewVC Help
Powered by ViewVC 1.1.21