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

Diff of /trunk/phylmd/ini_histins.f

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

trunk/Sources/phylmd/ini_histins.f revision 191 by guez, Mon May 9 19:56:28 2016 UTC trunk/phylmd/ini_histins.f revision 254 by guez, Mon Feb 5 10:39:38 2018 UTC
# Line 6  module ini_histins_m Line 6  module ini_histins_m
6    
7  contains  contains
8    
9    subroutine ini_histins(dtime)    subroutine ini_histins(dtime, ok_newmicro)
10    
11      ! From phylmd/ini_histins.h, v 1.2 2005/05/25 13:10:09      ! From phylmd/ini_histins.h, version 1.2, 2005/05/25 13:10:09
12    
13      use clesphys, only: ecrit_ins, ok_instan      use clesphys, only: ecrit_ins, ok_instan
14      use clesphys2, only: conv_emanuel      use clesphys2, only: conv_emanuel
15      use dimens_m, only: iim, jjm, llm, nqmx      use dimens_m, only: iim, jjm, llm, nqmx
     use dimphy, only: klon  
16      use disvert_m, only: presnivs      use disvert_m, only: presnivs
17      use dynetat0_m, only: day_ref, annee_ref, rlatu, rlonv      use dynetat0_m, only: day_ref, annee_ref, rlatu, rlonv
18      USE histbeg_totreg_m, ONLY : histbeg_totreg      USE histbeg_totreg_m, ONLY : histbeg_totreg
# Line 27  contains Line 26  contains
26      USE ymds2ju_m, only: ymds2ju      USE ymds2ju_m, only: ymds2ju
27    
28      REAL, intent(in):: dtime ! pas temporel de la physique (s)      REAL, intent(in):: dtime ! pas temporel de la physique (s)
29        logical, intent(in):: ok_newmicro
30        
31      ! Local:      ! Local:
32      real zjulian, zsto, zout      real zjulian, zsto, zout
33      integer nhori, nvert, nsrf, iq, it      integer nhori, nvert, nsrf, iq, it
34    
35      !-------------------------------------------------------------------      !-------------------------------------------------------------------
36    
37      IF (ok_instan) THEN      print *, 'Call sequence information: ini_histins'
38    
39        test_ok_instan: IF (ok_instan) THEN
40         zsto = dtime * ecrit_ins         zsto = dtime * ecrit_ins
41         zout = dtime * ecrit_ins         zout = dtime * ecrit_ins
42         CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)         CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
43         CALL histbeg_totreg("histins", rlonv(:iim) / pi * 180., &         CALL histbeg_totreg("histins", rlonv(:iim) / pi * 180., &
44              rlatu / pi * 180., 1, iim, &              rlatu / pi * 180., 1, iim, &
45              1, jjm + 1, itau_phy, zjulian, dtime, nhori, nid_ins)              1, jjm + 1, itau_phy, zjulian, dtime, nhori, nid_ins)
46         write(*, *)'Inst ', itau_phy, zjulian         print *, 'itau_phy = ', itau_phy
47           print *, "zjulian = ", zjulian
48         CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", &         CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", &
49              presnivs/100., nvert)              presnivs/100., nvert)
50         CALL histdef(nid_ins, "phis", "Surface geop. height", "-", &        
51           CALL histdef(nid_ins, "phis", "surface geopotential", "m2 s-2", &
52              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
53              "once", zsto, zout)              "once", zsto, zout)
54         CALL histdef(nid_ins, "aire", "Grid area", "-", &         CALL histdef(nid_ins, "aire", "Grid area", "-", &
# Line 56  contains Line 60  contains
60         CALL histdef(nid_ins, "tsol", "Surface Temperature", "K", &         CALL histdef(nid_ins, "tsol", "Surface Temperature", "K", &
61              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
62              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
63         CALL histdef(nid_ins, "t2m", "Temperature 2m", "K", &         CALL histdef(nid_ins, "t2m", "Temperature 2m", "K", &
64              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
65              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
66         CALL histdef(nid_ins, "q2m", "Specific humidity 2m", "Kg/Kg", &         CALL histdef(nid_ins, "q2m", "Specific humidity 2m", "Kg/Kg", &
67              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
68              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
69         CALL histdef(nid_ins, "u10m", "Vent zonal 10m", "m/s", &         CALL histdef(nid_ins, "u10m", "Vent zonal 10m", "m/s", &
70              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
71              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
72         CALL histdef(nid_ins, "v10m", "Vent meridien 10m", "m/s", &         CALL histdef(nid_ins, "v10m", "Vent meridien 10m", "m/s", &
73              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
74              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
75         CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa", &         CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa", &
76              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
77              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
78         CALL histdef(nid_ins, "plul", "Large-scale Precip.", "mm/day", &         CALL histdef(nid_ins, "plul", "Large-scale Precip.", "mm/day", &
79              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
80              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
81         CALL histdef(nid_ins, "pluc", "Convective Precip.", "mm/day", &         CALL histdef(nid_ins, "pluc", "Convective Precip.", "mm/day", &
82              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
83              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
84         CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-", &         CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-", &
85              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
86              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
87         CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-", &         CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-", &
88              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
89              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
90         CALL histdef(nid_ins, "precip", "Precipitation Totale liq+sol",  &         CALL histdef(nid_ins, "precip", "Precipitation Totale liq+sol",  &
91              "kg/(s*m2)", &              "kg/(s*m2)", &
92              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
93              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
94         CALL histdef(nid_ins, "snow", "Snow fall", "kg/(s*m2)", &         CALL histdef(nid_ins, "snow", "Snow fall", "kg/(s*m2)", &
95              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
96              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
        !        CALL histdef(nid_ins, "snow_mass", "Snow Mass", "kg/m2",  
        !    .                iim, (jjm + 1), nhori, 1, 1, 1, -99,  
        !    .                "inst(X)", zsto, zout)  
   
97         CALL histdef(nid_ins, "topl", "OLR", "W/m2", &         CALL histdef(nid_ins, "topl", "OLR", "W/m2", &
98              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
99              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
100         CALL histdef(nid_ins, "evap", "Evaporation", "kg/(s*m2)", &         CALL histdef(nid_ins, "evap", "Evaporation", "kg/(s*m2)", &
101              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
102              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
103         CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2", &         CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2", &
104              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
105              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
106         CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2", &         CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2", &
107              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
108              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
109         CALL histdef(nid_ins, "solldown", "Down. IR rad. at surface",  &         CALL histdef(nid_ins, "solldown", "Down. IR rad. at surface",  &
110              "W/m2", iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              "W/m2", iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
111              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
112         CALL histdef(nid_ins, "bils", "Surf. total heat flux", "W/m2", &         CALL histdef(nid_ins, "bils", "Surf. total heat flux", "W/m2", &
113              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
114              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
115         CALL histdef(nid_ins, "sens", "Sensible heat flux", "W/m2", &         CALL histdef(nid_ins, "sens", "Sensible heat flux", "W/m2", &
116              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
117              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
118         CALL histdef(nid_ins, "fder", "Heat flux derivation", "W/m2", &         CALL histdef(nid_ins, "fder", "Heat flux derivation", "W/m2", &
119              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
120              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
121         CALL histdef(nid_ins, "dtsvdfo", "Boundary-layer dTs(o)", "K/s", &         CALL histdef(nid_ins, "dtsvdfo", "Boundary-layer dTs(o)", "K/s", &
122              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
123              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
124         CALL histdef(nid_ins, "dtsvdft", "Boundary-layer dTs(t)", "K/s", &         CALL histdef(nid_ins, "dtsvdft", "Boundary-layer dTs(t)", "K/s", &
125              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
126              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
127         CALL histdef(nid_ins, "dtsvdfg", "Boundary-layer dTs(g)", "K/s", &         CALL histdef(nid_ins, "dtsvdfg", "Boundary-layer dTs(g)", "K/s", &
128              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
129              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
130         CALL histdef(nid_ins, "dtsvdfi", "Boundary-layer dTs(g)", "K/s", &         CALL histdef(nid_ins, "dtsvdfi", "Boundary-layer dTs(g)", "K/s", &
131              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
132              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
133           CALL histdef(nid_ins, "msnow", "surface snow amount", "kg/m2", &
134                iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout)
135    
136         DO nsrf = 1, nbsrf         DO nsrf = 1, nbsrf
   
137            call histdef(nid_ins, "pourc_"//clnsurf(nsrf),  &            call histdef(nid_ins, "pourc_"//clnsurf(nsrf),  &
138                 "% "//clnsurf(nsrf), "%",   &                 "% "//clnsurf(nsrf), "%",   &
139                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
140                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
   
141            call histdef(nid_ins, "fract_"//clnsurf(nsrf),  &            call histdef(nid_ins, "fract_"//clnsurf(nsrf),  &
142                 "Fraction "//clnsurf(nsrf), "1",   &                 "Fraction "//clnsurf(nsrf), "1",   &
143                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
144                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
   
145            call histdef(nid_ins, "sens_"//clnsurf(nsrf),  &            call histdef(nid_ins, "sens_"//clnsurf(nsrf),  &
146                 "Sensible heat flux "//clnsurf(nsrf), "W/m2",   &                 "Sensible heat flux "//clnsurf(nsrf), "W/m2",   &
147                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
148                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
   
149            call histdef(nid_ins, "tsol_"//clnsurf(nsrf),  &            call histdef(nid_ins, "tsol_"//clnsurf(nsrf),  &
150                 "Surface Temperature"//clnsurf(nsrf), "W/m2",   &                 "Surface Temperature"//clnsurf(nsrf), "W/m2",   &
151                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
152                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
   
153            call histdef(nid_ins, "lat_"//clnsurf(nsrf),  &            call histdef(nid_ins, "lat_"//clnsurf(nsrf),  &
154                 "Latent heat flux "//clnsurf(nsrf), "W/m2",   &                 "Latent heat flux "//clnsurf(nsrf), "W/m2",   &
155                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
156                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
   
157            call histdef(nid_ins, "taux_"//clnsurf(nsrf),  &            call histdef(nid_ins, "taux_"//clnsurf(nsrf),  &
158                 "Zonal wind stress"//clnsurf(nsrf), "Pa", &                 "Zonal wind stress"//clnsurf(nsrf), "Pa", &
159                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
160                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
   
161            call histdef(nid_ins, "tauy_"//clnsurf(nsrf),  &            call histdef(nid_ins, "tauy_"//clnsurf(nsrf),  &
162                 "Meridional xind stress "//clnsurf(nsrf), "Pa",   &                 "Meridional xind stress "//clnsurf(nsrf), "Pa",   &
163                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
164                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
   
165            call histdef(nid_ins, "albe_"//clnsurf(nsrf),  &            call histdef(nid_ins, "albe_"//clnsurf(nsrf),  &
166                 "Albedo "//clnsurf(nsrf), "-",   &                 "Albedo "//clnsurf(nsrf), "-",   &
167                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
168                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
   
169            call histdef(nid_ins, "rugs_"//clnsurf(nsrf),  &            call histdef(nid_ins, "rugs_"//clnsurf(nsrf),  &
170                 "rugosite "//clnsurf(nsrf), "-",   &                 "rugosite "//clnsurf(nsrf), "-",   &
171                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
172                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
173            !XXX            call histdef(nid_ins, "u10m_"//clnsurf(nsrf),  &
174                   "zonal wind 10 m "//clnsurf(nsrf), "m s-1",   &
175                   iim, (jjm + 1), nhori, 1, 1, 1, -99, &
176                   "inst(X)", zsto, zout)
177              call histdef(nid_ins, "v10m_"//clnsurf(nsrf),  &
178                   "meridional wind 10 m "//clnsurf(nsrf), "m s-1",   &
179                   iim, (jjm + 1), nhori, 1, 1, 1, -99, &
180                   "inst(X)", zsto, zout)
181         END DO         END DO
182    
183         CALL histdef(nid_ins, "rugs", "rugosity", "-", &         CALL histdef(nid_ins, "rugs", "rugosity", "-", &
184              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
185              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
186         CALL histdef(nid_ins, "albs", "Surface albedo", "-", &         CALL histdef(nid_ins, "albs", "Surface albedo", "-", &
187              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
188              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
189         CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m", &         CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m", &
190              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
191              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
192         CALL histdef(nid_ins, "s_pblt", "T at Boundary Layer Height",  &         CALL histdef(nid_ins, "s_pblt", "T at Boundary Layer Height",  &
193              "K", &              "K", &
194              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
195              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
196         CALL histdef(nid_ins, "s_lcl", "Condensation level", "m", &         CALL histdef(nid_ins, "s_lcl", "Condensation level", "m", &
197              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
198              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
199         CALL histdef(nid_ins, "s_capCL", "Conv avlbl pot ener for ABL", "J/m2", &         CALL histdef(nid_ins, "s_capCL", "Conv avlbl pot ener for ABL", "J/m2", &
200              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
201              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
202         CALL histdef(nid_ins, "s_oliqCL", "Liq Water in BL", "kg/m2", &         CALL histdef(nid_ins, "s_oliqCL", "Liq Water in BL", "kg/m2", &
203              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
204              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
205         CALL histdef(nid_ins, "s_cteiCL", "Instability criteria (ABL)", "K", &         CALL histdef(nid_ins, "s_cteiCL", "Instability criteria (ABL)", "K", &
206              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
207              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
208         CALL histdef(nid_ins, "s_therm", "Exces du thermique", "K", &         CALL histdef(nid_ins, "s_therm", "Exces du thermique", "K", &
209              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
210              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
211           CALL histdef(nid_ins, "qsurf", "Surface Air humidity", "", &
        CALL histdef(nid_ins, "s_trmb1", "deep_cape(HBTM2)", "J/m2", &  
             iim, (jjm + 1), nhori, 1, 1, 1, -99, &  
             "inst(X)", zsto, zout)  
   
        CALL histdef(nid_ins, "s_trmb2", "inhibition (HBTM2)", "J/m2", &  
             iim, (jjm + 1), nhori, 1, 1, 1, -99, &  
             "inst(X)", zsto, zout)  
   
        CALL histdef(nid_ins, "s_trmb3", "Point Omega (HBTM2)", "m", &  
212              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
213              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
214    
215         if (conv_emanuel) CALL histdef(nid_ins, "ptop", "cloud top pressure", &         if (conv_emanuel) then
216              "Pa", iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout)            CALL histdef(nid_ins, "ptop", "cloud top pressure", &
217                   "Pa", iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout)
218              CALL histdef(nid_ins, "dnwd0", "unsaturated downdraft", &
219                   "kg/m2/s", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
220                   zsto, zout)
221           end if
222    
223         ! Champs 3D:         ! Champs 3D:
224    
225           CALL histdef(nid_ins, "tro3", "ozone mole fraction", "-", &
226                iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", zsto, zout)
227         CALL histdef(nid_ins, "temp", "Temperature", "K", &         CALL histdef(nid_ins, "temp", "Temperature", "K", &
228              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
229              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
230         CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s", &         CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s", &
231              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
232              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
233         CALL histdef(nid_ins, "vitv", "Merid wind", "m/s", &         CALL histdef(nid_ins, "vitv", "Merid wind", "m/s", &
234              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
235              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
236         CALL histdef(nid_ins, "geop", "Geopotential height", "m", &         CALL histdef(nid_ins, "geop", "Geopotential height", "m", &
237              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
238              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
239         CALL histdef(nid_ins, "pres", "Air pressure", "Pa", &         CALL histdef(nid_ins, "pres", "Air pressure", "Pa", &
240              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
241              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
242         CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s", &         CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s", &
243              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
244              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
245         CALL histdef(nid_ins, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s", &         CALL histdef(nid_ins, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s", &
246              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
247              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
   
248         CALL histdef(nid_ins, "zmasse", "column density of air in cell", &         CALL histdef(nid_ins, "zmasse", "column density of air in cell", &
249              "kg m-2", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &              "kg m-2", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
250              zsto, zout)              zsto, zout)
   
251         CALL histdef(nid_ins, "rhum", "Relative humidity", &         CALL histdef(nid_ins, "rhum", "Relative humidity", &
252              "", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &              "", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
253              zsto, zout)              zsto, zout)
254           CALL histdef(nid_ins, "d_t_ec", "kinetic dissipation dT", &
255                "K/s", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
256                zsto, zout)
257           CALL histdef(nid_ins, "dtsw0", "CS SW radiation dT", &
258                "K/s", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
259                zsto, zout)
260           CALL histdef(nid_ins, "dtlw0", "CS LW radiation dT", &
261                "K/s", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
262                zsto, zout)
263    
264           if (ok_newmicro) then
265              CALL histdef(nid_ins, "re", "cloud droplet effective radius", &
266                   "micrometer", iim, jjm + 1, nhori, llm, 1, llm, nvert, &
267                   "inst(X)", zsto, zout)
268              CALL histdef(nid_ins, "fl", &
269                   "denominator of Cloud droplet effective radius", "", iim, &
270                   jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", zsto, zout)
271           end if
272    
273         DO it = 1, nqmx - 2         DO it = 1, nqmx - 2
274            ! champ 2D            ! champ 2D
# Line 300  contains Line 278  contains
278            CALL histdef(nid_ins, "fl"//tname(iq), "Flux "//ttext(iq), &            CALL histdef(nid_ins, "fl"//tname(iq), "Flux "//ttext(iq), &
279                 "U/m2/s", iim, jjm+1, nhori, llm, 1, llm, nvert, &                 "U/m2/s", iim, jjm+1, nhori, llm, 1, llm, nvert, &
280                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
   
           !---Ajout Olivia  
281            CALL histdef(nid_ins, "d_tr_th_"//tname(iq), &            CALL histdef(nid_ins, "d_tr_th_"//tname(iq), &
282                 "tendance thermique"// ttext(iq), "?", &                 "tendance thermique"// ttext(iq), "?", &
283                 iim, jjm+1, nhori, llm, 1, llm, nvert, &                 iim, jjm+1, nhori, llm, 1, llm, nvert, &
# Line 314  contains Line 290  contains
290                 "tendance couche limite"// ttext(iq), "?", &                 "tendance couche limite"// ttext(iq), "?", &
291                 iim, jjm+1, nhori, llm, 1, llm, nvert, &                 iim, jjm+1, nhori, llm, 1, llm, nvert, &
292                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
           !---fin Olivia  
   
293         ENDDO         ENDDO
294    
295         CALL histend(nid_ins)         CALL histend(nid_ins)
296      ENDIF      ENDIF test_ok_instan
297    
298    end subroutine ini_histins    end subroutine ini_histins
299    
Legend:
Removed from v.191  
changed lines
  Added in v.254
  ViewVC Help
Powered by ViewVC 1.1.21