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

Timestamp:

Jul 24, 2023, 4:50:57 PM (15 months ago)

Author:

acozic

Message:

add a flag for the dms cycle between ocean and atmosphere
remove cpp key CPLOCNINCA activated by option cplocninca
remove compilation option cplocninca

File:

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

LMDZ6/trunk/libf/phylmd/oasis.F90

@@ -24 +24 @@
   USE mod_prism_get_proto
   USE mod_prism_put_proto
-#endif
-#ifdef CPP_CPLOCNINCA
-  USE incaoasis, ONLY : inforcv
 #endif
 #endif
@@ -77 +74 @@
   INTEGER, PARAMETER :: idr_curenz = 7
   INTEGER, PARAMETER :: idr_oceco2 = 8
-
+  ! bulk salinity of the surface layer of the ocean, in ppt
   INTEGER, PARAMETER :: idr_sss = 9
-  ! bulk salinity of the surface layer of the ocean, in ppt
-
-  INTEGER, PARAMETER :: maxrecv    = 9  ! Maximum number of fields to receive
-  
-#ifdef CPP_CPLOCNINCA
-  INTEGER, PARAMETER :: idr_ocedms = 1
-  INTEGER, PARAMETER :: maxrcv = 1
-#endif
+  INTEGER, PARAMETER :: idr_ocedms = 10
+
+  INTEGER, PARAMETER :: maxrecv      = 10     ! 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)
+                                              ! and then send by routine in INCA model 
+  
 
   TYPE, PUBLIC ::   FLD_CPL            ! Type for coupling field information
@@ -118 +114 @@
     USE carbon_cycle_mod, ONLY : carbon_cycle_cpl
     USE wxios, ONLY : wxios_context_init 
+    USE chemistry_cycle_mod, ONLY : dms_cycle_cpl
     USE lmdz_xios  
     USE print_control_mod, ONLY: lunout
@@ -265 +262 @@
        inforecv(idr_oceco2)%action = .TRUE. ; inforecv(idr_oceco2)%name = 'SICO2FLX'
    ENDIF
-#ifdef CPP_CPLOCNINCA
-       inforcv(idr_ocedms)%action = .TRUE. ; inforcv(idr_ocedms)%name = 'SIDMSFLX'
-#endif
+   IF (dms_cycle_cpl) THEN 
+      inforecv(idr_ocedms)%action = .TRUE. ; inforecv(idr_ocedms)%name = 'SIDMSFLX'
+   ENDIF
+ 
 
 !************************************************************************************
@@ -353 +351 @@
     END DO
 
-! Now, if also coupling CPL with INCA, initialize here fields to be exchanged.
-#ifdef CPP_CPLOCNINCA
-    DO jf=1,maxrcv 
-       IF (inforcv(jf)%action) THEN
-          CALL prism_def_var_proto(inforcv(jf)%nid, inforcv(jf)%name, il_part_id, &
-               il_var_nodims, PRISM_In, il_var_actual_shape, il_var_type, &
-               ierror)
-          IF (ierror .NE. PRISM_Ok) THEN
-             WRITE(lunout,*) 'inicma : Problem with prism_def_var_proto for field : ',&
-                  inforcv(jf)%name
-             abort_message=' Problem in call to prism_def_var_proto for fields to receive'
-             CALL abort_physic(modname,abort_message,1)
-          ENDIF
-       ENDIF
-    END DO
+
+#ifdef INCA
+    IF (dms_cycle_cpl) THEN
+       CALL init_inca_oasis(inforecv(idr_ocedms))
+    ENDIF
 #endif
  
@@ -424 +412 @@
 ! Output arguments
 !************************************************************************************
-    REAL, DIMENSION(nbp_lon, jj_nb,maxrecv), INTENT(OUT) :: tab_get
+    REAL, DIMENSION(nbp_lon, jj_nb,maxrecv_phys), INTENT(OUT) :: tab_get
 
 ! Local variables
@@ -446 +434 @@
     ENDIF
     
-    DO i = 1, maxrecv
+    DO i = 1, maxrecv_phys
       IF (inforecv(i)%action .AND. inforecv(i)%nid .NE. -1) THEN
           field(:) = -99999.