| 1 | !IM Implemente en modes sequentiel et parallele |
|---|
| 2 | |
|---|
| 3 | CALL gather(latitude_deg,rlat_glo) |
|---|
| 4 | CALL bcast(rlat_glo) |
|---|
| 5 | CALL gather(longitude_deg,rlon_glo) |
|---|
| 6 | CALL bcast(rlon_glo) |
|---|
| 7 | |
|---|
| 8 | !$OMP MASTER |
|---|
| 9 | if (is_mpi_root) then |
|---|
| 10 | ! |
|---|
| 11 | zstophy = pdtphys |
|---|
| 12 | zout = -1 |
|---|
| 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(X)", zstophy,zout) |
|---|
| 36 | ! |
|---|
| 37 | CALL histdef(nid_ctesGCM, "R_peri", & |
|---|
| 38 | "Equinoxe","-", & |
|---|
| 39 | 1,1,nhori, 1,1,1, -99, 32, & |
|---|
| 40 | "ave(X)", zstophy,zout) |
|---|
| 41 | ! |
|---|
| 42 | CALL histdef(nid_ctesGCM, "R_incl", & |
|---|
| 43 | "Inclinaison","deg", & |
|---|
| 44 | 1,1,nhori, 1,1,1, -99, 32, & |
|---|
| 45 | "ave(X)", 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(X)", 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 | !================================================================= |
|---|
| 78 | ! |
|---|
| 79 | CALL histend(nid_ctesGCM) |
|---|
| 80 | |
|---|
| 81 | endif !(is_mpi_root) |
|---|
| 82 | !$OMP END MASTER |
|---|
| 83 | !================================================================= |
|---|
