Changeset 4076 for trunk

Timestamp:

Feb 18, 2026, 9:50:57 AM (7 days ago)

Author:

jbclement

Message:

PEM:
Fallback for 128-bit integers since this type is not supported by Intel Fortran.
JBC

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• changelog.txt (modified) (1 diff)
• display.F90 (modified) (3 diffs)
• numerics.F90 (modified) (1 diff)
• pem.F90 (modified) (2 diffs)

trunk/LMDZ.COMMON/libf/evolution/changelog.txt

@@ -883 +883 @@
 == 17/02/2026 == JBC
 Addition of runs ID for the current cycle in "pem_workflow.sts".
+
+== 18/02/2026 == JBC
+Fallback for 128-bit integers since this type is not supported by Intel Fortran.

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

@@ -91 +91 @@
 
 !=======================================================================
-SUBROUTINE print_end(n_yr_run,dur_secs,r_plnt2earth_yr,pem_ini_date,n_yr_sim)
+SUBROUTINE print_end(i_pem_run,n_yr_run,n_yr_sim,dur_secs,pem_ini_date,r_plnt2earth_yr)
 !-----------------------------------------------------------------------
 ! NAME
@@ -108 +108 @@
 ! DEPENDENCIES
 ! ------------
-use utility,         only: format_duration, int2str
-use workflow_status, only: i_pem_run
+use utility, only: format_duration, int2str
 
 ! DECLARATION
@@ -117 +116 @@
 ! ARGUMENTS
 ! ---------
-real(dp), intent(in) :: r_plnt2earth_yr, pem_ini_date, n_yr_sim, n_yr_run, dur_secs
+real(dp),    intent(in) :: n_yr_run, n_yr_sim, dur_secs, pem_ini_date, r_plnt2earth_yr
+integer(di), intent(in) :: i_pem_run
 
 ! CODE

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

@@ -28 +28 @@
 ! ----------
 ! Integers
-integer, parameter :: ti = int8                   ! Tiny integer      = 8 bits
-integer, parameter :: si = int16                  ! Short integer     = 16 bits
-integer, parameter :: di = int32                  ! Standard integer  = 32 bits (default)
-integer, parameter :: li = int64                  ! Long integer      = 64 bits
-integer, parameter :: lli = selected_int_kind(30) ! Long long integer = 128 bits
-integer, parameter :: wi = di                     ! Working integer
+integer, parameter :: lli_candidate = selected_int_kind(38)
+logical, parameter :: has_int128 = (lli_candidate > 0)
+integer, parameter :: ti = int8                                ! Tiny integer      = 8 bits
+integer, parameter :: si = int16                               ! Short integer     = 16 bits
+integer, parameter :: di = int32                               ! Standard integer  = 32 bits (default)
+integer, parameter :: li = int64                               ! Long integer      = 64 bits
+integer, parameter :: lli = merge(lli_candidate,li,has_int128) ! Long long integer = 128 bits (if available, otherwise fallback to int64)
+integer, parameter :: wi = di                                  ! Working integer
 
 ! Reals
-integer, parameter :: sp = real32  ! Simple precision    = 32  bits = 1-bit sign + 8-bits  exponent + 23-bits  significand
-integer, parameter :: dp = real64  ! Double precision    = 64  bits = 1-bit sign + 11-bits exponent + 52-bits  significand
-integer, parameter :: qp = real128 ! Quadruple precision = 128 bits = 1-bit sign + 15-bits exponent + 112-bits significand
+integer, parameter :: sp = real32  ! Simple precision    = 32  bits = 1-bit sign + 8-bit  exponent + 23-bit  significand
+integer, parameter :: dp = real64  ! Double precision    = 64  bits = 1-bit sign + 11-bit exponent + 52-bit  significand
+integer, parameter :: qp = real128 ! Quadruple precision = 128 bits = 1-bit sign + 15-bit exponent + 112-bit significand
 integer, parameter :: wp = dp      ! Working precision
 

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

@@ -50 +50 @@
 use tracers,            only: adapt_tracers2pressure, build4PCM_tracers, nq
 use utility,            only: real2str
-use workflow_status,    only: read_workflow_status, update_workflow_status
+use workflow_status,    only: i_pem_run, read_workflow_status, update_workflow_status
 use xios_data,          only: load_xios_data
 
@@ -533 +533 @@
 ! Footer
 call system_clock(c2)
-call print_end(n_yr_run,real((c2 - c1),dp)/real(cr,dp),r_plnt2earth_yr,pem_ini_date,n_yr_sim)
+call print_end(i_pem_run,n_yr_run,n_yr_sim,real((c2 - c1),dp)/real(cr,dp),pem_ini_date,r_plnt2earth_yr)
 
 END PROGRAM pem