source: LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/ini_histREGDYN.h @ 1373

Last change on this file since 1373 was 776, checked in by Laurent Fairhead, 17 years ago

Suite du merge entre la version et la HEAD: quelques modifications
de Yann sur le

LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.3 KB
Line 
1!
2! $Header$
3!
4
5      IF (ok_regdyn) THEN
6     
7        if (is_sequential) then
8c
9cIM      PRINT*, 'La frequence de sortie REGDYN est de ', ecrit_mth
10c
11         idayref = day_ref
12         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
13c
14c axe vertical pour les differents niveaux des histogrammes
15      DO iw=1, iwmax
16        zx_o500(iw)=wmin+(iw-1./2.)*pas_w
17      ENDDO
18
19         CALL histbeg("histREGDYN", kmaxm1,zx_tau, lmaxm1,zx_pc,
20     .                 1,kmaxm1,1,lmaxm1, itau_phy, zjulian, dtime, 
21     .                 nhoriRD, nid_regdyn)
22
23         CALL histvert(nid_regdyn, "omeganivs", "Omega levels", 
24     .                 "mb/day",
25     .                 iwmax, zx_o500, komega)
26c
27c   pour les champs instantannes, il faut mettre la meme valeur pour
28c   zout et zsto.
29c   dtime est passe par ailleurs a histbeg
30c
31c        zout = dtime * REAL(NINT(86400./dtime*ecrit_regdyn))
32c        zsto = zout
33c        print*,'zout,zsto=',zout,zsto
34c
35c stockage a chaque pas de temps de la physique
36c
37         zstophy = dtime
38cIM 020904      zstophy = dtime * nbapp_isccp
39
40c ecriture mensuelle
41c
42         zout = dtime * ecrit_mth
43cIM 020904     
44c        zout = dtime * ecrit_day
45c        zout = dtime * REAL(NINT(86400./dtime*ecrit_regdyn))
46
47c
48c Champs 3D:
49c
50c TROP
51         CALL histdef(nid_regdyn, "hw1", "Tropics Histogram ", "%",
52     &                kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, 
53     &                "ave(X)", zstophy,zout)
54
55         CALL histdef(nid_regdyn, "nh1", "Nb of pixels Tropics Histo",
56     &                "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
57     &                32,"ave(X)", zstophy,zout)
58c
59
60         CALL histdef(nid_regdyn, "nht1",
61     &                "Total Nb pixels Tropics Histo",
62     &                "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
63     &                32,"ave(X)", zstophy,zout)
64c
65c PAN
66         CALL histdef(nid_regdyn, "hw2", "North Pacific Histogram", "%",
67     &                kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, 
68     &                "ave(X)", zstophy,zout)
69
70         CALL histdef(nid_regdyn, "nh2", "Nb of pixels North Pacific",
71     &                "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
72     &                32,"ave(X)", zstophy,zout)
73c
74
75         CALL histdef(nid_regdyn, "nht2",
76     &                "Total Nb pixels North Pacific Histo",
77     &                "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
78     &                32,"ave(X)", zstophy,zout)
79c CAL
80         CALL histdef(nid_regdyn, "hw3", "California Histogram", "%",
81     &                kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, 
82     &                "ave(X)", zstophy,zout)
83
84         CALL histdef(nid_regdyn, "nh3", 
85     &                "Nb of pixels California Histo",
86     &                "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
87     &                32,"ave(X)", zstophy,zout)
88c
89
90         CALL histdef(nid_regdyn, "nht3",
91     &                "Total Nb pixels California Histo",
92     &                "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
93     &                32,"ave(X)", zstophy,zout)
94c HAW
95         CALL histdef(nid_regdyn, "hw4", "Hawai Histogram", "%",
96     &                kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, 
97     &                "ave(X)", zstophy,zout)
98
99         CALL histdef(nid_regdyn, "nh4", "Nb of pixels Hawai Histo",
100     &                "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
101     &                32,"ave(X)", zstophy,zout)
102c
103
104         CALL histdef(nid_regdyn, "nht4",
105     &                "Total Nb pixels Hawai Histo",
106     &                "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
107     &                32,"ave(X)", zstophy,zout)
108c WAP
109         CALL histdef(nid_regdyn, "hw5", "Warm Pool Histogram", "%",
110     &                kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, 
111     &                "ave(X)", zstophy,zout)
112
113         CALL histdef(nid_regdyn, "nh5", "Nb of pixels Warm Pool Histo",
114     &                "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
115     &                32,"ave(X)", zstophy,zout)
116c
117
118         CALL histdef(nid_regdyn, "nht5",
119     &                "Total Nb pixels Warm Pool Histo",
120     &                "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
121     &                32,"ave(X)", zstophy,zout)
122c
123         CALL histend(nid_regdyn)
124         
125         endif ! is_sequential
126
127      endif ! ok_regdyn
Note: See TracBrowser for help on using the repository browser.