/[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 159 by guez, Tue Jul 21 15:29:52 2015 UTC
# Line 8  contains Line 8  contains
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 dimens_m, only: iim, jjm, llm, nqmx
13      USE calendar, only: ymds2ju      use dimphy, only: klon
14        use disvert_m, only: presnivs
15        use dynetat0_m, only: day_ref, annee_ref
16      USE histbeg_totreg_m, ONLY : histbeg_totreg      USE histbeg_totreg_m, ONLY : histbeg_totreg
17      USE histdef_m, ONLY : histdef      USE histdef_m, ONLY : histdef
18      USE histend_m, ONLY : histend      USE histend_m, ONLY : histend
19      USE histvert_m, ONLY : histvert      USE histvert_m, ONLY : histvert
20        use indicesol, only: nbsrf, clnsurf
21        use iniadvtrac_m, only: tname, ttext
22      use phyetat0_m, only: rlon, rlat      use phyetat0_m, only: rlon, rlat
23      use disvert_m, only: presnivs      use temps, only: itau_phy
24        USE ymds2ju_m, only: ymds2ju
25    
26      REAL, intent(in):: dtime ! pas temporel de la physique (s)      REAL, intent(in):: dtime ! pas temporel de la physique (s)
27      logical, intent(in):: ok_instan      logical, intent(in):: ok_instan
# Line 27  contains Line 29  contains
29    
30      REAL zx_lon(iim, jjm + 1), zx_lat(iim, jjm + 1)      REAL zx_lon(iim, jjm + 1), zx_lat(iim, jjm + 1)
31      real zjulian, zsto, zout      real zjulian, zsto, zout
32      integer i, nhori, nvert, idayref, nsrf      integer i, nhori, nvert, nsrf, iq, it
33    
34      !-------------------------------------------------------------------      !-------------------------------------------------------------------
35    
36      IF (ok_instan) THEN      IF (ok_instan) THEN
   
37         zsto = dtime * ecrit_ins         zsto = dtime * ecrit_ins
38         zout = dtime * ecrit_ins         zout = dtime * ecrit_ins
39    
40         idayref = day_ref         CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
        CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)  
41    
42         CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlon, zx_lon)         CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlon, zx_lon)
43         DO i = 1, iim         DO i = 1, iim
# Line 217  contains Line 217  contains
217         CALL histdef(nid_ins, "albs", "Surface albedo", "-", &         CALL histdef(nid_ins, "albs", "Surface albedo", "-", &
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)
        CALL histdef(nid_ins, "albslw", "Surface albedo LW", "-", &  
             iim, (jjm + 1), nhori, 1, 1, 1, -99,  &  
             "inst(X)", zsto, zout)  
220    
        !IM cf. AM 081204 BEG  
        ! HBTM2  
221         CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m", &         CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m", &
222              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
223              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
# Line 296  contains Line 291  contains
291              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
292              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
293    
294           CALL histdef(nid_ins, "zmasse", "column density of air in cell", &
295                "kg m-2", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
296                zsto, zout)
297    
298           DO it = 1, nqmx - 2
299              ! champ 2D
300              iq=it+2
301              CALL histdef(nid_ins, tname(iq), ttext(iq), "U/kga", iim, jjm+1, &
302                   nhori, llm, 1, llm, nvert, "inst(X)", zsto, zout)
303              CALL histdef(nid_ins, "fl"//tname(iq), "Flux "//ttext(iq), &
304                   "U/m2/s", iim, jjm+1, nhori, llm, 1, llm, nvert, &
305                   "inst(X)", zsto, zout)
306    
307              !---Ajout Olivia
308              CALL histdef(nid_ins, "d_tr_th_"//tname(iq), &
309                   "tendance thermique"// ttext(iq), "?", &
310                   iim, jjm+1, nhori, llm, 1, llm, nvert, &
311                   "inst(X)", zsto, zout)
312              CALL histdef(nid_ins, "d_tr_cv_"//tname(iq), &
313                   "tendance convection"// ttext(iq), "?", &
314                   iim, jjm+1, nhori, llm, 1, llm, nvert, &
315                   "inst(X)", zsto, zout)
316              CALL histdef(nid_ins, "d_tr_cl_"//tname(iq), &
317                   "tendance couche limite"// ttext(iq), "?", &
318                   iim, jjm+1, nhori, llm, 1, llm, nvert, &
319                   "inst(X)", zsto, zout)
320              !---fin Olivia
321    
322           ENDDO
323    
324         CALL histend(nid_ins)         CALL histend(nid_ins)
325      ENDIF      ENDIF
326    
Legend:
Removed from v.82  
changed lines
  Added in v.159
  ViewVC Help
Powered by ViewVC 1.1.21