/[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

revision 155 by guez, Wed Jul 8 17:03:45 2015 UTC revision 190 by guez, Thu Apr 14 15:15:56 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 clesphys2, only: conv_emanuel
13        use dimens_m, only: iim, jjm, llm, nqmx
14      use dimphy, only: klon      use dimphy, only: klon
15      use disvert_m, only: presnivs      use disvert_m, only: presnivs
16      use dynetat0_m, only: day_ref, annee_ref      use dynetat0_m, only: day_ref, annee_ref, rlatu, rlonv
17      USE histbeg_totreg_m, ONLY : histbeg_totreg      USE histbeg_totreg_m, ONLY : histbeg_totreg
18      USE histdef_m, ONLY : histdef      USE histdef_m, ONLY : histdef
19      USE histend_m, ONLY : histend      USE histend_m, ONLY : histend
20      USE histvert_m, ONLY : histvert      USE histvert_m, ONLY : histvert
21      use indicesol, only: nbsrf, clnsurf      use indicesol, only: nbsrf, clnsurf
22      use phyetat0_m, only: rlon, rlat      use iniadvtrac_m, only: tname, ttext
23      use temps, only: itau_phy      use nr_util, only: pi
24      USE ymds2ju_m, only: ymds2ju      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
28      integer, intent(out):: nid_ins      integer, intent(out):: nid_ins
29        integer, intent(in):: itau_phy
30    
31      REAL zx_lon(iim, jjm + 1), zx_lat(iim, jjm + 1)      ! Local:
32      real zjulian, zsto, zout      real zjulian, zsto, zout
33      integer i, nhori, nvert, nsrf      integer 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         CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)         CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
41           CALL histbeg_totreg("histins", rlonv(:iim) / pi * 180., &
42         CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlon, zx_lon)              rlatu / pi * 180., 1, iim, &
43         DO i = 1, iim              1, jjm + 1, itau_phy, zjulian, dtime, nhori, nid_ins)
           zx_lon(i, 1) = rlon(i+1)  
           zx_lon(i, (jjm + 1)) = rlon(i+1)  
        ENDDO  
        CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlat, zx_lat)  
        CALL histbeg_totreg("histins", zx_lon(:, 1), zx_lat(1, :), 1, iim, 1, &  
             jjm + 1, itau_phy, zjulian, dtime, nhori, nid_ins)  
44         write(*, *)'Inst ', itau_phy, zjulian         write(*, *)'Inst ', itau_phy, zjulian
45         CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", &         CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", &
46              presnivs/100., nvert)              presnivs/100., nvert)
   
47         CALL histdef(nid_ins, "phis", "Surface geop. height", "-", &         CALL histdef(nid_ins, "phis", "Surface geop. height", "-", &
48              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
49              "once", zsto, zout)              "once", zsto, zout)
   
50         CALL histdef(nid_ins, "aire", "Grid area", "-", &         CALL histdef(nid_ins, "aire", "Grid area", "-", &
51              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
52              "once", zsto, zout)              "once", zsto, zout)
# Line 259  contains Line 251  contains
251              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
252              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
253    
254         !IM cf. AM 081204 END         if (conv_emanuel) CALL histdef(nid_ins, "ptop", "cloud top pressure", &
255                "Pa", iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout)
256    
257         ! Champs 3D:         ! Champs 3D:
258    
# Line 291  contains Line 284  contains
284              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
285              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
286    
287           CALL histdef(nid_ins, "zmasse", "column density of air in cell", &
288                "kg m-2", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
289                zsto, zout)
290    
291           CALL histdef(nid_ins, "rhum", "Relative humidity", &
292                "", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
293                zsto, zout)
294    
295           DO it = 1, nqmx - 2
296              ! champ 2D
297              iq=it+2
298              CALL histdef(nid_ins, tname(iq), ttext(iq), "U/kga", iim, jjm+1, &
299                   nhori, llm, 1, llm, nvert, "inst(X)", zsto, zout)
300              CALL histdef(nid_ins, "fl"//tname(iq), "Flux "//ttext(iq), &
301                   "U/m2/s", iim, jjm+1, nhori, llm, 1, llm, nvert, &
302                   "inst(X)", zsto, zout)
303    
304              !---Ajout Olivia
305              CALL histdef(nid_ins, "d_tr_th_"//tname(iq), &
306                   "tendance thermique"// ttext(iq), "?", &
307                   iim, jjm+1, nhori, llm, 1, llm, nvert, &
308                   "inst(X)", zsto, zout)
309              CALL histdef(nid_ins, "d_tr_cv_"//tname(iq), &
310                   "tendance convection"// ttext(iq), "?", &
311                   iim, jjm+1, nhori, llm, 1, llm, nvert, &
312                   "inst(X)", zsto, zout)
313              CALL histdef(nid_ins, "d_tr_cl_"//tname(iq), &
314                   "tendance couche limite"// ttext(iq), "?", &
315                   iim, jjm+1, nhori, llm, 1, llm, nvert, &
316                   "inst(X)", zsto, zout)
317              !---fin Olivia
318    
319           ENDDO
320    
321         CALL histend(nid_ins)         CALL histend(nid_ins)
322      ENDIF      ENDIF
323    
Legend:
Removed from v.155  
changed lines
  Added in v.190
  ViewVC Help
Powered by ViewVC 1.1.21