source: LMDZ4/trunk/libf/phylmd/ini_histREGDYN.h @ 684

Last change on this file since 684 was 644, checked in by Laurent Fairhead, 20 years ago

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