source: LMDZ5/branches/LMDZ5-DOFOCO/libf/phylmd/ini_paramLMDZ_phy.h @ 5005

Last change on this file since 5005 was 1538, checked in by musat, 13 years ago

Ajout fichier 1D mensuel paramLMDZ_phy.nc contenant

  • les parametres orbitaux
  • les taux des GES
  • les moyennes globales (bils, evap, evap_land, flat, nettop0, nettop, precip, tsol, t2m, prw)

IM
  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 4.6 KB
Line 
1cIM    Implemente en modes sequentiel et parallele
2
3       CALL gather(rlat,rlat_glo)
4       CALL bcast(rlat_glo)
5       CALL gather(rlon,rlon_glo)
6       CALL bcast(rlon_glo)
7
8c$OMP MASTER
9      if (is_mpi_root) then
10c
11       zstophy = dtime
12       zout = mth_len*un_jour
13c
14       idayref = day_ref
15       CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
16c
17       CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon_glo,zx_lon)
18       if (iim.gt.1) then
19       DO i = 1, iim
20         zx_lon(i,1) = rlon_glo(i+1)
21         zx_lon(i,jjmp1) = rlon_glo(i+1)
22       ENDDO
23       endif
24       CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat_glo,zx_lat)
25c
26       CALL histbeg("paramLMDZ_phy.nc", 
27     .                 np,zx_lon(np:np,1), np,zx_lat(1,np:np),
28     .                 1,1,1,1,
29     .                 itau_phy, zjulian, dtime,
30     .                 nhori, nid_ctesGCM)
31c
32       CALL histdef(nid_ctesGCM, "R_ecc", 
33     .                "Excentricite","-",
34     .                1,1,nhori, 1,1,1, -99, 32,
35     .                "ave", zstophy,zout)
36c
37       CALL histdef(nid_ctesGCM, "R_peri", 
38     .                "Equinoxe","-",
39     .                1,1,nhori, 1,1,1, -99, 32,
40     .                "ave", zstophy,zout)
41c
42       CALL histdef(nid_ctesGCM, "R_incl", 
43     .                "Inclinaison","deg",
44     .                1,1,nhori, 1,1,1, -99, 32,
45     .                "ave", zstophy,zout)
46c
47       CALL histdef(nid_ctesGCM, "solaire", 
48     .                "Constante solaire","W/m2",
49     .                1,1,nhori, 1,1,1, -99, 32,
50     .                "ave", zstophy,zout)
51c
52       CALL histdef(nid_ctesGCM, "co2_ppm", 
53     .                "Concentration du CO2", "ppm",
54     .                1,1,nhori, 1,1,1, -99, 32, 
55     .                "ave(X)", zstophy,zout)
56c
57       CALL histdef(nid_ctesGCM, "CH4_ppb", 
58     .                "Concentration du CH4", "ppb",
59     .                1,1,nhori, 1,1,1, -99, 32, 
60     .                "ave(X)", zstophy,zout)
61c
62       CALL histdef(nid_ctesGCM, "N2O_ppb",
63     .                "Concentration du N2O", "ppb",
64     .                1,1,nhori, 1,1,1, -99, 32,
65     .                "ave(X)", zstophy,zout)
66c
67       CALL histdef(nid_ctesGCM, "CFC11_ppt",
68     .                "Concentration du CFC11", "ppt",
69     .                1,1,nhori, 1,1,1, -99, 32,
70     .                "ave(X)", zstophy,zout)
71c
72       CALL histdef(nid_ctesGCM, "CFC12_ppt",
73     .                "Concentration du CFC12", "ppt",
74     .                1,1,nhori, 1,1,1, -99, 32,
75     .                "ave(X)", zstophy,zout)
76c
77       CALL histdef(nid_ctesGCM, "bils",
78     .                "Surface total heat flux", "W m-2",
79     .                1,1,nhori, 1,1,1, -99, 32,
80     .                "ave", zstophy,zout)
81c
82       CALL histdef(nid_ctesGCM, "evap",
83     .                "Evaporation", "kg m-2 s-1",
84     .                1,1,nhori, 1,1,1, -99, 32,
85     .                "ave", zstophy,zout)
86c
87       CALL histdef(nid_ctesGCM, "evap_land",
88     .                "Land evaporation", "kg m-2 s-1",
89     .                1,1,nhori, 1,1,1, -99, 32,
90     .                "ave", zstophy,zout)
91c
92       CALL histdef(nid_ctesGCM, "flat",
93     .                "Latent heat flux", "W m-2",
94     .                1,1,nhori, 1,1,1, -99, 32,
95     .                "ave", zstophy,zout)
96c
97       CALL histdef(nid_ctesGCM, "nettop0",
98     .                "Clear sky net downward radiatif flux at TOA", 
99     .                "W m-2",
100     .                1,1,nhori, 1,1,1, -99, 32,
101     .                "ave", zstophy,zout)
102c
103       CALL histdef(nid_ctesGCM, "nettop",
104     .                "Net downward radiatif flux at TOA", "W m-2",
105     .                1,1,nhori, 1,1,1, -99, 32,
106     .                "ave", zstophy,zout)
107c
108       CALL histdef(nid_ctesGCM, "precip",
109     .                "Total precipitation (liq+sol)", "kg m-2 s-1",
110     .                1,1,nhori, 1,1,1, -99, 32,
111     .                "ave", zstophy,zout)
112c
113       CALL histdef(nid_ctesGCM, "tsol",
114     .                "Surface temperature", "K",
115     .                1,1,nhori, 1,1,1, -99, 32,
116     .                "ave", zstophy,zout)
117c
118       CALL histdef(nid_ctesGCM, "t2m",
119     .                "Temperature at 2m", "K",
120     .                1,1,nhori, 1,1,1, -99, 32,
121     .                "ave", zstophy,zout)
122c
123       CALL histdef(nid_ctesGCM, "prw",
124     .                "Precipitable water", "kg m-2",
125     .                1,1,nhori, 1,1,1, -99, 32,
126     .                "ave", zstophy,zout) 
127c=================================================================
128c
129       CALL histend(nid_ctesGCM)
130       
131       endif !(is_mpi_root)
132c$OMP END MASTER
133c=================================================================
Note: See TracBrowser for help on using the repository browser.