source: LMDZ.3.3/trunk/libf/phylmd/ini_histREGDYN.h @ 1590

Last change on this file since 1590 was 484, checked in by lmdzadmin, 21 years ago

Initial import

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