Changeset 4755 for LMDZ6/trunk/libf/phylmd/StratAer

Timestamp:

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

Author:

lebasn

Message:

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

File:

• LMDZ6/trunk/libf/phylmd/StratAer/stratH2O_methox.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/StratAer/stratH2O_methox.F90

@@ -4 +4 @@
 SUBROUTINE stratH2O_methox(debutphy,paprs,dq_ch4mmr)
 !
-! output: CH4VMR in mmr (mass mixing ratio/sec: kg H2O/kg air)
+! output: CH4VMR in MMR/s (mass mixing ratio/s or kg H2O/kg air/s)
 
   USE netcdf95, ONLY: nf95_close, nf95_gw_var, nf95_inq_dimid, & 
@@ -26 +26 @@
   include "dimensions.h"
   
-! Variable input
+! Input variables
   REAL paprs(klon,klev+1)
-  LOGICAL, INTENT(IN) :: debutphy   ! le flag de l'initialisation de la physique
-! Variable output
-! tendance buffer pour appel de add_phys_tend
+  LOGICAL, INTENT(IN) :: debutphy   ! flag for first physiq step
+! Output variables
+! tendency buffer used in add_phys_tend subroutine (in physiq_mod)
   REAL, INTENT(INOUT), DIMENSION(klon,klev)  :: dq_ch4mmr
   
-! Variables locales
+! Local variables
   INTEGER n_lat   ! number of latitudes in the input data
   INTEGER n_lon   ! number of longitudes in the input data
@@ -42 +42 @@
 !$OMP THREADPRIVATE(mth_pre)
   
-! Champs reconstitues
+! Reconstitutes fields
   REAL paprs_glo(klon_glo,klev+1)
   
@@ -54 +54 @@
 ! levels of input data
   
-!stratospheric H2O source from CH4 oxidation
-! fixed climatos
-! H2O production in VMR/sec)
+! Stratospheric H2O source from CH4 oxidation from fixed climatos
+! (H2O production in VMR/sec)
   REAL, ALLOCATABLE :: CH4RVMR_in(:, :, :, :)
   REAL, ALLOCATABLE :: CH4RVMR_mth(:, :, :)
@@ -122 +121 @@
         !---correct latitudinal order,convert input from volume mixing ratio to mass mixing ratio
         DO j=1,n_lat
-           ! convert VMR/sec in mmr (mass mixing ratio/sec: kg H2O/kg air)
+           ! convert VMR/s in MMR/s (mass mixing ratio/s or kg H2O/kg air/s)
            ! x2 because CH4->2*H2O
            CH4RVMR_mth(:,j,:) = 2*CH4RVMR_in(:,n_lat+1-j,:,mth_cur)*mH2Omol/mAIRmol