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

Diff of /trunk/phylmd/ini_histins.f90

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

trunk/Sources/phylmd/ini_histins.f revision 225 by guez, Mon Oct 16 12:35:41 2017 UTC trunk/phylmd/ini_histins.f revision 321 by guez, Tue Dec 11 22:48:09 2018 UTC
# Line 6  module ini_histins_m Line 6  module ini_histins_m
6    
7  contains  contains
8    
9    subroutine ini_histins(dtime, ok_newmicro)    subroutine ini_histins(ok_newmicro)
10    
11      ! From phylmd/ini_histins.h, version 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 comconst, only: dtphys
16        use dimensions, only: iim, jjm, llm, nqmx
17      use disvert_m, only: presnivs      use disvert_m, only: presnivs
18      use dynetat0_m, only: day_ref, annee_ref, rlatu, rlonv      use dynetat0_m, only: rlatu, rlonv
19        use dynetat0_chosen_m, only: day_ref, annee_ref
20      USE histbeg_totreg_m, ONLY : histbeg_totreg      USE histbeg_totreg_m, ONLY : histbeg_totreg
21      USE histdef_m, ONLY : histdef      USE histdef_m, ONLY : histdef
22      USE histend_m, ONLY : histend      USE histend_m, ONLY : histend
23      USE histvert_m, ONLY : histvert      USE histvert_m, ONLY : histvert
24      use indicesol, only: nbsrf, clnsurf      use indicesol, only: nbsrf, clnsurf
25      use iniadvtrac_m, only: tname, ttext      use infotrac_init_m, only: tname, ttext
26      use nr_util, only: pi      use nr_util, only: pi
27      use phyetat0_m, only: itau_phy      use phyetat0_m, only: itau_phy
28      USE ymds2ju_m, only: ymds2ju      USE ymds2ju_m, only: ymds2ju
29    
     REAL, intent(in):: dtime ! pas temporel de la physique (s)  
