Skip to content
Snippets Groups Projects
Commit be21503a authored by Thomas Feuvrier's avatar Thomas Feuvrier
Browse files

Création de la fonction main de 6S et création du main qui appelle cette fonction.

parent dd306d5c
No related branches found
No related tags found
No related merge requests found
Expand all Hide whitespace changes
Inline Side-by-side
Showing
with 8434 additions and 1 deletion
Code/Visu/CMakeLists.txt
+ 3
− 0
View file @ be21503a
......@@ -5,6 +5,9 @@ FILE(GLOB OTBVisu_SRCS "*.cxx" )
ADD_LIBRARY(OTBVisu ${OTBVisu_SRCS})
TARGET_LINK_LIBRARIES (OTBVisu OTBGui ITKCommon OTBCommon ${FLTK_LIBRARIES} ${OPENGL_LIBRARIES})
MESSAGE("Warning: FLTK_LIBRARIES : ${FLTK_LIBRARIES} et OPENGL_LIBRARIES : ${OPENGL_LIBRARIES}")
INSTALL(TARGETS OTBVisu
RUNTIME DESTINATION ${OTB_INSTALL_BIN_DIR} COMPONENT RuntimeLibraries
LIBRARY DESTINATION ${OTB_INSTALL_LIB_DIR} COMPONENT RuntimeLibraries
......
Testing/Utilities/CMakeLists.txt
+ 1
− 0
View file @ be21503a
......@@ -104,6 +104,7 @@ svmGenericKernelFunctor.cxx
svmTest.cxx
svmGenericKernelTest.cxx
6SMainTest.cxx
6SMainOtbWithCallFunctionTest.cxx
)
INCLUDE_DIRECTORIES("${OTBTesting_BINARY_DIR}")
......
Testing/Utilities/otbUtilitiesTests.cxx
+ 1
− 0
View file @ be21503a
......@@ -34,4 +34,5 @@ REGISTER_TEST(svmGenericKernelFunctor);
REGISTER_TEST(svmTest);
REGISTER_TEST(svmGenericKernelTest);
REGISTER_TEST(SIXSMainTest);
REGISTER_TEST(SIXSMainOtbWithCallFunctionTest);
}
Utilities/otb6S/CMakeLists.txt
+ 8
− 1
View file @ be21503a
......@@ -49,6 +49,13 @@ LIST(REMOVE_ITEM LIBF2C_SRCS "${OTB6S_SOURCE_DIR}/libf2c/main.c" )
LIST(REMOVE_ITEM LIBF2C_SRCS "${OTB6S_SOURCE_DIR}/libf2c/getarg_.c" )
LIST(REMOVE_ITEM LIBF2C_SRCS "${OTB6S_SOURCE_DIR}/libf2c/iargc_.c" )
# Main code source
SET(OTB6S_MAIN_SRCS
otb_main6S_function.c # otb 6S main function
main.c # ORIGINAL 6S main
otb_main6S_program_call_function.c # OTB 6S main (which call otb_main6S_function function)
)
INCLUDE_DIRECTORIES(${OTB6S_SOURCE_DIR}/libf2c)
# Adjust the compiler flags to avoid problems with f2c-generated code.
......@@ -75,7 +82,7 @@ INSTALL(FILES ${__files1}
COMPONENT Development)
IF(NOT OTB_DISABLE_UTILITIES_COMPILATION)
ADD_LIBRARY(otb6S ${LIBF2C_SRCS} ${OTB6S_SRCS} )
ADD_LIBRARY(otb6S ${LIBF2C_SRCS} ${OTB6S_SRCS} ${OTB6S_MAIN_SRCS})
IF(UNIX)
TARGET_LINK_LIBRARIES (otb6S m)
......
Utilities/otb6S/main_6s.h
+ 1
− 0
View file @ be21503a
......@@ -6,6 +6,7 @@ extern "C" {
#endif
int ssssss_ ();
int call6s_();
#ifdef __cplusplus
}
......
Utilities/otb6S/otb_6S_mangle.h
+ 2
− 0
View file @ be21503a
......@@ -13,6 +13,8 @@ nm libotb_6s.a |grep " [TR] "
/* Functions */
#define varsol_ otb_6s_varsol_
/* otb 6S Main function */
#define ssssss_otb_main_function otb_6s_ssssss_otb_main_function
#define adaptquad_ otb_6s_adaptquad_
#define balanc_ otb_6s_balanc_
......
Utilities/otb6S/otb_6S_prototypes.h
+ 3
− 0
View file @ be21503a
......@@ -114,3 +114,6 @@
#include "6SV1.0B/WAVA4.h"
#include "6SV1.0B/WAVA5.h"
#include "6SV1.0B/WAVA6.h"
/* otb 6S Main function */
#include "otb_main6S_function.h"
Utilities/otb6S/otb_6S_unmangle.h
+ 2
− 0
View file @ be21503a
......@@ -5,6 +5,8 @@
/* Functions */
#undef varsol_
/* otb 6S Main function */
#undef ssssss_otb_main_function
#undef adaptquad_
#undef balanc_
......
Utilities/otb6S/otb_main6S_function.P 0 → 100644
+ 26
− 0
View file @ be21503a
extern int otb_6s__(real *otb_asol__, real *otb_phi0__, real *otb_avis__, real *otb_phiv__, integer *otb_month__, integer *otb_jday__, real *otb_pressure__, real *otb_uw__, real *otb_uo3__, integer *otb_iaer__, real *otb_taer55__, real *otb_wlinf__, real *otb_wlsup__, real *otb_s__, real *otb_ratm__, real *otb_sast__, real *otb_tgasm__, real *otb_sdtott__, real *otb_sutott__);
/* comlen sixs_ier__ 8 */
/* comlen mie_in__ 800 */
/* comlen sixs_planesim__ 680 */
/* comlen sixs_test__ 4 */
/* comlen sixs_sos__ 3984 */
/* comlen sixs_trunc__ 656 */
/* comlen sixs_ffu__ 6012 */
/* comlen sixs_del__ 8 */
/* comlen sixs_atm__ 680 */
/* comlen sixs_aer__ 160 */
/* comlen sixs_disc__ 840 */
/*:ref: gauss_ 14 5 6 6 6 6 4 */
/*:ref: varsol_ 14 3 4 4 6 */
/*:ref: us62_ 14 0 */
/*:ref: aeroso_ 14 6 4 6 6 6 13 124 */
/*:ref: pressure_ 14 3 6 6 6 */
/*:ref: presplane_ 14 4 6 6 6 6 */
/*:ref: equivwl_ 14 4 4 4 6 6 */
/*:ref: discom_ 14 18 4 4 6 6 6 6 6 6 6 4 4 4 6 6 6 6 6 6 */
/*:ref: specinterp_ 14 7 6 6 6 6 6 6 6 */
/*:ref: odrayl_ 14 2 6 6 */
/*:ref: solirr_ 14 2 6 6 */
/*:ref: abstra_ 14 34 4 6 6 6 6 6 6 6 4 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 */
/*:ref: interp_ 14 25 4 4 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 */
/*:ref: enviro_ 14 8 6 6 6 6 6 6 6 6 */
Utilities/otb6S/otb_main6S_function.c 0 → 100644
+ 5028
− 0
View file @ be21503a
This diff is collapsed.
Utilities/otb6S/otb_main6S_function.f 0 → 100755
+ 2698
− 0
View file @ be21503a
This diff is collapsed.
Utilities/otb6S/otb_main6S_function.h 0 → 100644
+ 20
− 0
View file @ be21503a
extern int otb_6s_ssssss_otb_main_function( otb_6s_real *otb_asol__,
otb_6s_real *otb_phi0__,
otb_6s_real *otb_avis__,
otb_6s_real *otb_phiv__,
otb_6s_integer *otb_month__,
otb_6s_integer *otb_jday__,
otb_6s_real *otb_pressure__,
otb_6s_real *otb_uw__,
otb_6s_real *otb_uo3__,
otb_6s_integer *otb_iaer__,
otb_6s_real *otb_taer55__,
otb_6s_real *otb_wlinf__,
otb_6s_real *otb_wlsup__,
otb_6s_real *otb_s__,
otb_6s_real *otb_ratm__,
otb_6s_real *otb_sast__,
otb_6s_real *otb_tgasm__,
otb_6s_real *otb_sdtott__,
otb_6s_real *otb_sutott__
);
Utilities/otb6S/otb_main6S_program_call_function.P 0 → 100644
+ 1
− 0
View file @ be21503a
/*:ref: otb_6s__ 14 19 6 6 6 6 4 4 6 6 6 4 6 6 6 6 6 6 6 6 6 */
Utilities/otb6S/otb_main6S_program_call_function.c 0 → 100644
+ 457
− 0
View file @ be21503a
/* otb_main6S_program_call_function.f -- translated by f2c (version 19970805).
You must link the resulting object file with the libraries:
-lf2c -lm (in that order)
*/
#ifdef __cplusplus
extern "C" {
#endif
#include "otb_6S.h"
/* Table of constant values */
static integer c__9 = 9;
static integer c__1 = 1;
static integer c__4 = 4;
static integer c__3 = 3;
/*< program call6S >*/
/* Main program */ OTB_MAIN__()
{
/* Format strings */
static char fmt_100[] = "(10x,\002 -----> atmospheric reflectance : \
\002,f6.5)";
static char fmt_200[] = "(10x,\002 -----> atmospheric spherical albedo \
: \002,f6.5)";
static char fmt_300[] = "(10x,\002 -----> total gaseous transmission : \
\002,f6.5)";
static char fmt_400[] = "(10x,\002 -----> downward transmittance : \
\002,f6.5)";
static char fmt_500[] = "(10x,\002 -----> upward transmittance : \
\002,f6.5)";
/* System generated locals */
integer i__1;
/* Builtin functions */
integer s_wsle(cilist *), do_lio(integer *, integer *, char *, ftnlen),
e_wsle(), s_rsle(cilist *), e_rsle(), s_wsfe(cilist *), do_fio(
integer *, char *, ftnlen), e_wsfe();
/* Local variables */
real otb_ratm__, pressure;
integer i__, l;
real s[1501];
integer ik;
real uw, uo3, phi0;
integer iaer, iinf, jday;
real asol, avis, phiv, sast;
integer isup, iread;
real taer55, tgasm, wlinf;
integer month;
real wlsup;
extern /* Subroutine */ int otb_6s__(real *, real *, real *, real *,
integer *, integer *, real *, real *, real *, integer *, real *,
real *, real *, real *, real *, real *, real *, real *, real *);
real sdtott, sutott;
/* Fortran I/O blocks */
static cilist io___4 = { 0, 6, 0, 0, 0 };
static cilist io___5 = { 0, 0, 0, 0, 0 };
static cilist io___12 = { 0, 6, 0, 0, 0 };
static cilist io___13 = { 0, 6, 0, 0, 0 };
static cilist io___14 = { 0, 6, 0, 0, 0 };
static cilist io___15 = { 0, 0, 0, 0, 0 };
static cilist io___19 = { 0, 6, 0, 0, 0 };
static cilist io___20 = { 0, 6, 0, 0, 0 };
static cilist io___21 = { 0, 6, 0, 0, 0 };
static cilist io___22 = { 0, 0, 0, 0, 0 };
static cilist io___24 = { 0, 6, 0, 0, 0 };
static cilist io___25 = { 0, 6, 0, 0, 0 };
static cilist io___26 = { 0, 6, 0, 0, 0 };
static cilist io___27 = { 0, 6, 0, 0, 0 };
static cilist io___28 = { 0, 0, 0, 0, 0 };
static cilist io___30 = { 0, 6, 0, 0, 0 };
static cilist io___31 = { 0, 6, 0, 0, 0 };
static cilist io___34 = { 0, 0, 0, 0, 0 };
static cilist io___37 = { 0, 6, 0, 0, 0 };
static cilist io___38 = { 0, 6, 0, 0, 0 };
static cilist io___40 = { 0, 0, 0, 0, 0 };
static cilist io___47 = { 0, 6, 0, 0, 0 };
static cilist io___48 = { 0, 6, 0, fmt_100, 0 };
static cilist io___49 = { 0, 6, 0, fmt_200, 0 };
static cilist io___50 = { 0, 6, 0, fmt_300, 0 };
static cilist io___51 = { 0, 6, 0, fmt_400, 0 };
static cilist io___52 = { 0, 6, 0, fmt_500, 0 };
/*< IMPLICIT NONE >*/
/*< real wlinf,wlsup,s(1501) >*/
/*< real asol,phi0,avis,phiv >*/
/*< integer month,jday >*/
/*< real uw,uo3,taer55 >*/
/*< real tgasm,sdtott,sutott,sast >*/
/*< integer iwave,iinf,isup,inhomo,idirec,iaer >*/
/*< integer iread,i,ik,l >*/
/* _otb Atmospheric reflectance */
/*< real otb_ratm !Atmospheric reflectance >*/
/*< real pressure !Atmospheric pressure >*/
/* ***********************************************************************
*/
/* Parameters initialization */
/* ***********************************************************************
*/
/*< iread=5 >*/
iread = 5;
/*< iinf=1 >*/
iinf = 1;
/*< isup=1501 >*/
isup = 1501;
/*< write(6,*) "------> Execution de l'interface d'appel de 6S " >*/
s_wsle(&io___4);
do_lio(&c__9, &c__1, "------> Execution de l'interface d'appel de 6S ",
47L);
e_wsle();
/* **********************************************************************c
*/
/* c
*/
/* * sun c
*/
/* \ * / c
*/
/* * * * * * c
*/
/* z / * \ c
*/
/* + /+ c
*/
/* satellite / + / c
*/
/* o/ + / c
*/
/* /.\ + /. c
*/
/* / . \ _avis-_+_-asol_/ . c
*/
/* . \- -+ / . north c
*/
/* . \ + / . + c
*/
/* . \ + / .+ c
*/
/* . \ + / +. c
*/
/* . \ + / + . c
*/
/* . \ + / + . c
*/
/* . \ +/ + . c
*/
/* west + + + + + + + . + + + + +\+ + + + + . + + + + + + + + east c
*/
/* . +.. . c
*/
/* . + . . . c
*/
/* . + . . . c
*/
/* . + . .'. . c
*/
/* . + .. . , ' .. c
*/
/* .+ . \ . c
*/
/* +. . \ . c
*/
/* + . . \ . c
*/
/* south . . (phiv-phi0) c
*/
/* c
*/
/* c
*/
/* c
*/
/* **********************************************************************c
*/
/* **********************************************************************c
*/
/* igeom geometrical conditions c
*/
/* -------------------------------------- c
*/
/* c
*/
/* you choose your own conditions; igeom=0 c
*/
/* 0 enter solar zenith angle (in degrees ) c
*/
/* solar azimuth angle " c
*/
/* satellite zenith angle " c
*/
/* satellite azimuth angle " c
*/
/* month c
*/
/* day of the month c
*/
/* c
*/
/* **********************************************************************c
*/
/*< read(iread,*) asol,phi0,avis,phiv,month,jday >*/
io___5.ciunit = iread;
s_rsle(&io___5);
do_lio(&c__4, &c__1, (char *)&asol, (ftnlen)sizeof(real));
do_lio(&c__4, &c__1, (char *)&phi0, (ftnlen)sizeof(real));
do_lio(&c__4, &c__1, (char *)&avis, (ftnlen)sizeof(real));
do_lio(&c__4, &c__1, (char *)&phiv, (ftnlen)sizeof(real));
do_lio(&c__3, &c__1, (char *)&month, (ftnlen)sizeof(integer));
do_lio(&c__3, &c__1, (char *)&jday, (ftnlen)sizeof(integer));
e_rsle();
/*< write(6,*) "asol,phi0,avis,phiv,month,jday : " >*/
s_wsle(&io___12);
do_lio(&c__9, &c__1, "asol,phi0,avis,phiv,month,jday : ", 33L);
e_wsle();
/*< write(6,*) asol,phi0,avis,phiv,month,jday >*/
s_wsle(&io___13);
do_lio(&c__4, &c__1, (char *)&asol, (ftnlen)sizeof(real));
do_lio(&c__4, &c__1, (char *)&phi0, (ftnlen)sizeof(real));
do_lio(&c__4, &c__1, (char *)&avis, (ftnlen)sizeof(real));
do_lio(&c__4, &c__1, (char *)&phiv, (ftnlen)sizeof(real));
do_lio(&c__3, &c__1, (char *)&month, (ftnlen)sizeof(integer));
do_lio(&c__3, &c__1, (char *)&jday, (ftnlen)sizeof(integer));
e_wsle();
/*< write(6,*) >*/
s_wsle(&io___14);
e_wsle();
/* **********************************************************************c
*/
/* idatm atmospheric model c
*/
/* -------------------- c
*/
/* pressure ( in mb ) c
*/
/* uw (in g/cm2 ) c
*/
/* uo3 (in cm-atm) c
*/
/* profil is taken from us62 c
*/
/* c
*/
/* **********************************************************************c
*/
/*< read(iread,*) pressure, uw, uo3 >*/
io___15.ciunit = iread;
s_rsle(&io___15);
do_lio(&c__4, &c__1, (char *)&pressure, (ftnlen)sizeof(real));
do_lio(&c__4, &c__1, (char *)&uw, (ftnlen)sizeof(real));
do_lio(&c__4, &c__1, (char *)&uo3, (ftnlen)sizeof(real));
e_rsle();
/*< write(6,*) "pressure, uw, uo3" >*/
s_wsle(&io___19);
do_lio(&c__9, &c__1, "pressure, uw, uo3", 17L);
e_wsle();
/*< write(6,*) pressure, uw, uo3 >*/
s_wsle(&io___20);
do_lio(&c__4, &c__1, (char *)&pressure, (ftnlen)sizeof(real));
do_lio(&c__4, &c__1, (char *)&uw, (ftnlen)sizeof(real));
do_lio(&c__4, &c__1, (char *)&uo3, (ftnlen)sizeof(real));
e_wsle();
/*< write(6,*) >*/
s_wsle(&io___21);
e_wsle();
/* **********************************************************************c
*/
/* c
*/
/* iaer aerosol model(type) c
*/
/* -------------- c
*/
/* c
*/
/* c
*/
/* you select one of the following standard aerosol models: c
*/
/* 0 no aerosols c
*/
/* 1 continental model ) c
*/
/* 2 maritime model ) according to sra models c
*/
/* 3 urban model ) c
*/
/* 5 shettle model for background desert aerosol c
*/
/* c
*/
/* **********************************************************************c
*/
/*< read(iread,*) iaer >*/
io___22.ciunit = iread;
s_rsle(&io___22);
do_lio(&c__3, &c__1, (char *)&iaer, (ftnlen)sizeof(integer));
e_rsle();
/*< write(6,*) "iaer = ", iaer >*/
s_wsle(&io___24);
do_lio(&c__9, &c__1, "iaer = ", 7L);
do_lio(&c__3, &c__1, (char *)&iaer, (ftnlen)sizeof(integer));
e_wsle();
/*< write(6,*) >*/
s_wsle(&io___25);
e_wsle();
/*< >*/
if (iaer != 0 && iaer != 1 && iaer != 2 && iaer != 3 && iaer != 5) {
/*< write(6,*) "Choix iaer non compatible : ", iaer >*/
s_wsle(&io___26);
do_lio(&c__9, &c__1, "Choix iaer non compatible : ", 28L);
do_lio(&c__3, &c__1, (char *)&iaer, (ftnlen)sizeof(integer));
e_wsle();
/*< write(6,*) "Choisir 0, 1, 2, 3 ou 5 uniquement" >*/
s_wsle(&io___27);
do_lio(&c__9, &c__1, "Choisir 0, 1, 2, 3 ou 5 uniquement", 34L);
e_wsle();
/*< goto 999 >*/
goto L999;
/*< endif >*/
}
/*< read(iread,*) taer55 >*/
io___28.ciunit = iread;
s_rsle(&io___28);
do_lio(&c__4, &c__1, (char *)&taer55, (ftnlen)sizeof(real));
e_rsle();
/*< write(6,*) "taer55 = ", taer55 >*/
s_wsle(&io___30);
do_lio(&c__9, &c__1, "taer55 = ", 9L);
do_lio(&c__4, &c__1, (char *)&taer55, (ftnlen)sizeof(real));
e_wsle();
/*< write(6,*) >*/
s_wsle(&io___31);
e_wsle();
/* **********************************************************************c
*/
/* iwave input of the spectral conditions c
*/
/* -------------------------------- c
*/
/* c
*/
/* 1 enter wlinf, wlsup and user's filter function s(lambda) c
*/
/* ( by step of 0.0025 micrometer). c
*/
/* c
*/
/* **********************************************************************c
*/
/*< do 38 l=iinf,isup >*/
i__1 = isup;
for (l = iinf; l <= i__1; ++l) {
/*< s(l)=0. >*/
s[l - 1] = (float)0.;
/*< 38 continue >*/
/* L38: */
}
/*< 110 read(iread,*) wlinf,wlsup >*/
/* L110: */
io___34.ciunit = iread;
s_rsle(&io___34);
do_lio(&c__4, &c__1, (char *)&wlinf, (ftnlen)sizeof(real));
do_lio(&c__4, &c__1, (char *)&wlsup, (ftnlen)sizeof(real));
e_rsle();
/*< write(6,*) " wlinf,wlsup : ", wlinf,wlsup >*/
s_wsle(&io___37);
do_lio(&c__9, &c__1, " wlinf,wlsup : ", 15L);
do_lio(&c__4, &c__1, (char *)&wlinf, (ftnlen)sizeof(real));
do_lio(&c__4, &c__1, (char *)&wlsup, (ftnlen)sizeof(real));
e_wsle();
/*< write(6,*) >*/
s_wsle(&io___38);
e_wsle();
/*< iinf=(wlinf-.25)/0.0025+1.5 >*/
iinf = (wlinf - (float).25) / (float).0025 + (float)1.5;
/*< isup=(wlsup-.25)/0.0025+1.5 >*/
isup = (wlsup - (float).25) / (float).0025 + (float)1.5;
/*< do 1113 ik=iinf,isup >*/
i__1 = isup;
for (ik = iinf; ik <= i__1; ++ik) {
/*< s(ik)=0. >*/
s[ik - 1] = (float)0.;
/*< 1113 continue >*/
/* L1113: */
}
/*< read(iread,*) (s(i),i=iinf,isup) >*/
io___40.ciunit = iread;
s_rsle(&io___40);
i__1 = isup;
for (i__ = iinf; i__ <= i__1; ++i__) {
do_lio(&c__4, &c__1, (char *)&s[i__ - 1], (ftnlen)sizeof(real));
}
e_rsle();
/* **********************************************************************c
*/
/* c
*/
/* c
*/
/* start of computations c
*/
/* c
*/
/* c
*/
/* c
*/
/* **********************************************************************c
*/
/* goto 800 */
/*< >*/
otb_6s_ssssss_otb_main_function(&asol, &phi0, &avis, &phiv, &month, &jday, &pressure, &uw, &uo3,
&iaer, &taer55, &wlinf, &wlsup, s, &otb_ratm__, &sast, &tgasm, &
sdtott, &sutott);
/*< write(6,*) >*/
s_wsle(&io___47);
e_wsle();
/*< 800 write(6,100) otb_ratm >*/
/* L800: */
s_wsfe(&io___48);
do_fio(&c__1, (char *)&otb_ratm__, (ftnlen)sizeof(real));
e_wsfe();
/*< write(6,200) sast >*/
s_wsfe(&io___49);
do_fio(&c__1, (char *)&sast, (ftnlen)sizeof(real));
e_wsfe();
/*< write(6,300) tgasm >*/
s_wsfe(&io___50);
do_fio(&c__1, (char *)&tgasm, (ftnlen)sizeof(real));
e_wsfe();
/*< write(6,400) sdtott >*/
s_wsfe(&io___51);
do_fio(&c__1, (char *)&sdtott, (ftnlen)sizeof(real));
e_wsfe();
/*< write(6,500) sutott >*/
s_wsfe(&io___52);
do_fio(&c__1, (char *)&sutott, (ftnlen)sizeof(real));
e_wsfe();
/*< 100 format(10x,40h -----> atmospheric reflectance : , f6.5) >*/
/*< 200 format(10x,40h -----> atmospheric spherical albedo : , f6.5) >*/
/*< 300 format(10x,40h -----> total gaseous transmission : , f6.5) >*/
/*< 400 format(10x,40h -----> downward transmittance : , f6.5) >*/
/*< 500 format(10x,40h -----> upward transmittance : , f6.5) >*/
/*< 999 continue >*/
L999:
/*< end >*/
return 0;
} /* OTB_MAIN__ */
/* Main program alias */ int call6s_ () { OTB_MAIN__ (); return 0; }
#ifdef __cplusplus
}
#endif
Utilities/otb6S/otb_main6S_program_call_function.f 0 → 100755
+ 183
− 0
View file @ be21503a
program call6S
IMPLICIT NONE
real wlinf,wlsup,s(1501)
real asol,phi0,avis,phiv
integer month,jday
real uw,uo3,taer55
real tgasm,sdtott,sutott,sast
integer iwave,iinf,isup,inhomo,idirec,iaer
integer iread,i,ik,l
C_otb Atmospheric reflectance
real otb_ratm !Atmospheric reflectance
real pressure !Atmospheric pressure
c***********************************************************************
c Parameters initialization
c***********************************************************************
iread=5
iinf=1
isup=1501
write(6,*) "------> Execution de l'interface d'appel de 6S "
c**********************************************************************c
c c
c * sun c
c \ * / c
c * * * * * c
c z / * \ c
c + /* c
c satellite / + / c
c o/ + / c
c /.\ + /. c
c / . \ _avis-_+_-asol_/ . c
c . \- -+ / . north c
c . \ + / . + c
c . \ + / .+ c
c . \ + / +. c
c . \ + / + . c
c . \ + / + . c
c . \ +/ + . c
c west + + + + + + + . + + + + +\+ + + + + . + + + + + + + + east c
c . +.. . c
c . + . . . c
c . + . . . c
c . + . .'. . c
c . + .. . , ' .. c
c .+ . \ . c
c +. . \ . c
c + . . \ . c
c south . . (phiv-phi0) c
c c
c c
c c
c**********************************************************************c
c**********************************************************************c
c igeom geometrical conditions c
c -------------------------------------- c
c c
c you choose your own conditions; igeom=0 c
c 0 enter solar zenith angle (in degrees ) c
c solar azimuth angle " c
c satellite zenith angle " c
c satellite azimuth angle " c
c month c
c day of the month c
c c
c**********************************************************************c
read(iread,*) asol,phi0,avis,phiv,month,jday
write(6,*) "asol,phi0,avis,phiv,month,jday : "
write(6,*) asol,phi0,avis,phiv,month,jday
write(6,*)
c**********************************************************************c
c idatm atmospheric model c
c -------------------- c
c pressure ( in mb ) c
c uw (in g/cm2 ) c
c uo3 (in cm-atm) c
c profil is taken from us62 c
c c
c**********************************************************************c
read(iread,*) pressure, uw, uo3
write(6,*) "pressure, uw, uo3"
write(6,*) pressure, uw, uo3
write(6,*)
c**********************************************************************c
c c
c iaer aerosol model(type) c
c -------------- c
c c
c c
c you select one of the following standard aerosol models: c
c 0 no aerosols c
c 1 continental model ) c
c 2 maritime model ) according to sra models c
c 3 urban model ) c
c 5 shettle model for background desert aerosol c
c c
c**********************************************************************c
read(iread,*) iaer
write(6,*) "iaer = ", iaer
write(6,*)
if ((iaer.ne.0).and.(iaer.ne.1).and.(iaer.ne.2).and.
c (iaer.ne.3).and.(iaer.ne.5)) then
write(6,*) "Choix iaer non compatible : ", iaer
write(6,*) "Choisir 0, 1, 2, 3 ou 5 uniquement"
goto 999
endif
read(iread,*) taer55
write(6,*) "taer55 = ", taer55
write(6,*)
c**********************************************************************c
c iwave input of the spectral conditions c
c -------------------------------- c
c c
c 1 enter wlinf, wlsup and user's filter function s(lambda) c
c ( by step of 0.0025 micrometer). c
c c
c**********************************************************************c
do 38 l=iinf,isup
s(l)=0.
38 continue
110 read(iread,*) wlinf,wlsup
write(6,*) " wlinf,wlsup : ", wlinf,wlsup
write(6,*)
iinf=(wlinf-.25)/0.0025+1.5
isup=(wlsup-.25)/0.0025+1.5
do 1113 ik=iinf,isup
s(ik)=0.
1113 continue
read(iread,*) (s(i),i=iinf,isup)
c**********************************************************************c
c c
c c
c start of computations c
c c
c c
c c
c**********************************************************************c
c goto 800
call otb_6S(asol,phi0,avis,phiv,month,jday,
I pressure, uw, uo3,
I iaer, taer55,
I wlinf, wlsup, s,
O otb_ratm,sast,tgasm,sdtott,sutott)
write(6,*)
800 write(6,100) otb_ratm
write(6,200) sast
write(6,300) tgasm
write(6,400) sdtott
write(6,500) sutott
100 format(10x,40h -----> atmospheric reflectance : , f6.5)
200 format(10x,40h -----> atmospheric spherical albedo : , f6.5)
300 format(10x,40h -----> total gaseous transmission : , f6.5)
400 format(10x,40h -----> downward transmittance : , f6.5)
500 format(10x,40h -----> upward transmittance : , f6.5)
999 continue
end
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment