Changeset 5654 for LMDZ6/trunk

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

Location:

LMDZ6/trunk/libf/phylmd

Files:

• cpl_mod.F90 (modified) (3 diffs)
• oasis.F90 (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/cpl_mod.F90

@@ -363 +363 @@
        CALL histdef(nidct, 'tmp_lat','tmp_lat', &
             "-",nbp_lon,nbp_lat, nhoridct, 1, 1, 1, -99, 32, "inst", dtime,dtime)
-       DO jf=1,maxsend
+       DO jf=1,maxsend_phys
          IF (infosend(i)%action) THEN
              CALL histdef(nidct, infosend(i)%name ,infosend(i)%name , &
@@ -1280 +1280 @@
     REAL, DIMENSION(nbp_lon, jj_nb)                          :: tmp_calv
 ! Table with all fields to send to coupler
-    REAL, DIMENSION(nbp_lon, jj_nb, maxsend)                 :: tab_flds
+    REAL, DIMENSION(nbp_lon, jj_nb, maxsend_phys)            :: tab_flds
     REAL, DIMENSION(klon_mpi)                                :: rlon_mpi, rlat_mpi
     REAL  :: calving(nb_zone_calving)
@@ -1510 +1510 @@
 !*************************************************************************************
     IF (is_sequential) THEN
-        DO j=1,maxsend
+        DO j=1,maxsend_phys
           IF (infosend(j)%action) CALL histwrite(nidct,infosend(j)%name, itau_w, &
              tab_flds(:,:,j),nbp_lon*(nbp_lat),ndexct)

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/))