Changeset 5791 for LMDZ6/branches/contrails/libf/phylmd/oasis.F90

Timestamp:

Jul 28, 2025, 7:23:15 PM (7 days ago)

Author:

aborella

Message:

Merge with trunk r5789

Location:

LMDZ6/branches/contrails

Files:

• . (modified) (1 prop)
• libf/phylmd/oasis.F90 (modified) (11 diffs)

LMDZ6/branches/contrails/libf/phylmd/oasis.F90

@@ -61 +61 @@
   INTEGER, PARAMETER :: ids_delta_sst = 30, ids_delta_sal = 31, ids_dter = 32, &
        ids_dser = 33, ids_dt_ds = 34
+  INTEGER, PARAMETER :: ids_atmn2o = 35
+  INTEGER, PARAMETER :: ids_atmndp = 36
+  INTEGER, PARAMETER :: ids_atmnh3 = 37
   
-  INTEGER, PARAMETER :: maxsend    = 34  ! Maximum number of fields to send
-  
+  INTEGER, PARAMETER :: maxsend    = 37  ! 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 
+
   ! Id for fields received from ocean
 
@@ -78 +82 @@
   INTEGER, PARAMETER :: idr_ocedms = 10
   INTEGER, PARAMETER :: idr_ocen2o = 11
-
-  INTEGER, PARAMETER :: maxrecv      = 11     ! Maximum number of fields to receive
+  INTEGER, PARAMETER :: idr_ocenh3 = 12
+
+  INTEGER, PARAMETER :: maxrecv      = 12     ! 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)
@@ -115 +120 @@
     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, nh3_cycle_cpl
     USE lmdz_xios  
     USE print_control_mod, ONLY: lunout
@@ -227 +232 @@
             infosend(ids_atmco2)%action = .TRUE. ; infosend(ids_atmco2)%name = 'COATMCO2'
         ENDIF
+        IF (n2o_cycle_cpl) THEN
+            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
+        IF (nh3_cycle_cpl) THEN
+            infosend(ids_atmnh3)%action = .TRUE. ; infosend(ids_atmnh3)%name = 'COATMNH3'
+        ENDIF
         infosend(ids_qraioc)%action = .TRUE. ; infosend(ids_qraioc)%name = 'COQRAIOC'
         infosend(ids_qsnooc)%action = .TRUE. ; infosend(ids_qsnooc)%name = 'COQSNOOC'
@@ -272 +286 @@
       inforecv(idr_ocen2o)%action = .TRUE. ; inforecv(idr_ocen2o)%name = 'SIN2OFLX'
    ENDIF
- 
+   IF (nh3_cycle_cpl) THEN
+      inforecv(idr_ocenh3)%action = .TRUE. ; inforecv(idr_ocenh3)%name = 'SINH3FLX'
+   ENDIF
 
 !************************************************************************************
@@ -358 +374 @@
     END DO
 
-
-IF (CPPKEY_INCA) THEN
-    IF (dms_cycle_cpl .OR. n2o_cycle_cpl) THEN
-       CALL init_inca_oasis(inforecv(idr_ocedms:idr_ocen2o))
-    ENDIF
-END IF
  
 !************************************************************************************
@@ -397 +407 @@
     ENDIF
 
+
+IF (CPPKEY_INCA) THEN
+    IF (dms_cycle_cpl .OR. n2o_cycle_cpl .OR. ndp_cycle_cpl .OR. nh3_cycle_cpl) THEN
+       CALL init_inca_oasis(inforecv(idr_ocedms:idr_ocenh3),infosend(ids_atmn2o:ids_atmnh3))
+    ENDIF
+END IF
+
 !$OMP END MASTER
     
@@ -451 +468 @@
              .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)
@@ -479 +496 @@
     INTEGER, INTENT(IN)                              :: ktime
     LOGICAL, INTENT(IN)                              :: last
-    REAL, DIMENSION(nbp_lon, jj_nb, maxsend), INTENT(IN) :: tab_put
+    REAL, DIMENSION(nbp_lon, jj_nb, maxsend_phys), INTENT(IN) :: tab_put
 
 ! Local variables
@@ -513 +530 @@
        IF (is_south_pole_dyn) wend=iend-nbp_lon+1
        
-       DO i = 1, maxsend
+       DO i = 1, maxsend_phys
           IF (infosend(i)%action) THEN
              field = RESHAPE(tab_put(:,:,i),(/nbp_lon*jj_nb/))
@@ -525 +542 @@
 !************************************************************************************
 
-    DO i = 1, maxsend
+    DO i = 1, maxsend_phys
       IF (infosend(i)%action .AND. infosend(i)%nid .NE. -1 ) THEN
           field = RESHAPE(tab_put(:,:,i),(/nbp_lon*jj_nb/))