Ignore:
Timestamp:
Jan 21, 2022, 4:25:19 PM (3 years ago)
Author:
oboucher
Message:

add changes made to physiq_mod in phylmd

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

    r4056 r4060  
    169169       t_seri,q_seri,ql_seri,qs_seri,u_seri,v_seri,tr_seri, &
    170170       ! Dynamic tendencies (diagnostics)
    171        d_t_dyn,d_q_dyn,d_ql_dyn,d_qs_dyn,d_u_dyn,d_v_dyn,d_tr_dyn, &
     171       d_t_dyn,d_q_dyn,d_ql_dyn,d_qs_dyn,d_u_dyn,d_v_dyn,d_tr_dyn,d_rneb_dyn, &
    172172       d_q_dyn2d,d_ql_dyn2d,d_qs_dyn2d, &
    173173       ! Physic tendencies
     
    13891389          WRITE (lunout, *) ' iflag_ice_thermo==1 requires 3 H2O tracers ', &
    13901390               '(H2Ov, H2Ol, H2Oi) but nqo=', nqo, '. Might as well stop here.'
     1391          abort_message='see above'
     1392          CALL abort_physic(modname,abort_message,1)
     1393       ENDIF
     1394
     1395       IF ((iflag_ice_sursat.GT.0).AND.(iflag_ice_thermo.EQ.0)) THEN
     1396          WRITE (lunout, *) ' iflag_ice_sursat=1 requires iflag_ice_thermo=1 as well'
     1397          abort_message='see above'
     1398          CALL abort_physic(modname,abort_message,1)
     1399       ENDIF
     1400
     1401       IF ((iflag_ice_sursat.GT.0).AND.(nqo.NE.4)) THEN
     1402          WRITE (lunout, *) ' iflag_ice_sursat=1 requires 4 H2O tracers ', &
     1403               '(H2Ov, H2Ol, H2Oi, rnebi) but nqo=', nqo, '. Might as well stop here.'
    13911404          abort_message='see above'
    13921405          CALL abort_physic(modname,abort_message,1)
     
    25522565       IF (nqtot > nqo) d_tr_dyn(:,:,:)=(tr_seri(:,:,:)-tr_ancien(:,:,:))/phys_tstep
    25532566       ! !! RomP <<<
     2567       !!d_rneb_dyn(:,:)=(rneb_seri(:,:)-rneb_ancien(:,:))/phys_tstep
     2568       d_rneb_dyn(:,:)=0.0
    25542569
    25552570#ifdef ISO
     
    26542669       ! !! RomP >>>   td dyn traceur
    26552670       IF (nqtot > nqo) d_tr_dyn(:,:,:)= 0.0
     2671       d_rneb_dyn(:,:)=0.0
    26562672       ! !! RomP <<<
    26572673       ancien_ok = .TRUE.
     
    46064622    CALL lscp(phys_tstep,paprs,pplay, &
    46074623         t_seri, q_seri,ptconv,ratqs, &
    4608          d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, cldliq, picefra, &
    4609          rain_lsc, snow_lsc, &
     4624         d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, rneb_seri, &
     4625         cldliq, picefra, rain_lsc, snow_lsc, &
    46104626         pfrac_impa, pfrac_nucl, pfrac_1nucl, &
    46114627         frac_impa, frac_nucl, beta_prec_fisrt, &
    46124628         prfl, psfl, rhcl,  &
    46134629         zqasc, fraca,ztv,zpspsk,ztla,zthl,iflag_cld_th, &
    4614          iflag_ice_thermo)
     4630         iflag_ice_thermo, iflag_ice_sursat)
    46154631
    46164632    ELSE
     
    64736489          d_qx(i,k,iliq) = ( ql_seri(i,k) - qx(i,k,iliq) ) / phys_tstep
    64746490          !CR: on ajoute le contenu en glace
    6475           IF (nqo.eq.3) THEN
     6491          IF (nqo.gt.3) THEN
    64766492             d_qx(i,k,isol) = ( qs_seri(i,k) - qx(i,k,isol) ) / phys_tstep
     6493          ENDIF
     6494          !--ice_sursat: nqo=4, on ajoute rneb
     6495          IF (nqo.eq.4) THEN
     6496             d_qx(i,k,irneb) = ( rneb_seri(i,k) - qx(i,k,irneb) ) / phys_tstep
    64776497          ENDIF
    64786498       ENDDO
Note: See TracChangeset for help on using the changeset viewer.