Changeset 5918 for LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

Timestamp:

Dec 8, 2025, 5:23:42 PM (4 weeks ago)

Author:

Sebastien Nguyen

Message:

report FC changes to phylmd/physiq_mod.F90 in phylmdiso/physiq_mod.F90

File:

• LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90 (modified) (14 diffs)

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -1 +1 @@
 
-! $Id: physiq_mod.F90 5842 2025-09-25 15:09:22Z rkazeroni $
+! $Id: physiq_mod.F90 5913 2025-12-04 16:00:50Z fairhead $
 !
 !#define IO_DEBUG
@@ -31 +31 @@
     USE change_srf_frac_mod
     USE conf_phys_m, ONLY: conf_phys
-    USE carbon_cycle_mod, ONLY : infocfields_init, RCO2_glo, carbon_cycle_rad
+    USE carbon_cycle_mod, ONLY : infocfields_init, RCO2_glo, carbon_cycle_rad   ! PC
+
     USE CFMIP_point_locations   ! IM stations CFMIP
     USE clouds_bigauss_mod, ONLY: clouds_bigauss
@@ -143 +144 @@
 #endif
 
+    USE ch4n2o_correction_mod, ONLY: ch4n2o_init, ch4n2o_compute, &                       ! PC
+  &                                  ok_CH4_eff_ppb, ok_N2O_eff_ppb, ghg_year, ghg_mode   ! PC
     USE phys_local_var_mod, ONLY: d_q_emiss
     USE strataer_local_var_mod, ONLY: strataer_init,flag_emit,flh2o,ok_qemiss,flag_verbose_strataer, &
@@ -1675 +1678 @@
        WRITE(lunout,*) 'Call to infocfields from physiq'
        CALL infocfields_init
+
+       !--PC:reads keys once: broadcasts config to all ranks
+       CALL ch4n2o_init()    ! reads getin_p, broadcasts configuration 
+
 
        !AI 08 2023
@@ -2176 +2183 @@
        IF (phys_tstep*REAL(radpas).GT.21600..AND.iflag_cycle_diurne.GE.1) THEN
           WRITE(lunout,*)'Nbre d appels au rayonnement insuffisant'
+          WRITE(lunout,*)'phys_tstep =', phys_tstep, radpas, iflag_cycle_diurne
           WRITE(lunout,*)"Au minimum 4 appels par jour si cycle diurne"
           abort_message='Nbre d appels au rayonnement insuffisant'
@@ -3310 +3318 @@
            do i=1,klon                                                       !
               t_env(i,k)=t_seri(i,k)                                         !
-              q_env(i,k)=q_seri(i,k)   
+              q_env(i,k)=q_seri(i,k)                                         !
 #ifdef ISO
               do ixt=1,ntraciso
@@ -3525 +3533 @@
              (d_deltat_vdf, d_deltaq_vdf, dsig0, dsig0, ddens0, ddens0, wkoccur1, 'vdf', abortphy &
 #ifdef ISO
-               ,d_deltaxt_vdf &
-#endif
-                ) 
+             ,d_deltaxt_vdf &
+#endif
+             ) 
        ELSE
           d_deltat_vdf(:,:) = 0.
@@ -3552 +3560 @@
                'vdf',abortphy,flag_inhib_tend,itap,0 &
 #ifdef ISO
-     &    ,d_xt_vdf,dxtl0,dxti0 &
+     &         ,d_xt_vdf,dxtl0,dxti0 &
 #endif      
-     &   )
+     &         )
        ENDIF
        CALL prt_enerbil('vdf',itap)
@@ -3928 +3936 @@
                  ,d_deltaxt_ajs_cv  &
 #endif
-                 ) 
+                 )
              ENDIF  ! (iflag_adjwk == 2 .AND. OK_bug_ajs_cv)
           ENDIF  ! (iflag_adjwk >= 1)
@@ -5043 +5051 @@
                   , d_t_ajsb, d_q_ajsb &
 #ifdef ISO
-     &      ,xt_seri,d_xt_ajsb      &
+     &            ,xt_seri,d_xt_ajsb      &
 #endif         
