Changeset 4754 for LMDZ6

Timestamp:

Nov 21, 2023, 2:51:40 PM (13 months ago)

Author:

acozic

Message:

Add possibility to coupled n2o between ocean and atmosphere

Location:

LMDZ6/trunk/libf/phylmd

Files:

• chemistry_cycle_mod.F90 (modified) (1 diff)
• conf_phys_m.F90 (modified) (6 diffs)
• oasis.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/chemistry_cycle_mod.F90

@@ -15 +15 @@
 
 ! Variables read from parmeter file physiq.def
-  LOGICAL, SAVE, PUBLIC :: dms_cycle_cpl       ! Coupling of CO2 fluxes between LMDZ/ORCHIDEE and LMDZ/OCEAN(PISCES) 
+  LOGICAL, SAVE, PUBLIC :: dms_cycle_cpl       ! Coupling of DMS fluxes between LMDZ/ORCHIDEE and LMDZ/OCEAN(PISCES) 
 !$OMP THREADPRIVATE(dms_cycle_cpl)   
+
+  LOGICAL, SAVE, PUBLIC :: n2o_cycle_cpl       ! Coupling of N2O fluxes between LMDZ/ORCHIDEE and LMDZ/OCEAN(PISCES) 
+!$OMP THREADPRIVATE(n2o_cycle_cpl)   
 
 

LMDZ6/trunk/libf/phylmd/conf_phys_m.F90

@@ -29 +29 @@
     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
+    USE chemistry_cycle_mod, ONLY: dms_cycle_cpl, n2o_cycle_cpl
     USE mod_grid_phy_lmdz, ONLY: klon_glo
     USE print_control_mod, ONLY: lunout
@@ -233 +233 @@
     LOGICAL, SAVE :: carbon_cycle_cpl_omp
     LOGICAL, SAVE :: carbon_cycle_rad_omp
-    LOGICAL, SAVE :: dms_cycle_cpl_omp
+    LOGICAL, SAVE :: dms_cycle_cpl_omp, n2o_cycle_cpl_omp
     INTEGER, SAVE :: level_coupling_esm_omp
     LOGICAL, SAVE :: read_fco2_ocean_cor_omp
@@ -2254 +2254 @@
     !Config Help = can be activated only if lmdz use inca for the atmospheric chemistry 
 
+    n2o_cycle_cpl_omp=.FALSE.
+    CALL getin('transm_n2o_oa',n2o_cycle_cpl_omp)
+    !
+    !Config Key  = n2o_cycle_cpl
+    !Config Desc = receive n2o from pisces via oasis
+    !Config Def  = .FALSE.
+    !Config Help = can be activated only if lmdz use inca for the atmospheric chemistry 
+    
     
     !$OMP END MASTER
@@ -2527 +2535 @@
     var_fco2_land_cor = var_fco2_land_cor_omp
     dms_cycle_cpl = dms_cycle_cpl_omp
+    n2o_cycle_cpl = n2o_cycle_cpl_omp
     !rajout Olivier Torres
     kz0=kz0_omp
@@ -2652 +2661 @@
 
     ! Test on chemistry cycle 
-    IF ((type_trac .ne. "inca" .AND. type_trac .ne. "inco") .AND. dms_cycle_cpl ) THEN 
-       CALL abort_physic('conf_phys', 'dms_cycle_cpl has to 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)  ) THEN 
+       CALL abort_physic('conf_phys', 'dms_cycle_cpl or n2o_cycle_cpl has to be TRUE only with INCA coupling model',1)
     ENDIF
     
@@ -2865 +2874 @@
     WRITE(lunout,*) ' var_fco2_land_cor = ', var_fco2_land_cor
     WRITE(lunout,*) ' dms_cycle_cpl = ', dms_cycle_cpl
+    WRITE(lunout,*) ' n2o_cycle_cpl = ', n2o_cycle_cpl
     WRITE(lunout,*) ' iflag_tsurf_inlandsis = ', iflag_tsurf_inlandsis
     WRITE(lunout,*) ' iflag_temp_inlandsis = ', iflag_temp_inlandsis

LMDZ6/trunk/libf/phylmd/oasis.F90

@@ -77 +77 @@
   INTEGER, PARAMETER :: idr_sss = 9
   INTEGER, PARAMETER :: idr_ocedms = 10
-
-  INTEGER, PARAMETER :: maxrecv      = 10     ! Maximum number of fields to receive
+  INTEGER, PARAMETER :: idr_ocen2o = 11
+
+  INTEGER, PARAMETER :: maxrecv      = 11     ! Maximum number of fields to receive
   INTEGER, PARAMETER :: maxrecv_phys = 9      ! Maximum number of fields to receive in physiq (without fields received in INCA model )
                                               ! will be changed in next version - INCA fields will be received in LMDZ (like for ORCHIDEE fields)
@@ -114 +115 @@
     USE carbon_cycle_mod, ONLY : carbon_cycle_cpl
     USE wxios, ONLY : wxios_context_init 
-    USE chemistry_cycle_mod, ONLY : dms_cycle_cpl
+    USE chemistry_cycle_mod, ONLY : dms_cycle_cpl, n2o_cycle_cpl
     USE lmdz_xios  
     USE print_control_mod, ONLY: lunout
@@ -265 +266 @@
       inforecv(idr_ocedms)%action = .TRUE. ; inforecv(idr_ocedms)%name = 'SIDMSFLX'
    ENDIF
+   IF (n2o_cycle_cpl) THEN 
+      inforecv(idr_ocen2o)%action = .TRUE. ; inforecv(idr_ocen2o)%name = 'SIN2OFLX'
+   ENDIF
  
 
@@ -353 +357 @@
 
 #ifdef INCA
-    IF (dms_cycle_cpl) THEN
-       CALL init_inca_oasis(inforecv(idr_ocedms))
+    IF (dms_cycle_cpl .OR. n2o_cycle_cpl) THEN
+       CALL init_inca_oasis(inforecv(idr_ocedms:idr_ocen2o))
     ENDIF
 #endif