| [1910] | 1 | !IM Implemente en modes sequentiel et parallele |
|---|
| [766] | 2 | |
|---|
| [2435] | 3 | CALL gather(latitude_deg,rlat_glo) |
|---|
| [1538] | 4 | CALL bcast(rlat_glo) |
|---|
| [2435] | 5 | CALL gather(longitude_deg,rlon_glo) |
|---|
| [1538] | 6 | CALL bcast(rlon_glo) |
|---|
| 7 | |
|---|
| [1910] | 8 | !$OMP MASTER |
|---|
| [1538] | 9 | if (is_mpi_root) then |
|---|
| [1910] | 10 | ! |
|---|
| [2408] | 11 | zstophy = pdtphys |
|---|
| [1538] | 12 | zout = mth_len*un_jour |
|---|
| [1910] | 13 | ! |
|---|
| [683] | 14 | idayref = day_ref |
|---|
| 15 | CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) |
|---|
| [1910] | 16 | ! |
|---|
| [2408] | 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 |
|---|
| [1538] | 20 | zx_lon(i,1) = rlon_glo(i+1) |
|---|
| [2408] | 21 | zx_lon(i,nbp_lat) = rlon_glo(i+1) |
|---|
| [683] | 22 | ENDDO |
|---|
| [879] | 23 | endif |
|---|
| [2408] | 24 | CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlat_glo,zx_lat) |
|---|
| [1910] | 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 | ! |
|---|
| [683] | 129 | CALL histend(nid_ctesGCM) |
|---|
| [766] | 130 | |
|---|
| [1538] | 131 | endif !(is_mpi_root) |
|---|
| [1910] | 132 | !$OMP END MASTER |
|---|
| 133 | !================================================================= |
|---|
