Sorry, a lot of errors in my last mail (the names in my code
are not P1 and P2) ...
You should read :
INTEGER, PARAMETER :: &
& SP=kind(1.0) &
& ,DP=kind(1.0D0)
! if precision(1.0) < 13 then P2_mask=0 --> P2=DP
! >=13 1 --> P2=SP
INTEGER, PARAMETER :: &
& SP_prec = precision(1.0_SP) &
& ,P2_prec = 13 &
& ,P2_mask = min( max( SP_prec-P2_prec+1, 0 ), 1 )
INTEGER, PARAMETER, PUBLIC :: &
& P1 = SP
& P2 = SP * P2_mask + DP * (1-P2_mask)
--
Pierre Hugonnet Seismic Data Processing R&D
COMPAGNIE GENERALE DE GEOPHYSIQUE
mail: 1 rue Leon Migaux phone.....(33) 164 47 45 59
91302 MASSY fax.......(33) 164 47 32 49
FRANCE [log in to unmask]
My opinions are not necessarily those of CGG
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|