Changeset 6005 for LMDZ6/trunk


Ignore:
Timestamp:
Jan 7, 2026, 8:20:00 PM (2 weeks ago)
Author:
evignon
Message:

poursuite du travail sur ratqs avec des efforts de replayisation

Location:
LMDZ6/trunk/libf
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/dynphy_lonlat/phylmd/etat0phys_netcdf.f90

    r5662 r6005  
    144144                   fact_cldcon, facttemps,ok_newmicro,iflag_radia,      &
    145145                   iflag_cldcon,                                        &
    146                    ratqsbas,ratqshaut,tau_ratqs,            &
    147              !GG      iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs,            &
    148146                   ok_ade, ok_aie, ok_alw, ok_cdnc, ok_volcan, flag_volc_surfstrat,     &
    149147                   aerosol_couple, chemistry_couple, flag_aerosol, flag_aerosol_strat,  &

  • LMDZ6/trunk/libf/phylmd/conf_phys_m.f90

    r5918 r6005  
    1515       callstats,&
    1616       solarlong0,seuil_inversion, &
    17        fact_cldcon, facttemps,ok_newmicro,iflag_radia,&
     17       fact_cldcon,facttemps,ok_newmicro,iflag_radia,&
    1818       iflag_cld_th, &
    19        ratqsbas,ratqshaut,tau_ratqs, &
    2019       ok_ade, ok_aie, ok_alw, ok_cdnc, ok_volcan, flag_volc_surfstrat, aerosol_couple, &
    2120       chemistry_couple, flag_aerosol, flag_aerosol_strat, flag_aer_feedback, &
    … …  
    8685    LOGICAL              :: flag_bc_internal_mixture
    8786    REAL                 :: bl95_b0, bl95_b1
    88     REAL                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut,tau_ratqs
     87    REAL                 :: fact_cldcon,facttemps
    8988    INTEGER              :: iflag_cld_th
    9089     
    … …  
    117116    REAL,SAVE           :: freq_ISCCP_omp, ecrit_ISCCP_omp
    118117    REAL,SAVE           :: freq_COSP_omp, freq_AIRS_omp
    119     REAL,SAVE           :: fact_cldcon_omp, facttemps_omp,ratqsbas_omp
     118    REAL,SAVE           :: fact_cldcon_omp,facttemps_omp
    120119    REAL,SAVE           :: tau_cld_cv_omp, coefw_cld_cv_omp
    121120    INTEGER, SAVE       :: iflag_cld_cv_omp
    122121
    123     REAL, SAVE          :: ratqshaut_omp
    124     REAL, SAVE          :: tau_ratqs_omp
    125122    REAL, SAVE          :: t_coupl_omp
    126123    INTEGER, SAVE       :: iflag_radia_omp
    … …  
    130127    INTEGER, SAVE       :: NSW_omp
    131128    INTEGER, SAVE       :: iflag_cld_th_omp, ip_ebil_phy_omp
    132     INTEGER, SAVE       :: iflag_ratqs_omp
    133129
    134130    REAL, SAVE          :: f_cdrag_ter_omp,f_cdrag_oce_omp
    … …  
    12651261    CALL getin('iflag_clw',iflag_clw_omp)
    12661262   
    1267 
    1268     !Config Key  = iflag_ratqs
    1269     !Config Desc =
    1270     !Config Def  = 1
    1271     !Config Help =
    1272     !
    1273     iflag_ratqs_omp = 1
    1274     CALL getin('iflag_ratqs',iflag_ratqs_omp)
    1275 
    12761263    !
    12771264    !Config Key  = iflag_radia
    … …  
    13761363    CALL getin('coefw_cld_cv',coefw_cld_cv_omp)
    13771364
    1378 
    1379 
    1380 
    13811365    !
    13821366    !Config Key  = fact_cldcon
    … …  
    14051389    ok_newmicro_omp = .TRUE.
    14061390    CALL getin('ok_newmicro',ok_newmicro_omp)
    1407     !
    1408     !Config Key  = ratqsbas
    1409     !Config Desc = 
    1410     !Config Def  = 0.01
    1411     !Config Help =
    1412     !
    1413     ratqsbas_omp = 0.01
    1414     CALL getin('ratqsbas',ratqsbas_omp)
    1415     !
    1416     !Config Key  = ratqshaut
    1417     !Config Desc = 
    1418     !Config Def  = 0.3
    1419     !Config Help =
    1420     !
    1421     ratqshaut_omp = 0.3
    1422     CALL getin('ratqshaut',ratqshaut_omp)
    1423 
    1424     !Config Key  = tau_ratqs
    1425     !Config Desc = 
    1426     !Config Def  = 1800.
    1427     !Config Help =
    1428     !
    1429     tau_ratqs_omp = 1800.
    1430     CALL getin('tau_ratqs',tau_ratqs_omp)
    1431 
    14321391    !
    14331392    !-----------------------------------------------------------------------
    … …  
    26292588    bl95_b0 = bl95_b0_omp
    26302589    bl95_b1 = bl95_b1_omp
     2590    facttemps = facttemps_omp
    26312591    fact_cldcon = fact_cldcon_omp
    2632     facttemps = facttemps_omp
    2633     ratqsbas = ratqsbas_omp
    2634     ratqshaut = ratqshaut_omp
    2635     tau_ratqs = tau_ratqs_omp
    2636 
    26372592    iflag_radia = iflag_radia_omp
    26382593    iflag_rrtm = iflag_rrtm_omp
    … …  
    26442599    tau_cld_cv = tau_cld_cv_omp
    26452600    coefw_cld_cv = coefw_cld_cv_omp
    2646     iflag_ratqs = iflag_ratqs_omp
    26472601    ip_ebil_phy = ip_ebil_phy_omp
    26482602    iflag_thermals = iflag_thermals_omp
    … …  
    29602914    WRITE(lunout,*) ' iflag_albedo = ', iflag_albedo !albedo SB
    29612915    WRITE(lunout,*) ' ok_chlorophyll =',ok_chlorophyll ! albedo SB
    2962     WRITE(lunout,*) ' iflag_ratqs = ', iflag_ratqs
    29632916    WRITE(lunout,*) ' seuil_inversion = ', seuil_inversion
    29642917    WRITE(lunout,*) ' fact_cldcon = ', fact_cldcon
    29652918    WRITE(lunout,*) ' facttemps = ', facttemps
    29662919    WRITE(lunout,*) ' ok_newmicro = ',ok_newmicro
    2967     WRITE(lunout,*) ' ratqsbas = ',ratqsbas
    2968     WRITE(lunout,*) ' ratqshaut = ',ratqshaut
    2969     WRITE(lunout,*) ' tau_ratqs = ',tau_ratqs
    29702920    WRITE(lunout,*) ' top_height = ',top_height
    29712921    WRITE(lunout,*) ' rad_froid = ',rad_froid

  • LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.f90

    r6002 r6005  
    1414  !$OMP THREADPRIVATE(RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RV, RG, RPI, EPS_W)
    1515 
    16   INTEGER, SAVE, PROTECTED :: iflag_ratqs        ! control of ratqs option
    17   !$OMP THREADPRIVATE(iflag_ratqs)
    18  
     16 
    1917  REAL, SAVE, PROTECTED :: seuil_neb=0.001       ! cloud fraction threshold: a cloud can precipitate when exceeded
    2018  !$OMP THREADPRIVATE(seuil_neb)
    … …  
    347345  ! parameters for subgrid variability (ratqs) routines
    348346
     347  INTEGER, SAVE, PROTECTED :: iflag_ratqs = 1               ! control of ratqs option
     348  !$OMP THREADPRIVATE(iflag_ratqs)
     349 
     350  REAL, SAVE, PROTECTED :: ratqshaut = 0.3                  ! ratqs [-] value at the top of the atmosphere
     351  !$OMP THREADPRIVATE(ratqshaut)
     352
     353  REAL, SAVE, PROTECTED :: ratqsbas = 0.01                  ! ratqs [-] value near the ground
     354  !$OMP THREADPRIVATE(ratqsbas)
     355 
     356  REAL, SAVE, PROTECTED :: ratqsdp = 20000.                 ! delta pressure [Pa] that controls the atan variation of ratqs with pressure
     357  !$OMP THREADPRIVATE(ratqsdp)
     358
     359  REAL, SAVE, PROTECTED :: ratqsp0 = 50000.                 ! pressure [Pa] where the atan variation of ratqs with pressure occurs
     360  !$OMP THREADPRIVATE(ratqsp0)
     361
     362  REAL, SAVE, PROTECTED :: tau_ratqs = 1800.                ! timescale [s] used when combining convective and stratiform ratqs
     363  !$OMP THREADPRIVATE(tau_ratqs)
     364
    349365  REAL, SAVE, PROTECTED :: a_ratqs_cv = 1.                  ! a_ratqs_cv and a_ratqs_wake control the contibution of cold pools to
    350366                                                            ! water subgrid variability
    … …  
    381397
    382398        SUBROUTINE lscp_ini(dtime, klon, klev, iflag_thermals, lunout_in,    &
    383                     prt_level_in, ok_ice_supersat_in, iflag_ratqs_in, fl_cor_ebil_in,                        &
    384                     RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in, RVTMP2_in,                                        &
     399                    prt_level_in, ok_ice_supersat_in,  fl_cor_ebil_in,       &
     400                    RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in, RVTMP2_in,        &
    385401                    RTT_in, RD_in, RV_in, RG_in, RPI_in, EPS_W_in)
    386402
    … …  
    390406   REAL, INTENT(IN)      :: dtime
    391407   INTEGER, INTENT(IN)   :: klon, klev, iflag_thermals
    392    INTEGER, INTENT(IN)   :: lunout_in,prt_level_in,iflag_ratqs_in,fl_cor_ebil_in
     408   INTEGER, INTENT(IN)   :: lunout_in,prt_level_in,fl_cor_ebil_in
    393409   LOGICAL, INTENT(IN)   :: ok_ice_supersat_in
    394410   REAL, INTENT(IN)      :: RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in
    … …  
    403419    prt_level=prt_level_in
    404420    fl_cor_ebil=fl_cor_ebil_in
    405     iflag_ratqs=iflag_ratqs_in
    406421    ok_ice_supersat=ok_ice_supersat_in
    407422
    … …  
    531546
    532547    ! for subgrid water variability (ratqs) routines
     548    CALL getin_p('iflag_ratqs', iflag_ratqs)
     549    CALL getin_p('ratqshaut', ratqshaut)
     550    CALL getin_p('ratqsbas', ratqsbas)
     551    CALL getin_p('ratqsdp', ratqshaut)
     552    CALL getin_p('ratqsp0', ratqsp0)
     553    CALL getin_p('tau_ratqs', tau_ratqs)
    533554    CALL getin_p('tau_var', tau_var)
    534555    CALL getin_p('fac_tau', fac_tau)
    … …  
    632653    WRITE(lunout,*) 'lscp_ini, iflag_cloudth_vert_noratqs:', iflag_cloudth_vert_noratqs
    633654    ! for subgrid water variability (ratqs) routines
     655    WRITE(lunout,*) 'lscp_ini, iflag_ratqs:', iflag_ratqs
     656    WRITE(lunout,*) 'lscp_ini, ratqshaut:', ratqshaut
     657    WRITE(lunout,*) 'lscp_ini, ratqsbas:', ratqsbas
     658    WRITE(lunout,*) 'lscp_ini, ratqsdp:', ratqsdp
     659    WRITE(lunout,*) 'lscp_ini, ratqsp0', ratqsp0
     660    WRITE(lunout,*) 'lscp_ini, tau_ratqs:', tau_ratqs
    634661    WRITE(lunout,*) 'lscp_ini, tau_var:', tau_var
    635662    WRITE(lunout,*) 'lscp_ini, fac_tau:', fac_tau
    … …  
    685712    ! of the model and an argument of the routines
    686713
     714    allocate(povariance(klon,klev))
     715    allocate(var_conv(klon,klev))
     716
    687717    IF (klon.EQ.1) THEN
    688718        DO k=1,klev

  • LMDZ6/trunk/libf/phylmd/lmdz_lscp_subgridvarq.f90

    r6002 r6005  
    3232!=======================================================================
    3333SUBROUTINE ratqs_main(klon,klev,nbsrf,is_ter,is_lic,          &
    34            iflag_ratqs,iflag_cld_th,pdtphys,                  &
    35            ratqsbas,ratqshaut,ratqsp0,ratqsdp,                &
    36            pctsrf,s_pblh,zstd,                                &
    37            tau_ratqs,fact_cldcon,wake_s, wake_deltaq,         &
     34           iflag_cld_th,fact_cldcon,pdtphys,                  &
     35           pctsrf,s_pblh,zstd, wake_s, wake_deltaq,           &
    3836           ptconv, clwcon0th, rnebcon0th,                     &
    3937           paprs,pplay,t_seri,q_seri,                         &
    … …  
    4543
    4644USE clouds_gno_mod,     ONLY: clouds_gno
    47 USE lmdz_lscp_ini,      ONLY: prt_level, lunout
     45USE lmdz_lscp_ini,      ONLY: prt_level, lunout, iflag_ratqs
     46USE lmdz_lscp_ini,      ONLY: ratqsbas, ratqshaut
     47USE lmdz_lscp_ini,      ONLY: tau_ratqs, ratqsp0, ratqsdp
    4848
    4949implicit none
    … …  
    7171integer, intent(in) :: nbsrf,is_ter,is_lic ! number of subgrid tiles and indices for land and landice
    7272integer, intent(in) :: iflag_cld_th        ! flag that controls cloud properties in presence of thermals
    73 integer, intent(in) :: iflag_ratqs         ! flag that controls ratqs options
    74 
     73real,intent(in)     :: fact_cldcon         ! factor for convective clouds [-]
    7574real,intent(in)     :: pdtphys             ! physics time step [s]
    76 
    77 real, intent(in)    :: ratqsbas,ratqshaut
    78 real, intent(in)    :: fact_cldcon,tau_ratqs
    79 real,intent(in)     :: ratqsp0, ratqsdp
    80 
    81 
    8275real, dimension(klon,klev+1), intent(in) :: paprs ! pressure at layer interfaces [Pa]
    8376real, dimension(klon,klev), intent(in)   :: pplay ! pressure at middle of layers [Pa]
    … …  
    8780real, dimension(klon,klev), intent(in)   :: entr_therm ! thermal plume entrainment rate * dz [kg/s/m2]
    8881real, dimension(klon,klev), intent(in)   :: detr_therm ! thermal plume detrainment rate * dz [kg/s/m2]
    89 real, dimension(klon,klev), intent(in)   :: qtc_cv     !
    90 real, dimension(klon,klev), intent(in)   :: sigt_cv    !
    91 real, dimension(klon,klev), intent(in)   :: detrain_cv ! deep convection detrainment
     82real, dimension(klon,klev), intent(in)   :: qtc_cv     ! total specific moisture in convective saturated draughts  [kg/kg]
     83real, dimension(klon,klev), intent(in)   :: sigt_cv    ! surface fraction of convective saturated draughts [-]
     84real, dimension(klon,klev), intent(in)   :: detrain_cv ! deep convection detrainment of specific humidity variance [kg/s/m2*(kg/kg)^2]
    9285real, dimension(klon,klev), intent(in)   :: fm_cv  ! deep convective mass flux [kg/s/m2]
    9386real, dimension(klon,klev), intent(in)   :: fqd    ! specific humidity tendency due to convective precip [kg/kg/s]
    … …  
    110103
    111104! Output
    112 real, dimension(klon,klev), intent(out) :: ratqs_inter_ ! interactive ratqs
    113 real, dimension(klon,klev), intent(out) :: sigma_qtherm ! standard deviation of humidity in thermals [kg/kg]
     105real, dimension(klon,klev), intent(out) :: ratqs_inter_  ! interactive ratqs
     106real, dimension(klon,klev), intent(out) :: sigma_qtherm  ! standard deviation of humidity in thermals [kg/kg]
    114107real, dimension(klon,klev), intent(out)   :: rnebcon0th  ! cloud fraction associated with thermal updrafts (old method) [-]
    115108
    … …  
    127120
    128121!   old-style convective ratqs computation as a function of q(z=0)-q / q
    129 !   on ecrase le tableau ratqsc calcule par clouds_gno
     122!   the ratqsc computed by clouds_gno is replaced
    130123      if (iflag_cld_th.eq.1) then
    131124         do k=1,klev
    … …  
    141134
    142135!  through log-normal distribution inversion
    143 !-----------------------------------------------------------------------
    144136      else if (iflag_cld_th.eq.4) then
    145137         ptconvthfalse(:,:)=.false.
    … …  
    153145
    154146!   stratiform ratqs
    155 !---------------------
    156 
    157147      if (iflag_ratqs.eq.0) then
    158148
    … …  
    261251
    262252      if (iflag_cld_th.eq.1 .or.iflag_cld_th.eq.2.or.iflag_cld_th.eq.4) then
    263 
    264 ! We add a small constant value to ratqsc*2 to account for small-scale fluctuations
     253         ! We add a small constant value to ratqsc*2 to account for small-scale fluctuations
    265254         do k=1,klev
    266255            do i=1,klon
    … …  
    272261
    273262!   ratqs are a combination of ratqss et ratqsc
    274        if(prt_level.ge.9) write(lunout,*)'PHYLMD NEW TAU_RATQS ',tau_ratqs
     263         if(prt_level.ge.9) write(lunout,*)'PHYLMD NEW TAU_RATQS ',tau_ratqs
    275264
    276265         if (tau_ratqs>1.e-10) then
    … …  
    282271         ratqs(:,:)=max(ratqs(:,:),ratqss(:,:))
    283272      else if (iflag_cld_th<=6) then
    284 !  we only keep the stable ratqs for lscp
     273        ! we only keep the stable ratqs for lscp
    285274         ratqs(:,:)=ratqss(:,:)
    286275      else
    … …  
    339328real, dimension(klon,klev), intent(in)   :: entr_therm ! thermal plume entrainment rate * dz [kg/s/m2]
    340329real, dimension(klon,klev), intent(in)   :: detr_therm ! thermal plume detrainment rate * dz [kg/s/m2]
    341 real, dimension(klon,klev), intent(in)   :: qtc_cv     !
    342 real, dimension(klon,klev), intent(in)   :: sigt_cv    !
    343 real, dimension(klon,klev), intent(in)   :: detrain_cv ! deep convection detrainment
     330real, dimension(klon,klev), intent(in)   :: qtc_cv     ! total specific moisture in convective saturated draughts  [kg/kg]
     331real, dimension(klon,klev), intent(in)   :: sigt_cv    ! surface fraction of convective saturated draughts [-]
     332real, dimension(klon,klev), intent(in)   :: detrain_cv ! deep convection detrainment of specific humidity variance [kg/s/m2*(kg/kg)^2]
    344333real, dimension(klon,klev), intent(in)   :: fm_cv  ! deep convective mass flux [kg/s/m2]
    345334real, dimension(klon,klev), intent(in)   :: fqd    ! specific humidity tendency due to convective precip [kg/kg/s]

  • LMDZ6/trunk/libf/phylmd/physiq_mod.F90

    r6003 r6005  
    10261026    REAL pde_u(klon,klev), pde_d(klon,klev)
    10271027    INTEGER kcbot(klon), kctop(klon), kdtop(klon)
    1028     !
    1029     REAL ratqsbas,ratqshaut,tau_ratqs
    1030     SAVE ratqsbas,ratqshaut,tau_ratqs
    1031     !$OMP THREADPRIVATE(ratqsbas,ratqshaut,tau_ratqs)
    1032     REAL, SAVE :: ratqsp0=50000., ratqsdp=20000.
    1033     !$OMP THREADPRIVATE(ratqsp0, ratqsdp)
    10341028
    10351029    ! Parametres lies au nouveau schema de nuages (SB, PDF)
    1036     REAL, SAVE :: fact_cldcon
    1037     REAL, SAVE :: facttemps
    1038     !$OMP THREADPRIVATE(fact_cldcon,facttemps)
     1030    REAL, SAVE :: facttemps, fact_cldcon
     1031    !$OMP THREADPRIVATE(facttemps,fact_cldcon)
    10391032    LOGICAL, SAVE :: ok_newmicro
    10401033    !$OMP THREADPRIVATE(ok_newmicro)
    1041 
    10421034    INTEGER, SAVE :: iflag_cld_th
    10431035    !$OMP THREADPRIVATE(iflag_cld_th)
    1044     !IM logical ptconv(klon,klev)  !passe dans phys_local_var_mod
    1045     !IM cf. AM 081204 BEG
    10461036    LOGICAL ptconvth(klon,klev)
    1047 
    10481037    REAL picefra(klon,klev)
    10491038    REAL nm_oro(klon)
    1050     !IM cf. AM 081204 END
    1051     !
    10521039    ! Variables liees a l'ecriture de la bande histoire physique
    10531040    !
    … …  
    13921379            callstats, &
    13931380            solarlong0,seuil_inversion, &
    1394             fact_cldcon, facttemps,ok_newmicro,iflag_radia, &
    1395             iflag_cld_th,ratqsbas,ratqshaut,tau_ratqs, &
     1381            fact_cldcon,facttemps,ok_newmicro,iflag_radia, &
     1382            iflag_cld_th, &
    13961383            ok_ade, ok_aie, ok_alw, ok_cdnc, ok_volcan, flag_volc_surfstrat, aerosol_couple, &
    13971384            chemistry_couple, flag_aerosol, flag_aerosol_strat, flag_aer_feedback, &
    … …  
    15161503       CALL getin_p('oliqmax',oliqmax)
    15171504       CALL getin_p('oicemax',oicemax)
    1518        CALL getin_p('ratqsp0',ratqsp0)
    1519        CALL getin_p('ratqsdp',ratqsdp)
    15201505       iflag_wake_tend = 0
    15211506       CALL getin_p('iflag_wake_tend',iflag_wake_tend)
    … …  
    15521537       WRITE(lunout,*) 'qtcon_multistep_max=',dtcon_multistep_max
    15531538       WRITE(lunout,*) 'qdcon_multistep_max=',dqcon_multistep_max
    1554        WRITE(lunout,*) 'ratqsp0=',            ratqsp0
    1555        WRITE(lunout,*) 'ratqsdp=',            ratqsdp
    15561539       WRITE(lunout,*) 'iflag_wake_tend=',    iflag_wake_tend
    15571540       WRITE(lunout,*) 'ok_bad_ecmwf_thermo=',ok_bad_ecmwf_thermo
    … …  
    18921875            &    RG,RD,RCPD,RKAPPA,RLVTT,RETV)
    18931876       CALL lscp_ini(pdtphys,klon,klev,iflag_thermals,lunout,prt_level,       &
    1894                      ok_ice_supersat,iflag_ratqs,fl_cor_ebil,                                     &
     1877                     ok_ice_supersat,fl_cor_ebil,                             &
    18951878                     RCPD,RLSTT,RLVTT,RLMLT,RVTMP2,RTT,RD,RV,RG,RPI,EPS_W)
    18961879       CALL blowing_snow_ini(RCPD, RLSTT, RLVTT, RLMLT, &
    … …  
    39493932    CALL ratqs_main_first(klon, cell_area)
    39503933    CALL ratqs_main(klon,klev,nbsrf,is_ter,is_lic,        &
    3951          iflag_ratqs,iflag_cld_th,pdtphys,  &
    3952          ratqsbas,ratqshaut,ratqsp0, ratqsdp, &
    3953          pctsrf,s_pblh,zstd, &
    3954          tau_ratqs,fact_cldcon,wake_s, wake_deltaq,   &
    3955          ptconv,clwcon0th, rnebcon0th,     &
     3934         iflag_cld_th,fact_cldcon, pdtphys,  &
     3935         pctsrf,s_pblh,zstd, wake_s, wake_deltaq,   &
     3936         ptconv, clwcon0th, rnebcon0th,     &
    39563937         paprs,pplay,t_seri,q_seri, &
    39573938         qtc_cv, sigt_cv,detrain_cv,fm_cv,fqd,fqcomp,sigd,zqsat, &

  • LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

    r6003 r6005  
    11201120    INTEGER kcbot(klon), kctop(klon), kdtop(klon)
    11211121    !
    1122     REAL ratqsbas,ratqshaut,tau_ratqs
    1123     SAVE ratqsbas,ratqshaut,tau_ratqs
    1124     !$OMP THREADPRIVATE(ratqsbas,ratqshaut,tau_ratqs)
    1125     REAL, SAVE :: ratqsp0=50000., ratqsdp=20000.
    1126     !$OMP THREADPRIVATE(ratqsp0, ratqsdp)
    11271122
    11281123    ! Parametres lies au nouveau schema de nuages (SB, PDF)
    1129     REAL, SAVE :: fact_cldcon
    1130     REAL, SAVE :: facttemps
    1131     !$OMP THREADPRIVATE(fact_cldcon,facttemps)
     1124    REAL, SAVE :: facttemps,fact_cldcon
     1125    !$OMP THREADPRIVATE(facttemps,fact_cldcon)
    11321126    LOGICAL, SAVE :: ok_newmicro
    11331127    !$OMP THREADPRIVATE(ok_newmicro)
    … …  
    14961490            callstats, &
    14971491            solarlong0,seuil_inversion, &
    1498             fact_cldcon, facttemps,ok_newmicro,iflag_radia, &
    1499             iflag_cld_th,ratqsbas,ratqshaut,tau_ratqs, &
     1492            fact_cldcon,facttemps,ok_newmicro,iflag_radia, &
     1493            iflag_cld_th, &
    15001494            ok_ade, ok_aie, ok_alw, ok_cdnc, ok_volcan, flag_volc_surfstrat, aerosol_couple, &
    15011495            chemistry_couple, flag_aerosol, flag_aerosol_strat, flag_aer_feedback, &
    … …  
    16261620       CALL getin_p('oliqmax',oliqmax)
    16271621       CALL getin_p('oicemax',oicemax)
    1628        CALL getin_p('ratqsp0',ratqsp0)
    1629        CALL getin_p('ratqsdp',ratqsdp)
    16301622       iflag_wake_tend = 0
    16311623       CALL getin_p('iflag_wake_tend',iflag_wake_tend)
    … …  
    16621654       WRITE(lunout,*) 'qtcon_multistep_max=',dtcon_multistep_max
    16631655       WRITE(lunout,*) 'qdcon_multistep_max=',dqcon_multistep_max
    1664        WRITE(lunout,*) 'ratqsp0=',            ratqsp0
    1665        WRITE(lunout,*) 'ratqsdp=',            ratqsdp
    16661656       WRITE(lunout,*) 'iflag_wake_tend=',    iflag_wake_tend
    16671657       WRITE(lunout,*) 'ok_bad_ecmwf_thermo=',ok_bad_ecmwf_thermo
    … …  
    20352025   &    RG,RD,RCPD,RKAPPA,RLVTT,RETV)
    20362026       CALL lscp_ini(pdtphys,klon,klev,iflag_thermals,lunout,prt_level,       &
    2037                      ok_ice_supersat,iflag_ratqs,fl_cor_ebil,                                     &
     2027                     ok_ice_supersat,fl_cor_ebil,                             &
    20382028                     RCPD,RLSTT,RLVTT,RLMLT,RVTMP2,RTT,RD,RV,RG,RPI,EPS_W)
    20392029       CALL blowing_snow_ini(RCPD, RLSTT, RLVTT, RLMLT, &
    … …  
    51055095
    51065096    CALL ratqs_main_first(klon, cell_area)
    5107     CALL ratqs_main(klon,klev,nbsrf,is_ter,is_lic,        &
    5108          iflag_ratqs,iflag_cld_th,pdtphys,  &
    5109          ratqsbas,ratqshaut,ratqsp0, ratqsdp, &
    5110          pctsrf,s_pblh,zstd, &
    5111          tau_ratqs,fact_cldcon,wake_s, wake_deltaq,   &
     5097    CALL ratqs_main(klon,klev,nbsrf,is_ter,is_lic,  &
     5098         iflag_cld_th,fact_cldcon,pdtphys,          &
     5099         pctsrf,s_pblh,zstd, wake_s, wake_deltaq,   &
    51125100         ptconv,clwcon0th, rnebcon0th,     &
    51135101         paprs,pplay,t_seri,q_seri, &
Note: See TracChangeset for help on using the changeset viewer.