Changeset 2841 for trunk/LMDZ.COMMON/libf/evolution

Timestamp:

Dec 7, 2022, 5:50:50 PM (2 years ago)

Author:

romain.vande

Message:

Mars PEM:
PEM is now adapted to run with XIOS diurnal averages (when they will work properly)
RV

File:

• trunk/LMDZ.COMMON/libf/evolution/nb_time_step_GCM.F90 (modified) (2 diffs)

trunk/LMDZ.COMMON/libf/evolution/nb_time_step_GCM.F90

@@ -38 +38 @@
       IF (ierr .NE. nf90_noerr) THEN
         write(*,*)"read_data_GCM: Le champ <temps> est absent"
-        write(*,*)"read_data_GCM: J essaie <Time>"
-        ierr = nf90_inq_varid (fID, "Time", vID)
+        write(*,*)"read_data_GCM: J essaie <time_counter>"
+        ierr = nf90_inq_varid (fID, "time_counter", vID)
         IF (ierr .NE. nf90_noerr) THEN
-           write(*,*)"read_data_GCM: Le champ <Time> est absent"
-           write(*,*)trim(nf90_strerror(ierr))
-           CALL ABORT_gcm("read_data_GCM", "", 1)
+          write(*,*)"read_data_GCM: Le champ <time_counter> est absent"
+          write(*,*)"read_data_GCM: J essaie <Time>"
+          IF (ierr .NE. nf90_noerr) THEN
+            write(*,*)"read_data_GCM: Le champ <Time> est absent"
+            write(*,*)trim(nf90_strerror(ierr))
+            CALL ABORT_gcm("read_data_GCM", "", 1)
+          ENDIF
+          ! Get the length of the "Time" dimension
+          ierr = nf90_inq_dimid(fID,"Time",vID)
+          ierr = nf90_inquire_dimension(fID,vID,len=timelen)          
         ENDIF
-        ! Get the length of the "Time" dimension
-        ierr = nf90_inq_dimid(fID,"Time",vID)
+        ! Get the length of the "time_counter" dimension
+        ierr = nf90_inq_dimid(fID,"time_counter",vID)
         ierr = nf90_inquire_dimension(fID,vID,len=timelen)
       ELSE   
@@ -59 +66 @@
 
   CONTAINS
-
 
 SUBROUTINE err(ierr,typ,nam)