Ignore:
Timestamp:
Feb 17, 2009, 3:44:06 PM (15 years ago)
Author:
lguez
Message:

"filtrez/filtreg.F" works with either "matmul" or a direct call to
BLAS.
"gcm.def", "physiq.def" and "run.def" adapted from
"http://forge.ipsl.jussieu.fr/igcmg/browser/CONFIG/LMDZOR/branches/LMDZ4OR_v3/LMDZOR/PARAM".

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4-dev/physiq.def

    r541 r1101  
     1### type_ocean = force / slab  /couple
     2type_ocean=force
     3### VEGET= y si ORCHIDEE, =n si bucket
     4VEGET=n
    15#
    2 # $Header$
     6# Parametres fichiers de sortie
    37#
    4 
    5 OCEAN=force
    6 VEGET=y
    7 OK_journe=n
     8### type_run = type run par rapport aux fichiers et variables de sortie
     9# - type_run = CLIM/ENSP (=1)
     10# - type_run = AMIP/CFMI (=2)
     11type_run=AMIP
     12### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc
     13OK_journe=y
     14### OK_mensuel= y sortir fichier mensuel histmth.nc, =n pas de fichier histmth.nc
    815OK_mensuel=y
     16### OK_instan=y, ecrire sorties "instantannees" (chaque pas de temps de la  physique)
    917OK_instan=n
     18### OK_hf=y, ecrire sorties hautes frequence histhf.nc, =n pas de fichier histhf.nc
     19ok_hf=n
     20#
     21# Parametres niveau de sorties differents fichiers
     22#
     23### lev_histhf=0-4, niveau de sortie fichier "histhf.nc"
     24# - lev_histhf=0 => pas de sorties histhf.nc
     25# - lev_histhf=2 => defaut
     26# - lev_histhf=3 => variables sur niveaux standards
     27# - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele => fichier. histhf3d.nc
     28lev_histhf=2
     29### lev_histday=0-5, niveau de sortie fichier "histday.nc"
     30# - lev_histday=0 => pas de sorties lev_histday.nc
     31# - lev_histday=2 => defaut
     32# - lev_histday=3 => + champs 3D => F. Lott
     33# - lev_histday=4 => + champs sous-surfaces
     34# - lev_histday=5 => + champs F. Aires
     35lev_histday=2
     36### lev_histmth=0-4, niveau de sortie fichier "histmth.nc"
     37# - lev_histmth=0 => pas de sorties lev_histmth.nc
     38# - lev_histmth=2 => defaut
     39# - lev_histmth=3 => albedo, rugosite sous-surfaces
     40# - lev_histmth=4 => champs tendances 3d
     41lev_histmth=2
     42### ecrit_hf = frequence ecriture fichier histhf.nc en jours
     43ecrit_hf=0.250
     44### ecrit_day = frequence ecriture fichier histday.nc en jours
     45ecrit_day=1.
     46### ecrit_mth = frequence ecriture fichier histmth.nc en jours
     47ecrit_mth=30.
     48### freqin_isccp = frequence input en secondes du simulateur ISCCP
     49freq_ISCCP=10800.
     50### freqout_isccp = frequence output en jours du simulateur ISCCP
     51ecrit_ISCCP=30.
     52### niveau du diagnostique de conservation d energie
    1053if_ebil=0
    1154#
    1255# parametres KE
    1356#
     57### epmax = Efficacite precipitation maximale
    1458epmax = .99
     59### ok_adj_ema = ?? pas utilise
    1560ok_adj_ema = n
     61### iflag_clw Flag calcul eau liquide
     62# - iflag_clw=0 : qcond_incld(i,l) = em_qcondc(l)
     63# - iflag_clw=1 : qcond_incld(i,l) = em_qcond(l)
     64# - iflag_clw=2 : eau liquide diagnostique en fonction de la Precip
    1665iflag_clw = 1
    1766#
    1867# parametres nuages
    1968#
    20 cld_lc_lsc = 0.00026
    21 cld_lc_con = 0.00026
     69### cld_lc_lsc  contenu en eau liquide des nuages large-scale (fisrtilp)
     70cld_lc_lsc = 2.6e-4
     71### cld_lc_con  contenu en eau liquide des nuages convectifs (fisrtilp)
     72cld_lc_con = 2.6e-4
     73### cld_tau_lsc  cte de temps utilisee pour eliminer l eau large-scale (fisrtilp)
    2274cld_tau_lsc = 3600.
     75### cld_tau_con  cte de temps utilisee pour eliminer l eau convective (fisrtilp)
    2376cld_tau_con = 3600.
     77### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace large-scale (fisrtilp)
    2478ffallv_lsc = 1.
     79### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace convectifs (fisrtilp)
    2580ffallv_con = 1.
    26 coef_eva = 0.00002
     81### coef_eva  coef evaporation precips eau/glace (fisrtilp/fisrtilp_tr?/conlmd?)
     82coef_eva = 2.e-5
     83### reevap_ice  reevaporation de toute la precip dans la couche du dessous pour la glace (fisrtilp)
    2784reevap_ice = y
     85### iflag_cldcon  flag pour calculer ratqsc=F(ratqsbas,fact_cldcon,q_seri) (physiq)
     86# - iflag_cldcon<=-1 diag. rain_Tiedtke
     87# - iflag_cldcon=1, ratqsc=ratqsbas+fact_cldcon*(q_seri(1)-q_seri(k))/q_seri(k)
     88# - iflag_cldcon=1/2, ratqs=max(ratqs,ratqsc)
     89# - iflag_cldcon=3,   ratqs=ratqss
    2890iflag_cldcon = 3
     91### iflag_pdf :  flag calcul distribution sous-maille de l eau et des nuages
     92# - iflag_pdf=0, version ratqs,
     93# - iflag_pdf=1, calcul eau condensee, fraction nuageuse, eau nuageuse a partir
     94# -              des PDFS de Sandrine Bony
    2995iflag_pdf = 1
     96### fact_cldcon  constante calcul ratqsc (voir iflag_cldcon) et proprietes nuages convectifs, clwcon0 (physiq.F)
    3097fact_cldcon = 1.
    31 #facttemps = 1.e-4
    32 facttemps = 0.0001
     98### facttemps=   facteur de relaxation de ratqs (iflag_cldcon=1/2) et rnebcon (iflag_cldcon=3)
     99facttemps = 1.e-4
     100## ok_newmicro   =y appel newmicro , =n appel nuage (calcul epaisseur optique et emmissivite des nuages)
    33101ok_newmicro = y
     102### ratqsbas     ratqs en bas si iflag_cldcon=1
    34103ratqsbas = 0.005
     104### ratqshaut    ratqs en haut pour ratqss "stables" 
    35105ratqshaut = 0.33
     106### iflag_ratqs=0 correspond a la version IPCC AR4
     107iflag_ratqs = 0
     108### rad_froid    rayon cristaux des nuages de glace "froids"
    36109rad_froid = 35
     110### rad_chau1    rayon goutelettes d eau chauds", en haut: k=4-klev
    37111rad_chau1=12
     112### rad_chau2    rayon goutelettes d eau chauds", en bas: k=1-3 
    38113rad_chau2=11
    39 #ksta_ter=1.e-7
    40 ksta_ter=0.0000001
     114## ksta_ter      coef.diffusion minimale sur terre/sic/lic
     115ksta_ter=1.e-7
     116### cdmmax =     cdrag maximum pour le moment
     117cdmmax = 2.5E-3
     118### cdhmax =     cdrag maximum pour l energie
     119cdhmax = 2.0E-3
    41120#
    42 # parametres outputs
     121# Parametres "orbitaux/ ere geologique"
    43122#
    44 #niveau de sortie "hf" lev_histhf avec
    45 # - lev_histhf=3 => defaut
    46 # - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele
    47 lev_histhf=3
    48 #niveau de sortie "day" lev_histday
    49 # - lev_histday=2 => defaut
    50 # - lev_histday=3 => + champs 3D => F. Lott
    51 # - lev_histday=4 => + champs sous-surfaces
    52 lev_histday=2
    53 #niveau de sortie "mth" lev_histmth avec
    54 # - lev_histmth=2 => defaut
    55 # - lev_histmth=3 => albedo, rugosite sous-surfaces
    56 # - lev_histmth=4 => champs tendances 3d
    57 lev_histmth=2
     123### R_ecc =      Excentricite
     124R_ecc = 0.016715
     125### R_peri =     Equinoxe
     126R_peri = 102.7
     127### R_incl =     Inclinaison
     128R_incl = 23.441
     129### solaire =    Constante solaire
     130solaire = 1365.
    58131#
    59 # parametres climatique
     132# Taux gaz a effet de serre
    60133#
    61 R_ecc = 0.016715
    62 R_peri = 102.7
    63 R_incl = 23.441
    64 solaire = 1365.
     134### co2_ppm =    taux CO2 en ppm
    65135co2_ppm = 348.
     136### CH4_ppb =    taux CH4 en ppb
    66137CH4_ppb = 1650.
     138### N2O_ppb =    taux N2O en ppb
    67139N2O_ppb = 306.
     140### CFC11_ppt =  taux CFC11 en ppt
    68141CFC11_ppt = 280.
     142### CFC12_ppt =  taux CFC12 en ppt
    69143CFC12_ppt = 484.
    70144#
    71 # parametres simulateur ISCCP
     145# Parametres effets directs/indirects des "aerosols"
    72146#
    73 top_height = 3
    74 #overlap = 1, 2 ou 3
     147### ok_ade=y/n   flag Aerosol direct effect
     148ok_ade=n
     149### ok_aie=y/n   flag Aerosol indirect effect
     150ok_aie=n
     151### aer_type =   Aerosol variation type : actuel / preind / scenario
     152aer_type=actuel
     153### bl95_b0 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
     154bl95_b0=1.7
     155### bl95_b1 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
     156bl95_b1=0.2
     157#
     158# Parametres simulateur ISCCP
     159#
     160### ok_isccp=y/n flag simulateur ISCCP
     161ok_isccp=n
     162### top_height = flag choix calcul nuages par le simulateur en utilisant
     163# -              les donnees IR et/ou VIS et l algorithme ISCCP-D1
     164# - top_height = 1 -> algo IR-VIS
     165# - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)"
     166# - top_height = 3 -> algo IR
     167top_height = 1
     168### overlap =    Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP
     169# - overlap=1    Max overlap
     170# - overlap=2    Random overlap
     171# - overlap=3    Max/Random overlap
    75172overlap = 3
    76 #cdmmax
    77 #cdmmax = 2.5E-3
    78 cdmmax = 0.0025
    79 #cdhmax
    80 #cdhmax = 2.0E-3
    81 cdhmax = 0.002
    82 iflag_pbl = 1
    83 iflag_thermals = 0
Note: See TracChangeset for help on using the changeset viewer.