Changeset 5654 for LMDZ6/trunk/libf/phylmd/oasis.F90

Timestamp:

May 15, 2025, 4:20:14 PM (3 weeks ago)

Author:

acozic

Message:

Add possibility to send N2O from atm chemistry (with Inca model - but definition for oasis is done in LMDZ) to NEMO

File:

• LMDZ6/trunk/libf/phylmd/oasis.F90 (modified) (7 diffs)

LMDZ6/trunk/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 :: maxsend    = 34  ! Maximum number of fields to send
-  
+  INTEGER, PARAMETER :: maxsend    = 35  ! 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
 
@@ -227 +229 @@
             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
         infosend(ids_qraioc)%action = .TRUE. ; infosend(ids_qraioc)%name = 'COQRAIOC'
         infosend(ids_qsnooc)%action = .TRUE. ; infosend(ids_qsnooc)%name = 'COQSNOOC'
@@ -358 +363 @@
     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
  
 !************************************************************************************
@@ -396 +395 @@
        WRITE(lunout,*) 'inicma : endef psmile ok '
     ENDIF
+
+
+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))
+    ENDIF
+END IF
 
 !$OMP END MASTER
@@ -479 +485 @@
     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 +519 @@
        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 +531 @@
 !************************************************************************************
 
-    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/))