diff --git a/source/calcul_Lstar_o.f b/source/calcul_Lstar_o.f index 99c11dc4..a91275c3 100644 --- a/source/calcul_Lstar_o.f +++ b/source/calcul_Lstar_o.f @@ -80,6 +80,7 @@ SUBROUTINE calcul_Lstar_opt(t_resol,r_resol, common /rconst/rad,pi C C + Ilflag = 0 Nder=Nder_def*r_resol Nreb=Nreb_def Ntet=Ntet_def*t_resol @@ -286,14 +287,10 @@ SUBROUTINE calcul_Lstar_opt(t_resol,r_resol, phi(I) = phi(I-1)+2.D0*pi/Nder Iflag_I = 0 c write(6,*)tet(I) - IF (Ilflag.EQ.0) THEN tetl = tet(I-1) IF (I.GT.2) tetl = 2.D0*tet(I-1)-tet(I-2) tet1 = tetl - ELSE - tetl = tet(I) - tet1 = tetl - ENDIF + c write(6,*)tetl c read(5,*) leI1 = baddata diff --git a/source/onera_desp_lib.f b/source/onera_desp_lib.f index 08a968d2..5f010325 100644 --- a/source/onera_desp_lib.f +++ b/source/onera_desp_lib.f @@ -64,7 +64,7 @@ SUBROUTINE make_lstar1(ntime,kext,options,sysaxes,iyearsat, c c Declare internal variables INTEGER*4 isat,iyear,kint,ifail - INTEGER*4 t_resol,r_resol,Ilflag,Ilflag_old + INTEGER*4 t_resol,r_resol REAL*8 mlon,mlon1 REAL*8 xGEO(3),xMAG(3),xSUN(3),rM,MLAT real*8 alti,lati,longi @@ -75,12 +75,9 @@ SUBROUTINE make_lstar1(ntime,kext,options,sysaxes,iyearsat, REAL*8 Lm(ntime_max),Lstar(ntime_max) C COMMON /magmod/k_ext,k_l,kint - COMMON /flag_L/Ilflag DATA xSUN /1.d0,0.d0,0.d0/ integer*4 int_field_select, ext_field_select C - Ilflag=0 - Ilflag_old=Ilflag if (options(3).lt.0 .or. options(3).gt.9) options(3)=0 t_resol=options(3)+1 r_resol=options(4)+1 @@ -132,12 +129,6 @@ SUBROUTINE make_lstar1(ntime,kext,options,sysaxes,iyearsat, c CALL calcul_Lstar_opt(t_resol,r_resol,XGeo & ,Lm(isat),Lstar(isat),XJ(isat),BLOCAL(isat),BMIN(isat)) - if (Ilflag_old .eq.1 .and. Lstar(isat).eq. Baddata) then - Ilflag=0 - CALL calcul_Lstar_opt(t_resol,r_resol,xGeo - & ,Lm(isat),Lstar(isat),XJ(isat),BLOCAL(isat),BMIN(isat)) - endif - Ilflag_old=Ilflag 99 continue