Changeset 5773 for LMDZ6/trunk/libf/phylmd/conf_phys_m.f90

Timestamp:

Jul 11, 2025, 11:11:36 AM (3 weeks ago)

Author:

acozic

Message:

Add possibility to couple nh3 between pisces and inca
use oasis initialization in lmdz
dev by T. Lurton

File:

• LMDZ6/trunk/libf/phylmd/conf_phys_m.f90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmd/conf_phys_m.f90

@@ -36 +36 @@
     USE carbon_cycle_mod,  ONLY: read_fco2_ocean_cor,var_fco2_ocean_cor
     USE carbon_cycle_mod,  ONLY: read_fco2_land_cor,var_fco2_land_cor
-    USE chemistry_cycle_mod, ONLY: dms_cycle_cpl, n2o_cycle_cpl, ndp_cycle_cpl
+    USE chemistry_cycle_mod, ONLY: dms_cycle_cpl, n2o_cycle_cpl, ndp_cycle_cpl, nh3_cycle_cpl
     USE mod_grid_phy_lmdz, ONLY: klon_glo
     USE print_control_mod, ONLY: lunout
@@ -240 +240 @@
     LOGICAL, SAVE :: carbon_cycle_cpl_omp
     LOGICAL, SAVE :: carbon_cycle_rad_omp
-    LOGICAL, SAVE :: dms_cycle_cpl_omp, n2o_cycle_cpl_omp, ndp_cycle_cpl_omp
+    LOGICAL, SAVE :: dms_cycle_cpl_omp, n2o_cycle_cpl_omp, ndp_cycle_cpl_omp, nh3_cycle_cpl_omp
     INTEGER, SAVE :: level_coupling_esm_omp
     LOGICAL, SAVE :: read_fco2_ocean_cor_omp
@@ -2411 +2411 @@
     !Config Desc = send ndep to pisces via oasis
     !Config Def  = .FALSE.
-    !Config Help = can be activated only if lmdz use inca for the atmospheric chemistry 
+    !Config Help = can be activated only if lmdz use inca for the atmospheric chemistry
+
+    
+    nh3_cycle_cpl_omp=.FALSE.
+    CALL getin('transm_nh3_oa',nh3_cycle_cpl_omp)
+    !
+    !Config Key  = nh3_cycle_cpl
+    !Config Desc = send nh3 from and to pisces via oasis
+    !Config Def  = .FALSE.
+    !Config Help = can be activated only if lmdz use inca for the atmospheric chemistry
     
     !$OMP END MASTER
@@ -2712 +2721 @@
     n2o_cycle_cpl = n2o_cycle_cpl_omp
     ndp_cycle_cpl = ndp_cycle_cpl_omp
+    nh3_cycle_cpl = nh3_cycle_cpl_omp
     !rajout Olivier Torres
     kz0=kz0_omp
@@ -2837 +2847 @@
 
     ! Test on chemistry cycle 
-    IF ((type_trac .ne. "inca" .AND. type_trac .ne. "inco") .AND. ( dms_cycle_cpl .OR. n2o_cycle_cpl .OR. ndp_cycle_cpl )  ) THEN 
-       CALL abort_physic('conf_phys', 'dms_cycle_cpl, n2o_cycle_cpl or ndp_cycle_cpl can be TRUE only with INCA coupling model',1)
+    IF ((type_trac .ne. "inca" .AND. type_trac .ne. "inco") .AND. ( dms_cycle_cpl .OR. n2o_cycle_cpl .OR. ndp_cycle_cpl .OR. nh3_cycle_cpl )  ) THEN 
+       CALL abort_physic('conf_phys', 'transm_dms_oa, transm_n2o_oa, transm_ndp_ao or transm_nh3_oa can be TRUE only with INCA coupling model',1)
     ENDIF
     
@@ -3075 +3085 @@
     WRITE(lunout,*) ' n2o_cycle_cpl = ', n2o_cycle_cpl
     WRITE(lunout,*) ' ndp_cycle_cpl = ', ndp_cycle_cpl
+    WRITE(lunout,*) ' nh3_cycle_cpl = ', nh3_cycle_cpl
     WRITE(lunout,*) ' iflag_tsurf_inlandsis = ', iflag_tsurf_inlandsis
     WRITE(lunout,*) ' iflag_temp_inlandsis = ', iflag_temp_inlandsis