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

Diff of /trunk/Sources/phylmd/ini_histins.f

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

trunk/phylmd/ini_histins.f revision 82 by guez, Wed Mar 5 14:57:53 2014 UTC trunk/Sources/phylmd/ini_histins.f revision 186 by guez, Mon Mar 21 15:36:26 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_fi_ecrit_m, only: gr_fi_ecrit
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)         CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlon, zx_lon)
46         DO i = 1, iim         DO i = 1, iim
# 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.82  
changed lines
  Added in v.186
  ViewVC Help
Powered by ViewVC 1.1.21