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)