Changeset 4247

Timestamp:

Sep 15, 2022, 9:54:51 AM (20 months ago)

Author:

tlurton

Message:

Add-ons in oasis.F90, bld.cfg and makelmdz_fcm to support coupling of species between PISCES and INCA in the IPSLESM/CO2AER configuration.

Location:

LMDZ6/trunk

Files:

• bld.cfg (modified) (1 diff)
• libf/phylmd/oasis.F90 (modified) (4 diffs)
• makelmdz_fcm (modified) (4 diffs)

LMDZ6/trunk/bld.cfg

@@ -100 +100 @@
 bld::excl_dep        use::ifile_attr
 bld::excl_dep        use::ixml_tree
+# The following for INCA used with CPL
+bld::excl_dep        use::incaoasis
 
 # Don't generate interface files

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
@@ -79 +82 @@
   INTEGER, PARAMETER :: maxrecv    = 9  ! Maximum number of fields to receive
   
+#ifdef CPP_CPLOCNINCA
+  INTEGER, PARAMETER :: idr_ocedms = 1
+  INTEGER, PARAMETER :: maxrcv = 1
+#endif
 
   TYPE, PUBLIC ::   FLD_CPL            ! Type for coupling field information
@@ -247 +254 @@
        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
 
 !************************************************************************************
@@ -331 +341 @@
        ENDIF
     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
+#endif
+ 
 !************************************************************************************
 ! Atmospheric Fields to send

LMDZ6/trunk/makelmdz_fcm

@@ -45 +45 @@
 isotrac=false
 force_compile=false
+cplocninca=false
 
 arch_defined="FALSE"
@@ -101 +102 @@
           true        : (obsolete; for backward compatibility) use ORCHIDEE tag 1.9-1.9.6
 [-chimie INCA/false]       : with INCA chemistry model or without (default: false)
+[-cplocninca ]        : allow transfer of species from PISCES to INCA (default: false)
 [-cosp true/false]    : compile with/without cosp package (default: false)
 [-cosp2 true/false]    : compile with/without cosp2 package (default: false)
@@ -179 +181 @@
           chimie="$2" ; shift ; shift ;;
 
+     "-cplocninca")
+          cplocninca=TRUE ; shift ;;
+
       "-isotopes")
           isotopes="$2" ; shift ; shift ;;
@@ -529 +534 @@
 fi
 
+if [[ "$cplocninca" == "TRUE" ]]
+then
+    CPP_KEY="${CPP_KEY} CPP_CPLOCNINCA"
+fi