source: LMDZ5/trunk/libf/phylmd/clesphys.h @ 1822

Last change on this file since 1822 was 1753, checked in by idelkadi, 12 years ago

Concerns energy conservation.


The source terms of the TKE prognostic equation are diagnosed
from the tendencies (du, dv, dT) associated with subrgrid scale
motions and treated as an additional heat source.
Controled by a new key, iflag_ener_conserv :

0 no conservation

-1 old adhoc correction for kinetic E only (used for CMIP5)

1 conservation
101 conversion from kinetic to heat only
110 conversion from potential to heat only

iflag_ener_conserv=-1 kept as default value for a test period
iflag_ener_conserv=1 is the advisable value
Concerns clesphys.h, and conf_phys.F90
New routine : ener_conserv.F90, called by physic.
New outputs :
bils_ec, contribution to the energy budget of the column of the

additional heat source (in W/m2)

bils_kinetic : change kinetic energy of the column in physics (W/m2)
bils_enthalp : idem for the total column enthalphy
bils_latent : idem for latent heat
Modified files : clesphys.h, conf_phys_m.F90, physiq.F,

phys_output_mod.F90, phys_output_var_mod.F90, phys_output_write.h,
ener_conserv.F90

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.7 KB
RevLine 
[1403]1
[524]2!
[1279]3! $Id: clesphys.h 1753 2013-05-13 12:53:28Z emillour $
[524]4!
[793]5!  ATTENTION!!!!: ce fichier include est compatible format fixe/format libre
6!                 veillez à n'utiliser que des ! pour les commentaires
7!                 et à bien positionner les & des lignes de continuation
8!                 (les placer en colonne 6 et en colonne 73)
9!
10!..include cles_phys.h
11!
[524]12       LOGICAL cycle_diurne,soil_model,new_oliq,ok_orodr,ok_orolf
13       LOGICAL ok_limitvrai
[1753]14       INTEGER nbapp_rad, iflag_con,iflag_ener_conserv
[1279]15       REAL co2_ppm, co2_ppm0, solaire
16       REAL(kind=8) RCO2, RCH4, RN2O, RCFC11, RCFC12 
[1539]17       REAL(kind=8) RCO2_act, RCH4_act, RN2O_act, RCFC11_act, RCFC12_act 
[1279]18       REAL(kind=8) CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt
[1539]19!IM ajout CFMIP2/CMIP5
20       REAL co2_ppm_per
[1737]21       LOGICAL ok_4xCO2atm
[1539]22       REAL(kind=8) RCO2_per,RCH4_per,RN2O_per,RCFC11_per,RCFC12_per
23       REAL(kind=8) CH4_ppb_per,N2O_ppb_per,CFC11_ppt_per,CFC12_ppt_per
[766]24
[793]25!OM ---> correction du bilan d'eau global
26!OM Correction sur precip KE
[766]27       REAL cvl_corr
[793]28!OM Fonte calotte dans bilan eau
[766]29       LOGICAL ok_lic_melt
30
[793]31!IM simulateur ISCCP
[524]32       INTEGER top_height, overlap
[793]33!IM seuils cdrm, cdrh
[524]34       REAL cdmmax, cdhmax
[793]35!IM param. stabilite s/ terres et en dehors
[524]36       REAL ksta, ksta_ter
[793]37!IM ok_kzmin : clef calcul Kzmin dans la CL de surface cf FH
[524]38       LOGICAL ok_kzmin
[900]39!IM, MAFo fmagic, pmagic : parametres - additionnel et multiplicatif - 
40!                          pour regler l albedo sur ocean
41       REAL fmagic, pmagic
[889]42! Hauteur (imposee) du contenu en eau du sol
43           REAL qsol0
[1279]44! Frottement au sol (Cdrag)
45       Real f_cdrag_ter,f_cdrag_oce
46! Rugoro
47       Real f_rugoro
48
[793]49!IM lev_histhf  : niveau sorties 6h
50!IM lev_histday : niveau sorties journalieres
51!IM lev_histmth : niveau sorties mensuelles
[1352]52!IM lev_histdayNMC : on peut sortir soit sur 8 (comme AR5) ou bien
53!                    sur 17 niveaux de pression
[524]54       INTEGER lev_histhf, lev_histday, lev_histmth
[1352]55       INTEGER lev_histdayNMC
[1279]56       Integer lev_histins, lev_histLES 
[1374]57!IM ok_histNMC  : sortie fichiers niveaux de pression (histmthNMC, histdayNMC, histhfNMC)
[1352]58!IM freq_outNMC : frequences de sortie fichiers niveaux de pression (histmthNMC, histdayNMC, histhfNMC)
59!IM freq_calNMC : frequences de calcul fis. hist*NMC.nc
60!IM pasphys : pas de temps de physique (secondes)
61       REAL pasphys
[1374]62       LOGICAL ok_histNMC(3)
[1352]63       REAL freq_outNMC(3) , freq_calNMC(3)
[1279]64       CHARACTER(len=4) type_run
65! aer_type: pour utiliser un fichier constant dans readaerosol
[1146]66       CHARACTER*8 :: aer_type
[644]67       LOGICAL ok_isccp, ok_regdyn
68       REAL lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
[1563]69       REAL ecrit_ins, ecrit_hf, ecrit_day
[684]70       REAL ecrit_mth, ecrit_tra, ecrit_reg
[1054]71       REAL ecrit_LES
[828]72       REAL freq_ISCCP, ecrit_ISCCP
[1279]73       REAL freq_COSP
[1368]74       LOGICAL :: ok_cosp,ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP
[998]75       INTEGER :: ip_ebil_phy, iflag_rrtm
[1001]76       LOGICAL :: ok_strato
77       LOGICAL :: ok_hines
[524]78
[793]79       COMMON/clesphys/cycle_diurne, soil_model, new_oliq,              &
80     &     ok_orodr, ok_orolf, ok_limitvrai, nbapp_rad, iflag_con       &
[1753]81     &     , iflag_ener_conserv                                         &
[1539]82     &     , co2_ppm, solaire                                           &
83     &     , RCO2, RCH4, RN2O, RCFC11, RCFC12                           &
84     &     , RCO2_act, RCH4_act, RN2O_act, RCFC11_act, RCFC12_act       &
[1737]85     &     , ok_4xCO2atm                                                & 
[1539]86     &     , RCO2_per, RCH4_per, RN2O_per, RCFC11_per, RCFC12_per       &
[793]87     &     , CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt                     &
[1539]88     &     , CH4_ppb_per, N2O_ppb_per, CFC11_ppt_per, CFC12_ppt_per     &
[793]89     &     , top_height, overlap, cdmmax, cdhmax, ksta, ksta_ter        &
[900]90     &     , ok_kzmin, fmagic, pmagic                                   &
[1279]91     &     , f_cdrag_ter,f_cdrag_oce,f_rugoro                           &
[900]92     &     , lev_histhf, lev_histday, lev_histmth                       &
[1352]93     &     , lev_histins, lev_histLES, lev_histdayNMC                   &
[1374]94     &     , pasphys, ok_histNMC, freq_outNMC, freq_calNMC              &
[1279]95     &     , type_run, ok_isccp, ok_regdyn, ok_cosp                     &
[1368]96     &     , ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP                     &
[793]97     &     , lonmin_ins, lonmax_ins, latmin_ins, latmax_ins             &
[1564]98     &     , ecrit_ins, ecrit_hf, ecrit_day                             &
[793]99     &     , ecrit_mth, ecrit_tra, ecrit_reg                            &
[1279]100     &     , freq_ISCCP, ecrit_ISCCP, freq_COSP, ip_ebil_phy            &
[1146]101     &     , ok_lic_melt, cvl_corr, aer_type                            &
[1279]102     &     , qsol0, iflag_rrtm, ok_strato,ok_hines,ecrit_LES            &
[1421]103     &     , co2_ppm0
[766]104     
[1572]105       save /clesphys/
[793]106!$OMP THREADPRIVATE(/clesphys/)
[766]107 
Note: See TracBrowser for help on using the repository browser.