Changeset 4530


Ignore:
Timestamp:
May 11, 2023, 8:46:15 PM (13 months ago)
Author:
evignon
Message:

commission suite a l'atelier nuages du 10 mai 2023:
ajouter les fractions precipitantes de ludo dans les sorties

Location:
LMDZ6/trunk
Files:
9 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/DefLists/field_def_lmdz.xml

    r4523 r4530  
    612612        <field id="rnebcon"    long_name="Convective Cloud Fraction"    unit="-" />
    613613        <field id="rnebls"    long_name="LS Cloud fraction"    unit="-" />
    614         <field id="rneblsvol" long_name="LS Cloud fraction by volume"    unit="-" />
     614        <field id="rneblsvol" long_name="LS Cloud fraction by volume"    unit="-" />
     615        <field id="pfraclr" long_name="LS precip fraction clear-sky part"    unit="-" />               
     616        <field id="pfracld" long_name="LS precip fraction cloudy part"    unit="-" />   
    615617        <field id="rhum"    long_name="Relative humidity"    unit="-" />
    616618        <field id="rhl" long_name="Relative humidity wrt liquid" unit="%" />

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

    r4425 r4530  
    99     paprs,pplay,t,q,ptconv,ratqs,                      &
    1010     d_t, d_q, d_ql, d_qi, rneb, rneblsvol, rneb_seri,  &
    11      radocond, radicefrac, rain, snow,                    &
     11     pfraclr,pfracld,                                   &
     12     radocond, radicefrac, rain, snow,                  &
    1213     frac_impa, frac_nucl, beta,                        &
    1314     prfl, psfl, rhcl, zqta, fraca,                     &
    … …  
    154155  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: rneb             ! cloud fraction [-]
    155156  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: rneblsvol        ! cloud fraction per unit volume [-] 
     157  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: pfraclr          ! precip fraction clear-sky part [-]
     158  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: pfracld          ! precip fraction cloudy part [-]
    156159  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: radocond         ! condensed water used in the radiation scheme [kg/kg]
    157160  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: radicefrac       ! ice fraction of condensed water for radiation scheme
    … …  
    297300d_qi(:,:) = 0.0
    298301rneb(:,:) = 0.0
     302pfraclr(:,:)=0.0
     303pfracld(:,:)=0.0
    299304radocond(:,:) = 0.0
    300305radicefrac(:,:) = 0.0
    … …  
    11411146    ! Outputs:
    11421147    ! Precipitation fluxes at layer interfaces
    1143     ! and temperature and water species tendencies
     1148    ! + precipitation fractions +
     1149    ! temperature and water species tendencies
    11441150    DO i = 1, klon
    11451151        psfl(i,k)=zifl(i)
    11461152        prfl(i,k)=zrfl(i)
     1153        pfraclr(i,k)=znebprecipclr(i)
     1154        pfracld(i,k)=znebprecipcld(i)
    11471155        d_ql(i,k) = (1-zfice(i))*zoliq(i)
    11481156        d_qi(i,k) = zfice(i)*zoliq(i)

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

    r4523 r4530  
    476476      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: rneb,rnebjn,rneblsvol
    477477!$OMP THREADPRIVATE(rneb,rnebjn,rneblsvol)
     478      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: pfraclr,pfracld
     479!$OMP THREADPRIVATE(pfraclr,pfracld)
    478480
    479481! variables de sorties MM
    … …  
    884886      ALLOCATE(beta_prec(klon,klev))
    885887      ALLOCATE(rneb(klon,klev),rnebjn(klon,klev),rneblsvol(klon,klev))
     888      ALLOCATE(pfraclr(klon,klev),pfracld(klon,klev))
     889      pfraclr(:,:)=0. ; pfracld(:,:)=0. ! because not always defined
    886890
    887891
    … …  
    11911195      DEALLOCATE(beta_prec)
    11921196      DEALLOCATE(rneb)
     1197      DEALLOCATE(pfraclr,pfracld)
    11931198      DEALLOCATE (zxsnow,snowhgt,qsnow,to_ice,sissnow,runoff,albsol3_lic)
    11941199

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

    r4523 r4530  
    14731473  TYPE(ctrl_out), SAVE :: o_rneblsvol = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11, 11/), &
    14741474    'rneblsvol', 'LS Cloud fraction by volume', '-', (/ ('', i=1, 10) /))
     1475  TYPE(ctrl_out), SAVE :: o_pfraclr = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
     1476    'pfraclr', 'LS precipitation fraction clear-sky part', '-', (/ ('', i=1, 10) /))
     1477  TYPE(ctrl_out), SAVE :: o_pfracld = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
     1478    'pfracld', 'LS precipitation fraction cloudy part', '-', (/ ('', i=1, 10) /))
    14751479  TYPE(ctrl_out), SAVE :: o_rhum = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11, 11/), &
    14761480    'rhum', 'Relative humidity', '-', (/ ('', i=1, 10) /))

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

    r4525 r4530  
    139139         o_zfull, o_zhalf, o_rneb, o_rnebjn, o_rnebcon, &
    140140         o_rnebls, o_rneblsvol, o_rhum, o_rhl, o_rhi, o_ozone, o_ozone_light, &
     141         o_pfraclr, o_pfracld, &
    141142         o_duphy, o_dtphy, o_dqphy, o_dqphy2d, o_dqlphy, o_dqlphy2d, &
    142143         o_dqsphy, o_dqsphy2d, o_dqbsphy, o_dqbsphy2d, o_albe_srf, o_z0m_srf, o_z0h_srf, &
    … …  
    330331         ql_seri, qs_seri, qbs_seri, tr_seri, qbs_seri,&
    331332         zphi, u_seri, v_seri, omega, cldfra, &
    332          rneb, rnebjn, rneblsvol, zx_rh, zx_rhl, zx_rhi, d_t_dyn,  &
     333         rneb, rnebjn, rneblsvol, zx_rh, zx_rhl, zx_rhi, &
     334         pfraclr, pfracld, d_t_dyn,  &
    333335         d_q_dyn,  d_ql_dyn, d_qs_dyn, d_qbs_dyn,  &
    334336         d_q_dyn2d,  d_ql_dyn2d, d_qs_dyn2d, d_qbs_dyn2d, &
    … …  
    18661868          CALL histwrite_phy(o_rhi, zx_tmp_fi3d)
    18671869       ENDIF
     1870     
     1871       IF (ok_new_lscp) THEN
     1872           CALL histwrite_phy(o_pfraclr, pfraclr)
     1873           CALL histwrite_phy(o_pfracld, pfracld)
     1874       ENDIF
    18681875
    18691876!--aviation & supersaturation

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

    r4527 r4530  
    293293       !
    294294       rneblsvol, &
     295       pfraclr,pfracld, &
    295296       zqsatl, zqsats, &
    296297       qclr, qcld, qss, qvc, rnebclr, rnebss, gamma_ss, &
    … …  
    36953696         t_seri, q_seri,ptconv,ratqs, &
    36963697         d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, rneblsvol, rneb_seri, &
     3698         pfraclr,pfracld, &
    36973699         radocond, picefra, rain_lsc, snow_lsc, &
    36983700         frac_impa, frac_nucl, beta_prec_fisrt, &

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

    r4523 r4530  
    591591      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: rneb,rnebjn,rneblsvol
    592592!$OMP THREADPRIVATE(rneb,rnebjn,rneblsvol)
     593      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: pfraclr,pfracld
     594!$OMP THREADPRIVATE(pfraclr,pfracld)
    593595
    594596! variables de sorties MM
    … …  
    10761078      ALLOCATE(beta_prec(klon,klev))
    10771079      ALLOCATE(rneb(klon,klev),rnebjn(klon,klev),rneblsvol(klon,klev))
     1080      ALLOCATE(pfraclr(klon,klev),pfracld(klon,klev))
     1081      pfraclr(:,:)=0. ; pfracld(:,:)=0. ! because not always defined
    10781082
    10791083
    … …  
    14471451      DEALLOCATE(beta_prec)
    14481452      DEALLOCATE(rneb)
     1453      DEALLOCATE(pfraclr,pfracld)
    14491454      DEALLOCATE (zxsnow,snowhgt,qsnow,to_ice,sissnow,runoff,albsol3_lic)
    14501455#ifdef ISO

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

    r4524 r4530  
    14731473  TYPE(ctrl_out), SAVE :: o_rneblsvol = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11, 11/), &
    14741474    'rneblsvol', 'LS Cloud fraction by volume', '-', (/ ('', i=1, 10) /))
     1475  TYPE(ctrl_out), SAVE :: o_pfraclr = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
     1476    'pfraclr', 'LS precipitation fraction clear-sky part', '-', (/ ('', i=1, 10) /))
     1477  TYPE(ctrl_out), SAVE :: o_pfracld = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
     1478    'pfracld', 'LS precipitation fraction cloudy part', '-', (/ ('', i=1, 10) /))
    14751479  TYPE(ctrl_out), SAVE :: o_rhum = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11, 11/), &
    14761480    'rhum', 'Relative humidity', '-', (/ ('', i=1, 10) /))

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

    r4523 r4530  
    331331       !
    332332       rneblsvol, &
     333       pfraclr,pfracld, &
    333334       zqsatl, zqsats, &
    334335       qclr, qcld, qss, qvc, rnebclr, rnebss, gamma_ss, &
    … …  
    47904791    CALL lscp(klon,klev,phys_tstep,missing_val,paprs,pplay, &
    47914792         t_seri, q_seri,ptconv,ratqs, &
    4792          d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, rneblsvol, rneb_seri, &
     4793         d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, rneblsvol, rneb_seri, &
     4794         pfraclr,pfracld, &
    47934795         radocond, picefra, rain_lsc, snow_lsc, &
    47944796         frac_impa, frac_nucl, beta_prec_fisrt, &
Note: See TracChangeset for help on using the changeset viewer.