Changeset 4098 for LMDZ6/trunk


Ignore:
Timestamp:
Mar 15, 2022, 2:00:42 PM (2 years ago)
Author:
oboucher
Message:

tidying up and correcting a few omissions for nqo > 3

File:
1 edited

Legend:

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

    r4089 r4098  
    13071307       IF (ok_ice_sursat.AND.(nqo.NE.4)) THEN
    13081308          WRITE (lunout, *) ' ok_ice_sursat=y requires 4 H2O tracers ', &
    1309                '(H2Ov, H2Ol, H2Oi, rnebi) but nqo=', nqo, '. Might as well stop here.'
     1309               '(H2Ov, H2Ol, H2Oi, H2Or) but nqo=', nqo, '. Might as well stop here.'
    13101310          abort_message='see above'
    13111311          CALL abort_physic(modname,abort_message,1)
     
    22522252          ql_seri(i,k) = qx(i,k,iliq)
    22532253          !CR: ATTENTION, on rajoute la variable glace
    2254           IF (nqo.eq.2) THEN
     2254          IF (nqo.EQ.2) THEN             !--vapour and liquid only
    22552255             qs_seri(i,k) = 0.
    22562256             rneb_seri(i,k) = 0.
    2257           ELSE IF (nqo.eq.3) THEN
     2257          ELSE IF (nqo.EQ.3) THEN        !--vapour, liquid and ice
    22582258             qs_seri(i,k) = qx(i,k,isol)
    22592259             rneb_seri(i,k) = 0.
    2260           ELSE IF (nqo.eq.4) THEN
     2260          ELSE IF (nqo.EQ.4) THEN        !--vapour, liquid, ice and rneb
    22612261             qs_seri(i,k) = qx(i,k,isol)
    22622262             rneb_seri(i,k) = qx(i,k,irneb)
     
    36333633       ENDDO
    36343634    ENDDO
    3635     IF (nqo==3) THEN
     3635    IF (nqo >= 3) THEN
    36363636    DO k = 1, klev
    36373637       DO i = 1, klon
     
    49244924    ELSE
    49254925       sh_in(:,:) = qx(:,:,ivap)
    4926        IF (nqo .EQ. 3) THEN
     4926       IF (nqo >= 3) THEN
    49274927          ch_in(:,:) = qx(:,:,iliq) + qx(:,:,isol)
    49284928       ELSE
     
    51145114          d_qx(i,k,iliq) = ( ql_seri(i,k) - qx(i,k,iliq) ) / phys_tstep
    51155115          !CR: on ajoute le contenu en glace
    5116           IF (nqo.ge.3) THEN
     5116          IF (nqo >= 3) THEN
    51175117             d_qx(i,k,isol) = ( qs_seri(i,k) - qx(i,k,isol) ) / phys_tstep
    51185118          ENDIF
    51195119          !--ice_sursat: nqo=4, on ajoute rneb
    5120           IF (nqo.eq.4) THEN
     5120          IF (nqo == 4) THEN
    51215121             d_qx(i,k,irneb) = ( rneb_seri(i,k) - qx(i,k,irneb) ) / phys_tstep
    51225122          ENDIF
Note: See TracChangeset for help on using the changeset viewer.