1 !IM Implemente en modes sequentiel et parallele
18 if (nbp_lon.gt.1) then
20 zx_lon(i,1) = rlon_glo(i+1)
21 zx_lon(i,nbp_lat) = rlon_glo(i+1)
24 CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlat_glo,
zx_lat)
26 CALL
histbeg("paramLMDZ_phy.nc", &
32 CALL histdef(nid_ctesGCM, "
R_ecc", &
34 1,1,nhori, 1,1,1, -99, 32, &
37 CALL histdef(nid_ctesGCM, "
R_peri", &
39 1,1,nhori, 1,1,1, -99, 32, &
42 CALL histdef(nid_ctesGCM, "R_incl", &
43 "Inclinaison","deg", &
44 1,1,nhori, 1,1,1, -99, 32, &
47 CALL histdef(nid_ctesGCM, "solaire", &
48 "Constante solaire","W/m2", &
49 1,1,nhori, 1,1,1, -99, 32, &
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)
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)
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)
68 "Concentration
du CFC11", "ppt", &
69 1,1,nhori, 1,1,1, -99, 32, &
70 "ave(X)", zstophy,
zout)
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)
77 CALL histdef(nid_ctesGCM, "
bils", &
78 "Surface total
heat flux", "W m-2", &
79 1,1,nhori, 1,1,1, -99, 32, &
82 CALL histdef(nid_ctesGCM, "
evap", &
83 "Evaporation", "kg m-2 s-1", &
84 1,1,nhori, 1,1,1, -99, 32, &
87 CALL histdef(nid_ctesGCM, "evap_land", &
88 "Land evaporation", "kg m-2 s-1", &
89 1,1,nhori, 1,1,1, -99, 32, &
92 CALL histdef(nid_ctesGCM, "
flat", &
93 "Latent
heat flux", "W m-2", &
94 1,1,nhori, 1,1,1, -99, 32, &
97 CALL histdef(nid_ctesGCM, "nettop0", &
98 "Clear sky net downward radiatif flux
at TOA", &
100 1,1,nhori, 1,1,1, -99, 32, &
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, &
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, &
113 CALL histdef(nid_ctesGCM, "
tsol", &
114 "Surface temperature", "K", &
115 1,1,nhori, 1,1,1, -99, 32, &
118 CALL histdef(nid_ctesGCM, "
t2m", &
119 "Temperature
at 2m", "K", &
120 1,1,nhori, 1,1,1, -99, 32, &
123 CALL histdef(nid_ctesGCM, "
prw", &
124 "Precipitable water", "kg m-2", &
125 1,1,nhori, 1,1,1, -99, 32, &
127 !=================================================================
129 CALL histend(nid_ctesGCM)
133 !=================================================================
!$Id Turb_fcg_gcssold if(prt_level.ge.1) then print *
INTERFACE SUBROUTINE RRTM_ECRT_140GP klon
!IM Implemente en modes sequentiel et parallele CALL rlon_glo CALL bcast(rlon_glo)!$OMP MASTER if(is_mpi_root) then!zstophy
real, dimension(:), allocatable, save tsol
!$Id zjulian!correction pour l heure initiale!jyg!jyg CALL ymds2ju(annee_ref, 1, day_ref, hour, zjulian)!jyg CALL histbeg_phy("histrac"
real, dimension(:), allocatable, save bils
c c zjulian c cym CALL iim cym klev cym zjulian
real, dimension(:,:), allocatable, save heat
!$Id ok_orolf LOGICAL ok_limitvrai LOGICAL ok_all_xml INTEGER iflag_ener_conserv REAL solaire RCFC12 RCFC12_act CFC11_ppt
!$Header!integer nvarmx dtime
!$Id ok_orolf LOGICAL ok_limitvrai LOGICAL ok_all_xml INTEGER iflag_ener_conserv REAL solaire RCFC12 RCFC12_act CFC12_ppt!IM ajout CFMIP2 CMIP5 LOGICAL ok_4xCO2atm RCFC12_per CFC12_ppt_per!OM correction du bilan d eau global!OM Correction sur precip KE REAL cvl_corr!OM Fonte calotte dans bilan eau LOGICAL ok_lic_melt!IM simulateur ISCCP INTEGER overlap!IM seuils cdrh REAL cdhmax!IM param stabilite s terres et en dehors REAL f_ri_cd_min!IM MAFo pmagic evap0!Frottement au sol(Cdrag) Real f_cdrag_ter
c c zjulian c cym CALL iim cym klev cym cym nhori
!$Id klon initialisation mois suivants day_rain itap ENDIF!Calcul fin de nday_rain calcul nday_rain itap DO i
real, dimension(:,:), allocatable, save t2m
!IM Implemente en modes sequentiel et parallele CALL gather(rlat, rlat_glo) CALL bcast(rlat_glo) CALL gather(rlon
real, dimension(:), allocatable, save evap
!$Id RNAVO!A1 Astronomical constants REAL ROMEGA!A1 bis Constantes concernant l orbite de la R_incl!A1 Geoide REAL R1SA!A1 Radiation!REAL RI0 REAL RSIGMA!A1 Thermodynamic gas phase REAL RCVV REAL RETV Thermodynamic solid phases REAL RCS!A1 Thermodynamic transition of phase REAL RATM!A1 Curve of saturation REAL RGAMS REAL RGAMD!COMMON YOMCST RNAVO ROMEGA R_ecc
!$Id ok_orolf LOGICAL ok_limitvrai LOGICAL ok_all_xml INTEGER iflag_ener_conserv REAL co2_ppm
!$Id La frequence de sortie REGDYN est ecrit_mth c iwmax il faut mettre la meme valeur pour c zout et zsto c dtime est passe par ailleurs a histbeg c c zout
c c $Id c nbregdyn DO klon c rlat(i) c ENDIF!lon c ENDIF!lat ENDIF!pctsrf ENDDO!klon ENDDO!nbregdyn cIM 190504 ENDIF!ok_regdyn cIM somme de toutes les nhistoW BEG IF(debut) THEN DO nreg
!$Id zjulian!correction pour l heure initiale!jyg!jyg CALL pdtphys
real, dimension(:,:), pointer, save du
c c zjulian c cym CALL iim cym klev * zx_lat
real(kind=8), dimension(8, 3), parameter at
subroutine gr_fi_ecrit(nfield, nlon, iim, jjmp1, fi, ecrit)
logical, save is_mpi_root
real, dimension(:), allocatable, save prw
!$Id RNAVO!A1 Astronomical constants REAL ROMEGA!A1 bis Constantes concernant l orbite de la R_peri
c c zjulian c cym CALL iim cym klev iim cym jjmp1 cym On stoke le fichier bilKP instantanne s jmax_ins print On stoke le fichier bilKP instantanne s s cym cym nid_bilKPins ENDIF c cIM BEG c cIM cf AM BEG region cym CALL histbeg("histbilKP_ins", iim, zx_lon(:, 1), cym.jjmp1, zx_lat(1,:), cym.imin_ins, imax_ins-imin_ins+1, cym.jmin_ins, jmax_ins-jmin_ins+1, cym.itau_phy, zjulian, dtime, cym.nhori, nid_bilKPins) CALL histbeg_phy("histbilKP_ins"
c c zjulian c cym CALL iim cym zx_lon(i, 1)
!$Id ok_orolf LOGICAL ok_limitvrai LOGICAL ok_all_xml INTEGER iflag_ener_conserv REAL solaire RCFC12 RCFC12_act N2O_ppb
c c $Id c nbregdyn DO klon c rlon(i)