LMDZ
ini_paramLMDZ_phy.h
Go to the documentation of this file.
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 !
11  zstophy = pdtphys
12  zout = mth_len*un_jour
13 !
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, &
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 !
129  CALL histend(nid_ctesGCM)
130 
131  endif !(is_mpi_root)
132 !$OMP END MASTER
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
Definition: ini_bilKP_ave.h:26
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
Definition: clesphys.h:19
!$Header!integer nvarmx dtime
Definition: gradsdef.h:20
!$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
Definition: ini_bilKP_ave.h:26
integer, save mth_len
!$Id itau_phy
Definition: temps.h:15
!$Id klon initialisation mois suivants day_rain itap ENDIF!Calcul fin de nday_rain calcul nday_rain itap DO i
Definition: calcul_divers.h:24
real, dimension(:,:), allocatable, save t2m
!IM Implemente en modes sequentiel et parallele CALL gather(rlat, rlat_glo) CALL bcast(rlat_glo) CALL gather(rlon
!$Id day_ref
Definition: temps.h:15
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
Definition: YOMCST.h:11
!$Id ok_orolf LOGICAL ok_limitvrai LOGICAL ok_all_xml INTEGER iflag_ener_conserv REAL co2_ppm
Definition: clesphys.h:12
!$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 idayref
Definition: ini_bilKP_ave.h:11
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
Definition: ini_histrac.h:11
real, dimension(:,:), pointer, save du
c c zjulian c cym CALL iim cym klev * zx_lat
Definition: ini_bilKP_ave.h:23
real(kind=8), dimension(8, 3), parameter at
subroutine gr_fi_ecrit(nfield, nlon, iim, jjmp1, fi, ecrit)
Definition: gr_fi_ecrit.F90:5
real, dimension(:), allocatable, save prw
!$Id RNAVO!A1 Astronomical constants REAL ROMEGA!A1 bis Constantes concernant l orbite de la R_peri
Definition: YOMCST.h:11
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"
!$Id flat
Definition: flux_arp.h:11
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
Definition: clesphys.h:19
c c $Id c nbregdyn DO klon c rlon(i)
!$Id annee_ref
Definition: temps.h:15