/[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 175 by guez, Fri Feb 5 16:02:34 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 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 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, idayref, 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    
41         idayref = day_ref         CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
        CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)  
42    
43         CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlon, zx_lon)         CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlon, zx_lon)
44         DO i = 1, iim         DO i = 1, iim
# Line 217  contains Line 218  contains
218         CALL histdef(nid_ins, "albs", "Surface albedo", "-", &         CALL histdef(nid_ins, "albs", "Surface albedo", "-", &
219              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
220              "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)  
221    
        !IM cf. AM 081204 BEG  
        ! HBTM2  
222         CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m", &         CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m", &
223              iim, (jjm + 1), nhori, 1, 1, 1, -99, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
224              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
# Line 296  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           DO it = 1, nqmx - 2
300              ! champ 2D
301              iq=it+2
302              CALL histdef(nid_ins, tname(iq), ttext(iq), "U/kga", iim, jjm+1, &
303                   nhori, llm, 1, llm, nvert, "inst(X)", zsto, zout)
304              CALL histdef(nid_ins, "fl"//tname(iq), "Flux "//ttext(iq), &
305                   "U/m2/s", iim, jjm+1, nhori, llm, 1, llm, nvert, &
306                   "inst(X)", zsto, zout)
307    
308              !---Ajout Olivia
309              CALL histdef(nid_ins, "d_tr_th_"//tname(iq), &
310                   "tendance thermique"// ttext(iq), "?", &
311                   iim, jjm+1, nhori, llm, 1, llm, nvert, &
312                   "inst(X)", zsto, zout)
313              CALL histdef(nid_ins, "d_tr_cv_"//tname(iq), &
314                   "tendance convection"// ttext(iq), "?", &
315                   iim, jjm+1, nhori, llm, 1, llm, nvert, &
316                   "inst(X)", zsto, zout)
317              CALL histdef(nid_ins, "d_tr_cl_"//tname(iq), &
318                   "tendance couche limite"// ttext(iq), "?", &
319                   iim, jjm+1, nhori, llm, 1, llm, nvert, &
320                   "inst(X)", zsto, zout)
321              !---fin Olivia
322    
323           ENDDO
324    
325         CALL histend(nid_ins)         CALL histend(nid_ins)
326      ENDIF      ENDIF
327    
Legend:
Removed from v.82  
changed lines
  Added in v.175
  ViewVC Help
Powered by ViewVC 1.1.21