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 | !c |
---|
11 | zstophy = dtime |
---|
12 | zout = mth_len*un_jour |
---|
13 | !c |
---|
14 | idayref = day_ref |
---|
15 | CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) |
---|
16 | !c |
---|
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) |
---|
25 | !c |
---|
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 | !c |
---|
32 | CALL histdef(nid_ctesGCM, "R_ecc", & |
---|
33 | & "Excentricite","-", & |
---|
34 | & 1,1,nhori, 1,1,1, -99, 32, & |
---|
35 | & "ave", zstophy,zout) |
---|
36 | !c |
---|
37 | CALL histdef(nid_ctesGCM, "R_peri", & |
---|
38 | & "Equinoxe","-", & |
---|
39 | & 1,1,nhori, 1,1,1, -99, 32, & |
---|
40 | & "ave", zstophy,zout) |
---|
41 | !c |
---|
42 | CALL histdef(nid_ctesGCM, "R_incl", & |
---|
43 | & "Inclinaison","deg", & |
---|
44 | & 1,1,nhori, 1,1,1, -99, 32, & |
---|
45 | & "ave", zstophy,zout) |
---|
46 | !c |
---|
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 | !c |
---|
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 | !c |
---|
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 | !c |
---|
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 | !c |
---|
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 | !c |
---|
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 | !c |
---|
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 | !c |
---|
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 | !c |
---|
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 | !c |
---|
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 | !c |
---|
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 | !c |
---|
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 | !c |
---|
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 | !c |
---|
113 | CALL histdef(nid_ctesGCM, "tsol", & |
---|
114 | & "Surface temperature", "K", & |
---|
115 | & 1,1,nhori, 1,1,1, -99, 32, & |
---|
116 | & "ave", zstophy,zout) |
---|
117 | !c |
---|
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 | !c |
---|
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 | !c================================================================= |
---|
128 | !c |
---|
129 | CALL histend(nid_ctesGCM) |
---|
130 | |
---|
131 | endif !(is_mpi_root) |
---|
132 | !$OMP END MASTER |
---|
133 | !c================================================================= |
---|