source: LMDZ4/branches/unlabeled-1.1.1/libf/phylmd/ini_histREGDYN.h @ 648

Last change on this file since 648 was 524, checked in by lmdzadmin, 21 years ago

Initial revision

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