Changeset 4755 for LMDZ6/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

Nov 22, 2023, 3:46:17 PM (11 months ago)

Author:

lebasn

Message:

Methox: Update comments to correct units + bugfix on output variable

File:

• LMDZ6/trunk/libf/phylmd/physiq_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -5025 +5025 @@
     !DC Calcul de la tendance due au methane
     IF (ok_qch4) THEN
-!      d_q_ch4: H2O source in ppmv/sec
+!      d_q_ch4: H2O source from CH4 in MMR/s (mass mixing ratio/s or kg H2O/kg air/s)
 #ifdef CPP_StratAer
        CALL stratH2O_methox(debut,paprs,d_q_ch4)
 #else
-!      ecmwf routine METHOX
+!      ECMWF routine METHOX
        CALL METHOX(1,klon,klon,klev,q_seri,d_q_ch4,pplay)
 #endif
-       ! ajout de la tendance d'humidite due au methane
+       ! add humidity tendency due to methane
        d_q_ch4_dtime(:,:) = d_q_ch4(:,:)*phys_tstep
        CALL add_phys_tend(du0, dv0, dt0, d_q_ch4_dtime, dql0, dqi0, dqbs0, paprs, &
             'q_ch4', abortphy,flag_inhib_tend,itap,0)
-       d_q_ch4(:,:) = d_q_ch4_dtime(:,:)/phys_tstep
+       d_q_ch4(:,:) = d_q_ch4_dtime(:,:)/phys_tstep ! update with H2O conserv done in add_phys_tend
     ENDIF
     !
@@ -5050 +5050 @@
        
        SELECT CASE(flag_emit)
-       CASE(1) ! emission volc H2O dans LMDZ
+       CASE(1) ! emission volc H2O in LMDZ
           DO ieru=1, nErupt
              IF (year_cur==year_emit_vol(ieru).AND.&
@@ -5058 +5058 @@
                 
                 IF(flag_verbose_strataer) print *,'IN physiq_mod: date=',year_cur,mth_cur,day_cur
-                ! initialisation tendance q emission
+                ! initialisation of q tendency emission
                 d_q_emiss(:,:)=0.
                 ! daily injection mass emission - NL