Changeset 8531 for branches/ORCHIDEE_2_2/ORCHIDEE/src_driver

Timestamp:

2024-04-15T11:23:37+02:00 (3 months ago)

Author:

josefine.ghattas

Message:

Integrated [7912], [7925], [8179] and [8412] done on the trunk to output forcing file that can be used in offline mode. See ticket #899

File:

• branches/ORCHIDEE_2_2/ORCHIDEE/src_driver/forcing_tools.f90 (modified) (6 diffs)

branches/ORCHIDEE_2_2/ORCHIDEE/src_driver/forcing_tools.f90

@@ -4484 +4484 @@
   REAL(r_std), ALLOCATABLE, DIMENSION(:,:,:) :: tmp_slab2d
   CHARACTER(LEN=80) :: name
+  CHARACTER(LEN=80) :: coordinates ! temporary variable for attribut coordinates
   LOGICAL :: windzero
   !
@@ -4544 +4545 @@
      iret = NF90_GET_ATT(force_id(fileindex), varid, 'cell_methods', cellmethod)
      IF (iret /= NF90_NOERR) THEN
-        ! If the attribute is not found then we set a reasonable default : instantaneous and centered.
-        cellmethod="time: instantaneous"
+        ! If the attribute is not found then we set undefined
+        cellmethod="undefined"
      ENDIF
+
+     ! Read also attribute coordinates that might contain information about the name of the time axis.
+     ! If the forcing file is produced by XIOS, this is the case.
+     ! Reading of the time axis name can fail if the attribute coordinates contains something else than the time axis name.
+     iret = NF90_GET_ATT(force_id(fileindex), varid, 'coordinates', coordinates)
+     IF (iret == NF90_NOERR) THEN
+        ! The attribute 'coordinates' was found. This is probably containg the name of the time axis. 
+        ! We add this information in the begining of the celllmethod variable.
+        cellmethod=TRIM(coordinates)//":"//TRIM(cellmethod)
+     ENDIF
+
      !
      !
@@ -4645 +4657 @@
   ENDDO
   !
-  ! Go through all the time axes we have to find the right one.
+  ! Go through all the time axes we have to find the right one based on the information
+  ! read from attribute cell_methods or coordinates
   !
   timeindex=0
@@ -4682 +4695 @@
            ENDDO
            !
-           ! If there is no "(" then we have to find the centered axis.
         ELSE 
+           ! There is no "(" so we set the centered axis by defaut
            DO im=1,nbtmethods
               IF ( INDEX(time_cellmethod(itbase+im), "cent") > 0 ) THEN
@@ -4691 +4704 @@
         ENDIF
         !
-        ! The name of the time axis was found bu no method could be identified
+        ! The name of the time axis was found but no method could be identified
         !
         IF ( timeindex < 1 ) THEN
@@ -4702 +4715 @@
      ENDIF
   ENDDO
-  !
-  ! Should no corresponding time axis name be found, 
-  ! then we use the first centered one.
-  !
-  itax=1
-  DO WHILE ( timeindex < 1 ) 
-     IF ( INDEX(time_cellmethod(itax), "cent") > 0 ) THEN
-        timeindex = itax
-     ELSE
-        itax = itax + 1
-     ENDIF
-  ENDDO
-  !
+
+
+  !
+  ! If still no corresponding time axis name has be found, then we use the first centered one.
+  ! This is the case if none of the attributes cell_methods or coordinates were found or if they 
+  ! didn't include the name of the corresponding time axis.
+  IF ( timeindex < 1 ) THEN
+     CALL ipslerr(2,'forcing_attributetimeaxe',&
+          'No information about the name of the time axis was found in the attributes cell_method or coordinantes', &
+          'The first centered time axis will be used as defaut.', &
+          'Attribute coordinates is optional but must not contain something else than the name of the time axis.')
+     itax=1
+     DO WHILE ( timeindex < 1 ) 
+        IF ( INDEX(time_cellmethod(itax), "cent") > 0 ) THEN
+           timeindex = itax
+        ELSE
+           itax = itax + 1
+        ENDIF
+     ENDDO
+  END IF
+  
 END SUBROUTINE forcing_attributetimeaxe