Ignore:
Timestamp:
Jul 1, 2010, 11:02:53 AM (14 years ago)
Author:
Laurent Fairhead
Message:

Merged LMDZ4V5.0-dev branch changes r1292:r1399 to trunk.

Validation:
Validation consisted in compiling the HEAD revision of the trunk,
LMDZ4V5.0-dev branch and the merged sources and running different
configurations on local and SX8 machines comparing results.

Local machine: bench configuration, 32x24x11, gfortran

  • IPSLCM5A configuration (comparison between trunk and merged sources):
    • numerical convergence on dynamical fields over 3 days
    • start files are equivalent (except for RN and PB fields)
    • daily history files equivalent
  • MH07 configuration, new physics package (comparison between LMDZ4V5.0-dev branch and merged sources):
    • numerical convergence on dynamical fields over 3 days
    • start files are equivalent (except for RN and PB fields)
    • daily history files equivalent

SX8 machine (brodie), 96x95x39 on 4 processors:

  • IPSLCM5A configuration:
    • start files are equivalent (except for RN and PB fields)
    • monthly history files equivalent
  • MH07 configuration:
    • start files are equivalent (except for RN and PB fields)
    • monthly history files equivalent

Changes to the makegcm and create_make_gcm scripts to take into account
main programs in F90 files

Fusion de la branche LMDZ4V5.0-dev (r1292:r1399) au tronc principal

Validation:
La validation a consisté à compiler la HEAD de le trunk et de la banche
LMDZ4V5.0-dev et les sources fusionnées et de faire tourner le modéle selon
différentes configurations en local et sur SX8 et de comparer les résultats

En local: 32x24x11, config bench/gfortran

  • pour une config IPSLCM5A (comparaison tronc/fusion):
    • convergence numérique sur les champs dynamiques après 3 jours
    • restart et restartphy égaux (à part sur RN et Pb)
    • fichiers histoire égaux
  • pour une config nlle physique (MH07) (comparaison LMDZ4v5.0-dev/fusion):
    • convergence numérique sur les champs dynamiques après 3 jours
    • restart et restartphy égaux
    • fichiers histoire équivalents

Sur brodie, 96x95x39 sur 4 proc:

  • pour une config IPSLCM5A:
    • restart et restartphy égaux (à part sur RN et PB)
    • pas de différence dans les fichiers histmth.nc
  • pour une config MH07
    • restart et restartphy égaux (à part sur RN et PB)
    • pas de différence dans les fichiers histmth.nc

Changement sur makegcm et create_make-gcm pour pouvoir prendre en compte des
programmes principaux en *F90

