/[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/libf/phylmd/ini_histins.f90 revision 66 by guez, Thu Sep 20 13:00:41 2012 UTC trunk/Sources/phylmd/ini_histins.f revision 189 by guez, Tue Mar 29 15:20:23 2016 UTC
# Line 4  module ini_histins_m Line 4  module ini_histins_m
4    
5  contains  contains
6    
7    subroutine ini_histins(dtime, ok_instan, nid_ins)    subroutine ini_histins(dtime, ok_instan, nid_ins, itau_phy)
8    
9      ! From phylmd/ini_histins.h, v 1.2 2005/05/25 13:10:09      ! From phylmd/ini_histins.h, v 1.2 2005/05/25 13:10:09
10    
     use dimens_m, only: iim, jjm, llm  
     use dimphy, only: klon  
     use temps, only: itau_phy, day_ref, annee_ref  
11      use clesphys, only: ecrit_ins      use clesphys, only: ecrit_ins
12      use indicesol, only: nbsrf, clnsurf      use clesphys2, only: conv_emanuel
13      USE calendar, only: ymds2ju      use dimens_m, only: iim, jjm, llm, nqmx
14        use dimphy, only: klon
15        use disvert_m, only: presnivs
16        use dynetat0_m, only: day_ref, annee_ref
17        use gr_phy_write_m, only: gr_phy_write
18      USE histbeg_totreg_m, ONLY : histbeg_totreg      USE histbeg_totreg_m, ONLY : histbeg_totreg
19      USE histdef_m, ONLY : histdef      USE histdef_m, ONLY : histdef
20      USE histend_m, ONLY : histend      USE histend_m, ONLY : histend
21      USE histvert_m, ONLY : histvert      USE histvert_m, ONLY : histvert
22        use indicesol, only: nbsrf, clnsurf
23        use iniadvtrac_m, only: tname, ttext
24      use phyetat0_m, only: rlon, rlat      use phyetat0_m, only: rlon, rlat
25      use disvert_m, only: presnivs      USE ymds2ju_m, only: ymds2ju
26    
27      REAL, intent(in):: dtime ! pas temporel de la physique (s)      REAL, intent(in):: dtime ! pas temporel de la physique (s)
28      logical, intent(in):: ok_instan      logical, intent(in):: ok_instan
29      integer, intent(out):: nid_ins      integer, intent(out):: nid_ins
30        integer, intent(in):: itau_phy
31    
32        ! Local:
33      REAL zx_lon(iim, jjm + 1), zx_lat(iim, jjm + 1)      REAL zx_lon(iim, jjm + 1), zx_lat(iim, jjm + 1)
34      real zjulian, zsto, zout      real zjulian, zsto, zout
35      integer i, nhori, nvert, idayref, nsrf      integer i, nhori, nvert, nsrf, iq, it
36    
37      !-------------------------------------------------------------------      !-------------------------------------------------------------------
38    
39      IF (ok_instan) THEN      IF (ok_instan) THEN
   
40         zsto = dtime * ecrit_ins         zsto = dtime * ecrit_ins
41         zout = dtime * ecrit_ins         zout = dtime * ecrit_ins
42    
43         idayref = day_ref         CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
        CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)  
44    
45         CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlon, zx_lon)         zx_lon = gr_phy_write(rlon)
46         DO i = 1, iim         DO i = 1, iim
47            zx_lon(i, 1) = rlon(i+1)            zx_lon(i, 1) = rlon(i+1)
48            zx_lon(i, (jjm + 1)) = rlon(i+1)            zx_lon(i, (jjm + 1)) = rlon(i+1)
49         ENDDO         ENDDO
50         CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlat, zx_lat)         zx_lat = gr_phy_write(rlat)
51         CALL histbeg_totreg("histins", zx_lon(:, 1), zx_lat(1, :), 1, iim, 1, &         CALL histbeg_totreg("histins", zx_lon(:, 1), zx_lat(1, :), 1, iim, 1, &
52              jjm + 1, itau_phy, zjulian, dtime, nhori, nid_ins)              jjm + 1, itau_phy, zjulian, dtime, nhori, nid_ins)
53         write(*, *)'Inst ', itau_phy, zjulian         write(*, *)'Inst ', itau_phy, zjulian
54         CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", &         CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", &
55              llm, presnivs/100., nvert)              presnivs/100., nvert)
56    
57         CALL histdef(nid_ins, "phis", "Surface geop. height", "-", &         CALL histdef(nid_ins, "phis", "Surface geop. height", "-", &
58              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
# Line 217  contains Line 220  contains
220         CALL histdef(nid_ins, "albs", "Surface albedo", "-", &         CALL histdef(nid_ins, "albs", "Surface albedo", "-", &
221              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
222              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
        CALL histdef(nid_ins, "albslw", "Surface albedo LW", "-", &  
             iim, (jjm + 1), nhori, 1, 1, 1, -99,  &  
             "inst(X)", zsto, zout)  
223    
        !IM cf. AM 081204 BEG  
        ! HBTM2  
224         CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m", &         CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m", &
225              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
226              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
# Line 264  contains Line 262  contains
262              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
263              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
264    
265         !IM cf. AM 081204 END         if (conv_emanuel) CALL histdef(nid_ins, "ptop", "cloud top pressure", &
266                "Pa", iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout)
267    
268         ! Champs 3D:         ! Champs 3D:
269    
# Line 296  contains Line 295  contains
295              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
296              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
297    
298           CALL histdef(nid_ins, "zmasse", "column density of air in cell", &
299                "kg m-2", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
300                zsto, zout)
301    
302           CALL histdef(nid_ins, "rhum", "Relative humidity", &
303                "", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
304                zsto, zout)
305    
306           DO it = 1, nqmx - 2
307              ! champ 2D
308              iq=it+2
309              CALL histdef(nid_ins, tname(iq), ttext(iq), "U/kga", iim, jjm+1, &
310                   nhori, llm, 1, llm, nvert, "inst(X)", zsto, zout)
311              CALL histdef(nid_ins, "fl"//tname(iq), "Flux "//ttext(iq), &
312                   "U/m2/s", iim, jjm+1, nhori, llm, 1, llm, nvert, &
313                   "inst(X)", zsto, zout)
314    
315              !---Ajout Olivia
316              CALL histdef(nid_ins, "d_tr_th_"//tname(iq), &
317                   "tendance thermique"// ttext(iq), "?", &
318                   iim, jjm+1, nhori, llm, 1, llm, nvert, &
319                   "inst(X)", zsto, zout)
320              CALL histdef(nid_ins, "d_tr_cv_"//tname(iq), &
321                   "tendance convection"// ttext(iq), "?", &
322                   iim, jjm+1, nhori, llm, 1, llm, nvert, &
323                   "inst(X)", zsto, zout)
324              CALL histdef(nid_ins, "d_tr_cl_"//tname(iq), &
325                   "tendance couche limite"// ttext(iq), "?", &
326                   iim, jjm+1, nhori, llm, 1, llm, nvert, &
327                   "inst(X)", zsto, zout)
328              !---fin Olivia
329    
330           ENDDO
331    
332         CALL histend(nid_ins)         CALL histend(nid_ins)
333      ENDIF      ENDIF
334    
Legend:
Removed from v.66  
changed lines
  Added in v.189
  ViewVC Help
Powered by ViewVC 1.1.21