Changeset 10214 for branches

Timestamp:

2018-10-23T16:41:23+02:00 (6 years ago)

Author:

frrh

Message:

Working version with OASIS3-MCT vn3.0 and 4.0

File:

• branches/UKMO/dev_r5518_GO6_runoff_O4test/NEMOGCM/NEMO/OPA_SRC/SBC/cpl_oasis3.F90 (modified) (3 diffs)

branches/UKMO/dev_r5518_GO6_runoff_O4test/NEMOGCM/NEMO/OPA_SRC/SBC/cpl_oasis3.F90

@@ -163 +163 @@
       INTEGER :: id_part_rnf_1d ! Partition for 1d river outflow fields
       INTEGER :: id_part_temp   ! Temperary partition used to choose either 0d or 1d partitions
-      INTEGER :: vector_length  ! Length of 0d or 1d variables (0d variables will have vector_length=1)
       INTEGER :: paral(5)       ! OASIS3 box partition
-      INTEGER :: ishape(4)      ! Shape of arrays passed to PSMILe. 
+      INTEGER :: ishape(4)      ! Shape of 2D arrays passed to PSMILe. 
                                 ! Redundant from OASIS3-MCT vn4.0 onwards but required
                                 ! to satisfy interface and for backward compatibility. 
+      INTEGER :: ishape0d1d(2)  ! Shape of 0D or 1D arrays passed to PSMILe. 
       INTEGER :: var_nodims(2)  ! Number of coupling field dimensions.
                                 ! var_nodims(1) is redundant from OASIS3-MCT vn4.0 onwards
@@ -209 +209 @@
       ishape(4) = nlej-nldj+1
 
+      ishape0d1d(1) = 0
+      ishape0d1d(2) = 0
 
       !
@@ -358 +360 @@
                     IF (nproc == 0) THEN                    
                        
+                       ishape0d1d(1) = 1
                        IF (srcv(ji)%dimensions == 0) THEN
                        
                           ! If 0D then set temporary variables to 0D components
                           id_part_temp = id_part_0d
-                          vector_length = 1
+                          ishape0d1d(2) = 1
                        ELSE
                        
                           ! If 1D then set temporary variables to river outflow components
                           id_part_temp = id_part_rnf_1d
-                          vector_length = nn_cpl_river
+                          ishape0d1d(2)= nn_cpl_river
                           
                        END IF
-                       
+
                        CALL oasis_def_var (srcv(ji)%nid(jc,jm), zclname, id_part_temp   , var_nodims,   &
-                                   OASIS_In           , (/ 1, vector_length /) , OASIS_REAL, nerror )
-                    ELSE
+                                   OASIS_In           , ishape0d1d(1:2) , OASIS_REAL, nerror )
+
+                    ELSE                       
+                       ishape0d1d(1) = 0
+                       ishape0d1d(2) = 0
                        ! Dummy call to keep OASIS3-MCT happy. 
                        CALL oasis_def_var (srcv(ji)%nid(jc,jm), zclname, id_part_0d   , var_nodims,   &
-                                   OASIS_In           , (/ 0, 0 /) , OASIS_REAL, nerror )
-                    END IF
+                                   OASIS_In           , ishape0d1d(1:2) , OASIS_REAL, nerror )
+                    
+                    ENDIF
                   ELSE 
                     ! It's a "normal" 2D (or pseudo 3D) coupling field.