Hi Zhiqiang,
Just in case you might want to run the HL DOS executable in Linux
through wine, I figured out that you don't need to create a FAT32
partition. Now there is an EXT3 (Linux) partition driver for Windows :
http://ext2fsd.sourceforge.net/
With this tool, you can mount and read/write Linux EXT2 and EXT3
partitions from Windows. After installing ext2fsd in Windows, you only
have to install HL in windows into your EXT3 Linux partition and voila!
You can switch to Linux and use the HL executable (and the GUI) through
Wine or run your own executable you compiled for Linux.
a minor detail but just in case...
best
Servet
Hello Zhiqiang,
Are you sure you don't want to run HL through "wine" in Linux and avoid
compiling? It's very easy. First install "wine".
Then using the GUI, in windows, create your HL run (let's name it ac9_013).
Then, in Linux, open a shell prompt and type :
cd /h42/maincode
wine maincode.exe < ../run/batch/Iac9_013.txt
If you don't have the new Linux drivers that allow read/write to a
windows NTFS partition, you might have to install HL on a FAT32
partition so that Linux can write to it.
***************************************************************************
If you really want to compile the code, here is how I did...
As I told you, I used Intel Fortran compiler under Linux, which is a
commercial software. You might want to try out other compilers. I'd be
interested in knowing the results.
*in IOshorten.f90, comment out all the windows module declarations
and calls to the function GetShortPathName
*In IOshorten.f90 (near line 47), modify the open statement
IF(lreturn.gt.0) then
open(iounit, file=shortfn, status='old', ACTION='read')
to read
IF(1) then
open(iounit, file=longfn, status='old', ACTION='read')
* Rename the file Dimens.inc to DIMENS.INC
* Remove the last lines of
SLAcom.f
BLAS.f
* Modify Initial.f line 693 for
write(10,fmt="(//' ERROR:',/' depths are not monotonically incre
1asing, or are too close together (deltazK =',f6.4,')',)") deltazK
to read
write(10,fmt="(//' ERROR:',/' depths are not monotonically incre
1asing, or are too close together (deltazK =',f6.4,')')") deltazK
* Modify irradat.f line 13 for
INCLUDE 'dimens.inc'
to read
INCLUDE 'DIMENS.INC'
* Modify setdflts.f around lines 64-68
comment lines 64-65 to read
c sl='\ '
c sysl=sl(1:1)
and comment-out lines 67-68 to read
bsl='/ '
sysl=bsl(1:1)
*IN gcirrad.f, main.f and setdflts.f:
replace all occurrences of
character*24 gcifile
to
character*30 gcifile
*Run the GUI to create a specific HL run (for example, let's name this
run ac9_St_013). The GUI will create the file batch/ac9_St_013.for
* To compile hydrolight:
cd maincode
ifort -c *.f
(this will create all the object files)
ifort batch/ac9_St_013.for *.o -I. -o hdrl_exec_ac9
(this will compile the hdrl_exec_ac9.for file that was created by the
GUI for your specific HL run and link it into the executable file called
hdrl_exec_ac9 ).
To run HL, type in the shell prompt:
hdrl_exec_ac9 < ../run/batch/Iac9_St_013.txt
Hope it works
good luck
servet
Zhiqiang Chen wrote:
> Hi Servet:
>
> Thank you very much for your message. It is exactly my purpose to avoid
> GUI to run HL because i am trying to run quite a lot simulations.
>
> Could you please fill me up more information about how to compile those
> code? I am having HL version 4.2 and using Linux FC4 OS. I am not
> familar with compilation of codes, thus you are really appreciated if
> you can provide information as detailed as possible.
> such as
> what compiler??
> what switch should be turn on?
> is a makefile available?
>
> thanks;
>
> zhiqiang
>
> */Servet Ahmet Cizmeli <sa.cizmeli@...>/* wrote:
>
> Hi Zhiquang
>
> You can always compile the source code under Linux. I succeeded under
> Linux with the (commercial) Intel fortran compiler. But this means that
> you won't have access to the graphical user interface (GUI) and you'll
> have to modifiy the input files and run the program by hand. In
> order to
> compile, you might have to modify a couple of fortran files but the
> modifications are minor. Let me know if you are willing to take the
> challenge.
>
> I also tried to run the GUI under Linux through "wine" (a software for
> Linux to run windows programs). I nearly succeeded (went through
> many of
> the input windows) but had an error at the end. I did not follow it up
> because I did not have time at that moment but I feel that the GUI can
> be made functional with "wine".
>
> Whether or not you succeed to run the GUI in "wine", if you are under
> Linux, you can always easily run the DOS executable maincode.exe under
> "wine". I am sure it would work like a charm.
>
> good luck
> servet
>
> Zhiqiang Chen wrote:
> >
> >
> > Hi all:
> >
> > Can I ask a question about HL?
> >
> > I read a message talking about a Linux version of HL? Has someone
> > suceeded in this effort? If so, how can I get a copy of it?
> >
> > thanks;
> >
> > zhiqiang
> >
>
Hi Zhiquang
You can always compile the source code under Linux. I succeeded under
Linux with the (commercial) Intel fortran compiler. But this means that
you won't have access to the graphical user interface (GUI) and you'll
have to modifiy the input files and run the program by hand. In order to
compile, you might have to modify a couple of fortran files but the
modifications are minor. Let me know if you are willing to take the
challenge.
I also tried to run the GUI under Linux through "wine" (a software for
Linux to run windows programs). I nearly succeeded (went through many of
the input windows) but had an error at the end. I did not follow it up
because I did not have time at that moment but I feel that the GUI can
be made functional with "wine".
Whether or not you succeed to run the GUI in "wine", if you are under
Linux, you can always easily run the DOS executable maincode.exe under
"wine". I am sure it would work like a charm.
good luck
servet
Zhiqiang Chen wrote:
>
>
> Hi all:
>
> Can I ask a question about HL?
>
> I read a message talking about a Linux version of HL? Has someone
> suceeded in this effort? If so, how can I get a copy of it?
>
> thanks;
>
> zhiqiang
>
> ------------------------------------------------------------------------
> Fussy? Opinionated? Impossible to please? Perfect. Join Yahoo!'s user
> panel
>
<http://us.rd.yahoo.com/evt=48516/*http://surveylink.yahoo.com/gmrs/yahoo_panel_\
invite.asp?a=7
> > and lay it on us.
>
You can change several limits, including the maximum number of wavebands within the UI.
From the RUN ID form at the beginning of the UI, select "change limits" to access the PARAMETER LIMIT form. This is discussed in the Hydrolight Users' Guide in Section 5.1, page 34.
Regards,
Lydia Sundman
On Aug 2, 2007, at 4:42 AM, Elena Torrecilla wrote:
Hello everyone,
I am a new Hydrolight's user. I would like to know if there is any way to increase the maximum number of bands in the multi-wavelength simulations? Is that parameter easy to be extended?
Hello everyone,
I am a new Hydrolight's user. I would like to know if there is any way
to increase the maximum number of bands in the multi-wavelength
simulations? Is that parameter easy to be extended?
Thanks,
Elena
It really looks like the most important information is Eddif and eddir
if you look at what the skyirrad subroutine returns. Since those are the
only things returned, that is all you really need.
-M
Servet Ahmet Cizmeli wrote:
>Hi Marcos,
>
>Thanks a lot for your precious comments. I'll try to write a fortran
>routine that will read the diffuse spectral irradiance along with Edtot.
>
>best
>Servet
>
>
>
>Yahoo! Groups Links
>
>
>
>
>
>
>
--
Marcos Montes, Ph.D.
Research Physicist
Code 7232, Naval Research Laboratory, Washington, DC 20375
marcos.montes@... \W: 202-767-7308 \ Fax: 202-404-5689
Hi Marcos,
Thanks a lot for your precious comments. I'll try to write a fortran
routine that will read the diffuse spectral irradiance along with Edtot.
best
Servet
Servet,
Just to clarify a few items. On Friday in an email to the HLUG, I said:
> From my reading of the code, there is no way for you to put in the
> diffuse and direct components without changing it. The code needs to
> be edited in order to do that. The simplest approach is probably:
> 1) Make your input irradiance file have two additional columns, one
> for diffuse and one for direct irradiance
> 2) Edit irradat.f to correctly read, store , and return this information
> 3) Edit incfiles.for to not call gcirrad
Note that the above steps do not perform a compilation. They merely
consist of editing files.
Really, I need to add a few steps:
First, you should really completely read the Hydrolight 4.2 Technical
Documentation (h42techdoc.pdf, in the Documents folder), especially
portions of Appendix A (especially "Other Troubleshooting" and "General
Information") and Appendix D.1.
4) At this point, you need to compile the files to make an executable.
This may be done as described in the Hydrolight Technical Manual - see
Appendix A "General Information" section. What you need to do is to open
a windows command prompt (or dos command prompt) window. Change drives
to the drive where the Hydrolight directory lives on. Change directories
to the Hydrolight's maincode directory. Follow the compilation
instructions given in Appendix A. This creates an executable known as
maincode.exe.
Then change to the Hydrolight run directory. Then you need to issue the
following command:
..\maincode\maincode.exe < .\batch\Iroot.txt
but instead of "Iroot.txt" use the name of the input file that you want
to use, which usually exists in the HL's run\batch folder.
Good luck,
-M
Servet Ahmet Cizmeli wrote:
>I understand it much better now, especially after Marcos' comments.
>
>Thanks Marcos, thanks Lydia.
>
>One more question: Looking at my newly compiled incfiles.for, I see that
>the routine GCIRRAD is still being called to estimate the diffuse and
>direct components of Edtot (which is now read off disk from my data
>file). Is there a way to incorporate into H42 our spectral measurements
>of eddir and eddif? The only option I see in the GUI is the parameter
>"cloud" which is a scalar between [0-1] and does not fully account for
>the spectral shape of the components.
>
>
>
>Yahoo! Groups Links
>
>
>
>
>
>
>
--
Marcos Montes, Ph.D.
Research Physicist
Code 7232, Naval Research Laboratory, Washington, DC 20375
marcos.montes@... \W: 202-767-7308 \ Fax: 202-404-5689
I understand it much better now, especially after Marcos' comments.
Thanks Marcos, thanks Lydia.
One more question: Looking at my newly compiled incfiles.for, I see that
the routine GCIRRAD is still being called to estimate the diffuse and
direct components of Edtot (which is now read off disk from my data
file). Is there a way to incorporate into H42 our spectral measurements
of eddir and eddif? The only option I see in the GUI is the parameter
"cloud" which is a scalar between [0-1] and does not fully account for
the spectral shape of the components.
Servet,
Ah, good, you found incfiles.for. When I wrote the email, I had
forgotten the incfiles.for is the actual one that gets compiled. There
is also a copy in the batch folder.
From my reading of the code, there is no way for you to put in the
diffuse and direct components without changing it. The code needs to be
edited in order to do that. The simplest approach is probably:
1) Make your input irradiance file have two additional columns, one for
diffuse and one for direct irradiance
2) Edit irradat.f to correctly read, store , and return this information
3) Edit incfiles.for to not call gcirrad
-M
Servet Ahmet Cizmeli wrote:
>I understand it much better now, especially after Marcos' comments.
>
>Thanks Marcos, thanks Lydia.
>
>One more question: Looking at my newly compiled incfiles.for, I see that
>the routine GCIRRAD is still being called to estimate the diffuse and
>direct components of Edtot (which is now read off disk from my data
>file). Is there a way to incorporate into H42 our spectral measurements
>of eddir and eddif? The only option I see in the GUI is the parameter
>"cloud" which is a scalar between [0-1] and does not fully account for
>the spectral shape of the components.
>
>
>
>Yahoo! Groups Links
>
>
>
>
>
>
>
--
Marcos Montes, Ph.D.
Research Physicist
Code 7232, Naval Research Laboratory, Washington, DC 20375
marcos.montes@... \W: 202-767-7308 \ Fax: 202-404-5689
I understand it much better now, especially after Marcos' comments.
Thanks Marcos, thanks Lydia.
One more question: Looking at my newly compiled incfiles.for, I see that
the routine GCIRRAD is still being called to estimate the diffuse and
direct components of Edtot (which is now read off disk from my data
file). Is there a way to incorporate into H42 our spectral measurements
of eddir and eddif? The only option I see in the GUI is the parameter
"cloud" which is a scalar between [0-1] and does not fully account for
the spectral shape of the components.
Folks,
Marcos is correct that modifying runs outside of the UI is a task
that requires a lot of courage and should not be done cavalierly.
I also wanted to add that when we [finally] ,release v.5 of the code
(HE5, a free upgrade to all H42 users), the code will need to
recompile only very RARELY, in cases where the user has written their
own subroutine(s). All input parameters will be stored solely in
the IROOT.TXT file and the vast majority of runs can be made using a
standard, pre-compiled (or compile once, if dimens.inc is changed or
you port to a different OS) code. We did this to (1) make the
Iroot.txt files 'stand-alone' so runs can be repeated later, and (2)
to make it easier for users to call HE5 recursively from other apps
(though since it is not a DLL will still require some work on said
user's side). Of course, this required some changes to the IROOT.TXT
format.
Other improvements include better ozone cllimatology (based on the
TOMS v8), the ability to run either Hydrolight or Ecolight (screaming
fast, azimuthally-averaged version), extending data files farther
into the near IR and UV, and surface files that account for the
change in index of refraction with wavelength.
Anyway, just a glimpse into the (hopefully very near) future!
Regards,
Lydia
In Hydrolight (HL), not everything is controlled by the input file.
Some parameters are specified in the input file. Some are specified in
the source code. Anytime that any piece of the source code changes, HL
needs to be recompiled. The user interface (UI) edits and write the
input file and it also edits some of the source code.
So, when certain
parameters are changed for calculations, what ends up happening is that
both the source code and the input file are modified. This is why HL
needs to be recompiled. And this is why you should always run HL from
the UI unless you understand all the ins and outs about the how the
source code gets modified (and you should also know Fortran, and be
comfortable compiling HL by hand).
The piece of source code
that usually gets modified is in the maincode\batch\ folder
(directory); maincode is somewhere in the the main Hydrolight
Directory. [Note that the file dimens.inc can also be
modified through
the user interface, but it is of no concern for this question.] In your
case the file name is probably
maincode\batch\irrad.for (sort by date and see the newest
file). If you
go completely through the UI, you'll see that this file changes
depending on the specifications of the run. This file has several
subroutines. The one we are concerned with currently is the one called subroutine
skyirrad. When you choose to input the irradiance, then yours
probably looks something like: c -----------------------------------------
subroutine skyirrad(suntheta,sunphi,eddif,eddir)
c
c This file selects the sky irradiance model to be used in the run,
c and is called by routine qasky.
c
real suntheta,sunphi,eddif,eddir
c
INCLUDE 'dimens.inc'
c
Character surfname*120,pfname*120,
1 Drootname*120,Srootname*120,Mrootname*120,
2 datafiles*120
COMMON /Cfilenames/ surfname,pfname(mxcomp),
1 Drootname,Srootname,Mrootname,
2 datafiles(0:7+mxcomp)
real dirfrac, Edtot
c call irradat to read the specified datafile
call irradat( Edtot, datafiles(11) )
c
c Routine GCIRRAD provides the RADTRAN direct:diffuse ratio
write(10,531)
call gcirrad(suntheta,sunphi,eddif,eddir)
dirfrac = eddir/(eddir+eddif)
c
c Use dirfrac from RADTRAN and the total Ed from the data file
c to calculate the direct and diffuse parts
eddir = dirfrac * Edtot
eddif = Edtot - eddir
c
write(10,533) dirfrac
write(10,532) Eddif,Eddir,Edtot
531 format(//2x,'Gregg and Carder values used to get the ratio of ',
1'diffuse to direct irradiance'/)
532 format(//2x,'Sky spectral irradiances used by Hydrolight for ',
1'this run (from datafile and G&C ratio):'/10x,'Ed(diffuse) =',
2 1pe11.3,5x,'Ed(direct) =',e11.3,5x,'Ed(total) =',e11.3)
533 format(5x,'Ratio of direct to total Irradiance ',
1'as calculated from the Gregg and Carder model:',f6.3)
return
end
c -----------------------------------------
The line I highlighted in red is the one that calls subroutine irradat (in the
maincode directory) that reads the specified file (whose name is stored
in the variable datafiles(11), but it is
not always 11 - it depends on the number of components in your model).
Note that folks who use the UI and seleclt Radtran to set the
irradiance would typically have a skyirrad subroutine that looks something like: c
-----------------------------------------
subroutine skyirrad(suntheta,sunphi,eddif,eddir)
c
c This file selects the sky irradiance model to be used in the run,
c and is called by routine qasky.
c
real suntheta,sunphi,eddif,eddir
c
INCLUDE 'dimens.inc'
c
call gcirrad(suntheta,sunphi,eddif,eddir)
return
end
c ----------------------------------------- There is
no call to the subroutine irradat and the input file is not read in
this case, as expected.
Oh, it is also informative to look at the rest of the file in the batch
directory. If you are using fixed component concentrations, you'll see
that they are listed in the same file (again, in your case,
maincode\batch\irrad.for) and are in fact hard wired
into the source
code in the subroutine
abmodel.
So, in order to do what you want using the current generation of the
code:
1) If you are changing only the input irradiance
file name, change that name in the
input file (the last line). Make sure you've compiled a version that
includes a
subroutine
skyirrad
that has a call
irradat
line. Then you can make
several input files each with a different irradiance file name, and run
to your heart's content without recompiling.
2) If you are changing multiple parameters in set of HydroLight runs,
then you need to be more
careful because, as mentioned at the top of this message, some items
are compiled into the code (and so not affected by the input file) and
some are not. Any time one of the items that is compiled into the code
is changed, then HydroLight needs to be recompiled.
-Marcos
Servet Ahmet Cizmeli wrote:
Hi everybody,
I would like to use the spectral downwelling irradiance measurements
we performed above water surface instead of the Radtran sky model.
With the GUI (H42), I created two runs:
1. a run that uses the option RADTRAN
2. a second run that is exactly similar to the first one except that
it uses the option "USER-SUPPLIED DATA FILE" (E:\H42\data\irradata.txt)
When I compare the two Input files associated with these two runs, I
see that they are exactly the same (except the first-two lines). To my understanding, the two runs should create output files that
would exactly be the same... However in the output files, Ed above
surface are slightly different...
Am I missing something here? What is the parameter in the Input file
that makes Hydrolight use the sky model we want and to what RECORD
GROUP it does belong? I need to perfectly understand how it works as I
am building an automated system that will bypass the GUI and modify
the Input file for each run.
thanks from advance for your input
best
servet
My Input file:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
irrad
irrad
0 1 2 1 0 2 4 1 440 0.1 0.014 3 440 0.1 0.014 -2 440 1 0.014 3 440 0 1 ..\data\pfh2oab.txt
..\data\defaults\apstarchl.txt
dummyastar.txt
..\data\defaults\apstarchl.txt
0 -999 -999 -999 -999 -999 1 -999 -999 -999 -999 -999 bstarDummy.txt
dummybstar.txt
0 0 0 0.0005 0 0 0 0.0005 pureh2o.dpf
avgpart.dpf
27 405, 415, 427, 440, 450, 460, 470, 485, 495, 505, 515, 525, 535, 545, 560, 570, 580, 590, 600, 615, 625, 635, 645, 660, 670, 680, 685, 700, 0, 1, 0, 0, 4
2, 3, 30, 0, 0
0.00
0, 0
0, 11, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, ..\data\pfh2oab.txt
1
E:\H42\data\a_phystar_2001_085.txt
dummyFilteredAc9.txt
dummyHscat.txt
E:\H42\data\CHL_2001_085.txt
dummyCDOMdata.txt
dummyR.bot
dummydata.txt
dummyComp.txt
E:\H42\data\irradata.txt
Yahoo! Groups Links
<*> To visit your group on the web, go to:
http://groups.yahoo.com/group/HydroLightUsers/
<*> Your email settings:
Individual Email | Traditional
<*> To change settings online go to:
http://groups.yahoo.com/group/HydroLightUsers/join
(Yahoo! ID required)
<*> To change settings via email:
mailto:HydroLightUsers-digest@yahoogroups.commailto:HydroLightUsers-fullfeatured@yahoogroups.com
<*> To unsubscribe from this group, send an email to:
HydroLightUsers-unsubscribe@yahoogroups.com
<*> Your use of Yahoo! Groups is subject to:
http://docs.yahoo.com/info/terms/
-- Marcos Montes, Ph.D.
Research Physicist
Code 7232, Naval Research Laboratory, Washington, DC 20375
marcos.montes@... \W: 202-767-7308 \ Fax: 202-404-5689
Even though the input file for your two runs look the same, the aren't the same run. The reason is that when you select RADTRAN, it calculates on the fly the irradiances based on parameters such as the day of year (for earth-sun distance), latitude, longitude, and atmospheric parameters (including ozone, atmos type, water vapor content, etc). In H42 those parameters are read from gcirrad.dat (written to the file by the UI). When you select RADTRAN as your atmosphere option, the input line that has the filename "..data/irradata.txt" is never used.
However, when you select the data file option, none of those parameters are used, no model is called -- the code simply reads the values from the data file and sets the downwellign Irradiance just above the surface to the values in the file.
It is purely coincidence that the two runs even have similar surface Eds -- it just means that the conditions you put into RADTRAN were similar to the ones when the irradata.txt file was generated.
Regards,
Lydia
On Jul 5, 2007, at 10:46 PM, Servet Ahmet Cizmeli wrote:
Hi everybody,
I would like to use the spectral downwelling irradiance measurements we performed above water surface instead of the Radtran sky model. With the GUI (H42), I created two runs:
1. a run that uses the option RADTRAN 2. a second run that is exactly similar to the first one except that it uses the option "USER-SUPPLIED DATA FILE" (E:\H42\data\irradata.txt)
When I compare the two Input files associated with these two runs, I see that they are exactly the same (except the first-two lines).
To my understanding, the two runs should create output files that would exactly be the same... However in the output files, Ed above surface are slightly different...
Am I missing something here? What is the parameter in the Input file that makes Hydrolight use the sky model we want and to what RECORD GROUP it does belong? I need to perfectly understand how it works as I am building an automated system that will bypass the GUI and modify the Input file for each run.
Hi everybody,
I would like to use the spectral downwelling irradiance measurements
we performed above water surface instead of the Radtran sky model.
With the GUI (H42), I created two runs:
1. a run that uses the option RADTRAN
2. a second run that is exactly similar to the first one except that
it uses the option "USER-SUPPLIED DATA FILE" (E:\H42\data\irradata.txt)
When I compare the two Input files associated with these two runs, I
see that they are exactly the same (except the first-two lines).
To my understanding, the two runs should create output files that
would exactly be the same... However in the output files, Ed above
surface are slightly different...
Am I missing something here? What is the parameter in the Input file
that makes Hydrolight use the sky model we want and to what RECORD
GROUP it does belong? I need to perfectly understand how it works as I
am building an automated system that will bypass the GUI and modify
the Input file for each run.
thanks from advance for your input
best
servet
My Input file:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
irrad
irrad
0 1 2 1 0
2 4
1 440 0.1 0.014
3 440 0.1 0.014
-2 440 1 0.014
3 440 0 1
..\data\pfh2oab.txt
..\data\defaults\apstarchl.txt
dummyastar.txt
..\data\defaults\apstarchl.txt
0 -999 -999 -999 -999
-999
1 -999 -999 -999 -999
-999
bstarDummy.txt
dummybstar.txt
0 0 0 0.0005
0 0 0 0.0005
pureh2o.dpf
avgpart.dpf
27
405, 415, 427, 440, 450, 460, 470, 485, 495, 505,
515, 525, 535, 545, 560, 570, 580, 590, 600, 615,
625, 635, 645, 660, 670, 680, 685, 700,
0, 1, 0, 0, 4
2, 3, 30, 0, 0
0.00
0, 0
0, 11, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
10,
..\data\pfh2oab.txt
1
E:\H42\data\a_phystar_2001_085.txt
dummyFilteredAc9.txt
dummyHscat.txt
E:\H42\data\CHL_2001_085.txt
dummyCDOMdata.txt
dummyR.bot
dummydata.txt
dummyComp.txt
E:\H42\data\irradata.txt
Correct. Thus a direction of theta = 45 deg, phi = 90 deg, is the same
direction either above or below the surface. The theta,phi coordinate
system does not change depending on where you are. If you're in the air
(depth flagged as -1) looking down in a given (theta,phi) direction,
then you're seeing the total upwelling radiance (in air) in that
direction. If you're in the water looking in the same downward
direction, you're seeing the upwelling (underwater) radiance in that
direction.
Refraction at the surface is accounted for by the fact that radiance
from one quad in the air maps into one (or more) different quads underwater.
Curt
cindyzju wrote:
>
> So in the full radiance output file, if the radiance is just below
> the water surface (like say 0), the angle in that radiance table
> (theda and phi) should be referring to under water angle while the
> angle for the above surface (like say -1) radiance should be in the
> air. Am I right?
>
> Thanks
> --- In HydroLightUsers@yahoogroups.com
> <mailto:HydroLightUsers%40yahoogroups.com>, Curtis Mobley
> <curtis.mobley@...> wrote:
> >
> > The directions in H are defined by the partitioning of all
> directions
> > into "quads", as shown in Fig 2 of the Users Guide, and as
> described in
> > section 2.4 of the Users Guide. These directions/quads are the
> same
> > above and below the sea surface. Thus, for example, if the sun is
> in a
> > particular quad above the surface, its refracted beam will be in a
> > different quad below the surface (for a level surface). Let me
> know if
> > this doesn't answer your question.
> >
> > Curt
> >
> > cindyzju wrote:
> >
> > >
> > > If I use the full radiance output setting, the angle I get is
> below
> > > the water or in the air or depends on the which depth I am
> gettting
> > > from?
> > >
> > > thanks
> > >
> > >
> >
> >
> > --
> > -----------------------------------
> > Curtis D. Mobley, Ph.D.
> > Vice President and Senior Scientist
> > Sequoia Scientific, Inc.
> > 2700 Richards Road, Suite 107
> > Bellevue, WA 98005
> > voice: 425-641-0944 ext 109
> > fax: 425-643-0595
> > email: curtis.mobley@...
> > WWW: www.sequoiasci.com
> > personal: www.curtismobley.com
> > -----------------------------------
> >
>
>
--
-----------------------------------
Curtis D. Mobley, Ph.D.
Vice President and Senior Scientist
Sequoia Scientific, Inc.
2700 Richards Road, Suite 107
Bellevue, WA 98005
voice: 425-641-0944 ext 109
fax: 425-643-0595
email: curtis.mobley@...
WWW: www.sequoiasci.com
personal: www.curtismobley.com
-----------------------------------
So in the full radiance output file, if the radiance is just below
the water surface (like say 0), the angle in that radiance table
(theda and phi) should be referring to under water angle while the
angle for the above surface (like say -1) radiance should be in the
air. Am I right?
Thanks
--- In HydroLightUsers@yahoogroups.com, Curtis Mobley
<curtis.mobley@...> wrote:
>
> The directions in H are defined by the partitioning of all
directions
> into "quads", as shown in Fig 2 of the Users Guide, and as
described in
> section 2.4 of the Users Guide. These directions/quads are the
same
> above and below the sea surface. Thus, for example, if the sun is
in a
> particular quad above the surface, its refracted beam will be in a
> different quad below the surface (for a level surface). Let me
know if
> this doesn't answer your question.
>
> Curt
>
> cindyzju wrote:
>
> >
> > If I use the full radiance output setting, the angle I get is
below
> > the water or in the air or depends on the which depth I am
gettting
> > from?
> >
> > thanks
> >
> >
>
>
> --
> -----------------------------------
> Curtis D. Mobley, Ph.D.
> Vice President and Senior Scientist
> Sequoia Scientific, Inc.
> 2700 Richards Road, Suite 107
> Bellevue, WA 98005
> voice: 425-641-0944 ext 109
> fax: 425-643-0595
> email: curtis.mobley@...
> WWW: www.sequoiasci.com
> personal: www.curtismobley.com
> -----------------------------------
>
The directions in H are defined by the partitioning of all directions
into "quads", as shown in Fig 2 of the Users Guide, and as described in
section 2.4 of the Users Guide. These directions/quads are the same
above and below the sea surface. Thus, for example, if the sun is in a
particular quad above the surface, its refracted beam will be in a
different quad below the surface (for a level surface). Let me know if
this doesn't answer your question.
Curt
cindyzju wrote:
>
> If I use the full radiance output setting, the angle I get is below
> the water or in the air or depends on the which depth I am gettting
> from?
>
> thanks
>
>
--
-----------------------------------
Curtis D. Mobley, Ph.D.
Vice President and Senior Scientist
Sequoia Scientific, Inc.
2700 Richards Road, Suite 107
Bellevue, WA 98005
voice: 425-641-0944 ext 109
fax: 425-643-0595
email: curtis.mobley@...
WWW: www.sequoiasci.com
personal: www.curtismobley.com
-----------------------------------
I have added a HydroLight Technical Note HTN5_Color.pdf discussing the
conversion of spectra to CIE (x,y,Y) and RGB values. The current H4.3
version outputs CIE (x,y,Y) values as described here. If you need
this and your version doesn't have the new code, let me know and I'll
send you the routines.
Hi everyone,
I've submitted a file at that describes a faster version of the
selpfbb subroutine. This subroutine is mostly used when there are
measured values of a, b, and bb. In this case, the new version of the
subroutine provides much faster calculations than versions provided
with Hydrolight. This is important because accurate calculations need
a fairly small value for the tolerance in bb/b at the ABACBB screen of
the user interface.
The improvement in the speed of the calculations and the results of
the calculations will vary depending on the IOPs, the physical
specifications (how deep, how many wavelengths, etc.) as well as
depending on the particular version of Hydrolight that you use. This
subroutine has been tested with a few different revisions of
Hydrolight 4.2 (an early version of Hydrolight 4.2, as well as the
with the latest revision). I have not yet tested it with Hydolight
4.1, some versions of which also have the selpfbb.f subroutine. I'd be
interested in hearing from anyone who has Hydrolight 4.1 about how
well it works.
The pdf file describing the subroutine is in the "Files" section at:
http://groups.yahoo.com/group/HydroLightUsers/files/selpfbb_faster_MJM.pdf
The new subroutine is also in the "Files" section, and is at:
http://groups.yahoo.com/group/HydroLightUsers/files/selpfbb.f
-Marcos Montes
--
The principal reason is as stated in the HL User's Guide.
For your particular example, choosing a maximum depth implies that
there are no sources below that depth (fluorescence, Raman scattering,
bioluminescence) whose light can reach the surface and also implies
that the inherent optical properties (absorption, scattering, phase
function) do not vary as a function of depth below the lowest depth.
Fluorescence is a source term, and by stopping at 1m, you are not
including the complete vertical region from which light from the
fluorescing materials can reach the surface. If the total attenuation
at the fluorescing wavelengths is c(wavelength) and if there is
no depth structure, you should go to a depth so that any fluorescence
at any wavelength makes a negligible contribution to the measured
light. Basically, you probably want to go to some depth so that
exp(-c * depth) << 1
(here, exp(x) represents exponentiation, i.e., e raised to the
power x).
Not knowing anything about your IOPs, and being just a bit
conservative, and since exp(-5) < 0.01, you should try a depth so
that c(wavelength) * max_depth = 5
or
max_depth = 5/c(wavelength) ,
where c(wavelength) is the depth independent total attenuation
at the fluorescing wavelength; if c is measured in inverse
meters, then max_depth will be in meters.
Considerations about relevant depths are listed in Section 8 of
Hydrolight User's Guide; also, online you should read the article at
the YahooGroups "Depth Selection for Inelastic Scattering Calculations"
(there is both an email and a file in the file section). The latter
article discusses some concerns about the proper output depth gird that
one needs to use for calculations involving fluorescence and Raman
scattering.
-Marcos Montes
cindyzju wrote:
Using the same IOP value as input with homogeneous and infinite deep water body assumption, the maximum output depth is set 1m for one run and 5m for another run. The output reflectance is the same when chlorophyll fluorescence is not included in the calculation. However if chlorophyll fluorescence is considered in each run, the output reflectance is significantly different in the fluorescence region. The resulted fluorescence with 5m output depth is much larger than that of 1m. Anybody can tell me what is the reason for that?
Yahoo! Groups Links
<*> To visit your group on the web, go to:
http://groups.yahoo.com/group/HydroLightUsers/
<*> To unsubscribe from this group, send an email to:
HydroLightUsers-unsubscribe@yahoogroups.com
<*> Your use of Yahoo! Groups is subject to:
http://docs.yahoo.com/info/terms/
-- Marcos Montes, Ph.D.
Code 7232, Naval Research Laboratory, Washington, DC 20375
marcos.montes@...
W: 202-767-7308, Fax: 202-404-5689
Using the same IOP value as input with homogeneous and infinite deep
water body assumption, the maximum output depth is set 1m for one run
and 5m for another run. The output reflectance is the same when
chlorophyll fluorescence is not included in the calculation. However if
chlorophyll fluorescence is considered in each run, the output
reflectance is significantly different in the fluorescence region. The
resulted fluorescence with 5m output depth is much larger than that of
1m. Anybody can tell me what is the reason for that?
Allow me to answer...
When you choose the Petzold average particle phase function (named
"avgpart.dpf") in the pull-down memu, that phase function is used at
every wavelength and depth for scattering by particles.
When you choose a particular backscatter fraction (bb/b value), then H
generates a Fournier-Forand (not Petzold) phase function having that
bb/b value at every wavelength and depth (generated as described in
Mobley et al, 2002, Phase function effects on oceanic light fields, Appl
Optics 41(6), 1035-1050).
The only way to have H use different phase functions for different
depths and/or wavelengths is to input a file of, say, ac9 data with
b(z,lambda) and a file of, say, HydroScat 6 data with bb(z,lambda).
Then H will used the bb(z,lambda)/b(z,lamba) ratio to determine a
Fournier-Forand phase function with a different backscatter fraction at
each depth and wavelength. When spectral b and bb data are inputted, H
linearly interpolates in depth and wavelength to get values at depths
and wavelengths not in the data sets. For depths/wavelengths outside
the range covered by the data, H uses the values at the nearest
depth/wavelength. E.g., if your b and bb data starts at 410 nm, but the
H run starts at 400 nm, H will use the IOPs at 410 for the wavelength at
400. Thus H interpolates within data ranges but it does not extrapolate
outside the data range.
Let me know if this doesn't answer your question.
Curt
sa.cizmeli wrote:
> Hello everybody,
>
> I would like to have more details about the way Hydrolight extrapolates
> the Petzold bb at other wavelengths than that of Petzold's measurements.
> I am using our ac-9 measurements as input into the Hydrolight and
> playing with the backscattering fraction that is used in the choice of
> the VSF.
>
> I looked throught the documentation but I don't seem to find anything
> about it.
>
> thanks from advance
> servet
>
>
--
-----------------------------------
Curtis D. Mobley, Ph.D.
Vice President and Senior Scientist
Sequoia Scientific, Inc.
2700 Richards Road, Suite 107
Bellevue, WA 98005
voice: 425-641-0944 ext 109
fax: 425-643-0595
email: curtis.mobley@...
WWW: www.sequoiasci.com
personal: www.curtismobley.com
-----------------------------------
Hello everybody,
I would like to have more details about the way Hydrolight extrapolates
the Petzold bb at other wavelengths than that of Petzold's measurements.
I am using our ac-9 measurements as input into the Hydrolight and
playing with the backscattering fraction that is used in the choice of
the VSF.
I looked throught the documentation but I don't seem to find anything
about it.
thanks from advance
servet