Changeset 1101


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".

Location:
LMDZ4/branches/LMDZ4-dev
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/LMDZ4-dev/gcm.def

    r524 r1101  
    1 #
    2 # $Header$
    3 #
    4 ## nombre de pas par jour (multiple de iperiod) ( ici pour  dt = 1 min )     
     1## nombre de pas par jour (multiple de iperiod)
    52day_step=480
    63## periode pour le pas Matsuno (en pas)
    74iperiod=5
    85## periode de la dissipation (en pas)
    9 idissip=15
     6idissip=5
    107## choix de l'operateur de dissipation (star ou  non star )
    118lstardis=y
    129## nombre d'iterations de l'operateur de dissipation   gradiv
    13 nitergdiv=1
     10nitergdiv=2
    1411## nombre d'iterations de l'operateur de dissipation  nxgradrot
    1512nitergrot=2
    … …  
    1714niterh=2
    1815## temps de dissipation des plus petites long.d ondes pour u,v (gradiv) 
    19 tetagdiv=36000.
     16tetagdiv=10800.
    2017## temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
    2118tetagrot=18000.
    … …  
    2724purmats=n
    2825## avec ou sans physique                                                 
    29 ##physic=n
    30 iflag_phys=1
     26physic=y
    3127## periode de la physique (en pas)                                       
    3228iphysiq=10
    3329## frequence (en  jours ) de l'ecriture du fichier histphy               
    34 ecritphy=1
     30ecritphy=30
    3531##  Cycle diurne  ou non                 
    3632cycle_diurne=y
    … …  
    4743## Nombre  d'appels des routines de rayonnements ( par jour)                 
    4844nbapp_rad=12
    49 ##  Flag  pour la convection (1 pour LMD, 2 pour Tiedtke, 3 KE, 4 KE vect)
    50 iflag_con=3
     45##  Flag  pour la convection : 1 pour LMD, 2 pour Tiedtke, 3 KE(nvlle version JYG), 30 KE(version IPCC AR4), 4 KE vect
     46iflag_con=30
     47## Facteur multiplication des precip convectives dans KE
     48cvl_corr=1.0
    5149## longitude en degres du centre du zoom                                 
    5250clon=0.

  • LMDZ4/branches/LMDZ4-dev/libf/filtrez/filtreg.F

    r1086 r1101  
    184184               
    185185               DO j = jdfil,jffil
     186#ifdef BLAS
    186187                  CALL DGEMM("N", "N", iim, nbniv, iim, 1.0,
    187188     &                 matrinvn(1,1,j),
    188189     &                 iim, champ(1,j,1), iip1*nlat, 0.0,
    189190     &                 eignq(1,j-jdfil+1,1), iim*nlat)
     191#else
     192                  eignq(:,j-jdfil+1,:)
     193     $                 = matmul(matrinvn(:,:,j), champ(:iim,j,:))
     194#endif
    190195               END DO
    191196               
    … …  
    193198               
    194199               DO j = jdfil,jffil
     200#ifdef BLAS
    195201                  CALL DGEMM("N", "N", iim, nbniv, iim, 1.0,
    196202     &                 matriceun(1,1,j),
    197203     &                 iim, champ(1,j,1), iip1*nlat, 0.0,
    198204     &                 eignq(1,j-jdfil+1,1), iim*nlat)
     205#else
     206                  eignq(:,j-jdfil+1,:)
     207     $                 = matmul(matriceun(:,:,j), champ(:iim,j,:))
     208#endif
    199209               END DO
    200210               
    … …  
    202212               
    203213               DO j = jdfil,jffil
     214#ifdef BLAS
    204215                  CALL DGEMM("N", "N", iim, nbniv, iim, 1.0,
    205216     &                 matricevn(1,1,j),
    206217     &                 iim, champ(1,j,1), iip1*nlat, 0.0,
    207218     &                 eignq(1,j-jdfil+1,1), iim*nlat)
     219#else
     220                  eignq(:,j-jdfil+1,:)
     221     $                 = matmul(matricevn(:,:,j), champ(:iim,j,:))
     222#endif
    208223               END DO
    209224               
    … …  
    215230               
    216231               DO j = jdfil,jffil
     232#ifdef BLAS
    217233                  CALL DGEMM("N", "N", iim, nbniv, iim, 1.0,
    218234     &                 matrinvs(1,1,j-jfiltsu+1),
    219235     &                 iim, champ(1,j,1), iip1*nlat, 0.0,
    220236     &                 eignq(1,j-jdfil+1,1), iim*nlat)
     237#else
     238                  eignq(:,j-jdfil+1,:)
     239     $                 = matmul(matrinvs(:,:,j-jfiltsu+1),
     240     $                 champ(:iim,j,:))
     241#endif
    221242               END DO
    222243               
    … …  
    225246               
    226247               DO j = jdfil,jffil
     248#ifdef BLAS
    227249                  CALL DGEMM("N", "N", iim, nbniv, iim, 1.0,
    228250     &                 matriceus(1,1,j-jfiltsu+1),
    229251     &                 iim, champ(1,j,1), iip1*nlat, 0.0,
    230252     &                 eignq(1,j-jdfil+1,1), iim*nlat)
     253#else
     254                  eignq(:,j-jdfil+1,:)
     255     $                 = matmul(matriceus(:,:,j-jfiltsu+1),
     256     $                 champ(:iim,j,:))
     257#endif
    231258               END DO
    232259                             
    … …  
    234261               
    235262               DO j = jdfil,jffil
     263#ifdef BLAS
    236264                  CALL DGEMM("N", "N", iim, nbniv, iim, 1.0,
    237265     &                 matricevs(1,1,j-jfiltsv+1),
    238266     &                 iim, champ(1,j,1), iip1*nlat, 0.0,
    239267     &                 eignq(1,j-jdfil+1,1), iim*nlat)
     268#else
     269                  eignq(:,j-jdfil+1,:)
     270     $                 = matmul(matricevs(:,:,j-jfiltsv+1),
     271     $                 champ(:iim,j,:))
     272#endif
    240273               END DO
    241274                             

  • 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

  • LMDZ4/branches/LMDZ4-dev/run.def

    r1083 r1101  
    1 #
    2 # $Header: /home/cvsroot/LMDZ4/run.def,v 1.1.1.1 2004-05-19 12:53:05 lmdzadmin Exp $
    3 #
    41INCLUDEDEF=physiq.def
    52INCLUDEDEF=gcm.def
    … …  
    96dayref=1
    107##  Annee de l'etat  initial (   avec  4  chiffres   )
    11 anneeref=1979
     8anneeref=1980
    129## Nombre de jours d'integration
    1310nday=1
    1411## periode de sortie des variables de controle (en pas)
    15 iconser=5
     12iconser=240
    1613## periode d'ecriture du fichier histoire (en jour)
    1714iecri=1
    18 ## periode de stockage fichier histmoy (en jour)
    19 periodav=1.
    20 ## unite de sortie des impressions
    21 lunout=7
     15## flag de sortie dynzon
     16ok_dynzon=n
     17## periode de stockage fichier dynzon (en jour)
     18periodav=30.
Note: See TracChangeset for help on using the changeset viewer.