source: LMDZ6/trunk/libf/phylmd/clesphys.h @ 3816

Last change on this file since 3816 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
RevLine 
[1279]1! $Id: clesphys.h 3435 2019-01-22 15:21:59Z lguez $
[524]2!
[793]3!  ATTENTION!!!!: ce fichier include est compatible format fixe/format libre
[2110]4!                 veillez \`a n'utiliser que des ! pour les commentaires
5!                 et \`a bien positionner les & des lignes de continuation
[793]6!                 (les placer en colonne 6 et en colonne 73)
7!
8!..include cles_phys.h
9!
[3317]10       INTEGER iflag_cycle_diurne
11       LOGICAL soil_model,new_oliq,ok_orodr,ok_orolf
[524]12       LOGICAL ok_limitvrai
[2114]13       LOGICAL ok_all_xml
[3048]14       LOGICAL ok_lwoff
[2730]15       INTEGER nbapp_rad, iflag_con, nbapp_cv, nbapp_wk, iflag_ener_conserv
[1279]16       REAL co2_ppm, co2_ppm0, solaire
[2952]17!FC
18       REAL Cd_frein
[2524]19       LOGICAL ok_suntime_rrtm
[1279]20       REAL(kind=8) RCO2, RCH4, RN2O, RCFC11, RCFC12 
[1539]21       REAL(kind=8) RCO2_act, RCH4_act, RN2O_act, RCFC11_act, RCFC12_act 
[1279]22       REAL(kind=8) CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt
[1539]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
[766]26
[793]27!OM ---> correction du bilan d'eau global
28!OM Correction sur precip KE
[766]29       REAL cvl_corr
[793]30!OM Fonte calotte dans bilan eau
[766]31       LOGICAL ok_lic_melt
[2946]32!OB Depot de vapeur d eau sur la calotte pour le bilan eau
33       LOGICAL ok_lic_cond
[766]34
[793]35!IM simulateur ISCCP
[524]36       INTEGER top_height, overlap
[793]37!IM seuils cdrm, cdrh
[524]38       REAL cdmmax, cdhmax
[793]39!IM param. stabilite s/ terres et en dehors
[2126]40       REAL ksta, ksta_ter, f_ri_cd_min
[793]41!IM ok_kzmin : clef calcul Kzmin dans la CL de surface cf FH
[524]42       LOGICAL ok_kzmin
[900]43!IM, MAFo fmagic, pmagic : parametres - additionnel et multiplicatif - 
44!                          pour regler l albedo sur ocean
[2561]45       REAL pbl_lmixmin_alpha
[900]46       REAL fmagic, pmagic
[889]47! Hauteur (imposee) du contenu en eau du sol
[1894]48           REAL qsol0,albsno0,evap0
[1279]49! Frottement au sol (Cdrag)
50       Real f_cdrag_ter,f_cdrag_oce
[2254]51       REAL min_wind_speed,f_gust_wk,f_gust_bl,f_qsat_oce,f_z0qh_oce
[2243]52       REAL z0m_seaice,z0h_seaice
53       INTEGER iflag_gusts,iflag_z0_oce
[2240]54
[1279]55! Rugoro
[2243]56       Real f_rugoro,z0min
[1279]57
[3327]58! tau_gl : constante de rappel de la temperature a la surface de la glace
59       REAL tau_gl
60
[793]61!IM lev_histhf  : niveau sorties 6h
62!IM lev_histday : niveau sorties journalieres
63!IM lev_histmth : niveau sorties mensuelles
[1352]64!IM lev_histdayNMC : on peut sortir soit sur 8 (comme AR5) ou bien
65!                    sur 17 niveaux de pression
[524]66       INTEGER lev_histhf, lev_histday, lev_histmth
[1352]67       INTEGER lev_histdayNMC
[1279]68       Integer lev_histins, lev_histLES 
[1374]69!IM ok_histNMC  : sortie fichiers niveaux de pression (histmthNMC, histdayNMC, histhfNMC)
[1352]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
[1374]72       LOGICAL ok_histNMC(3)
[1828]73       INTEGER levout_histNMC(3)
[1352]74       REAL freq_outNMC(3) , freq_calNMC(3)
[1279]75       CHARACTER(len=4) type_run
76! aer_type: pour utiliser un fichier constant dans readaerosol
[1938]77       CHARACTER(len=8) :: aer_type
[1877]78       LOGICAL ok_regdyn
[644]79       REAL lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
[1563]80       REAL ecrit_ins, ecrit_hf, ecrit_day
[684]81       REAL ecrit_mth, ecrit_tra, ecrit_reg
[1054]82       REAL ecrit_LES
[828]83       REAL freq_ISCCP, ecrit_ISCCP
[2580]84       REAL freq_COSP, freq_AIRS
[1368]85       LOGICAL :: ok_cosp,ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP
[2580]86       LOGICAL :: ok_airs
[2227]87       INTEGER :: ip_ebil_phy, iflag_rrtm, iflag_ice_thermo, NSW, iflag_albedo
88       LOGICAL :: ok_chlorophyll
[1001]89       LOGICAL :: ok_strato
[1938]90       LOGICAL :: ok_hines, ok_gwd_rando
[2136]91       LOGICAL :: ok_qch4
[2007]92       LOGICAL :: ok_conserv_q
[2788]93       LOGICAL :: adjust_tropopause
94       LOGICAL :: ok_daily_climoz
[2973]95! flag to bypass or not the phytrac module
96       INTEGER :: iflag_phytrac
[524]97
[1961]98       COMMON/clesphys/                                                 &
99! REAL FIRST
[2007]100     &       co2_ppm, solaire                                           &
[1539]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       &
[793]104     &     , CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt                     &
[1539]105     &     , CH4_ppb_per, N2O_ppb_per, CFC11_ppt_per, CFC12_ppt_per     &
[2561]106     &     , cdmmax,cdhmax,ksta,ksta_ter,f_ri_cd_min,pbl_lmixmin_alpha  &
[1961]107     &     , fmagic, pmagic                                             &
[3327]108     &     , f_cdrag_ter,f_cdrag_oce,f_rugoro,z0min,tau_gl              &
[2254]109     &     , min_wind_speed,f_gust_wk,f_gust_bl,f_qsat_oce,f_z0qh_oce   &
[2243]110     &     , z0m_seaice,z0h_seaice                                      &
[3435]111     &     , freq_outNMC, freq_calNMC                                   &
[1961]112     &     , lonmin_ins, lonmax_ins, latmin_ins, latmax_ins             &
[2580]113     &     , freq_ISCCP, ecrit_ISCCP, freq_COSP, freq_AIRS              &
[1961]114     &     , cvl_corr                                                   &
115     &     , qsol0,albsno0,evap0                                        &
116     &     , co2_ppm0                                                   &
[2952]117!FC
118     &     , Cd_frein                                                   &
[1961]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                                                 &
[3317]124     &     , iflag_cycle_diurne, soil_model, new_oliq                         &
[2707]125     &     , ok_orodr, ok_orolf, ok_limitvrai, nbapp_rad                &
[2730]126     &     , iflag_con, nbapp_cv, nbapp_wk                              &
[1961]127     &     , iflag_ener_conserv                                         &
[2524]128     &     , ok_suntime_rrtm                                            & 
[1961]129     &     , overlap                                                    &
130     &     , ok_kzmin                                                   &
[900]131     &     , lev_histhf, lev_histday, lev_histmth                       &
[1828]132     &     , lev_histins, lev_histLES, lev_histdayNMC, levout_histNMC   &
[1961]133     &     , ok_histNMC                                                 &
[2580]134     &     , type_run, ok_regdyn, ok_cosp, ok_airs                      &
[1368]135     &     , ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP                     &
[1961]136     &     , ip_ebil_phy                                                &
[2243]137     &     , iflag_gusts ,iflag_z0_oce                                  &
[2946]138     &     , ok_lic_melt, ok_lic_cond, aer_type                         &
[2136]139     &     , iflag_rrtm, ok_strato,ok_hines, ok_qch4                    &
[2227]140     &     , iflag_ice_thermo, ok_gwd_rando, NSW, iflag_albedo          &
[2788]141     &     , ok_chlorophyll,ok_conserv_q, adjust_tropopause             &
[3048]142     &     , ok_daily_climoz, ok_all_xml, ok_lwoff                      &
[2973]143     &     , iflag_phytrac
[766]144     
[1572]145       save /clesphys/
[793]146!$OMP THREADPRIVATE(/clesphys/)
Note: See TracBrowser for help on using the repository browser.