/[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 155 by guez, Wed Jul 8 17:03:45 2015 UTC revision 178 by guez, Fri Mar 11 18:47: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    
11      use clesphys, only: ecrit_ins      use clesphys, only: ecrit_ins
12      use dimens_m, only: iim, jjm, llm      use dimens_m, only: iim, jjm, llm, nqmx
13      use dimphy, only: klon      use dimphy, only: klon
14      use disvert_m, only: presnivs      use disvert_m, only: presnivs
15      use dynetat0_m, only: day_ref, annee_ref      use dynetat0_m, only: day_ref, annee_ref
# Line 18  contains Line 18  contains
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      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
     use temps, only: itau_phy  
23      USE ymds2ju_m, only: ymds2ju      USE ymds2ju_m, only: ymds2ju
24    
25      REAL, intent(in):: dtime ! pas temporel de la physique (s)      REAL, intent(in):: dtime ! pas temporel de la physique (s)
26      logical, intent(in):: ok_instan      logical, intent(in):: ok_instan
27      integer, intent(out):: nid_ins      integer, intent(out):: nid_ins
28        integer, intent(in):: itau_phy
29    
30        ! Local:
31      REAL zx_lon(iim, jjm + 1), zx_lat(iim, jjm + 1)      REAL zx_lon(iim, jjm + 1), zx_lat(iim, jjm + 1)
32      real zjulian, zsto, zout      real zjulian, zsto, zout
33      integer i, nhori, nvert, nsrf      integer i, nhori, nvert, nsrf, iq, it
34    
35      !-------------------------------------------------------------------      !-------------------------------------------------------------------
36    
37      IF (ok_instan) THEN      IF (ok_instan) THEN
   
38         zsto = dtime * ecrit_ins         zsto = dtime * ecrit_ins
39         zout = dtime * ecrit_ins         zout = dtime * ecrit_ins
40    
# Line 291  contains Line 292  contains
292              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
293              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
294    
295           CALL histdef(nid_ins, "zmasse", "column density of air in cell", &
296                "kg m-2", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
297                zsto, zout)
298    
299           CALL histdef(nid_ins, "rhum", "Relative humidity", &
300                "", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
301                zsto, zout)
302    
303           DO it = 1, nqmx - 2
304              ! champ 2D
305              iq=it+2
306              CALL histdef(nid_ins, tname(iq), ttext(iq), "U/kga", iim, jjm+1, &
307                   nhori, llm, 1, llm, nvert, "inst(X)", zsto, zout)
308              CALL histdef(nid_ins, "fl"//tname(iq), "Flux "//ttext(iq), &
309                   "U/m2/s", iim, jjm+1, nhori, llm, 1, llm, nvert, &
310                   "inst(X)", zsto, zout)
311    
312              !---Ajout Olivia
313              CALL histdef(nid_ins, "d_tr_th_"//tname(iq), &
314                   "tendance thermique"// ttext(iq), "?", &
315                   iim, jjm+1, nhori, llm, 1, llm, nvert, &
316                   "inst(X)", zsto, zout)
317              CALL histdef(nid_ins, "d_tr_cv_"//tname(iq), &
318                   "tendance convection"// ttext(iq), "?", &
319                   iim, jjm+1, nhori, llm, 1, llm, nvert, &
320                   "inst(X)", zsto, zout)
321              CALL histdef(nid_ins, "d_tr_cl_"//tname(iq), &
322                   "tendance couche limite"// ttext(iq), "?", &
323                   iim, jjm+1, nhori, llm, 1, llm, nvert, &
324                   "inst(X)", zsto, zout)
325              !---fin Olivia
326    
327           ENDDO
328    
329         CALL histend(nid_ins)         CALL histend(nid_ins)
330      ENDIF      ENDIF
331    
Legend:
Removed from v.155  
changed lines
  Added in v.178
  ViewVC Help
Powered by ViewVC 1.1.21