/[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 159 by guez, Tue Jul 21 15:29:52 2015 UTC 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    
11      use clesphys, only: ecrit_ins      use clesphys, only: ecrit_ins
12        use clesphys2, only: conv_emanuel
13      use dimens_m, only: iim, jjm, llm, nqmx      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
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
# Line 20  contains Line 22  contains
22      use indicesol, only: nbsrf, clnsurf      use indicesol, only: nbsrf, clnsurf
23      use iniadvtrac_m, only: tname, ttext      use iniadvtrac_m, only: tname, ttext
24      use phyetat0_m, only: rlon, rlat      use phyetat0_m, only: rlon, rlat
     use temps, only: itau_phy  
25      USE ymds2ju_m, only: ymds2ju      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, nsrf, iq, it      integer i, nhori, nvert, nsrf, iq, it
# Line 39  contains Line 42  contains
42    
43         CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)         CALL ymds2ju(annee_ref, 1, day_ref, 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
# Line 259  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 295  contains Line 299  contains
299              "kg m-2", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &              "kg m-2", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
300              zsto, zout)              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         DO it = 1, nqmx - 2
307            ! champ 2D            ! champ 2D
308            iq=it+2            iq=it+2
Legend:
Removed from v.159  
changed lines
  Added in v.189
  ViewVC Help
Powered by ViewVC 1.1.21