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

Last change on this file since 5456 was 776, checked in by Laurent Fairhead, 18 years ago

Suite du merge entre la version et la HEAD: quelques modifications
de Yann sur le

LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.3 KB
RevLine 
[524]1!
2! $Header$
3!
4
5      IF (ok_regdyn) THEN
[766]6     
[776]7        if (is_sequential) then
[524]8c
[644]9cIM      PRINT*, 'La frequence de sortie REGDYN est de ', ecrit_mth
10c
[524]11         idayref = day_ref
12         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
13c
14c 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)
[644]26c
[524]27c   pour les champs instantannes, il faut mettre la meme valeur pour
[644]28c   zout et zsto.
[524]29c   dtime est passe par ailleurs a histbeg
[644]30c
[524]31c        zout = dtime * REAL(NINT(86400./dtime*ecrit_regdyn))
32c        zsto = zout
33c        print*,'zout,zsto=',zout,zsto
[644]34c
35c stockage a chaque pas de temps de la physique
36c
37         zstophy = dtime
38cIM 020904      zstophy = dtime * nbapp_isccp
39
[524]40c ecriture mensuelle
41c
42         zout = dtime * ecrit_mth
[644]43cIM 020904     
44c        zout = dtime * ecrit_day
[524]45c        zout = dtime * REAL(NINT(86400./dtime*ecrit_regdyn))
46
47c
48c Champs 3D:
49c
50c TROP
51         CALL histdef(nid_regdyn, "hw1", "Tropics Histogram ", "%",
52     &                kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, 
[644]53     &                "ave(X)", zstophy,zout)
[524]54
55         CALL histdef(nid_regdyn, "nh1", "Nb of pixels Tropics Histo",
56     &                "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
[644]57     &                32,"ave(X)", zstophy,zout)
[524]58c
59
[644]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)
[524]64c
65c PAN
66         CALL histdef(nid_regdyn, "hw2", "North Pacific Histogram", "%",
67     &                kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, 
[644]68     &                "ave(X)", zstophy,zout)
[524]69
70         CALL histdef(nid_regdyn, "nh2", "Nb of pixels North Pacific",
71     &                "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
[644]72     &                32,"ave(X)", zstophy,zout)
[524]73c
74
[644]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)
[524]79c CAL
80         CALL histdef(nid_regdyn, "hw3", "California Histogram", "%",
81     &                kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, 
[644]82     &                "ave(X)", zstophy,zout)
[524]83
[644]84         CALL histdef(nid_regdyn, "nh3", 
85     &                "Nb of pixels California Histo",
[524]86     &                "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
[644]87     &                32,"ave(X)", zstophy,zout)
[524]88c
89
[644]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)
[524]94c HAW
95         CALL histdef(nid_regdyn, "hw4", "Hawai Histogram", "%",
96     &                kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, 
[644]97     &                "ave(X)", zstophy,zout)
[524]98
99         CALL histdef(nid_regdyn, "nh4", "Nb of pixels Hawai Histo",
100     &                "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
[644]101     &                32,"ave(X)", zstophy,zout)
[524]102c
103
[644]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)
[524]108c WAP
109         CALL histdef(nid_regdyn, "hw5", "Warm Pool Histogram", "%",
110     &                kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, 
[644]111     &                "ave(X)", zstophy,zout)
[524]112
113         CALL histdef(nid_regdyn, "nh5", "Nb of pixels Warm Pool Histo",
114     &                "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega,
[644]115     &                32,"ave(X)", zstophy,zout)
[524]116c
117
[644]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)
[524]122c
123         CALL histend(nid_regdyn)
[766]124         
[776]125         endif ! is_sequential
[524]126
127      endif ! ok_regdyn
Note: See TracBrowser for help on using the repository browser.