source: LMDZ6/branches/SETHET_DECOUPLE/libf/phylmd/clesphys.h @ 5441

Last change on this file since 5441 was 3435, checked in by Laurent Fairhead, 6 years ago

"Historic" :-) commit merging the physics branch used for DYNAMICO with the LMDZ trunk.
The same physics branch can now be used seamlessly with the traditional lon-lat LMDZ
dynamical core and DYNAMICO.
Testing consisted in running a lon-lat LMDZ bucket simulation with the NPv6.1 physics package
with the original trunk sources and the merged sources. Tests were succesful in the sense that
numeric continuity was preserved in the restart files from both simulation. Further tests
included running both versions of the physics codes for one year in a LMDZOR setting in which
the restart files also came out identical.

Caution:

  • as the physics package now manages unstructured grids, grid information needs to be transmitted

to the surface scheme ORCHIDEE. This means that the interface defined in surf_land_orchidee_mod.F90
is only compatible with ORCHIDEE version orchidee2.1 and later versions. If previous versions of
ORCHIDEE need to be used, the CPP key ORCHIDEE_NOUNSTRUCT needs to be set at compilation time.
This is done automatically if makelmdz/makelmdz_fcm are called with the veget orchidee2.0 switch

  • due to a limitation in XIOS, the time at which limit conditions will be read in by DYNAMICO will be

delayed by one physic timestep with respect to the time it is read in by the lon-lat model. This is caused
by the line

IF (MOD(itime-1, lmt_pas) == 0 .OR. (jour_lu /= jour .AND. grid_type /= unstructured)) THEN ! time to read

in limit_read_mod.F90

Work still needed on COSP integration and XML files for DYNAMICO

