1 | !IM 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 | |
---|
8 | !$OMP MASTER |
---|
9 | if (is_mpi_root) then |
---|
10 | ! |
---|
11 | zstophy = pdtphys |
---|
12 | zout = mth_len*un_jour |
---|
13 | ! |
---|
14 | idayref = day_ref |
---|
15 | CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) |
---|
16 | ! |
---|
17 | CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlon_glo,zx_lon) |
---|
18 | if (nbp_lon.gt.1) then |
---|
19 | DO i = 1, nbp_lon |
---|
20 | zx_lon(i,1) = rlon_glo(i+1) |
---|
21 | zx_lon(i,nbp_lat) = rlon_glo(i+1) |
---|
22 | ENDDO |
---|
23 | endif |
---|
24 | CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlat_glo,zx_lat) |
---|
25 | ! |
---|
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) |
---|
31 | ! |
---|
32 | CALL histdef(nid_ctesGCM, "R_ecc", & |
---|
33 | "Excentricite","-", & |
---|
34 | 1,1,nhori, 1,1,1, -99, 32, & |
---|
35 | "ave", zstophy,zout) |
---|
36 | ! |
---|
37 | CALL histdef(nid_ctesGCM, "R_peri", & |
---|
38 | "Equinoxe","-", & |
---|
39 | 1,1,nhori, 1,1,1, -99, 32, & |
---|
40 | "ave", zstophy,zout) |
---|
41 | ! |
---|
42 | CALL histdef(nid_ctesGCM, "R_incl", & |
---|
43 | "Inclinaison","deg", & |
---|
44 | 1,1,nhori, 1,1,1, -99, 32, & |
---|
45 | "ave", zstophy,zout) |
---|
46 | ! |
---|
47 | CALL histdef(nid_ctesGCM, "solaire", & |
---|
48 | "Constante solaire","W/m2", & |
---|
49 | 1,1,nhori, 1,1,1, -99, 32, & |
---|
50 | "ave", zstophy,zout) |
---|
51 | ! |
---|
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) |
---|
56 | ! |
---|
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) |
---|
61 | ! |
---|
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) |
---|
66 | ! |
---|
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) |
---|
71 | ! |
---|
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) |
---|
76 | ! |
---|
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) |
---|
81 | ! |
---|
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) |
---|
86 | ! |
---|
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) |
---|
91 | ! |
---|
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) |
---|
96 | ! |
---|
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) |
---|
102 | ! |
---|
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) |
---|
107 | ! |
---|
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) |
---|
112 | ! |
---|
113 | CALL histdef(nid_ctesGCM, "tsol", & |
---|
114 | "Surface temperature", "K", & |
---|
115 | 1,1,nhori, 1,1,1, -99, 32, & |
---|
116 | "ave", zstophy,zout) |
---|
117 | ! |
---|
118 | CALL histdef(nid_ctesGCM, "t2m", & |
---|
119 | "Temperature at 2m", "K", & |
---|
120 | 1,1,nhori, 1,1,1, -99, 32, & |
---|
121 | "ave", zstophy,zout) |
---|
122 | ! |
---|
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) |
---|
127 | !================================================================= |
---|
128 | ! |
---|
129 | CALL histend(nid_ctesGCM) |
---|
130 | |
---|
131 | endif !(is_mpi_root) |
---|
132 | !$OMP END MASTER |
---|
133 | !================================================================= |
---|