Changeset 5687

Timestamp:

May 28, 2025, 2:07:56 PM (5 weeks ago)

Author:

acozic

Message:

Add new flag to coupled N deposition between atm and ocean
coupling initialized in lmdz but send is in inca code
can be use only with inca chemistry

Location:

LMDZ6/trunk/libf/phylmd

Files:

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

LMDZ6/trunk/libf/phylmd/chemistry_cycle_mod.f90

@@ -21 +21 @@
 !$OMP THREADPRIVATE(n2o_cycle_cpl)   
 
+  LOGICAL, SAVE, PUBLIC :: ndp_cycle_cpl       ! Coupling of N deposition between LMDZ/ORCHIDEE and LMDZ/OCEAN(PISCES)
+!$OMP THREADPRIVATE(ndp_cycle_cpl)
 
   

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
+    USE chemistry_cycle_mod, ONLY: dms_cycle_cpl, n2o_cycle_cpl, ndp_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
+    LOGICAL, SAVE :: dms_cycle_cpl_omp, n2o_cycle_cpl_omp, ndp_cycle_cpl_omp
     INTEGER, SAVE :: level_coupling_esm_omp
     LOGICAL, SAVE :: read_fco2_ocean_cor_omp
@@ -2401 +2401 @@
     !
     !Config Key  = n2o_cycle_cpl
-    !Config Desc = receive n2o from pisces via oasis
+    !Config Desc = receive and send n2o from and to pisces via oasis
     !Config Def  = .FALSE.
     !Config Help = can be activated only if lmdz use inca for the atmospheric chemistry 
-    
+   
+    ndp_cycle_cpl_omp=.FALSE.
+    CALL getin('transm_ndp_ao',ndp_cycle_cpl_omp)
+    !
+    !Config Key  = ndp_cycle_cpl
+    !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 
     
     !$OMP END MASTER
@@ -2704 +2711 @@
     dms_cycle_cpl = dms_cycle_cpl_omp
     n2o_cycle_cpl = n2o_cycle_cpl_omp
+    ndp_cycle_cpl = ndp_cycle_cpl_omp
     !rajout Olivier Torres
     kz0=kz0_omp
@@ -2829 +2837 @@
 
     ! Test on chemistry cycle 
-    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)
+    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)
     ENDIF
     
@@ -3066 +3074 @@
 !GG
     WRITE(lunout,*) ' n2o_cycle_cpl = ', n2o_cycle_cpl
+    WRITE(lunout,*) ' ndp_cycle_cpl = ', ndp_cycle_cpl
     WRITE(lunout,*) ' iflag_tsurf_inlandsis = ', iflag_tsurf_inlandsis
     WRITE(lunout,*) ' iflag_temp_inlandsis = ', iflag_temp_inlandsis

LMDZ6/trunk/libf/phylmd/oasis.F90

@@ -62 +62 @@
        ids_dser = 33, ids_dt_ds = 34
   INTEGER, PARAMETER :: ids_atmn2o = 35
+  INTEGER, PARAMETER :: ids_atmndp = 36
   
-  INTEGER, PARAMETER :: maxsend    = 35  ! Maximum number of fields to send
+  INTEGER, PARAMETER :: maxsend    = 36  ! Maximum number of fields to send
   INTEGER, PARAMETER :: maxsend_phys = 34 ! Maximum number of fields to send in LMDZ phys - the last one will be send by Inca 
 
@@ -117 +118 @@
     USE carbon_cycle_mod, ONLY : carbon_cycle_cpl
     use wxios_mod, ONLY : wxios_context_init
-    USE chemistry_cycle_mod, ONLY : dms_cycle_cpl, n2o_cycle_cpl
+    USE chemistry_cycle_mod, ONLY : dms_cycle_cpl, n2o_cycle_cpl, ndp_cycle_cpl
     USE lmdz_xios  
     USE print_control_mod, ONLY: lunout
@@ -232 +233 @@
             infosend(ids_atmn2o)%action = .TRUE. ; infosend(ids_atmn2o)%name = 'COATMN2O'
         ENDIF
+        IF (ndp_cycle_cpl) THEN
+            infosend(ids_atmndp)%action = .TRUE. ; infosend(ids_atmndp)%name = 'COATMNDP'
+        ENDIF
         infosend(ids_qraioc)%action = .TRUE. ; infosend(ids_qraioc)%name = 'COQRAIOC'
         infosend(ids_qsnooc)%action = .TRUE. ; infosend(ids_qsnooc)%name = 'COQSNOOC'
@@ -398 +402 @@
 
 IF (CPPKEY_INCA) THEN
-    IF (dms_cycle_cpl .OR. n2o_cycle_cpl) THEN
-       CALL init_inca_oasis(inforecv(idr_ocedms:idr_ocen2o),infosend(ids_atmn2o:ids_atmn2o))
+    IF (dms_cycle_cpl .OR. n2o_cycle_cpl .OR. ndp_cycle_cpl) THEN
+       CALL init_inca_oasis(inforecv(idr_ocedms:idr_ocen2o),infosend(ids_atmn2o:ids_atmndp))
     ENDIF
 END IF
@@ -457 +461 @@
              .AND. ierror.NE.PRISM_Input .AND. ierror.NE.PRISM_RecvOut &
              .AND. ierror.NE.PRISM_FromRestOut) THEN
-              WRITE (lunout,*)  'Error with receiving filed : ', inforecv(i)%name, ktime   
+              WRITE (lunout,*)  'Error with receiving field: ', inforecv(i)%name, ktime   
               abort_message=' Problem in prism_get_proto '
               CALL abort_physic(modname,abort_message,1)