Updated documentation to reflect the possibility to use variable length arrays

Author: AntoineBrunet

docs/frames/frames.html

@@ -23,8 +23,7 @@
 !               SUBROUTINE coord_trans1: Generic coordinate transformation from one Earth or Heliospheric coordinate
 !                                        system to another one
 !               SUBROUTINE coord_trans_vec1: Generic coordinate transformation from one Earth or Heliospheric coordinate
-!                                        system to another one (handle up to
-!                                        ntime_max positions)
+!                                        system to another one 
 !
 !***************************************************************************************************
 C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
@@ -2706,8 +2705,6 @@ SUBROUTINE coord_trans1(sysaxesIN,sysaxesOUT,iyr,idoy,
 C      REAL*8 xINV(3,nmax), xOUTV(3,nmax)
 C      INTEGER*4 numpoints
 C
-C     As with all (most?) onera library calls, the maximum array size
-C     is limited to ntime_max elements
 C                              Contributed by Timothy Guild, 2.2.07
 C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 C
@@ -2716,7 +2713,6 @@ SUBROUTINE coord_trans_vec1(ntime,sysaxesIN,sysaxesOUT,
      &   iyear,idoy,secs,xINV,xOUTV)
 
       IMPLICIT NONE
-      INCLUDE 'ntime_max.inc'
       INCLUDE 'variables.inc'
 
       INTEGER*4 nmax,i,ntime, sysaxesIN, sysaxesOUT

source/CoordTrans.f

@@ -20669,7 +20669,6 @@ SUBROUTINE EmpiricalLstar1(ntime,kext,options,iyearsat,idoy,
 c
       IMPLICIT NONE
       INCLUDE 'variables.inc'
-      INCLUDE 'ntime_max.inc'
 c
 c declare inputs
       INTEGER*4    kext,options(5)

source/LAndI2Lstar.f

@@ -30,7 +30,6 @@ SUBROUTINE GET_Bderivs(ntime,kext,options,sysaxes,dX,
 C        diffB(i,j,t) = dB_i/dx_j for point t (t=1 to ntime)
 
       IMPLICIT NONE
-      INCLUDE 'ntime_max.inc'   ! include file created by make, defines ntime_max
       INCLUDE 'variables.inc'
 C
       COMMON /magmod/k_ext,k_l,kint
@@ -133,7 +132,6 @@ SUBROUTINE compute_grad_curv_curl(ntime,Bgeo,Bmag,gradBmag,diffB,
 C     computes gradient factors, curvature factors, and curl of B
 
       IMPLICIT NONE
-      INCLUDE 'ntime_max.inc'   ! include file created by make, defines ntime_max
       INCLUDE 'variables.inc'
 
 C     all coordinates are in GEO reference frame

source/get_bderivs.f

@@ -94,7 +94,6 @@ SUBROUTINE GET_HEMI_MULTI(ntime,kext,options,sysaxes,iyearsat,
 c     calls get_hemi1 multiple times (ntime, <= ntime)
       IMPLICIT NONE
       INCLUDE 'variables.inc'
-      INCLUDE 'ntime_max.inc'   ! include file created by make, defines ntime_max
 
 c     declare inputs
       INTEGER*4    ntime,kext,options(5)

source/get_hemi.f

@@ -172,12 +172,12 @@ SUBROUTINE make_lstar_shell_splitting1(ntime,Nipa,kext,options,
         INTEGER*4    kext,k_ext,k_l,options(5),Nalp,Nipa
       PARAMETER (Nalp=25)
         INTEGER*4    ntime,sysaxes
-      INTEGER*4    iyearsat(ntime_max)
-      integer*4    idoy(ntime_max)
-      real*8     UT(ntime_max)
-      real*8     xIN1(ntime_max),xIN2(ntime_max),xIN3(ntime_max)
+      INTEGER*4    iyearsat(ntime)
+      integer*4    idoy(ntime)
+      real*8     UT(ntime)
+      real*8     xIN1(ntime),xIN2(ntime),xIN3(ntime)
       real*8     alpha(Nalp)
-      real*8     maginput(25,ntime_max)
+      real*8     maginput(25,ntime)
 c
 c
 c Declare internal variables
@@ -190,8 +190,8 @@ SUBROUTINE make_lstar_shell_splitting1(ntime,Nipa,kext,options,
 c
 c Declare output variables
         REAL*8     BLOCAL(ntime_max,Nalp),BMIN(ntime_max)
-      REAL*8     XJ(ntime_max,Nalp),MLT(ntime_max)
-        REAL*8     Lm(ntime_max,Nalp),Lstar(ntime_max,Nalp)
+      REAL*8     XJ(ntime_max,Nipa),MLT(ntime_max)
+        REAL*8     Lm(ntime_max,Nipa),Lstar(ntime_max,Nipa)
 C
       COMMON /magmod/k_ext,k_l,kint
         COMMON /flag_L/Ilflag
@@ -1561,11 +1561,11 @@ subroutine msis86(ntime,whichAp,DOY,UT,ALT,LAT,LONG,F107A,
      &F107,AP,Dens,Temp)
 c
       IMPLICIT NONE
-      INTEGER*4 ISW,ntime,whichAp,DOY(100000),I,J
+      INTEGER*4 ISW,ntime,whichAp,DOY(ntime),I,J
       REAL*8 SV(25),STL
-      REAL*8 UT(100000),ALT(100000),LAT(100000),LONG(100000)
-      REAL*8 F107A(100000),F107(100000),AP(7,100000)
-      REAL*8 Dens(8,100000),Temp(2,100000),APin(7),D(8),T(2)
+      REAL*8 UT(ntime),ALT(ntime),LAT(ntime),LONG(ntime)
+      REAL*8 F107A(ntime),F107(ntime),AP(7,ntime)
+      REAL*8 Dens(8,ntime),Temp(2,ntime),APin(7),D(8),T(2)
 c
       COMMON/CSWI/ISW
       DO I=1,25
@@ -1605,11 +1605,11 @@ subroutine msise90(ntime,whichAp,DOY,UT,ALT,LAT,LONG,F107A,
      &F107,AP,Dens,Temp)
 c
       IMPLICIT NONE
-      INTEGER*4 ISW,ntime,whichAp,DOY(100000),I,J
+      INTEGER*4 ISW,ntime,whichAp,DOY(ntime),I,J
       REAL*8 SV(25),STL
-      REAL*8 UT(100000),ALT(100000),LAT(100000),LONG(100000)
-      REAL*8 F107A(100000),F107(100000),AP(7,100000)
-      REAL*8 Dens(8,100000),Temp(2,100000),APin(7),D(8),T(2)
+      REAL*8 UT(ntime),ALT(ntime),LAT(ntime),LONG(ntime)
+      REAL*8 F107A(ntime),F107(ntime),AP(7,ntime)
+      REAL*8 Dens(8,ntime),Temp(2,ntime),APin(7),D(8),T(2)
 c
       COMMON/CSWI/ISW
       DO I=1,25
@@ -1639,11 +1639,11 @@ subroutine nrlmsise00(ntime,whichAp,DOY,UT,ALT,LAT,LONG,F107A,
      &F107,AP,Dens,Temp)
 c
       IMPLICIT NONE
-      INTEGER*4 ISW,ntime,whichAp,DOY(100000),I,J
+      INTEGER*4 ISW,ntime,whichAp,DOY(ntime),I,J
       REAL*8 SV(25),STL
-      REAL*8 UT(100000),ALT(100000),LAT(100000),LONG(100000)
-      REAL*8 F107A(100000),F107(100000),AP(7,100000)
-      REAL*8 Dens(9,100000),Temp(2,100000),APin(7),D(8),T(2)
+      REAL*8 UT(ntime),ALT(ntime),LAT(ntime),LONG(ntime)
+      REAL*8 F107A(ntime),F107(ntime),AP(7,ntime)
+      REAL*8 Dens(9,ntime),Temp(2,ntime),APin(7),D(8),T(2)
 c
       COMMON/CSWI/ISW
       DO I=1,25