/[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

revision 215 by guez, Tue Mar 28 12:46:28 2017 UTC revision 227 by guez, Thu Nov 2 15:47:03 2017 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, version 1.2, 2005/05/25 13:10:09      ! From phylmd/ini_histins.h, version 1.2, 2005/05/25 13:10:09
12    
# Line 26  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
# Line 47  contains Line 48  contains
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                
51         CALL histdef(nid_ins, "phis", "Surface geop. height", "-", &         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 169  contains Line 170  contains
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              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", "-", &
# Line 208  contains Line 217  contains
217         CALL histdef(nid_ins, "s_trmb3", "Point Omega (HBTM2)", "m", &         CALL histdef(nid_ins, "s_trmb3", "Point Omega (HBTM2)", "m", &
218              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
219              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
220           CALL histdef(nid_ins, "qsurf", "Surface Air humidity", "", &
221                iim, (jjm + 1), nhori, 1, 1, 1, -99, &
222                "inst(X)", zsto, zout)
223    
224         if (conv_emanuel) then         if (conv_emanuel) then
225            CALL histdef(nid_ins, "ptop", "cloud top pressure", &            CALL histdef(nid_ins, "ptop", "cloud top pressure", &
# Line 258  contains Line 270  contains
270              "K/s", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &              "K/s", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
271              zsto, zout)              zsto, zout)
272    
273           if (ok_newmicro) then
274              CALL histdef(nid_ins, "re", "cloud droplet effective radius", &
275                   "micrometer", iim, jjm + 1, nhori, llm, 1, llm, nvert, &
276                   "inst(X)", zsto, zout)
277              CALL histdef(nid_ins, "fl", &
278                   "denominator of Cloud droplet effective radius", "", iim, &
279                   jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", zsto, zout)
280           end if
281    
282         DO it = 1, nqmx - 2         DO it = 1, nqmx - 2
283            ! champ 2D            ! champ 2D
284            iq=it+2            iq=it+2
Legend:
Removed from v.215  
changed lines
  Added in v.227
  ViewVC Help
Powered by ViewVC 1.1.21