30      logical, intent(in):: ok_newmicro      logical, intent(in):: ok_newmicro
31            
32      ! Local:      ! Local:
# Line 37  contains Line 38  contains
38      print *, 'Call sequence information: ini_histins'      print *, 'Call sequence information: ini_histins'
39    
40      test_ok_instan: IF (ok_instan) THEN      test_ok_instan: IF (ok_instan) THEN
41         zsto = dtime * ecrit_ins         zsto = dtphys * ecrit_ins
42         zout = dtime * ecrit_ins         zout = dtphys * ecrit_ins
43         CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)         CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
44         CALL histbeg_totreg("histins", rlonv(:iim) / pi * 180., &         CALL histbeg_totreg("histins", rlonv(:iim) / pi * 180., &
45              rlatu / pi * 180., 1, iim, &              rlatu / pi * 180., 1, iim, &
46              1, jjm + 1, itau_phy, zjulian, dtime, nhori, nid_ins)              1, jjm + 1, itau_phy, zjulian, dtphys, nhori, nid_ins)
47         print *, 'itau_phy = ', itau_phy         print *, 'itau_phy = ', itau_phy
48         print *, "zjulian = ", zjulian         print *, "zjulian = ", zjulian
49         CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", &         CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", &
50              presnivs/100., nvert)              presnivs/100., nvert)
51                
52         CALL histdef(nid_ins, "phis", "Surface geop. height", "-", &          ! Once:
53              iim, (jjm + 1), nhori, 1, 1, 1, -99, &         CALL histdef(nid_ins, "phis", "surface geopotential", "m2 s-2", &
54                iim, jjm + 1, nhori, 1, 1, 1, -99, &
55              "once", zsto, zout)              "once", zsto, zout)
56         CALL histdef(nid_ins, "aire", "Grid area", "-", &         CALL histdef(nid_ins, "aire", "Grid area", "-", &
57              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, jjm + 1, nhori, 1, 1, 1, -99, &
58              "once", zsto, zout)              "once", zsto, zout)
59    
60         ! Champs 2D:         ! Champs 2D:
61    
62         CALL histdef(nid_ins, "tsol", "Surface Temperature", "K", &         CALL histdef(nid_ins, "tsol", "Surface Temperature", "K", &
63              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, jjm + 1, nhori, 1, 1, 1, -99,  &
64              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
65         CALL histdef(nid_ins, "t2m", "Temperature 2m", "K", &         CALL histdef(nid_ins, "t2m", "Temperature 2m", "K", &
66              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, jjm + 1, nhori, 1, 1, 1, -99, &
67              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
68         CALL histdef(nid_ins, "q2m", "Specific humidity 2m", "Kg/Kg", &         CALL histdef(nid_ins, "q2m", "Specific humidity 2m", "Kg/Kg", &
69              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, jjm + 1, nhori, 1, 1, 1, -99, &
70              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
71         CALL histdef(nid_ins, "u10m", "Vent zonal 10m", "m/s", &         CALL histdef(nid_ins, "u10m", "Vent zonal 10m", "m/s", &
72              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, jjm + 1, nhori, 1, 1, 1, -99, &
73              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
74         CALL histdef(nid_ins, "v10m", "Vent meridien 10m", "m/s", &         CALL histdef(nid_ins, "v10m", "Vent meridien 10m", "m/s", &
75              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, jjm + 1, nhori, 1, 1, 1, -99, &
76              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
77         CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa", &         CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa", &
78              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, jjm + 1, nhori, 1, 1, 1, -99, &
79              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
80         CALL histdef(nid_ins, "plul", "Large-scale Precip.", "mm/day", &         CALL histdef(nid_ins, "plul", "Large-scale Precip.", "mm/day", &
81              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, jjm + 1, nhori, 1, 1, 1, -99, &
82              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
83         CALL histdef(nid_ins, "pluc", "Convective Precip.", "mm/day", &         CALL histdef(nid_ins, "pluc", "Convective Precip.", "mm/day", &
84              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, jjm + 1, nhori, 1, 1, 1, -99, &
85              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
86         CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-", &         CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-", &
87              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, jjm + 1, nhori, 1, 1, 1, -99,  &
88              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
89         CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-", &         CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-", &
90              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, jjm + 1, nhori, 1, 1, 1, -99,  &
91              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
92         CALL histdef(nid_ins, "precip", "Precipitation Totale liq+sol",  &         CALL histdef(nid_ins, "precip", "Precipitation Totale liq+sol",  &
93              "kg/(s*m2)", &              "kg/(s*m2)", &
94              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, jjm + 1, nhori, 1, 1, 1, -99,  &
95              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
96         CALL histdef(nid_ins, "snow", "Snow fall", "kg/(s*m2)", &         CALL histdef(nid_ins, "snow", "Snow fall", "kg/(s*m2)", &
97              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, jjm + 1, nhori, 1, 1, 1, -99,  &
98              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
99         CALL histdef(nid_ins, "topl", "OLR", "W/m2", &         CALL histdef(nid_ins, "topl", "OLR", "W/m2", &
100              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, jjm + 1, nhori, 1, 1, 1, -99, &
101              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
102         CALL histdef(nid_ins, "evap", "Evaporation", "kg/(s*m2)", &         CALL histdef(nid_ins, "evap", "Evaporation", "kg/(s*m2)", &
103              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, jjm + 1, nhori, 1, 1, 1, -99,  &
104              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
105         CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2", &         CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2", &
106              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, jjm + 1, nhori, 1, 1, 1, -99,  &
             "inst(X)", zsto, zout)  
        CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2", &  
             iim, (jjm + 1), nhori, 1, 1, 1, -99,  &  
107              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
108           CALL histdef(nid_ins, "rls", "surface net downward longwave flux", &
109                "W/m2", iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout)
110         CALL histdef(nid_ins, "solldown", "Down. IR rad. at surface",  &         CALL histdef(nid_ins, "solldown", "Down. IR rad. at surface",  &
111              "W/m2", iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              "W/m2", iim, jjm + 1, nhori, 1, 1, 1, -99,  &
112              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
113         CALL histdef(nid_ins, "bils", "Surf. total heat flux", "W/m2", &         CALL histdef(nid_ins, "bils", "Surf. total heat flux", "W/m2", &
114              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, jjm + 1, nhori, 1, 1, 1, -99,  &
115              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
116         CALL histdef(nid_ins, "sens", "Sensible heat flux", "W/m2", &         CALL histdef(nid_ins, "sens", "Sensible heat flux", "W/m2", &
117              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, jjm + 1, nhori, 1, 1, 1, -99,  &
118              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
119         CALL histdef(nid_ins, "fder", "Heat flux derivation", "W/m2", &         CALL histdef(nid_ins, "fder", "Heat flux derivation", "W/m2", &
120              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, jjm + 1, nhori, 1, 1, 1, -99,  &
121              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
122         CALL histdef(nid_ins, "dtsvdfo", "Boundary-layer dTs(o)", "K/s", &         CALL histdef(nid_ins, "dtsvdfo", "Boundary-layer dTs(o)", "K/s", &
123              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, jjm + 1, nhori, 1, 1, 1, -99,  &
124              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
125         CALL histdef(nid_ins, "dtsvdft", "Boundary-layer dTs(t)", "K/s", &         CALL histdef(nid_ins, "dtsvdft", "Boundary-layer dTs(t)", "K/s", &
126              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, jjm + 1, nhori, 1, 1, 1, -99,  &
127              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
128         CALL histdef(nid_ins, "dtsvdfg", "Boundary-layer dTs(g)", "K/s", &         CALL histdef(nid_ins, "dtsvdfg", "Boundary-layer dTs(g)", "K/s", &
129              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, jjm + 1, nhori, 1, 1, 1, -99,  &
130              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
131         CALL histdef(nid_ins, "dtsvdfi", "Boundary-layer dTs(g)", "K/s", &         CALL histdef(nid_ins, "dtsvdfi", "Boundary-layer dTs(g)", "K/s", &
132              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, jjm + 1, nhori, 1, 1, 1, -99,  &
133              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
134         CALL histdef(nid_ins, "msnow", "surface snow amount", "kg/m2", &         CALL histdef(nid_ins, "msnow", "surface snow amount", "kg/m2", &
135              iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout)              iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout)
136           CALL histdef(nid_ins, "zxfqcalving", "ice calving", "kg m-2 s-1", &
137                iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout)
138           CALL histdef(nid_ins, "run_off_lic", "land ice melt to ocean", &
139                "kg m-2 s-1", iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, &
140                zout)
141           CALL histdef(nid_ins, "flat", "latent heat flux", "W m-2", iim, &
142                jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout)
143    
144         DO nsrf = 1, nbsrf         DO nsrf = 1, nbsrf
           call histdef(nid_ins, "pourc_"//clnsurf(nsrf),  &  
                "% "//clnsurf(nsrf), "%",   &  
                iim, (jjm + 1), nhori, 1, 1, 1, -99, &  
                "inst(X)", zsto, zout)  
145            call histdef(nid_ins, "fract_"//clnsurf(nsrf),  &            call histdef(nid_ins, "fract_"//clnsurf(nsrf),  &
146                 "Fraction "//clnsurf(nsrf), "1",   &                 "Fraction "//clnsurf(nsrf), "1",   &
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, "sens_"//clnsurf(nsrf),  &            call histdef(nid_ins, "sens_"//clnsurf(nsrf),  &
150                 "Sensible heat flux "//clnsurf(nsrf), "W/m2",   &                 "Sensible heat flux "//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, "tsol_"//clnsurf(nsrf),  &            call histdef(nid_ins, "tsol_"//clnsurf(nsrf),  &
154                 "Surface Temperature"//clnsurf(nsrf), "W/m2",   &                 "Surface Temperature"//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, "lat_"//clnsurf(nsrf),  &            call histdef(nid_ins, "lat_"//clnsurf(nsrf),  &
158                 "Latent heat flux "//clnsurf(nsrf), "W/m2",   &                 "Latent heat flux "//clnsurf(nsrf), "W/m2",   &
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, "taux_"//clnsurf(nsrf),  &            call histdef(nid_ins, "taux_"//clnsurf(nsrf),  &
162                 "Zonal wind stress"//clnsurf(nsrf), "Pa", &                 "Zonal wind 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, "tauy_"//clnsurf(nsrf),  &            call histdef(nid_ins, "tauy_"//clnsurf(nsrf),  &
166                 "Meridional xind stress "//clnsurf(nsrf), "Pa",   &                 "Meridional xind stress "//clnsurf(nsrf), "Pa",   &
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, "albe_"//clnsurf(nsrf),  &            call histdef(nid_ins, "albe_"//clnsurf(nsrf),  &
170                 "Albedo "//clnsurf(nsrf), "-",   &                 "Albedo "//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            call histdef(nid_ins, "rugs_"//clnsurf(nsrf),  &            call histdef(nid_ins, "rugs_"//clnsurf(nsrf),  &
174                 "rugosite "//clnsurf(nsrf), "-",   &                 "rugosite "//clnsurf(nsrf), "-",   &
175                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &                 iim, jjm + 1, nhori, 1, 1, 1, -99, &
176                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
177            call histdef(nid_ins, "u10m_"//clnsurf(nsrf),  &            call histdef(nid_ins, "u10m_"//clnsurf(nsrf),  &
178                 "zonal wind 10 m "//clnsurf(nsrf), "m s-1",   &                 "zonal wind 10 m "//clnsurf(nsrf), "m s-1",   &
179                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &                 iim, jjm + 1, nhori, 1, 1, 1, -99, &
180                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
181            call histdef(nid_ins, "v10m_"//clnsurf(nsrf),  &            call histdef(nid_ins, "v10m_"//clnsurf(nsrf),  &
182                 "meridional wind 10 m "//clnsurf(nsrf), "m s-1",   &                 "meridional wind 10 m "//clnsurf(nsrf), "m s-1",   &
183                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &                 iim, jjm + 1, nhori, 1, 1, 1, -99, &
184                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
185         END DO         END DO
186    
187         CALL histdef(nid_ins, "rugs", "rugosity", "-", &         CALL histdef(nid_ins, "rugs", "rugosity", "-", &
188              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, jjm + 1, nhori, 1, 1, 1, -99,  &
189              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
190         CALL histdef(nid_ins, "albs", "Surface albedo", "-", &         CALL histdef(nid_ins, "albs", "Surface albedo", "-", &
191              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, jjm + 1, nhori, 1, 1, 1, -99,  &
192              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
193         CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m", &         CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m", &
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_pblt", "T at Boundary Layer Height",  &         CALL histdef(nid_ins, "s_pblt", "T at Boundary Layer Height",  &
197              "K", &              "K", &
198              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, jjm + 1, nhori, 1, 1, 1, -99, &
199              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
200         CALL histdef(nid_ins, "s_lcl", "Condensation level", "m", &         CALL histdef(nid_ins, "s_lcl", "Condensation level", "m", &
201              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, jjm + 1, nhori, 1, 1, 1, -99, &
202              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
203         CALL histdef(nid_ins, "s_capCL", "Conv avlbl pot ener for ABL", "J/m2", &         CALL histdef(nid_ins, "s_capCL", &
204              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              "Convective available potential energy for atmospheric boundary " &
205                // "layer", "J/m2", iim, jjm + 1, nhori, 1, 1, 1, -99, &
206              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
207         CALL histdef(nid_ins, "s_oliqCL", "Liq Water in BL", "kg/m2", &         CALL histdef(nid_ins, "s_oliqCL", "Liq Water in BL", "kg/m2", &
208              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, jjm + 1, nhori, 1, 1, 1, -99, &
209              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
210         CALL histdef(nid_ins, "s_cteiCL", "Instability criteria (ABL)", "K", &         CALL histdef(nid_ins, "s_cteiCL", "Instability criteria (ABL)", "K", &
211              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, jjm + 1, nhori, 1, 1, 1, -99, &
212              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
213         CALL histdef(nid_ins, "s_therm", "Exces du thermique", "K", &         CALL histdef(nid_ins, "s_therm", "Exces du thermique", "K", &
214              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, jjm + 1, nhori, 1, 1, 1, -99, &
             "inst(X)", zsto, zout)  
        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", &  
             iim, (jjm + 1), nhori, 1, 1, 1, -99, &  
215              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
216         CALL histdef(nid_ins, "qsurf", "Surface Air humidity", "", &         CALL histdef(nid_ins, "qsurf", "Surface Air humidity", "", &
217              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, jjm + 1, nhori, 1, 1, 1, -99, &
218              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
219    
220         if (conv_emanuel) then         if (conv_emanuel) then
# Line 234  contains Line 230  contains
230         CALL histdef(nid_ins, "tro3", "ozone mole fraction", "-", &         CALL histdef(nid_ins, "tro3", "ozone mole fraction", "-", &
231              iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", zsto, zout)              iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", zsto, zout)
232         CALL histdef(nid_ins, "temp", "Temperature", "K", &         CALL histdef(nid_ins, "temp", "Temperature", "K", &
233              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, jjm + 1, nhori, llm, 1, llm, nvert, &
234              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
235         CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s", &         CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s", &
236              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, jjm + 1, nhori, llm, 1, llm, nvert, &
237              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
238         CALL histdef(nid_ins, "vitv", "Merid wind", "m/s", &         CALL histdef(nid_ins, "vitv", "Merid wind", "m/s", &
239              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, jjm + 1, nhori, llm, 1, llm, nvert, &
240              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
241         CALL histdef(nid_ins, "geop", "Geopotential height", "m", &         CALL histdef(nid_ins, "geop", "Geopotential height", "m", &
242              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, jjm + 1, nhori, llm, 1, llm, nvert, &
243              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
244         CALL histdef(nid_ins, "pres", "Air pressure", "Pa", &         CALL histdef(nid_ins, "pres", "Air pressure", "Pa", &
245              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, jjm + 1, nhori, llm, 1, llm, nvert, &
246              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
247         CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s", &         CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s", &
248              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, jjm + 1, nhori, llm, 1, llm, nvert, &
249              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
250         CALL histdef(nid_ins, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s", &         CALL histdef(nid_ins, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s", &
251              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, jjm + 1, nhori, llm, 1, llm, nvert, &
252              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
253         CALL histdef(nid_ins, "zmasse", "column density of air in cell", &         CALL histdef(nid_ins, "zmasse", "column density of air in cell", &
254              "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)", &
Legend:
Removed from v.225  
changed lines
  Added in v.321
  ViewVC Help
Powered by ViewVC 1.1.21