source: LMDZ6/branches/Amaury_dev/libf/phylmd/ini_histREGDYN.h @ 5133

Last change on this file since 5133 was 5128, checked in by abarral, 5 months ago

Correct bug in vlspltqs_loc.f90 from r2270 where we call SSUM with incorrect arguments.
Merge the three different versions of abort_gcm into one
Fix seq, para 3D compilation broken from r5107 onwards
(lint) usual + Remove uneeded fixed-form continuations
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • 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      IF (ok_regdyn) THEN
5     
6        IF (is_sequential) THEN
7c
8cIM      PRINT*, 'La frequence de sortie REGDYN est de ', ecrit_mth
9c
10         idayref = day_ref
11         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
12c
13c axe vertical pour les differents niveaux des histogrammes
14      DO iw=1, iwmax
15        zx_o500(iw)=wmin+(iw-1./2.)*pas_w
16      ENDDO
17
18         CALL histbeg("histREGDYN", kmaxm1,zx_tau, lmaxm1,zx_pc,
19     .                 1,kmaxm1,1,lmaxm1, itau_phy, zjulian, dtime, 
20     .                 nhoriRD, nid_regdyn)
21
22         CALL histvert(nid_regdyn, "omeganivs", "Omega levels", 
23     .                 "mb/day",
24     .                 iwmax, zx_o500, komega)
25c
26c   pour les champs instantannes, il faut mettre la meme valeur pour
27c   zout et zsto.
28c   dtime est passe par ailleurs a histbeg
29c
30c        zout = dtime * REAL(NINT(86400./dtime*ecrit_regdyn))
31c        zsto = zout
32c        PRINT*,'zout,zsto=',zout,zsto
33c
34c stockage a chaque pas de temps de la physique
35c
36         zstophy = dtime
37cIM 020904      zstophy = dtime * nbapp_isccp
38
39c ecriture mensuelle
40c
41         zout = dtime * ecrit_mth
42cIM 020904     
43c        zout = dtime * ecrit_day
44c        zout = dtime * REAL(NINT(86400./dtime*ecrit_regdyn))
45
46c
47c Champs 3D:
48c
49c TROP
50         CALL histdef(nid_regdyn, "hw1", "Tropics Histogram ", "%",
51                      kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32,
52                      "ave(X)", zstophy,zout)
53
54         CALL histdef(nid_regdyn, "nh1", "Nb of pixels Tropics Histo",
55                      "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
56                      32,"ave(X)", zstophy,zout)
57c
58
59         CALL histdef(nid_regdyn, "nht1",
60                      "Total Nb pixels Tropics Histo",
61                      "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
62                      32,"ave(X)", zstophy,zout)
63c
64c PAN
65         CALL histdef(nid_regdyn, "hw2", "North Pacific Histogram", "%",
66                      kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32,
67                      "ave(X)", zstophy,zout)
68
69         CALL histdef(nid_regdyn, "nh2", "Nb of pixels North Pacific",
70                      "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
71                      32,"ave(X)", zstophy,zout)
72c
73
74         CALL histdef(nid_regdyn, "nht2",
75                      "Total Nb pixels North Pacific Histo",
76                      "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
77                      32,"ave(X)", zstophy,zout)
78c CAL
79         CALL histdef(nid_regdyn, "hw3", "California Histogram", "%",
80                      kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32,
81                      "ave(X)", zstophy,zout)
82
83         CALL histdef(nid_regdyn, "nh3", 
84                      "Nb of pixels California Histo",
85                      "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
86                      32,"ave(X)", zstophy,zout)
87c
88
89         CALL histdef(nid_regdyn, "nht3",
90                      "Total Nb pixels California Histo",
91                      "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
92                      32,"ave(X)", zstophy,zout)
93c HAW
94         CALL histdef(nid_regdyn, "hw4", "Hawai Histogram", "%",
95                      kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32,
96                      "ave(X)", zstophy,zout)
97
98         CALL histdef(nid_regdyn, "nh4", "Nb of pixels Hawai Histo",
99                      "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
100                      32,"ave(X)", zstophy,zout)
101c
102
103         CALL histdef(nid_regdyn, "nht4",
104                      "Total Nb pixels Hawai Histo",
105                      "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
106                      32,"ave(X)", zstophy,zout)
107c WAP
108         CALL histdef(nid_regdyn, "hw5", "Warm Pool Histogram", "%",
109                      kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32,
110                      "ave(X)", zstophy,zout)
111
112         CALL histdef(nid_regdyn, "nh5", "Nb of pixels Warm Pool Histo",
113                      "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
114                      32,"ave(X)", zstophy,zout)
115c
116
117         CALL histdef(nid_regdyn, "nht5",
118                      "Total Nb pixels Warm Pool Histo",
119                      "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
120                      32,"ave(X)", zstophy,zout)
121c
122         CALL histend(nid_regdyn)
123         
124         END IF ! is_sequential
125
126      endif ! ok_regdyn
Note: See TracBrowser for help on using the repository browser.