Changeset 255


Ignore:
Timestamp:
02/19/08 10:47:01 (16 years ago)
Author:
sdipsl
Message:

SD : up to date commented physiq.def regarding the LMD4_V3_? IO
handling. Set ISCCP input frequency and isccp average frequency
according to http://cfmip.metoffice.com/ISCCP.html

ISCCP is off by default. Wait for bugfix in LMDZ4_V3_? concerning ISCCP
simulator.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • CONFIG/trunk/IPSLCM4_v2/EXP00/PARAM/physiq.def

    r217 r255  
    22## $Id$ 
    33# 
     4### OCEAN = force / slab  /couple 
     5OCEAN=couple 
     6### VEGET= y si ORCHIDEE, =n si bucket 
     7VEGET=y 
    48# 
    5 # Automatically generated make config: don't edit 
     9# Parametres fichiers de sortie 
    610# 
    7  
    8 OCEAN=couple 
    9 VEGET=y 
     11### type_run = type run par rapport aux fichiers et variables de sortie 
     12# - type_run = CLIM/ENSP (=1) 
     13# - type_run = AMIP/CFMI (=2) 
     14type_run=CLIM 
     15### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc 
    1016OK_journe=n 
     17### OK_mensuel= y sortir fichier mensuel histmth.nc, =n pas de fichier histmth.nc 
    1118OK_mensuel=y 
     19### OK_instan=y, ecrire sorties "instantannees" (chaque pas de temps de la  physique) 
    1220OK_instan=n 
     21### OK_hf=y, ecrire sorties hautes frequence histhf.nc, =n pas de fichier histhf.nc 
     22ok_hf=n 
     23# 
     24# Parametres niveau de sorties differents fichiers  
     25# 
     26### lev_histhf=0-4, niveau de sortie fichier "histhf.nc"  
     27# - lev_histhf=0 => pas de sorties histhf.nc 
     28# - lev_histhf=2 => defaut 
     29# - lev_histhf=3 => variables sur niveaux standards 
     30# - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele => fichier. histhf3d.nc 
     31lev_histhf=2 
     32### lev_histday=0-5, niveau de sortie fichier "histday.nc" 
     33# - lev_histday=0 => pas de sorties lev_histday.nc 
     34# - lev_histday=2 => defaut 
     35# - lev_histday=3 => + champs 3D => F. Lott 
     36# - lev_histday=4 => + champs sous-surfaces 
     37# - lev_histday=5 => + champs F. Aires 
     38lev_histday=2 
     39### lev_histmth=0-4, niveau de sortie fichier "histmth.nc" 
     40# - lev_histmth=0 => pas de sorties lev_histmth.nc 
     41# - lev_histmth=2 => defaut 
     42# - lev_histmth=3 => albedo, rugosite sous-surfaces 
     43# - lev_histmth=4 => champs tendances 3d 
     44lev_histmth=2 
     45### ecrit_hf = frequence ecriture fichier histhf.nc en jours 
     46ecrit_hf=0.250 
     47### ecrit_day = frequence ecriture fichier histday.nc en jours 
     48ecrit_day=1. 
     49### ecrit_mth = frequence ecriture fichier histmth.nc en jours 
     50ecrit_mth=30. 
     51### freqin_isccp = frequence input en secondes du simulateur ISCCP 
     52freq_ISCCP=10800. 
     53### freqout_isccp = frequence output en jours du simulateur ISCCP 
     54ecrit_ISCCP=30. 
     55### niveau du diagnostique de conservation d energie 
    1356if_ebil=0 
    1457# 
    1558# parametres KE 
    1659# 
     60### epmax = Efficacite precipitation maximale 
    1761epmax = .99 
     62### ok_adj_ema = ?? pas utilise 
    1863ok_adj_ema = n 
     64### iflag_clw Flag calcul eau liquide 
     65# - iflag_clw=0 : qcond_incld(i,l) = em_qcondc(l) 
     66# - iflag_clw=1 : qcond_incld(i,l) = em_qcond(l) 
     67# - iflag_clw=2 : eau liquide diagnostique en fonction de la Precip 
    1968iflag_clw = 1 
    2069#  
    2170# parametres nuages 
    2271# 
     72### cld_lc_lsc  contenu en eau liquide des nuages large-scale (fisrtilp)  
    2373cld_lc_lsc = 2.6e-4 
     74### cld_lc_con  contenu en eau liquide des nuages convectifs (fisrtilp)  
    2475cld_lc_con = 2.6e-4 
     76### cld_tau_lsc  cte de temps utilisee pour eliminer l eau large-scale (fisrtilp)  
    2577cld_tau_lsc = 3600. 
     78### cld_tau_con  cte de temps utilisee pour eliminer l eau convective (fisrtilp)  
    2679cld_tau_con = 3600. 
     80### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace large-scale (fisrtilp)  
    2781ffallv_lsc = 1. 
     82### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace convectifs (fisrtilp)  
    2883ffallv_con = 1. 
     84### coef_eva  coef evaporation precips eau/glace (fisrtilp/fisrtilp_tr?/conlmd?)  
    2985coef_eva = 2.e-5 
     86### reevap_ice  reevaporation de toute la precip dans la couche du dessous pour la glace (fisrtilp)  
    3087reevap_ice = y 
     88### iflag_cldcon  flag pour calculer ratqsc=F(ratqsbas,fact_cldcon,q_seri) (physiq)  
     89# - iflag_cldcon<=-1 diag. rain_Tiedtke  
     90# - iflag_cldcon=1, ratqsc=ratqsbas+fact_cldcon*(q_seri(1)-q_seri(k))/q_seri(k) 
     91# - iflag_cldcon=1/2, ratqs=max(ratqs,ratqsc) 
     92# - iflag_cldcon=3,   ratqs=ratqss 
    3193iflag_cldcon = 3 
     94### iflag_pdf :  flag calcul distribution sous-maille de l eau et des nuages  
     95# - iflag_pdf=0, version ratqs,  
     96# - iflag_pdf=1, calcul eau condensee, fraction nuageuse, eau nuageuse a partir 
     97# -              des PDFS de Sandrine Bony  
    3298iflag_pdf = 1 
     99### fact_cldcon  constante calcul ratqsc (voir iflag_cldcon) et proprietes nuages convectifs, clwcon0 (physiq.F)  
    33100fact_cldcon = 1. 
     101### facttemps=   facteur de relaxation de ratqs (iflag_cldcon=1/2) et rnebcon (iflag_cldcon=3) 
    34102facttemps = 1.e-4 
     103## ok_newmicro   =y appel newmicro , =n appel nuage (calcul epaisseur optique et emmissivite des nuages)  
    35104ok_newmicro = y 
     105### ratqsbas     ratqs en bas si iflag_cldcon=1  
    36106ratqsbas = 0.005 
     107### ratqshaut    ratqs en haut pour ratqss "stables"   
    37108ratqshaut = 0.33 
     109### rad_froid    rayon cristaux des nuages de glace "froids"  
    38110rad_froid = 35 
     111### rad_chau1    rayon goutelettes d eau chauds", en haut: k=4-klev  
    39112rad_chau1=12 
     113### rad_chau2    rayon goutelettes d eau chauds", en bas: k=1-3   
    40114rad_chau2=11 
     115## ksta_ter      coef.diffusion minimale sur terre/sic/lic  
    41116ksta_ter=1.e-7 
     117### cdmmax =     cdrag maximum pour le moment 
     118cdmmax = 2.5E-3 
     119### cdhmax =     cdrag maximum pour l energie 
     120cdhmax = 2.0E-3 
    42121# 
    43 # parametres outputs 
     122# Parametres "orbitaux/ ere geologique" 
    44123# 
    45 # Flag pour hf : ok_hf=y/n 
    46 ok_hf=n 
    47 # Flag pour les variables de sorties : type_run=CLIM/AMIP 
    48 # Pour ENSEMBLES mettre type_run=CLIM 
    49 type_run=AMIP 
    50 #niveau de sortie "hf" lev_histhf avec 
    51 # - lev_histhf=3 => defaut 
    52 # - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele 
    53 lev_histhf=3 
    54 #niveau de sortie "day" lev_histday 
    55 # - lev_histday=2 => defaut 
    56 # - lev_histday=3 => + champs 3D => F. Lott 
    57 # - lev_histday=4 => + champs sous-surfaces 
    58 lev_histday=2 
    59 #niveau de sortie "mth" lev_histmth avec 
    60 # - lev_histmth=2 => defaut 
    61 # - lev_histmth=3 => albedo, rugosite sous-surfaces 
    62 # - lev_histmth=4 => champs tendances 3d 
    63 lev_histmth=2 
     124### R_ecc =      Excentricite 
     125R_ecc = 0.016715 
     126### R_peri =     Equinoxe 
     127R_peri = 102.7 
     128### R_incl =     Inclinaison 
     129R_incl = 23.441 
     130### solaire =    Constante solaire 
     131solaire = 1365. 
    64132# 
    65 # parametres climatique 
     133# Taux gaz a effet de serre 
    66134# 
    67 R_ecc = 0.016715 
    68 R_peri = 102.7 
    69 R_incl = 23.441 
    70 solaire = 1365. 
     135### co2_ppm =    taux CO2 en ppm 
    71136co2_ppm = 348. 
    72 #RCO2 = co2_ppm * 1.0e-06  * 44.011/28.97 
    73 #RCO2 = 348. * 1.0e-06  * 44.011/28.97 
    74 #RCO2 =   5.286789092164308E-04 
    75 #RCO2 = 425.43e-06 
     137### CH4_ppb =    taux CH4 en ppb 
    76138CH4_ppb = 1650. 
    77 #RCH4 = 1.65E-06* 16.043/28.97 
    78 #RCH4 =   9.137366240938903E-07 
     139### N2O_ppb =    taux N2O en ppb 
    79140N2O_ppb = 306. 
    80 #RN2O = 306.E-09* 44.013/28.97 
    81 #RN2O =    4.648939592682085E-07 
     141### CFC11_ppt =  taux CFC11 en ppt 
    82142CFC11_ppt = 280. 
    83 #RCFC11 = 280.E-12* 137.3686/28.97 
    84 #RCFC11 =    1.327690990680013E-09 
     143### CFC12_ppt =  taux CFC12 en ppt 
    85144CFC12_ppt = 484. 
    86 #RCFC12 = 484.E-12* 120.9140/28.97 
    87 #RCFC12 =    2.020102726958923E-09 
    88145# 
    89 # parametres simulateur ISCCP 
     146# Parametres effets directs/indirects des "aerosols" 
    90147# 
     148### ok_ade=y/n   flag Aerosol direct effect 
     149ok_ade=y 
     150### ok_aie=y/n   flag Aerosol indirect effect 
     151ok_aie=y 
     152### aer_type =   Aerosol variation type : actuel / preind / scenario 
     153aer_type=actuel 
     154### bl95_b0 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995) 
     155bl95_b0=1.7 
     156### bl95_b1 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995) 
     157bl95_b1=0.2 
     158# 
     159# Parametres simulateur ISCCP 
     160# 
     161### ok_isccp=y/n flag simulateur ISCCP 
    91162ok_isccp=n 
    92 #top_height = 1 ou 3 
    93 top_height = 3 
    94 #overlap = 1, 2 ou 3 
     163### top_height = flag choix calcul nuages par le simulateur en utilisant  
     164# -              les donnees IR et/ou VIS et l algorithme ISCCP-D1 
     165# - top_height = 1 -> algo IR-VIS 
     166# - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)" 
     167# - top_height = 3 -> algo IR 
     168top_height = 1 
     169### overlap =    Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP 
     170# - overlap=1    Max overlap 
     171# - overlap=2    Random overlap 
     172# - overlap=3    Max/Random overlap 
    95173overlap = 3 
    96 #cdmmax 
    97 cdmmax = 2.5E-3 
    98 #cdhmax 
    99 cdhmax = 2.0E-3 
    100 # effets direct et indirect des aérosols 
    101 ok_ade=y 
    102 ok_aie=y 
    103 aer_type=actuel 
    104 bl95_b0=1.7 
    105 bl95_b1=0.2 
    106  
Note: See TracChangeset for help on using the changeset viewer.