Ignore:
Timestamp:
Aug 17, 2006, 5:41:51 PM (18 years ago)
Author:
Laurent Fairhead
Message:

Inclusion des modifs de Y. Meurdesoif pour la version V3
LF

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/V3_test/libf/phylmd/conf_phys.F90

    r702 r704  
    2323include "clesphys.inc"
    2424include "compbl.h"
    25 #include "control.inc"
     25include "control.inc"
    2626!
    2727! Configuration de la "physique" de LMDZ a l'aide de la fonction
    … …  
    4949  REAL                 :: bl95_b0, bl95_b1
    5050  real                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut
    51 !IM integer            :: iflag_cldcon, if_ebil
    5251  integer              :: iflag_cldcon
     52
     53  character (len = 6),SAVE  :: ocean_omp
     54  logical,SAVE              :: ok_veget_omp, ok_newmicro_omp
     55  logical,SAVE        :: ok_journe_omp, ok_mensuel_omp, ok_instan_omp       
     56  LOGICAL,SAVE        :: ok_ade_omp, ok_aie_omp
     57  REAL,SAVE           :: bl95_b0_omp, bl95_b1_omp
     58  real,SAVE           :: fact_cldcon_omp, facttemps_omp,ratqsbas_omp
     59  real,SAVE           :: ratqshaut_omp
     60  integer,SAVE        :: iflag_cldcon_omp, ip_ebil_phy_omp
    5361
    5462! Local
    … …  
    5765
    5866  integer :: iflag_thermals,nsplit_thermals
    59 !
    60 !
    61 !
    62 
    63 
     67  integer,SAVE :: iflag_thermals_omp,nsplit_thermals_omp
     68
     69  REAL,SAVE :: R_ecc_omp,R_peri_omp,R_incl_omp,solaire_omp,co2_ppm_omp
     70  REAL,SAVE :: RCO2_omp,CH4_ppb_omp,RCH4_omp,N2O_ppb_omp,RN2O_omp,CFC11_ppt_omp
     71  REAL,SAVE :: RCFC11_omp,CFC12_ppt_omp,RCFC12_omp,epmax_omp
     72  LOGICAL,SAVE :: ok_adj_ema_omp
     73  INTEGER,SAVE :: iflag_clw_omp
     74  REAL,SAVE :: cld_lc_lsc_omp,cld_lc_con_omp,cld_tau_lsc_omp,cld_tau_con_omp
     75  REAL,SAVE :: ffallv_lsc_omp, ffallv_con_omp,coef_eva_omp
     76  LOGICAL,SAVE :: reevap_ice_omp
     77  INTEGER,SAVE :: iflag_pdf_omp
     78  REAL,SAVE :: rad_froid_omp, rad_chau1_omp, rad_chau2_omp
     79  INTEGER,SAVE :: top_height_omp,overlap_omp
     80  REAL,SAVE :: cdmmax_omp,cdhmax_omp,ksta_omp,ksta_ter_omp
     81  LOGICAL,SAVE :: ok_kzmin_omp
     82  INTEGER,SAVE :: iflag_pbl_omp,lev_histhf_omp,lev_histday_omp,lev_histmth_omp
     83  CHARACTER*4, SAVE :: type_run_omp
     84  LOGICAL,SAVE :: ok_isccp_omp
     85  REAL,SAVE :: lonmin_ins_omp, lonmax_ins_omp, latmin_ins_omp, latmax_ins_omp
     86  REAL,SAVE :: ecrit_hf_omp, ecrit_day_omp, ecrit_mth_omp, ecrit_reg_omp
     87  REAL,SAVE :: ecrit_tra_omp
     88  LOGICAL, SAVE :: ok_slab_sicOBS_omp
     89  REAL,SAVE :: cvl_corr_omp
     90  LOGICAL,SAVE :: ok_lic_melt_omp
     91!
     92!
     93!
     94
     95!$OMP MASTER
    6496!Config Key  = OCEAN
    6597!Config Desc = Type d'ocean
    … …  
    6799!Config Help = Type d'ocean utilise: force, slab,couple
    68100!
    69   ocean = 'force '
    70   call getin('OCEAN', ocean)
     101  ocean_omp = 'force '
     102  call getin('OCEAN', ocean_omp)
    71103!
    72104!Config Key  = VEGET
    … …  
    75107!Config Help = Type de modele de vegetation utilise
    76108!
    77   ok_veget = .false.
    78   call getin('VEGET', ok_veget)
     109  ok_veget_omp = .false.
     110  call getin('VEGET', ok_veget_omp)
    79111!
    80112!Config Key  = OK_journe
    … …  
    84116!              journalieres
    85117!
    86   ok_journe = .false.
    87   call getin('OK_journe', ok_journe)
     118  ok_journe_omp = .false.
     119  call getin('OK_journe', ok_journe_omp)
    88120!
    89121!Config Key  = OK_mensuel
    … …  
    93125!              mensuelles
    94126!
    95   ok_mensuel = .true.
    96   call getin('OK_mensuel', ok_mensuel)
     127  ok_mensuel_omp = .true.
     128  call getin('OK_mensuel', ok_mensuel_omp)
    97129!
    98130!Config Key  = OK_instan
    … …  
    102134!              instantanees
    103135!
    104   ok_instan = .false.
    105   call getin('OK_instan', ok_instan)
     136  ok_instan_omp = .false.
     137  call getin('OK_instan', ok_instan_omp)
    106138!
    107139!Config Key  = ok_ade
    … …  
    110142!Config Help = Used in radlwsw.F
    111143!
    112   ok_ade = .false.
    113   call getin('ok_ade', ok_ade)
     144  ok_ade_omp = .false.
     145  call getin('ok_ade', ok_ade_omp)
    114146
    115147!
    … …  
    119151!Config Help = Used in nuage.F and radlwsw.F
    120152!
    121   ok_aie = .false.
    122   call getin('ok_aie', ok_aie)
     153  ok_aie_omp = .false.
     154  call getin('ok_aie', ok_aie_omp)
    123155
    124156!
    … …  
    128160!Config Help = Used in nuage.F
    129161!
    130   bl95_b0 = 2.
    131   call getin('bl95_b0', bl95_b0)
     162  bl95_b0_omp = 2.
     163  call getin('bl95_b0', bl95_b0_omp)
    132164
    133165!Config Key  = bl95_b1
    … …  
    136168!Config Help = Used in nuage.F
    137169!
    138   bl95_b1 = 0.2
    139   call getin('bl95_b1', bl95_b1)
     170  bl95_b1_omp = 0.2
     171  call getin('bl95_b1', bl95_b1_omp)
    140172
    141173!
    … …  
    147179!               
    148180!
    149   ip_ebil_phy = 0
    150   call getin('ip_ebil_phy', ip_ebil_phy)
     181  ip_ebil_phy_omp = 0
     182  call getin('ip_ebil_phy', ip_ebil_phy_omp)
    151183!!
    152184!! Constante solaire & Parametres orbitaux & taux gaz effet de serre BEG
    … …  
    158190!               
    159191!valeur AMIP II
    160   R_ecc = 0.016715
    161   call getin('R_ecc', R_ecc)
     192  R_ecc_omp = 0.016715
     193  call getin('R_ecc', R_ecc_omp)
    162194!!
    163195!Config Key  = R_peri
    … …  
    168200!
    169201!valeur AMIP II
    170   R_peri = 102.7
    171   call getin('R_peri', R_peri)
     202  R_peri_omp = 102.7
     203  call getin('R_peri', R_peri_omp)
    172204!!
    173205!Config Key  = R_incl
    … …  
    178210!
    179211!valeur AMIP II
    180   R_incl = 23.441
    181   call getin('R_incl', R_incl)
     212  R_incl_omp = 23.441
     213  call getin('R_incl', R_incl_omp)
    182214!!
    183215!Config Key  = solaire
    … …  
    188220!
    189221!valeur AMIP II
    190   solaire = 1365.
    191   call getin('solaire', solaire)
     222  solaire_omp = 1365.
     223  call getin('solaire', solaire_omp)
    192224!!
    193225!Config Key  = co2_ppm
    … …  
    198230!
    199231!valeur AMIP II
    200   co2_ppm = 348.
    201   call getin('co2_ppm', co2_ppm)
     232  co2_ppm_omp = 348.
     233  call getin('co2_ppm', co2_ppm_omp)
    202234!!
    203235!Config Key  = RCO2
    … …  
    209241! RCO2 = 5.286789092164308E-04
    210242!ancienne valeur
    211   RCO2 = co2_ppm * 1.0e-06  * 44.011/28.97 ! pour co2_ppm=348.
     243  RCO2_omp = co2_ppm_omp * 1.0e-06  * 44.011/28.97 ! pour co2_ppm=348.
    212244
    213245!!  call getin('RCO2', RCO2)
    … …  
    228260  zzz = 1650.
    229261  call getin('CH4_ppb', zzz)
    230   CH4_ppb = zzz
    231   RCH4 = CH4_ppb * 1.0E-09 * 16.043/28.97
     262  CH4_ppb_omp = zzz
     263  RCH4_omp = CH4_ppb_omp * 1.0E-09 * 16.043/28.97
    232264!!
    233265!Config Key  = RN2O
    … …  
    246278  zzz=306.
    247279  call getin('N2O_ppb', zzz)
    248   N2O_ppb = zzz
    249   RN2O = N2O_ppb * 1.0E-09 * 44.013/28.97
     280  N2O_ppb_omp = zzz
     281  RN2O_omp = N2O_ppb_omp * 1.0E-09 * 44.013/28.97
    250282!!
    251283!Config Key  = RCFC11
    … …  
    258290  zzz = 280.
    259291  call getin('CFC11_ppt',zzz)
    260   CFC11_ppt = zzz
    261   RCFC11=CFC11_ppt* 1.0E-12 * 137.3686/28.97
     292  CFC11_ppt_omp = zzz
     293  RCFC11_omp=CFC11_ppt_omp* 1.0E-12 * 137.3686/28.97
    262294! RCFC11 = 1.327690990680013E-09
    263295!OK call getin('RCFC11', RCFC11)
    … …  
    272304  zzz = 484.
    273305  call getin('CFC12_ppt',zzz)
    274   CFC12_ppt = zzz
    275   RCFC12 = CFC12_ppt * 1.0E-12 * 120.9140/28.97
     306  CFC12_ppt_omp = zzz
     307  RCFC12_omp = CFC12_ppt_omp * 1.0E-12 * 120.9140/28.97
    276308! RCFC12 = 2.020102726958923E-09
    277309!OK call getin('RCFC12', RCFC12)
    … …  
    281313!! KE
    282314!
     315
     316!Config key  = cvl_corr
     317!Config Desc = Facteur multiplication des precip convectives dans KE
     318!Config Def  = 1.00
     319!Config Help = 1.02 pour un moderne ou un pre-ind. A ajuster pour un glaciaire
     320  cvl_corr_omp = 1.00
     321  CALL getin('cvl_corr', cvl_corr_omp)
     322
     323
    283324!Config Key  = epmax
    284325!Config Desc = Efficacite precip
    … …  
    286327!Config Help =
    287328!
    288   epmax = .993
    289   call getin('epmax', epmax)
     329  epmax_omp = .993
     330  call getin('epmax', epmax_omp)
    290331!
    291332!Config Key  = ok_adj_ema
    … …  
    294335!Config Help =
    295336!
    296   ok_adj_ema = .false.
    297   call getin('ok_adj_ema',ok_adj_ema)
     337  ok_adj_ema_omp = .false.
     338  call getin('ok_adj_ema',ok_adj_ema_omp)
    298339!
    299340!Config Key  = iflag_clw
    … …  
    302343!Config Help =
    303344!
    304   iflag_clw = 0
    305   call getin('iflag_clw',iflag_clw)
     345  iflag_clw_omp = 0
     346  call getin('iflag_clw',iflag_clw_omp)
    306347!
    307348!Config Key  = cld_lc_lsc
    … …  
    310351!Config Help =
    311352!
    312   cld_lc_lsc = 2.6e-4
    313   call getin('cld_lc_lsc',cld_lc_lsc)
     353  cld_lc_lsc_omp = 2.6e-4
     354  call getin('cld_lc_lsc',cld_lc_lsc_omp)
    314355!
    315356!Config Key  = cld_lc_con
    … …  
    318359!Config Help =
    319360!
    320   cld_lc_con = 2.6e-4
    321   call getin('cld_lc_con',cld_lc_con)
     361  cld_lc_con_omp = 2.6e-4
     362  call getin('cld_lc_con',cld_lc_con_omp)
    322363!
    323364!Config Key  = cld_tau_lsc
    … …  
    326367!Config Help =
    327368!
    328   cld_tau_lsc = 3600.
    329   call getin('cld_tau_lsc',cld_tau_lsc)
     369  cld_tau_lsc_omp = 3600.
     370  call getin('cld_tau_lsc',cld_tau_lsc_omp)
    330371!
    331372!Config Key  = cld_tau_con
    … …  
    334375!Config Help =
    335376!
    336   cld_tau_con = 3600.
    337   call getin('cld_tau_con',cld_tau_con)
     377  cld_tau_con_omp = 3600.
     378  call getin('cld_tau_con',cld_tau_con_omp)
    338379!
    339380!Config Key  = ffallv_lsc
    … …  
    342383!Config Help =
    343384!
    344   ffallv_lsc = 1.
    345   call getin('ffallv_lsc',ffallv_lsc)
     385  ffallv_lsc_omp = 1.
     386  call getin('ffallv_lsc',ffallv_lsc_omp)
    346387!
    347388!Config Key  = ffallv_con
    … …  
    350391!Config Help =
    351392!
    352   ffallv_con = 1.
    353   call getin('ffallv_con',ffallv_con)
     393  ffallv_con_omp = 1.
     394  call getin('ffallv_con',ffallv_con_omp)
    354395!
    355396!Config Key  = coef_eva
    … …  
    358399!Config Help =
    359400!
    360   coef_eva = 2.e-5
    361   call getin('coef_eva',coef_eva)
     401  coef_eva_omp = 2.e-5
     402  call getin('coef_eva',coef_eva_omp)
    362403!
    363404!Config Key  = reevap_ice
    … …  
    366407!Config Help =
    367408!
    368   reevap_ice = .false.
    369   call getin('reevap_ice',reevap_ice)
     409  reevap_ice_omp = .false.
     410  call getin('reevap_ice',reevap_ice_omp)
    370411!
    371412!Config Key  = iflag_cldcon
    … …  
    374415!Config Help =
    375416!
    376   iflag_cldcon = 1
    377   call getin('iflag_cldcon',iflag_cldcon)
     417  iflag_cldcon_omp = 1
     418  call getin('iflag_cldcon',iflag_cldcon_omp)
    378419
    379420!
    … …  
    383424!Config Help =
    384425!
    385   iflag_pdf = 0
    386   call getin('iflag_pdf',iflag_pdf)
     426  iflag_pdf_omp = 0
     427  call getin('iflag_pdf',iflag_pdf_omp)
    387428!
    388429!Config Key  = fact_cldcon
    … …  
    391432!Config Help =
    392433!
    393   fact_cldcon = 0.375
    394   call getin('fact_cldcon',fact_cldcon)
     434  fact_cldcon_omp = 0.375
     435  call getin('fact_cldcon',fact_cldcon_omp)
    395436
    396437!
    … …  
    400441!Config Help =
    401442!
    402   facttemps = 1.e-4
    403   call getin('facttemps',facttemps)
     443  facttemps_omp = 1.e-4
     444  call getin('facttemps',facttemps_omp)
    404445
    405446!
    … …  
    409450!Config Help =
    410451!
    411   ok_newmicro = .true.
    412   call getin('ok_newmicro',ok_newmicro)
     452  ok_newmicro_omp = .true.
     453  call getin('ok_newmicro',ok_newmicro_omp)
    413454!
    414455!Config Key  = ratqsbas
    … …  
    417458!Config Help =
    418459!
    419   ratqsbas = 0.01
    420   call getin('ratqsbas',ratqsbas)
     460  ratqsbas_omp = 0.01
     461  call getin('ratqsbas',ratqsbas_omp)
    421462!
    422463!Config Key  = ratqshaut
    … …  
    425466!Config Help =
    426467!
    427   ratqshaut = 0.3
    428   call getin('ratqshaut',ratqshaut)
     468  ratqshaut_omp = 0.3
     469  call getin('ratqshaut',ratqshaut_omp)
    429470
    430471!
    … …  
    434475!Config Help =
    435476!
    436   rad_froid = 35.0
    437   call getin('rad_froid',rad_froid)
     477  rad_froid_omp = 35.0
     478  call getin('rad_froid',rad_froid_omp)
    438479
    439480!
    … …  
    443484!Config Help =
    444485!
    445   rad_chau1 = 13.0
    446   call getin('rad_chau1',rad_chau1)
     486  rad_chau1_omp = 13.0
     487  call getin('rad_chau1',rad_chau1_omp)
    447488
    448489!
    … …  
    452493!Config Help =
    453494!
    454   rad_chau2 = 9.0
    455   call getin('rad_chau2',rad_chau2)
     495  rad_chau2_omp = 9.0
     496  call getin('rad_chau2',rad_chau2_omp)
    456497
    457498!
    … …  
    461502!Config Help =
    462503!
    463   top_height = 3
    464   call getin('top_height',top_height)
     504  top_height_omp = 3
     505  call getin('top_height',top_height_omp)
    465506
    466507!
    … …  
    470511!Config Help =
    471512!
    472   overlap = 3
    473   call getin('overlap',overlap)
     513  overlap_omp = 3
     514  call getin('overlap',overlap_omp)
    474515
    475516
    … …  
    481522!Config Help =
    482523!
    483   cdmmax = 1.3E-3
    484   call getin('cdmmax',cdmmax)
     524  cdmmax_omp = 1.3E-3
     525  call getin('cdmmax',cdmmax_omp)
    485526
    486527!
    … …  
    490531!Config Help =
    491532!
    492   cdhmax = 1.1E-3
    493   call getin('cdhmax',cdhmax)
     533  cdhmax_omp = 1.1E-3
     534  call getin('cdhmax',cdhmax_omp)
    494535
    495536!261103
    … …  
    500541!Config Help =
    501542!
    502   ksta = 1.0e-10
    503   call getin('ksta',ksta)
     543  ksta_omp = 1.0e-10
     544  call getin('ksta',ksta_omp)
    504545
    505546!
    … …  
    509550!Config Help =
    510551!
    511   ksta_ter = 1.0e-10
    512   call getin('ksta_ter',ksta_ter)
     552  ksta_ter_omp = 1.0e-10
     553  call getin('ksta_ter',ksta_ter_omp)
    513554
    514555!
    … …  
    518559!Config Help =
    519560!
    520   ok_kzmin = .true.
    521   call getin('ok_kzmin',ok_kzmin)
     561  ok_kzmin_omp = .true.
     562  call getin('ok_kzmin',ok_kzmin_omp)
     563
     564!Config Key = ok_lic_melt
     565!Config Desc = Prise en compte de la fonte de la calotte dans le bilan d'eau
     566!Config Def  = .false.
     567!Config Help = mettre a .false. pour assurer la conservation en eau
     568  ok_lic_melt_omp = .false.
     569  call getin('ok_lic_melt', ok_lic_melt_omp)
    522570
    523571!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    … …  
    530578!Config Help =
    531579!
    532   iflag_pbl = 1
    533   call getin('iflag_pbl',iflag_pbl)
     580  iflag_pbl_omp = 1
     581  call getin('iflag_pbl',iflag_pbl_omp)
    534582!
    535583!Config Key  = iflag_thermals
    … …  
    538586!Config Help =
    539587!
    540   iflag_thermals = 0
    541   call getin('iflag_thermals',iflag_thermals)
     588  iflag_thermals_omp = 0
     589  call getin('iflag_thermals',iflag_thermals_omp)
    542590!
    543591!
    … …  
    547595!Config Help =
    548596!
    549   nsplit_thermals = 1
    550   call getin('nsplit_thermals',nsplit_thermals)
     597  nsplit_thermals_omp = 1
     598  call getin('nsplit_thermals',nsplit_thermals_omp)
    551599
    552600
    … …  
    558606!Config Help =
    559607!
    560   lev_histhf = 1
    561   call getin('lev_histhf',lev_histhf)
     608  lev_histhf_omp = 1
     609  call getin('lev_histhf',lev_histhf_omp)
    562610
    563611!
    … …  
    567615!Config Help =
    568616!
    569   lev_histday = 1
    570   call getin('lev_histday',lev_histday)
     617  lev_histday_omp = 1
     618  call getin('lev_histday',lev_histday_omp)
    571619
    572620!
    … …  
    576624!Config Help =
    577625!
    578   lev_histmth = 2
    579   call getin('lev_histmth',lev_histmth)
    580 
    581 !
     626  lev_histmth_omp = 2
     627  call getin('lev_histmth',lev_histmth_omp)
     628
     629  !
    582630!Config Key  = type_run
    583631!Config Desc =
    … …  
    585633!Config Help =
    586634!
    587   type_run = 'AMIP'
    588   call getin('type_run',type_run)
    589 
    590 !
     635  type_run_omp = 'AMIP'
     636  call getin('type_run',type_run_omp)
     637
     638  !
    591639!Config Key  = ok_isccp
    592640!Config Desc =
    … …  
    595643!
    596644! ok_isccp = .true.
    597   ok_isccp = .false.
    598   call getin('ok_isccp',ok_isccp)
     645  ok_isccp_omp = .false.
     646  call getin('ok_isccp',ok_isccp_omp)
    599647
    600648!
    … …  
    607655!Config Help =
    608656!
    609    lonmin_ins = 100.
    610    call getin('lonmin_ins',lonmin_ins)
     657   lonmin_ins_omp = 100.
     658   call getin('lonmin_ins',lonmin_ins_omp)
    611659!
    612660!Config Key  = lonmax_ins
    … …  
    616664!
    617665   lonmax_ins = 130.
    618    call getin('lonmax_ins',lonmax_ins)
     666   call getin('lonmax_ins',lonmax_ins_omp)
    619667!
    620668!Config Key  = latmin_ins
    … …  
    623671!Config Help =
    624672!
    625    latmin_ins = -20.
    626    call getin('latmin_ins',latmin_ins)
     673   latmin_ins_omp = -20.
     674   call getin('latmin_ins',latmin_ins_omp)
    627675!
    628676!Config Key  = latmax_ins
    … …  
    631679!Config Help =
    632680!
    633    latmax_ins = 20.
    634    call getin('latmax_ins',latmax_ins)
     681   latmax_ins_omp = 20.
     682   call getin('latmax_ins',latmax_ins_omp)
    635683!
    636684!Config Key  = ecrit_hf
    … …  
    639687!Config Help =
    640688!
    641   ecrit_hf = 1./8.
    642   call getin('ecrit_hf',ecrit_hf)
     689  ecrit_hf_omp = 1./8.
     690  call getin('ecrit_hf',ecrit_hf_omp)
    643691!
    644692!Config Key  = ecrit_day
    … …  
    647695!Config Help =
    648696!
    649   ecrit_day = 1.0
    650   call getin('ecrit_day',ecrit_day)
     697  ecrit_day_omp = 1.0
     698  call getin('ecrit_day',ecrit_day_omp)
    651699!
    652700!Config Key  = ecrit_mth
    … …  
    655703!Config Help =
    656704!
    657   ecrit_mth = 30.
    658   call getin('ecrit_mth',ecrit_mth)
     705  ecrit_mth_omp = 30.
     706  call getin('ecrit_mth',ecrit_mth_omp)
    659707!
    660708!Config Key  = ecrit_tra
    … …  
    663711!Config Help =
    664712!
    665   ecrit_tra = 30.
    666   call getin('ecrit_tra',ecrit_tra)
     713  ecrit_tra_omp = 30.
     714  call getin('ecrit_tra',ecrit_tra_omp)
    667715!
    668716!Config Key  = ecrit_reg
    … …  
    671719!Config Help =
    672720!
    673   ecrit_reg = 0.25   !4 fois par jour
    674   call getin('ecrit_reg',ecrit_reg)
     721  ecrit_reg_omp = 0.25   !4 fois par jour
     722  call getin('ecrit_reg',ecrit_reg_omp)
    675723!
    676724!
    … …  
    682730!              de glace de mer Observee
    683731!
    684   ok_slab_sicOBS = .true.
    685   call getin('ok_slab_sicOBS', ok_slab_sicOBS)
     732  ok_slab_sicOBS_omp = .true.
     733  call getin('ok_slab_sicOBS', ok_slab_sicOBS_omp)
     734
     735
     736!$OMP END MASTER
     737!$OMP BARRIER
     738
     739    R_ecc = R_ecc_omp
     740    R_peri = R_peri_omp
     741    R_incl = R_incl_omp
     742    solaire = solaire_omp
     743    co2_ppm = co2_ppm_omp
     744    RCO2 = RCO2_omp
     745    CH4_ppb = CH4_ppb_omp
     746    RCH4 = RCH4_omp
     747    N2O_ppb = N2O_ppb_omp
     748    RN2O = RN2O_omp
     749    CFC11_ppt = CFC11_ppt_omp
     750    RCFC11 = RCFC11_omp
     751    CFC12_ppt = CFC12_ppt_omp
     752    RCFC12 = RCFC12_omp
     753    epmax = epmax_omp
     754    ok_adj_ema = ok_adj_ema_omp
     755    iflag_clw = iflag_clw_omp
     756    cld_lc_lsc = cld_lc_lsc_omp
     757    cld_lc_con = cld_lc_con_omp
     758    cld_tau_lsc = cld_tau_lsc_omp
     759    cld_tau_con = cld_tau_con_omp
     760    ffallv_lsc = ffallv_lsc_omp
     761    ffallv_con = ffallv_con_omp
     762    coef_eva = coef_eva_omp
     763    reevap_ice = reevap_ice_omp
     764    iflag_pdf = iflag_pdf_omp
     765    rad_froid = rad_froid_omp
     766    rad_chau1 = rad_chau1_omp
     767    rad_chau2 = rad_chau2_omp
     768    top_height = top_height_omp
     769    overlap = overlap_omp
     770    cdmmax = cdmmax_omp
     771    cdhmax = cdhmax_omp
     772    ksta = ksta_omp
     773    ksta_ter = ksta_ter_omp
     774    ok_kzmin = ok_kzmin_omp
     775    iflag_pbl = iflag_pbl_omp
     776    lev_histhf = lev_histhf_omp
     777    lev_histday = lev_histday_omp
     778    lev_histmth = lev_histmth_omp
     779
     780    ocean = ocean_omp
     781    ok_veget = ok_veget_omp
     782    ok_newmicro = ok_newmicro_omp
     783    ok_journe = ok_journe_omp
     784    ok_mensuel = ok_mensuel_omp
     785    ok_instan = ok_instan_omp
     786    ok_ade = ok_ade_omp
     787    ok_aie = ok_aie_omp
     788    bl95_b0 = bl95_b0_omp
     789    bl95_b1 = bl95_b1_omp
     790    fact_cldcon = fact_cldcon_omp
     791    facttemps = facttemps_omp
     792    ratqsbas = ratqsbas_omp
     793    ratqshaut = ratqshaut_omp
     794    iflag_cldcon = iflag_cldcon_omp
     795    ip_ebil_phy = ip_ebil_phy_omp
     796    iflag_thermals = iflag_thermals_omp
     797    nsplit_thermals = nsplit_thermals_omp
     798    type_run = type_run_omp
     799    ok_isccp = ok_isccp_omp
     800    lonmin_ins = lonmin_ins_omp
     801    lonmax_ins = lonmax_ins_omp
     802    latmin_ins = latmin_ins_omp
     803    latmax_ins = latmax_ins_omp
     804    ecrit_hf   = ecrit_hf_omp
     805    ecrit_day = ecrit_day_omp
     806    ecrit_mth = ecrit_mth_omp
     807    ecrit_tra = ecrit_tra_omp
     808    ecrit_reg = ecrit_reg_omp
     809    ok_slab_sicOBS = ok_slab_sicOBS_omp
     810    cvl_corr = cvl_corr_omp
     811    ok_lic_melt = ok_lic_melt_omp
     812
     813!$OMP MASTER
    686814
    687815  write(numout,*)' ##############################################'
    … …  
    703831  write(numout,*)' CFC11_ppt=',CFC11_ppt,' RCFC11 =  ',RCFC11
    704832  write(numout,*)' CFC12_ppt=',CFC12_ppt,' RCFC12 =  ',RCFC12
     833  write(numout,*)' cvl_corr=', cvl_corr
     834  write(numout,*)'ok_lic_melt=', ok_lic_melt
    705835  write(numout,*)' epmax = ', epmax
    706836  write(numout,*)' ok_adj_ema = ', ok_adj_ema
    … …  
    741871  write(numout,*)' lonmin lonmax latmin latmax bilKP_ins =',&
    742872 & lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
    743   write(numout,*)' ecrit_ hf, day, mth, reg, tra',&
     873  write(numout,*)' ecrit_ hf, day, mth, reg',&
    744874 & ecrit_hf, ecrit_day, ecrit_mth, ecrit_reg, ecrit_tra
     875
     876!$OMP_END_MASTER
     877
    745878  return
    746 
     879 
    747880  end subroutine conf_phys
    748881
    … …  
    761894
    762895  REAL          :: tau_calv
     896  REAL,SAVE     :: tau_calv_omp
    763897
    764898! Local
    … …  
    770904!Config Help =
    771905!
    772   tau_calv = 360.
    773   call getin('tau_calv',tau_calv)
    774 
     906  tau_calv_omp = 360.*10.
     907!$OMP MASTER
     908  call getin('tau_calv',tau_calv_omp)
     909!$OMP END MASTER
     910!$OMP BARRIER
     911
     912  tau_calv=tau_calv_omp
     913 
     914!$OMP MASTER
    775915  write(numout,*)' ##############################################'
    776916  WRITE(numout,*)' Configuration de l''interface atm/surfaces  : '
    777917  WRITE(numout,*)' tau_calv = ',tau_calv
     918!$OMP END MASTER
     919
    778920  return
    779921
Note: See TracChangeset for help on using the changeset viewer.