/[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 213 by guez, Mon Feb 27 15:44:55 2017 UTC trunk/phylmd/ini_histins.f revision 301 by guez, Thu Aug 2 17:23:07 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(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: day_ref, annee_ref, rlatu, rlonv
19      USE histbeg_totreg_m, ONLY : histbeg_totreg      USE histbeg_totreg_m, ONLY : histbeg_totreg
# Line 25  contains Line 26  contains
26      use phyetat0_m, only: itau_phy      use phyetat0_m, only: itau_phy
27      USE ymds2ju_m, only: ymds2ju      USE ymds2ju_m, only: ymds2ju
28    
29      REAL, intent(in):: dtime ! pas temporel de la physique (s)      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 35  contains Line 36  contains
36    
37      print *, 'Call sequence information: ini_histins'      print *, 'Call sequence information: ini_histins'
38    
39      IF (ok_instan) THEN      test_ok_instan: IF (ok_instan) THEN
40         zsto = dtime * ecrit_ins         zsto = dtphys * ecrit_ins
41         zout = dtime * ecrit_ins         zout = dtphys * 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, dtphys, nhori, nid_ins)
46         print *, 'itau_phy = ', itau_phy         print *, 'itau_phy = ', itau_phy
47         print *, "zjulian = ", zjulian         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                
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 129  contains Line 130  contains
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           CALL histdef(nid_ins, "zxfqcalving", "ice calving", "kg m-2 s-1", &
136                iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout)
137           CALL histdef(nid_ins, "run_off_lic", "land ice melt to ocean", &
138                "kg m-2 s-1", iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, &
139                zout)
140    
141         DO nsrf = 1, nbsrf         DO nsrf = 1, nbsrf
142            call histdef(nid_ins, "pourc_"//clnsurf(nsrf),  &            call histdef(nid_ins, "pourc_"//clnsurf(nsrf),  &
# Line 167  contains Line 175  contains
175                 "rugosite "//clnsurf(nsrf), "-",   &                 "rugosite "//clnsurf(nsrf), "-",   &
176                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
177                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
178              call histdef(nid_ins, "u10m_"//clnsurf(nsrf),  &
179                   "zonal wind 10 m "//clnsurf(nsrf), "m s-1",   &
180                   iim, (jjm + 1), nhori, 1, 1, 1, -99, &
181                   "inst(X)", zsto, zout)
182              call histdef(nid_ins, "v10m_"//clnsurf(nsrf),  &
183                   "meridional wind 10 m "//clnsurf(nsrf), "m s-1",   &
184                   iim, (jjm + 1), nhori, 1, 1, 1, -99, &
185                   "inst(X)", zsto, zout)
186         END DO         END DO
187    
188         CALL histdef(nid_ins, "rugs", "rugosity", "-", &         CALL histdef(nid_ins, "rugs", "rugosity", "-", &
# Line 185  contains Line 201  contains
201         CALL histdef(nid_ins, "s_lcl", "Condensation level", "m", &         CALL histdef(nid_ins, "s_lcl", "Condensation level", "m", &
202              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
203              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
204         CALL histdef(nid_ins, "s_capCL", "Conv avlbl pot ener for ABL", "J/m2", &         CALL histdef(nid_ins, "s_capCL", &
205              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              "Convective available potential energy for atmospheric boundary " &
206                // "layer", "J/m2", iim, (jjm + 1), nhori, 1, 1, 1, -99, &
207              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
208         CALL histdef(nid_ins, "s_oliqCL", "Liq Water in BL", "kg/m2", &         CALL histdef(nid_ins, "s_oliqCL", "Liq Water in BL", "kg/m2", &
209              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
# Line 197  contains Line 214  contains
214         CALL histdef(nid_ins, "s_therm", "Exces du thermique", "K", &         CALL histdef(nid_ins, "s_therm", "Exces du thermique", "K", &
215              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
216              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
217         CALL histdef(nid_ins, "s_trmb1", "deep_cape(HBTM2)", "J/m2", &         CALL histdef(nid_ins, "qsurf", "Surface Air humidity", "", &
             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", &  
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    
# Line 256  contains Line 267  contains
267              "K/s", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &              "K/s", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
268              zsto, zout)              zsto, zout)
269    
270           if (ok_newmicro) then
271              CALL histdef(nid_ins, "re", "cloud droplet effective radius", &
272                   "micrometer", iim, jjm + 1, nhori, llm, 1, llm, nvert, &
273                   "inst(X)", zsto, zout)
274              CALL histdef(nid_ins, "fl", &
275                   "denominator of Cloud droplet effective radius", "", iim, &
276                   jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", zsto, zout)
277           end if
278    
279         DO it = 1, nqmx - 2         DO it = 1, nqmx - 2
280            ! champ 2D            ! champ 2D
281            iq=it+2            iq=it+2
# Line 279  contains Line 299  contains
299         ENDDO         ENDDO
300    
301         CALL histend(nid_ins)         CALL histend(nid_ins)
302      ENDIF      ENDIF test_ok_instan
303    
304    end subroutine ini_histins    end subroutine ini_histins
305    
Legend:
Removed from v.213  
changed lines
  Added in v.301
  ViewVC Help
Powered by ViewVC 1.1.21