My Project
 All Classes Files Functions Variables Macros
ini_histREGDYN.h
Go to the documentation of this file.
1 !
2 ! $Header$
3 !
4 
5  IF (ok_regdyn) THEN
6 
7  if (is_sequential) then
8 c
9 cIM PRINT*, 'La frequence de sortie REGDYN est de ', ecrit_mth
10 c
11  idayref = day_ref
12  CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
13 c
14 c 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)
26 c
27 c pour les champs instantannes, il faut mettre la meme valeur pour
28 c zout et zsto.
29 c dtime est passe par ailleurs a histbeg
30 c
31 c zout = dtime * REAL(NINT(86400./dtime*ecrit_regdyn))
32 c zsto = zout
33 c print*,'zout,zsto=',zout,zsto
34 c
35 c stockage a chaque pas de temps de la physique
36 c
37  zstophy = dtime
38 cIM 020904 zstophy = dtime * nbapp_isccp
39 
40 c ecriture mensuelle
41 c
42  zout = dtime * ecrit_mth
43 cIM 020904
44 c zout = dtime * ecrit_day
45 c zout = dtime * REAL(NINT(86400./dtime*ecrit_regdyn))
46 
47 c
48 c Champs 3D:
49 c
50 c 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)
58 c
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)
64 c
65 c 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)
73 c
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)
79 c 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)
88 c
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)
94 c 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)
102 c
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)
108 c 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)
116 c
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)
122 c
123  CALL histend(nid_regdyn)
124 
125  endif ! is_sequential
126 
127  endif ! ok_regdyn