-     &   )
+     &            )
           ELSE IF (iflag_thermals>0) THEN
              ! Calling adjustment above the top of thermal plumes
@@ -5051 +5059 @@
                   , d_t_ajsb, d_q_ajsb &
 #ifdef ISO
-     &      ,xt_seri,d_xt_ajsb &      
+     &            ,xt_seri,d_xt_ajsb &      
 #endif         
-     &   )
+     &            )
           ENDIF
 
@@ -5061 +5069 @@
                'ajsb',abortphy,flag_inhib_tend,itap,0 &
 #ifdef ISO
-     &    ,d_xt_ajsb,dxtl0,dxti0 &
+     &        ,d_xt_ajsb,dxtl0,dxti0 &
 #endif      
-     &   )
+     &        )
           CALL prt_enerbil('ajsb',itap)
           d_t_ajs(:,:)=d_t_ajs(:,:)+d_t_ajsb(:,:)
@@ -5302 +5310 @@
          'lsc',abortphy,flag_inhib_tend,itap,0 &
 #ifdef ISO
-     &    ,d_xt_lsc,d_xtl_lsc,d_xti_lsc &
+     &   ,d_xt_lsc,d_xtl_lsc,d_xti_lsc &
 #endif      
      &   )
@@ -6136 +6144 @@
           !
           RCO2 = RCO2_act
+!          RCH4 = RCH4_act
+!          RN2O = RN2O_act
+#ifdef CPP_RRTM
+          ! PC
+          ! Only attempt to compute effective values if at least one
+          ! gas is actually routed as "effective"
+          IF (ok_CH4_eff_ppb .OR. ok_N2O_eff_ppb) THEN
+             ! ghg_year is provided by the driver (libIGCM); the transfer
+             ! functions currently ignore it and use it only for logging
+             CALL ch4n2o_compute(ghg_year, CH4_ppb, N2O_ppb, CH4_eff_ppb, N2O_eff_ppb)
+          END IF
+
+          ! CH4: choose which mass mixing ratio to feed RRTM
+          ! if "ok_CH4_eff_ppb" = .TRUE., we convert the effective
+          ! concentration (ppb) into a mass mixing ratio
+          ! Otherwise we use the "actual" RCH4_act 
+          IF (ok_CH4_eff_ppb) THEN
+             RCH4 = CH4_eff_ppb * 1.0E-09 * RMCH4 / RMD
+             !RCH4_eff = CH4_eff_ppb * 1.0E-09 * RMCH4 / RMD
+          ELSE
+          RCH4 = RCH4_act
+          END IF
+
+          ! N2O: choose which mass mixing ratio to feed RRTM
+          ! if "ok_N2O_eff_ppb" = .TRUE., we convert the effective
+          ! concentration (ppb) into a mass mixing ratio
+          ! Othewise we use the "actual" RN2O_act           
+          IF (ok_N2O_eff_ppb) THEN
+             RN2O = N2O_eff_ppb * 1.0E-09 * RMN2O / RMD
+             !RN2O_eff = N2O_eff_ppb * 1.0E-09 * RMN2O / RMD
+          ELSE
+          RN2O = RN2O_act
+          END IF
+          ! PC
+#else
+          IF (ok_CH4_eff_ppb .OR. ok_N2O_eff_ppb .OR. &
+                  TRIM(ghg_mode) /= "identity") THEN
+
+                  abort_message = 'CH4/N2O effective mapping (ghg_mode=' // &
+                          TRIM(ghg_mode) // ') is only supported ' // &
+                          ' for RRTM (CPP_RRTM). Please recompile ' // &
+                          ' with -rrtm or set ok_CH4_eff_ppb=n, ' // &
+                          'ok_N2O_eff_ppb=n and ghg_mode="identity".'
+                  CALL abort_physic(modname,abort_message, 1)
+          ENDIF        
+
+          ! In all cases we feed the unmodified "actual" mixing ratios
+          ! to the radiation scheme when not using RRTM.
           RCH4 = RCH4_act
           RN2O = RN2O_act
+#endif
+
+
           RCFC11 = RCFC11_act
           RCFC12 = RCFC12_act