KJ135
/
ISCE_INSAR
mirror of http://172.16.0.12:5000/WBJY/ISCE_INSAR.git
4
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

upgraded syntax of ALOS_pre_process/readOrbitPulse.f and ISSI/src/igrf2005_sub.f to compile with new gfortran (#435) 

* upgraded syntax of ALOS_pre_process/readOrbitPulse.f to compile with new gfortran

* updated subroutines in /ISSI/src/igrf2005_sub.f to be consistent
LT1AB
Eric Fielding 2022-02-15 16:07:32 -08:00 committed by GitHub
parent 1eaa62a459
commit 2717095323
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
2 changed files with 10 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
components/isceobj/Sensor/src/ALOS_pre_process/readOrbitPulse.f
View File

@ -144,7 +144,8 @@ c print *,val
end
integer*4 function unpackBytes(i1, i2, i3, i4)
integer*4 i1, i2, i3, i4
unpackBytes = iand(i1, 255)*256*256*256 + iand(i2, 255)*256*256 +
$ iand(i3, 255)*256 + iand(i4, 255)
integer*1 i1, i2, i3, i4, i255
i255 = 255 ! needed to keep gfortran >9 happy
unpackBytes = iand(i1, i255)*256*256*256 + iand(i2, i255)*256*256 +
$ iand(i3, i255)*256 + iand(i4, i255)
end function

8
contrib/issi/components/ISSI/src/igrf2005_sub.f
View File

@ -17,6 +17,7 @@ c----------------------------------------------------------------
INTEGER EGNR,AGNR,OGNR
REAL LATI,LONGI
COMMON/GENER/ UMR,ERA,AQUAD,BQUAD
character*80 PATH
C
CALL INITIZE
ibbb=0
@ -53,11 +54,12 @@ C 2005.00 04/25/05 CALL FELDI and DO 1111 I=1,7 (Alexey Petrov)
C 2005.01 11/10/05 added igrf_dip and geodip (MLAT)
C 2005.02 11/10/05 updated to IGRF-10 version
C 2006.00 12/21/06 GH2(120) -> GH2(144)
c 2022/01/23 added path to be consistent with igrf_bvector (Eric Fielding)
C
C*********************************************************************
subroutine igrf_sub(xlat,xlong,year,height,
& xl,icode,dip,dec)
& xl,icode,dip,dec,path)
c----------------------------------------------------------------
c INPUT:
c xlat geodatic latitude in degrees
@ -65,6 +67,7 @@ c xlong geodatic longitude in degrees
c year decimal year (year+month/12.0-0.5 or year+day-of-year/365
c or 366 if leap year)
c height height in km
c path the path to the data files
c OUTPUT:
c xl L value
c icode =1 L is correct; =2 L is not correct;
@ -76,6 +79,7 @@ c----------------------------------------------------------------
INTEGER EGNR,AGNR,OGNR
REAL LATI,LONGI
COMMON/GENER/ UMR,ERA,AQUAD,BQUAD
character*80 PATH
C
CALL INITIZE
ibbb=0
@ -86,7 +90,7 @@ C
c
C----------------CALCULATE PROFILES-----------------------------------
c
CALL FELDCOF(YEAR,DIMO)
CALL FELDCOF(YEAR,DIMO,PATH)
CALL FELDG(LATI,LONGI,HEIGHT,BNORTH,BEAST,BDOWN,BABS)
CALL SHELLG(LATI,LONGI,HEIGHT,DIMO,XL,ICODE,BAB1)
DIP=ASIN(BDOWN/BABS)/UMR