EM, YM, LF

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 6.8 KB
Line 
1! $Id: clesphys.h 3435 2019-01-22 15:21:59Z fhourdin $
2!
3!  ATTENTION!!!!: ce fichier include est compatible format fixe/format libre
4!                 veillez \`a n'utiliser que des ! pour les commentaires
5!                 et \`a bien positionner les & des lignes de continuation
6!                 (les placer en colonne 6 et en colonne 73)
7!
8!..include cles_phys.h
9!
10       INTEGER iflag_cycle_diurne
11       LOGICAL soil_model,new_oliq,ok_orodr,ok_orolf
12       LOGICAL ok_limitvrai
13       LOGICAL ok_all_xml
14       LOGICAL ok_lwoff
15       INTEGER nbapp_rad, iflag_con, nbapp_cv, nbapp_wk, iflag_ener_conserv
16       REAL co2_ppm, co2_ppm0, solaire
17!FC
18       REAL Cd_frein
19       LOGICAL ok_suntime_rrtm
20       REAL(kind=8) RCO2, RCH4, RN2O, RCFC11, RCFC12 
21       REAL(kind=8) RCO2_act, RCH4_act, RN2O_act, RCFC11_act, RCFC12_act 
22       REAL(kind=8) CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt
23!IM ajout CFMIP2/CMIP5
24       REAL(kind=8) RCO2_per,RCH4_per,RN2O_per,RCFC11_per,RCFC12_per
25       REAL(kind=8) CH4_ppb_per,N2O_ppb_per,CFC11_ppt_per,CFC12_ppt_per
26
27!OM ---> correction du bilan d'eau global
28!OM Correction sur precip KE
29       REAL cvl_corr
30!OM Fonte calotte dans bilan eau
31       LOGICAL ok_lic_melt
32!OB Depot de vapeur d eau sur la calotte pour le bilan eau
33       LOGICAL ok_lic_cond
34
35!IM simulateur ISCCP
36       INTEGER top_height, overlap
37!IM seuils cdrm, cdrh
38       REAL cdmmax, cdhmax
39!IM param. stabilite s/ terres et en dehors
40       REAL ksta, ksta_ter, f_ri_cd_min
41!IM ok_kzmin : clef calcul Kzmin dans la CL de surface cf FH
42       LOGICAL ok_kzmin
43!IM, MAFo fmagic, pmagic : parametres - additionnel et multiplicatif - 
44!                          pour regler l albedo sur ocean
45       REAL pbl_lmixmin_alpha
46       REAL fmagic, pmagic
47! Hauteur (imposee) du contenu en eau du sol
48           REAL qsol0,albsno0,evap0
49! Frottement au sol (Cdrag)
50       Real f_cdrag_ter,f_cdrag_oce
51       REAL min_wind_speed,f_gust_wk,f_gust_bl,f_qsat_oce,f_z0qh_oce
52       REAL z0m_seaice,z0h_seaice
53       INTEGER iflag_gusts,iflag_z0_oce
54
55! Rugoro
56       Real f_rugoro,z0min
57
58! tau_gl : constante de rappel de la temperature a la surface de la glace
59       REAL tau_gl
60
61!IM lev_histhf  : niveau sorties 6h
62!IM lev_histday : niveau sorties journalieres
63!IM lev_histmth : niveau sorties mensuelles
64!IM lev_histdayNMC : on peut sortir soit sur 8 (comme AR5) ou bien
65!                    sur 17 niveaux de pression
66       INTEGER lev_histhf, lev_histday, lev_histmth
67       INTEGER lev_histdayNMC
68       Integer lev_histins, lev_histLES 
69!IM ok_histNMC  : sortie fichiers niveaux de pression (histmthNMC, histdayNMC, histhfNMC)
70!IM freq_outNMC : frequences de sortie fichiers niveaux de pression (histmthNMC, histdayNMC, histhfNMC)
71!IM freq_calNMC : frequences de calcul fis. hist*NMC.nc
72       LOGICAL ok_histNMC(3)
73       INTEGER levout_histNMC(3)
74       REAL freq_outNMC(3) , freq_calNMC(3)
75       CHARACTER(len=4) type_run
76! aer_type: pour utiliser un fichier constant dans readaerosol
77       CHARACTER(len=8) :: aer_type
78       LOGICAL ok_regdyn
79       REAL lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
80       REAL ecrit_ins, ecrit_hf, ecrit_day
81       REAL ecrit_mth, ecrit_tra, ecrit_reg
82       REAL ecrit_LES
83       REAL freq_ISCCP, ecrit_ISCCP
84       REAL freq_COSP, freq_AIRS
85       LOGICAL :: ok_cosp,ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP
86       LOGICAL :: ok_airs
87       INTEGER :: ip_ebil_phy, iflag_rrtm, iflag_ice_thermo, NSW, iflag_albedo
88       LOGICAL :: ok_chlorophyll
89       LOGICAL :: ok_strato
90       LOGICAL :: ok_hines, ok_gwd_rando
91       LOGICAL :: ok_qch4
92       LOGICAL :: ok_conserv_q
93       LOGICAL :: adjust_tropopause
94       LOGICAL :: ok_daily_climoz
95! flag to bypass or not the phytrac module
96       INTEGER :: iflag_phytrac
97
98       COMMON/clesphys/                                                 &
99! REAL FIRST
100     &       co2_ppm, solaire                                           &
101     &     , RCO2, RCH4, RN2O, RCFC11, RCFC12                           &
102     &     , RCO2_act, RCH4_act, RN2O_act, RCFC11_act, RCFC12_act       &
103     &     , RCO2_per, RCH4_per, RN2O_per, RCFC11_per, RCFC12_per       &
104     &     , CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt                     &
105     &     , CH4_ppb_per, N2O_ppb_per, CFC11_ppt_per, CFC12_ppt_per     &
106     &     , cdmmax,cdhmax,ksta,ksta_ter,f_ri_cd_min,pbl_lmixmin_alpha  &
107     &     , fmagic, pmagic                                             &
108     &     , f_cdrag_ter,f_cdrag_oce,f_rugoro,z0min,tau_gl              &
109     &     , min_wind_speed,f_gust_wk,f_gust_bl,f_qsat_oce,f_z0qh_oce   &
110     &     , z0m_seaice,z0h_seaice                                      &
111     &     , freq_outNMC, freq_calNMC                                   &
112     &     , lonmin_ins, lonmax_ins, latmin_ins, latmax_ins             &
113     &     , freq_ISCCP, ecrit_ISCCP, freq_COSP, freq_AIRS              &
114     &     , cvl_corr                                                   &
115     &     , qsol0,albsno0,evap0                                        &
116     &     , co2_ppm0                                                   &
117!FC
118     &     , Cd_frein                                                   &
119     &     , ecrit_LES                                                  &
120     &     , ecrit_ins, ecrit_hf, ecrit_day                             &
121     &     , ecrit_mth, ecrit_tra, ecrit_reg                            &
122! THEN INTEGER AND LOGICALS
123     &     , top_height                                                 &
124     &     , iflag_cycle_diurne, soil_model, new_oliq                         &
125     &     , ok_orodr, ok_orolf, ok_limitvrai, nbapp_rad                &
126     &     , iflag_con, nbapp_cv, nbapp_wk                              &
127     &     , iflag_ener_conserv                                         &
128     &     , ok_suntime_rrtm                                            & 
129     &     , overlap                                                    &
130     &     , ok_kzmin                                                   &
131     &     , lev_histhf, lev_histday, lev_histmth                       &
132     &     , lev_histins, lev_histLES, lev_histdayNMC, levout_histNMC   &
133     &     , ok_histNMC                                                 &
134     &     , type_run, ok_regdyn, ok_cosp, ok_airs                      &
135     &     , ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP                     &
136     &     , ip_ebil_phy                                                &
137     &     , iflag_gusts ,iflag_z0_oce                                  &
138     &     , ok_lic_melt, ok_lic_cond, aer_type                         &
139     &     , iflag_rrtm, ok_strato,ok_hines, ok_qch4                    &
140     &     , iflag_ice_thermo, ok_gwd_rando, NSW, iflag_albedo          &
141     &     , ok_chlorophyll,ok_conserv_q, adjust_tropopause             &
142     &     , ok_daily_climoz, ok_all_xml, ok_lwoff                      &
143     &     , iflag_phytrac
144     
145       save /clesphys/
146!$OMP THREADPRIVATE(/clesphys/)
Note: See TracBrowser for help on using the repository browser.