Location:
LMDZ4/trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/trunk

  • LMDZ4/trunk/physiq.def

    r1279 r1403  
    1 #
    21## $Id$
    32#
    4 #
    5 # Automatically generated make config: don t edit
    6 #
    7 type_ocean=force
    8 # avec ou sans orchidee
    9 VEGET=n
    10 #type_run = AMIP, ENSP, clim
    11 type_run=AMIP
    12 #
    13 # Controle des sorties
    14 # sorties moyennees tous les jours  dans histday.nc
    15 OK_journe=y
    16 # sorties moyennees tous les mois  dans histmth.nc
    17 OK_mensuel=y
    18 # sorties moyennees toutes les 6 ou bien 3h dans histhf.nc
    19 ok_hf=n
    20 # sorties moyennees tous les pas de temps de la physique dans histins.nc
    21 OK_instan=n
    22 #
    23 ecrit_mth=30.
    24 ecrit_day=1.
    25 ecrit_hf=0.25
    26 #
    27 #niveau de sortie "hf" lev_histhf
    28 lev_histhf=4
    29 #niveau de sortie "day" lev_histday
    30 lev_histday=5
    31 #niveau de sortie "mth" lev_histmth
    32 lev_histmth=4
    33 
    34 # parametres KE
    35 if_ebil=0
    36 epmax = .99
    37 ok_adj_ema = n
    38 iflag_clw = 1
    39 #
    40 # parametres nuages
    41 cld_lc_lsc = 2.6e-4
    42 cld_lc_con = 2.6e-4
    43 cld_tau_lsc = 3600.
    44 cld_tau_con = 3600.
    45 ffallv_lsc = 1.
    46 ffallv_con = 1.
    47 coef_eva = 2.e-5
    48 reevap_ice = y
    49 iflag_cldcon = 3
    50 iflag_pdf = 1
    51 fact_cldcon = 1.
    52 facttemps = 1.e-4
    53 ok_newmicro = y
    54 iflag_ratqs=0
    55 ratqsbas = 0.005
    56 ratqshaut = 0.33
    57 rad_froid = 35
    58 rad_chau1=12
    59 rad_chau2=11
    60 ksta_ter=1.e-7
    61 ksta=1.e-10
    62 #ok_kzmin : calcul Kzmin dans la CL de surface
    63 ok_kzmin=y
    64 #
    65 # parametres climatique
    66 R_ecc = 0.016715
    67 R_peri = 102.7
    68 R_incl = 23.441
    69 solaire = 1365.
    70 co2_ppm = 348.
    71 #RCO2 = co2_ppm * 1.0e-06  * 44.011/28.97
    72 #RCO2 = 348. * 1.0e-06  * 44.011/28.97
    73 #RCO2 =   5.286789092164308E-04
    74 #RCO2 = 425.43e-06
    75 CH4_ppb = 1650.
    76 #RCH4 = 1.65E-06* 16.043/28.97
    77 #RCH4 =   9.137366240938903E-07
    78 N2O_ppb = 306.
    79 #RN2O = 306.E-09* 44.013/28.97
    80 #RN2O =    4.648939592682085E-07
    81 CFC11_ppt = 280.
    82 #RCFC11 = 280.E-12* 137.3686/28.97
    83 #RCFC11 =    1.327690990680013E-09
    84 CFC12_ppt = 484.
    85 #RCFC12 = 484.E-12* 120.9140/28.97
    86 #RCFC12 =    2.020102726958923E-09
    87 #
    88 # effets direct et indirect des aerosols
    89 ok_ade=n
    90 ok_aie=n
    91 bl95_b0=1.7
    92 bl95_b1=0.2
    93 #
    94 # parametres simulateur ISCCP
    95 #ok_isccp : y/n avec/sans simulateur ISCCP
    96 ok_isccp=n
    97 #top_height = 1 ou 3
    98 top_height = 1
    99 #overlap = 1, 2 ou 3
    100 overlap = 3
    101 #cdmmax
    102 cdmmax = 2.5E-3
    103 #cdhmax
    104 cdhmax = 2.0E-3
    105 #
    106 #ok_regdyn : y/n calcul/non des regymes dynamiques sur regions pre-definies
    107 ok_regdyn=y
    108 #
    109 # Flag  pour la convection (1 pour LMD, 2 pour Tiedtke, 3 KE nouvelle physique, 30 KE IPCC)
    110 iflag_con=30
    111 #
    112 # activation thermiques wake, ...
    113 iflag_thermals = 0
    114 nsplit_thermals =1
    115 tau_thermals=1800.
    116 iflag_pbl = 1
    117 iflag_coupl=0
    118 iflag_wake=0
    119 iflag_clos=0
    120 iflag_mix=1
    121 qqa1=0.
    122 qqa2=1.
    123 ## frequence (en  jours ) de l'ecriture du fichier histphy               
    124 ecritphy=30
     3# PARAMETRES ANCIENNEMENT DANS gcm.def
    1254##  Cycle diurne  ou non                 
    1265cycle_diurne=y
    … …  
    1298##  Choix ou non  de  New oliq               
    1309new_oliq=y
     10##  Activation ou non de la parametrisation de Hines pour la strato
     11ok_hines=n
    13112##  Orodr  ou  non   pour l orographie             
    13213ok_orodr=y
    … …  
    13617ok_limitvrai=n
    13718## Nombre  d'appels des routines de rayonnements ( par jour)                 
    138 nbapp_rad=12
     19nbapp_rad=24
     20##  Flag  pour la convection : 1 pour LMD, 2 pour Tiedtke, 3 KE(nvlle version JYG), 30 KE(version IPCC AR4), 4 KE vect
     21iflag_con=30
    13922## Facteur multiplication des precip convectives dans KE
    14023cvl_corr=1.0
     24##  Facteur additif pour l'albedo
     25pmagic=0.008
     26#
     27#
     28#
     29# Parametres fichiers de sortie
     30#
     31### type_run = type run par rapport aux fichiers et variables de sortie
     32# - type_run = CLIM/ENSP (=1)
     33# - type_run = AMIP/CFMI (=2)
     34type_run=AMIP
     35### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc
     36OK_journe=n
     37### OK_mensuel= y sortir fichier mensuel histmth.nc, =n pas de fichier histmth.nc
     38OK_mensuel=y
     39### OK_instan=y, ecrire sorties "instantannees" (chaque pas de temps de la  physique)
     40OK_instan=n
     41### OK_hf=y, ecrire sorties hautes frequence histhf.nc, =n pas de fichier histhf.nc
     42ok_hf=n
     43#
     44# Parametres niveau de sorties differents fichiers
     45#
     46### lev_histhf=0-4, niveau de sortie fichier "histhf.nc"
     47# - lev_histhf=0 => pas de sorties histhf.nc
     48# - lev_histhf=2 => defaut
     49# - lev_histhf=3 => variables sur niveaux standards
     50# - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele => fichier. histhf3d.nc
     51lev_histhf=2
     52### lev_histday=0-5, niveau de sortie fichier "histday.nc"
     53# - lev_histday=0 => pas de sorties lev_histday.nc
     54# - lev_histday=2 => defaut
     55# - lev_histday=3 => + champs 3D => F. Lott
     56# - lev_histday=4 => + champs sous-surfaces
     57# - lev_histday=5 => + champs F. Aires
     58lev_histday=2
     59### lev_histmth=0-4, niveau de sortie fichier "histmth.nc"
     60# - lev_histmth=0 => pas de sorties lev_histmth.nc
     61# - lev_histmth=2 => defaut
     62# - lev_histmth=3 => albedo, rugosite sous-surfaces
     63# - lev_histmth=4 => champs tendances 3d
     64lev_histmth=2
     65### ecrit_hf = frequence ecriture fichier histhf.nc en jours
     66ecrit_hf=0.250
     67### ecrit_day = frequence ecriture fichier histday.nc en jours
     68ecrit_day=1.
     69### ecrit_mth = frequence ecriture fichier histmth.nc en jours
     70ecrit_mth=30
     71### freqin_isccp = frequence input en secondes du simulateur ISCCP
     72freq_ISCCP=10800.
     73### freqout_isccp = frequence output en jours du simulateur ISCCP
     74ecrit_ISCCP=30
     75### niveau du diagnostique de conservation d energie
     76if_ebil=0
     77#
     78# parametres KE
     79#
     80### epmax = Efficacite precipitation maximale
     81epmax = .999
     82### ok_adj_ema = ?? pas utilise
     83ok_adj_ema = n
     84### iflag_clw Flag calcul eau liquide
     85# - iflag_clw=0 : qcond_incld(i,l) = em_qcondc(l)
     86# - iflag_clw=1 : qcond_incld(i,l) = em_qcond(l)
     87# - iflag_clw=2 : eau liquide diagnostique en fonction de la Precip
     88iflag_clw = 1
     89#
     90# parametres nuages
     91#
     92### cld_lc_lsc  contenu en eau liquide des nuages large-scale (fisrtilp)
     93cld_lc_lsc = 4.16e-4
     94### cld_lc_con  contenu en eau liquide des nuages convectifs (fisrtilp)
     95cld_lc_con = 4.16e-4
     96### cld_tau_lsc cte de temps utilisee pour eliminer l eau large-scale (fisrtilp)
     97cld_tau_lsc = 1800.
     98### cld_tau_con cte de temps utilisee pour eliminer l eau convective (fisrtilp)
     99cld_tau_con = 1800.
     100### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace large-scale (fisrtilp)
     101ffallv_lsc = 0.5
     102### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace convectifs (fisrtilp)
     103ffallv_con = 0.5
     104### coef_eva coef evaporation precips eau/glace (fisrtilp/fisrtilp_tr?/conlmd?)
     105coef_eva = 2.e-5
     106### reevap_ice  reevaporation de toute la precip dans la couche du dessous pour la glace (fisrtilp)
     107reevap_ice = y
     108### iflag_cldcon  flag pour calculer ratqsc=F(ratqsbas,fact_cldcon,q_seri) (physiq)
     109# - iflag_cldcon<=-1 diag. rain_Tiedtke
     110# - iflag_cldcon=1, ratqsc=ratqsbas+fact_cldcon*(q_seri(1)-q_seri(k))/q_seri(k)
     111# - iflag_cldcon=1/2, ratqs=max(ratqs,ratqsc)
     112# - iflag_cldcon=3,   ratqs=ratqss
     113iflag_cldcon = 3
     114### iflag_pdf :  flag calcul distribution sous-maille de l eau et des nuages
     115# - iflag_pdf=0, version ratqs,
     116# - iflag_pdf=1, calcul eau condensee, fraction nuageuse, eau nuageuse a partir
     117# -              des PDFS de Sandrine Bony
     118iflag_pdf = 1
     119### fact_cldcon  constante calcul ratqsc (voir iflag_cldcon) et proprietes nuages convectifs, clwcon0 (physiq.F)
     120fact_cldcon = 1.
     121### facttemps=   facteur de relaxation de ratqs (iflag_cldcon=1/2) et rnebcon (iflag_cldcon=3)
     122facttemps = 0.
     123## ok_newmicro   =y appel newmicro , =n appel nuage (calcul epaisseur optique et emmissivite des nuages)
     124ok_newmicro = y
     125### iflag_ratqs=0 correspond a la version IPCC AR4
     126iflag_ratqs=0
     127### ratqsbas     ratqs en bas si iflag_cldcon=1
     128ratqsbas = 0.005
     129### ratqshaut    ratqs en haut pour ratqss "stables" 
     130ratqshaut = 0.33
     131### rad_froid    rayon cristaux des nuages de glace "froids"
     132rad_froid = 35
     133### rad_chau1    rayon goutelettes d eau chauds", en haut: k=4-klev
     134rad_chau1=12
     135### rad_chau2    rayon goutelettes d eau chauds", en bas: k=1-3 
     136rad_chau2=11
     137#
     138# Coefficient et parametres sur les drags
     139#
     140f_cdrag_ter=1.
     141f_cdrag_stable=1.
     142f_cdrag_oce=0.8
     143f_rugoro=0.
     144### ksta_ter      coef.diffusion minimale sur terre/sic/lic
     145ksta_ter=1.e-7
     146### cdmmax =     cdrag maximum pour le moment
     147cdmmax = 2.5E-3
     148### cdhmax =     cdrag maximum pour l energie
     149cdhmax = 2.0E-3
     150#
     151# Parametres "orbitaux/ ere geologique"
     152#
     153### R_ecc =      Excentricite
     154R_ecc = 0.016715
     155### R_peri =     Equinoxe
     156R_peri = 102.7
     157### R_incl =     Inclinaison
     158R_incl = 23.441
     159### solaire =    Constante solaire
     160solaire = 1361.
     161#
     162# Taux gaz a effet de serre
     163#
     164### co2_ppm =    taux CO2 en ppm
     165co2_ppm = 367.
     166### CH4_ppb =    taux CH4 en ppb
     167CH4_ppb = 1760.
     168### N2O_ppb =    taux N2O en ppb
     169N2O_ppb = 316.
     170### CFC11_ppt =  taux CFC11 en ppt
     171CFC11_ppt = 741.2
     172### CFC12_ppt =  taux CFC12 en ppt
     173CFC12_ppt = 535.
     174#
     175# Parametres effets directs/indirects des "aerosols"
     176#
     177### ok_ade=y/n   flag Aerosol direct effect
     178ok_ade=n
     179### ok_aie=y/n   flag Aerosol indirect effect
     180ok_aie=n
     181### aer_type =   Aerosol variation type : actuel / preind / scenario / annuel
     182aer_type=actuel
     183###  type of coupled aerosol =1 (default) =2 => bc  only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol
     184flag_aerosol=6
     185### bl95_b0 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
     186bl95_b0=1.7
     187### bl95_b1 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
     188bl95_b1=0.2
     189#
     190# Parametre de lecture de l'ozone
     191#
     192# Allowed values are 0, 1 and 2
     193# 0: do not read an ozone climatology
     194# 1: read a single ozone climatology that will be used day and night
     195# 2: read two ozone climatologies, the average day and night climatology and the daylight climatology
     196read_climoz=0
     197#
     198# Parametres simulateur COSP (CFMIP Observational Simulator Package)
     199#
     200### ok_cosp=y/n flag simulateur COSP
     201ok_cosp=n
     202## freq_COSP = frequence d'appel de COSP en secondes
     203freq_COSP=10800.
     204#
     205# Parametres simulateur ISCCP
     206#
     207### ok_isccp=y/n flag simulateur ISCCP
     208ok_isccp=n
     209### top_height = flag choix calcul nuages par le simulateur en utilisant
     210# -              les donnees IR et/ou VIS et l algorithme ISCCP-D1
     211# - top_height = 1 -> algo IR-VIS
     212# - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)"
     213# - top_height = 3 -> algo IR
     214top_height = 1
     215### overlap =    Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP
     216# - overlap=1    Max overlap
     217# - overlap=2    Random overlap
     218# - overlap=3    Max/Random overlap
     219overlap = 3
Note: See TracChangeset for help on using the changeset viewer.