ccrmwiki - User contributions [en]

Tsunami simulation

2013-07-24T11:18:40Z

Jan:

A sample tsunami run (impact of Alaska tsunami waves on Cannon Beach, OR) can be downloaded [http://www.ccrm.vims.edu/yinglong/wiki_files/tsunami_ex15.tgz here]. Note that the files are compatible with the serial version of SELFE, but the idea for parallel version is similar.

Typically you need to follow these 2 steps in tsunami simulations:
<OL>
<LI> Deformation run (EX15/Def/ in the sample run): this simulates the earthquake and the set-up of the initial surface
waves. For this you need:
*bdef.gr3 (refer to the user manual) which specifies the total seafloor deformation.
*turn on hotout, dt=1s, 10 steps, rnday=1.158e-4
*We typically run this stage for 10 sec duration and at the end of the run you'll find an output files called 10_hotstart which is then used as hotstart.in for the next (propagation and inundation) stage.
*In addition, use mod_depth.f (inside the bundle), which takes info in bdef.gr3 and hgrid.gr3 (pre-earthquake depths) to generate hgrid.new (post-earthquake depths). The latter is used in the next stage. For completeness you need to attach the boundary condition part (b.tmp in the bundle) of hgrid.gr3 to the end of hgrid.new;

<LI> Propagation and inundation stage (EX15/ in the sample run): this run continues from the deformation
run above with no further seafloor movement. You'll notice that many input files are identical
to the Def/ run, but be careful of differences in hgrid.gr3 (linked to Def/hgrid.new) and param.in
(imm=0 etc).
*For this run we also turn off the advection in the open ocean (up to say 20m isobath) because it's negligible there (adv.gr3) - see for info on CFL number.
*Also note that hotstart.in is from the previous stage. A relatively small time step of 1 sec (by SELFE standard) is used because of the constraint from inundation stage.
The operational time step for tsunami applications is generally in the range of a few seconds.
The inundation results may also be sensitive to the min. depth used in the run (1cm in this example).
Also since only 2 S layers are used in the vertical (quasi-2D mode) we set the bottom friction to 0.
If you have to use non-zero friction you'll need to use more layers and make sure the bottom boundary
layer is resolved. If you use the recent parallel version, you can simply use a 2D model with a proper Manning coefficient.
<LI> After the run is done you can look at global outputs (elevation, depth-averaged velocity etc).
The maximum elevation (maxelev.gr3) and depth-averaged velocity (maxdahv.gr3) are also part of the
outputs (for parallel versions, use trunk/src/Utility/Combining_Scripts/combine_gr3.f90). The maximum inundation can be computed from maxelev.gr3.
</OL>

Tsunami simulation

2013-07-24T08:33:33Z

Jan:

A sample tsunami run (impact of Alaska tsunami waves on Cannon Beach, OR) can be downloaded [http://www.ccrm.vims.edu/yinglong/wiki_files/tsunami_ex15.tgz here]. Note that the files are compatible with the serial version of SELFE, but the idea for parallel version is similar.

Typically you need to follow these 2 steps in tsunami simulations:
<OL>
<LI> Deformation run (EX15/Def/ in the sample run): this simulates the earthquake and the set-up of the initial surface
waves. For this you need:
*bdef.gr3 (refer to the user manual) which specifies the total seafloor deformation.
*turn on hotout, dt=1s, 10 steps, rnday=1.158e-4
*We typically run this stage for 10 sec duration and at the end of the run you'll find an output files called 10_hotstart which is then used as hotstart.in for the next (propagation and inundation) stage.
*In addition, use mod_depth.f (inside the bundle), which takes info in bdef.gr3 and hgrid.gr3 (pre-earthquake depths) to generate hgrid.new (post-earthquake depths). The latter is used in the next stage. For completeness you need to attach the boundary condition part (b.tmp in the bundle) of hgrid.gr3 to the end of hgrid.new;

<LI> Propagation and inundation stage (EX15/ in the sample run): this run continues from the deformation
run above with no further seafloor movement. You'll notice that many input files are identical
to the Def/ run, but be careful of differences in hgrid.gr3 (linked to Def/hgrid.new) and param.in
(imm=0 etc). For this run we also turn off the advection in the open ocean (up to say 20m isobath) because it's negligible
there (adv.gr3) - see for info on CFL number. Also note that hotstart.in is from the previous stage. A relatively small
time step of 1 sec (by SELFE standard) is used because of the constraint from inundation stage.
The operational time step for tsunami applications is generally in the range of a few seconds.
The inundation results may also be sensitive to the min. depth used in the run (1cm in this example).
Also since only 2 S layers are used in the vertical (quasi-2D mode) we set the bottom friction to 0.
If you have to use non-zero friction you'll need to use more layers and make sure the bottom boundary
layer is resolved. If you use the recent parallel version, you can simply use a 2D model with a proper Manning coefficient.
<LI> After the run is done you can look at global outputs (elevation, depth-averaged velocity etc).
The maximum elevation (maxelev.gr3) and depth-averaged velocity (maxdahv.gr3) are also part of the
outputs (for parallel versions, use trunk/src/Utility/Combining_Scripts/combine_gr3.f90). The maximum inundation can be computed from maxelev.gr3.
</OL>

Atmospheric forcing

2013-01-31T13:28:05Z

Jan: /* prc-files (displayed using the ncdump utility) */

==Theory==

[http://ccrm.vims.edu/yinglong/wiki_files/HeatExchange-Theory.pdf Here] is some info about the heat exchange module inside SELFE, based on Dr. Mike Zuluaf's slides.

==File description==
In SELFE, atmoserpic forcings include wind, air pressure and temperature, precipitation, humidity and longwave and shortwave fluxes. These are specified in the netcdf files inside sflux/ dir, and conform to the [http://cf-pcmdi.llnl.gov/documents/cf-conventions/1.0/cf-conventions.html NetCDF Climate and Forecast (CF) Metadata Convention 1.0].

There are 4 types of files in sflux/dir; see [http://www.stccmop.org/CORIE/modeling/selfe/NARR/narr.html this site] for sample files.

<OL>
<LI>sflux_inputs.txt: 
$cat sflux_inputs.txt 
&sflux_inputs ! file name 
start_year = 1998, ! start year 
start_month = 09, ! start month 
start_day = 01, ! start day 
start_hour = 0.0, ! start hour 
utc_start = 0.0, ! offset from UTC in hours, with positive numbers for western hemishpere 
/ 
<LI>sflux_air_1.0*.nc: netcdf files that have time, wind speed, air temperature and pressure and specific humidity;
<LI>sflux_prc_1.0*.nc: netcdf files that have precipitation data;
<LI>sflux_rad_1.0*.nc: netcdf files that have downward long and short (solar) wave radiation fluxes.
</OL>

==Using NARR files for your simulation (North America only)==

First, make sure the [http://ccrm.vims.edu/yinglong/wiki_files/narr.ll NARR grid] covers your hgrid.ll.

In your run directory, mkdir sflux and inside it, create symbolic links to the NARR files. e.g., if you run starts from June 10, 2004 and ends June 20, 2004, then

sflux_air_1.001.nc --> narr_air.2004_06_10.nc 
sflux_air_1.002.nc --> narr_air.2004_06_11.nc 
... 
sflux_air_1.011.nc --> narr_air.2004_06_20.nc 
sflux_air_1.012.nc --> narr_air.2004_06_21.nc (extra day to account for time zone difference) 

Similarly for sflux_rad_*.nc and sflux_prc_*.nc. The number "1" after "air_" denotes first data set used; you can use up to 2 sets in SELFE (which combines them with some given weights set in sflux_subs.F90); here you only use 1 set.

Don't forget to update sflux_inputs.txt also.

==Preparing your own sflux inputs==

After familiarize yourself with the NARR files and their format, you may embark on creating your own nc files. The best way is to modify existing matlab scripts (src/Utility/Sflux_nc/readnc*.m) included in the source code bundle, which have extensively in-line comments to guide you along the way.

'''Conventions for .nc-files in sflux directory and additional files needed for sflux run:'''

'''wind:''' u-component is eastward, v-comp. is northward (normal math convention, not compass convention)

'''windrot_geo2proj.gr3:''' rotates winds in case they do not align with coordinate axes,
i.e. lat/lon

'''watertype.gr3:''' 6 is clear water, 7 is muddiest. Search in selfe_init.F90 for watertype.
===air-files (header displayed using the ncdump utility)===
<code>
netcdf sflux_air_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float uwind(time, ny_grid, nx_grid) ;
uwind:long_name = "Surface Eastward Air Velocity (10m AGL)" ;
uwind:standard_name = "eastward_wind" ;
uwind:units = "m/s" ;
float vwind(time, ny_grid, nx_grid) ;
vwind:long_name = "Surface Northward Air Velocity (10m AGL)" ;
vwind:standard_name = "northward_wind" ;
vwind:units = "m/s" ;
float prmsl(time, ny_grid, nx_grid) ;
prmsl:long_name = "Pressure reduced to MSL" ;
prmsl:standard_name = "air_pressure_at_sea_level" ;
prmsl:units = "Pa" ;
float stmp(time, ny_grid, nx_grid) ;
stmp:long_name = "Surface Air Temperature (2m AGL)" ;
stmp:standard_name = "air_temperature" ;
stmp:units = "K" ;
float spfh(time, ny_grid, nx_grid) ;
spfh:long_name = "Surface Specific Humidity (2m AGL)" ;
spfh:standard_name = "specific_humidity" ;
spfh:units = "1" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

===rad-files (header displayed using the ncdump utility)===
<code>
netcdf sflux_rad_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float dlwrf(time, ny_grid, nx_grid) ;
dlwrf:long_name = "Downward Long Wave Radiation Flux" ;
dlwrf:standard_name = "surface_downwelling_longwave_flux_in_air" ;
dlwrf:units = "W/m^2" ;
float dswrf(time, ny_grid, nx_grid) ;
dswrf:long_name = "Downward Short Wave Radiation Flux" ;
dswrf:standard_name = "surface_downwelling_shortwave_flux_in_air" ;
dswrf:units = "W/m^2" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

===prc-files (header displayed using the ncdump utility)===
<code>
netcdf sflux_prc_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float prate(time, ny_grid, nx_grid) ;
prate:long_name = "Surface Precipitation Rate" ;
prate:standard_name = "precipitation_flux" ;
prate:units = "kg/m^2/s" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

Atmospheric forcing

2013-01-31T13:27:58Z

Jan: /* rad-files (displayed using the ncdump utility) */

==Theory==

[http://ccrm.vims.edu/yinglong/wiki_files/HeatExchange-Theory.pdf Here] is some info about the heat exchange module inside SELFE, based on Dr. Mike Zuluaf's slides.

==File description==
In SELFE, atmoserpic forcings include wind, air pressure and temperature, precipitation, humidity and longwave and shortwave fluxes. These are specified in the netcdf files inside sflux/ dir, and conform to the [http://cf-pcmdi.llnl.gov/documents/cf-conventions/1.0/cf-conventions.html NetCDF Climate and Forecast (CF) Metadata Convention 1.0].

There are 4 types of files in sflux/dir; see [http://www.stccmop.org/CORIE/modeling/selfe/NARR/narr.html this site] for sample files.

<OL>
<LI>sflux_inputs.txt: 
$cat sflux_inputs.txt 
&sflux_inputs ! file name 
start_year = 1998, ! start year 
start_month = 09, ! start month 
start_day = 01, ! start day 
start_hour = 0.0, ! start hour 
utc_start = 0.0, ! offset from UTC in hours, with positive numbers for western hemishpere 
/ 
<LI>sflux_air_1.0*.nc: netcdf files that have time, wind speed, air temperature and pressure and specific humidity;
<LI>sflux_prc_1.0*.nc: netcdf files that have precipitation data;
<LI>sflux_rad_1.0*.nc: netcdf files that have downward long and short (solar) wave radiation fluxes.
</OL>

==Using NARR files for your simulation (North America only)==

First, make sure the [http://ccrm.vims.edu/yinglong/wiki_files/narr.ll NARR grid] covers your hgrid.ll.

In your run directory, mkdir sflux and inside it, create symbolic links to the NARR files. e.g., if you run starts from June 10, 2004 and ends June 20, 2004, then

sflux_air_1.001.nc --> narr_air.2004_06_10.nc 
sflux_air_1.002.nc --> narr_air.2004_06_11.nc 
... 
sflux_air_1.011.nc --> narr_air.2004_06_20.nc 
sflux_air_1.012.nc --> narr_air.2004_06_21.nc (extra day to account for time zone difference) 

Similarly for sflux_rad_*.nc and sflux_prc_*.nc. The number "1" after "air_" denotes first data set used; you can use up to 2 sets in SELFE (which combines them with some given weights set in sflux_subs.F90); here you only use 1 set.

Don't forget to update sflux_inputs.txt also.

==Preparing your own sflux inputs==

After familiarize yourself with the NARR files and their format, you may embark on creating your own nc files. The best way is to modify existing matlab scripts (src/Utility/Sflux_nc/readnc*.m) included in the source code bundle, which have extensively in-line comments to guide you along the way.

'''Conventions for .nc-files in sflux directory and additional files needed for sflux run:'''

'''wind:''' u-component is eastward, v-comp. is northward (normal math convention, not compass convention)

'''windrot_geo2proj.gr3:''' rotates winds in case they do not align with coordinate axes,
i.e. lat/lon

'''watertype.gr3:''' 6 is clear water, 7 is muddiest. Search in selfe_init.F90 for watertype.
===air-files (header displayed using the ncdump utility)===
<code>
netcdf sflux_air_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float uwind(time, ny_grid, nx_grid) ;
uwind:long_name = "Surface Eastward Air Velocity (10m AGL)" ;
uwind:standard_name = "eastward_wind" ;
uwind:units = "m/s" ;
float vwind(time, ny_grid, nx_grid) ;
vwind:long_name = "Surface Northward Air Velocity (10m AGL)" ;
vwind:standard_name = "northward_wind" ;
vwind:units = "m/s" ;
float prmsl(time, ny_grid, nx_grid) ;
prmsl:long_name = "Pressure reduced to MSL" ;
prmsl:standard_name = "air_pressure_at_sea_level" ;
prmsl:units = "Pa" ;
float stmp(time, ny_grid, nx_grid) ;
stmp:long_name = "Surface Air Temperature (2m AGL)" ;
stmp:standard_name = "air_temperature" ;
stmp:units = "K" ;
float spfh(time, ny_grid, nx_grid) ;
spfh:long_name = "Surface Specific Humidity (2m AGL)" ;
spfh:standard_name = "specific_humidity" ;
spfh:units = "1" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

===rad-files (header displayed using the ncdump utility)===
<code>
netcdf sflux_rad_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float dlwrf(time, ny_grid, nx_grid) ;
dlwrf:long_name = "Downward Long Wave Radiation Flux" ;
dlwrf:standard_name = "surface_downwelling_longwave_flux_in_air" ;
dlwrf:units = "W/m^2" ;
float dswrf(time, ny_grid, nx_grid) ;
dswrf:long_name = "Downward Short Wave Radiation Flux" ;
dswrf:standard_name = "surface_downwelling_shortwave_flux_in_air" ;
dswrf:units = "W/m^2" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

===prc-files (displayed using the ncdump utility)===
<code>
netcdf sflux_prc_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float prate(time, ny_grid, nx_grid) ;
prate:long_name = "Surface Precipitation Rate" ;
prate:standard_name = "precipitation_flux" ;
prate:units = "kg/m^2/s" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

Atmospheric forcing

2013-01-31T13:27:50Z

Jan: /* air-files (displayed using the ncdump utility) */

==Theory==

[http://ccrm.vims.edu/yinglong/wiki_files/HeatExchange-Theory.pdf Here] is some info about the heat exchange module inside SELFE, based on Dr. Mike Zuluaf's slides.

==File description==
In SELFE, atmoserpic forcings include wind, air pressure and temperature, precipitation, humidity and longwave and shortwave fluxes. These are specified in the netcdf files inside sflux/ dir, and conform to the [http://cf-pcmdi.llnl.gov/documents/cf-conventions/1.0/cf-conventions.html NetCDF Climate and Forecast (CF) Metadata Convention 1.0].

There are 4 types of files in sflux/dir; see [http://www.stccmop.org/CORIE/modeling/selfe/NARR/narr.html this site] for sample files.

<OL>
<LI>sflux_inputs.txt: 
$cat sflux_inputs.txt 
&sflux_inputs ! file name 
start_year = 1998, ! start year 
start_month = 09, ! start month 
start_day = 01, ! start day 
start_hour = 0.0, ! start hour 
utc_start = 0.0, ! offset from UTC in hours, with positive numbers for western hemishpere 
/ 
<LI>sflux_air_1.0*.nc: netcdf files that have time, wind speed, air temperature and pressure and specific humidity;
<LI>sflux_prc_1.0*.nc: netcdf files that have precipitation data;
<LI>sflux_rad_1.0*.nc: netcdf files that have downward long and short (solar) wave radiation fluxes.
</OL>

==Using NARR files for your simulation (North America only)==

First, make sure the [http://ccrm.vims.edu/yinglong/wiki_files/narr.ll NARR grid] covers your hgrid.ll.

In your run directory, mkdir sflux and inside it, create symbolic links to the NARR files. e.g., if you run starts from June 10, 2004 and ends June 20, 2004, then

sflux_air_1.001.nc --> narr_air.2004_06_10.nc 
sflux_air_1.002.nc --> narr_air.2004_06_11.nc 
... 
sflux_air_1.011.nc --> narr_air.2004_06_20.nc 
sflux_air_1.012.nc --> narr_air.2004_06_21.nc (extra day to account for time zone difference) 

Similarly for sflux_rad_*.nc and sflux_prc_*.nc. The number "1" after "air_" denotes first data set used; you can use up to 2 sets in SELFE (which combines them with some given weights set in sflux_subs.F90); here you only use 1 set.

Don't forget to update sflux_inputs.txt also.

==Preparing your own sflux inputs==

After familiarize yourself with the NARR files and their format, you may embark on creating your own nc files. The best way is to modify existing matlab scripts (src/Utility/Sflux_nc/readnc*.m) included in the source code bundle, which have extensively in-line comments to guide you along the way.

'''Conventions for .nc-files in sflux directory and additional files needed for sflux run:'''

'''wind:''' u-component is eastward, v-comp. is northward (normal math convention, not compass convention)

'''windrot_geo2proj.gr3:''' rotates winds in case they do not align with coordinate axes,
i.e. lat/lon

'''watertype.gr3:''' 6 is clear water, 7 is muddiest. Search in selfe_init.F90 for watertype.
===air-files (header displayed using the ncdump utility)===
<code>
netcdf sflux_air_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float uwind(time, ny_grid, nx_grid) ;
uwind:long_name = "Surface Eastward Air Velocity (10m AGL)" ;
uwind:standard_name = "eastward_wind" ;
uwind:units = "m/s" ;
float vwind(time, ny_grid, nx_grid) ;
vwind:long_name = "Surface Northward Air Velocity (10m AGL)" ;
vwind:standard_name = "northward_wind" ;
vwind:units = "m/s" ;
float prmsl(time, ny_grid, nx_grid) ;
prmsl:long_name = "Pressure reduced to MSL" ;
prmsl:standard_name = "air_pressure_at_sea_level" ;
prmsl:units = "Pa" ;
float stmp(time, ny_grid, nx_grid) ;
stmp:long_name = "Surface Air Temperature (2m AGL)" ;
stmp:standard_name = "air_temperature" ;
stmp:units = "K" ;
float spfh(time, ny_grid, nx_grid) ;
spfh:long_name = "Surface Specific Humidity (2m AGL)" ;
spfh:standard_name = "specific_humidity" ;
spfh:units = "1" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

===rad-files (displayed using the ncdump utility)===
<code>
netcdf sflux_rad_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float dlwrf(time, ny_grid, nx_grid) ;
dlwrf:long_name = "Downward Long Wave Radiation Flux" ;
dlwrf:standard_name = "surface_downwelling_longwave_flux_in_air" ;
dlwrf:units = "W/m^2" ;
float dswrf(time, ny_grid, nx_grid) ;
dswrf:long_name = "Downward Short Wave Radiation Flux" ;
dswrf:standard_name = "surface_downwelling_shortwave_flux_in_air" ;
dswrf:units = "W/m^2" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

===prc-files (displayed using the ncdump utility)===
<code>
netcdf sflux_prc_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float prate(time, ny_grid, nx_grid) ;
prate:long_name = "Surface Precipitation Rate" ;
prate:standard_name = "precipitation_flux" ;
prate:units = "kg/m^2/s" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

Atmospheric forcing

2013-01-31T13:27:26Z

Jan: /* prc-files */

==Theory==

[http://ccrm.vims.edu/yinglong/wiki_files/HeatExchange-Theory.pdf Here] is some info about the heat exchange module inside SELFE, based on Dr. Mike Zuluaf's slides.

==File description==
In SELFE, atmoserpic forcings include wind, air pressure and temperature, precipitation, humidity and longwave and shortwave fluxes. These are specified in the netcdf files inside sflux/ dir, and conform to the [http://cf-pcmdi.llnl.gov/documents/cf-conventions/1.0/cf-conventions.html NetCDF Climate and Forecast (CF) Metadata Convention 1.0].

There are 4 types of files in sflux/dir; see [http://www.stccmop.org/CORIE/modeling/selfe/NARR/narr.html this site] for sample files.

<OL>
<LI>sflux_inputs.txt: 
$cat sflux_inputs.txt 
&sflux_inputs ! file name 
start_year = 1998, ! start year 
start_month = 09, ! start month 
start_day = 01, ! start day 
start_hour = 0.0, ! start hour 
utc_start = 0.0, ! offset from UTC in hours, with positive numbers for western hemishpere 
/ 
<LI>sflux_air_1.0*.nc: netcdf files that have time, wind speed, air temperature and pressure and specific humidity;
<LI>sflux_prc_1.0*.nc: netcdf files that have precipitation data;
<LI>sflux_rad_1.0*.nc: netcdf files that have downward long and short (solar) wave radiation fluxes.
</OL>

==Using NARR files for your simulation (North America only)==

First, make sure the [http://ccrm.vims.edu/yinglong/wiki_files/narr.ll NARR grid] covers your hgrid.ll.

In your run directory, mkdir sflux and inside it, create symbolic links to the NARR files. e.g., if you run starts from June 10, 2004 and ends June 20, 2004, then

sflux_air_1.001.nc --> narr_air.2004_06_10.nc 
sflux_air_1.002.nc --> narr_air.2004_06_11.nc 
... 
sflux_air_1.011.nc --> narr_air.2004_06_20.nc 
sflux_air_1.012.nc --> narr_air.2004_06_21.nc (extra day to account for time zone difference) 

Similarly for sflux_rad_*.nc and sflux_prc_*.nc. The number "1" after "air_" denotes first data set used; you can use up to 2 sets in SELFE (which combines them with some given weights set in sflux_subs.F90); here you only use 1 set.

Don't forget to update sflux_inputs.txt also.

==Preparing your own sflux inputs==

After familiarize yourself with the NARR files and their format, you may embark on creating your own nc files. The best way is to modify existing matlab scripts (src/Utility/Sflux_nc/readnc*.m) included in the source code bundle, which have extensively in-line comments to guide you along the way.

'''Conventions for .nc-files in sflux directory and additional files needed for sflux run:'''

'''wind:''' u-component is eastward, v-comp. is northward (normal math convention, not compass convention)

'''windrot_geo2proj.gr3:''' rotates winds in case they do not align with coordinate axes,
i.e. lat/lon

'''watertype.gr3:''' 6 is clear water, 7 is muddiest. Search in selfe_init.F90 for watertype.
===air-files (displayed using the ncdump utility)===
<code>
netcdf sflux_air_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float uwind(time, ny_grid, nx_grid) ;
uwind:long_name = "Surface Eastward Air Velocity (10m AGL)" ;
uwind:standard_name = "eastward_wind" ;
uwind:units = "m/s" ;
float vwind(time, ny_grid, nx_grid) ;
vwind:long_name = "Surface Northward Air Velocity (10m AGL)" ;
vwind:standard_name = "northward_wind" ;
vwind:units = "m/s" ;
float prmsl(time, ny_grid, nx_grid) ;
prmsl:long_name = "Pressure reduced to MSL" ;
prmsl:standard_name = "air_pressure_at_sea_level" ;
prmsl:units = "Pa" ;
float stmp(time, ny_grid, nx_grid) ;
stmp:long_name = "Surface Air Temperature (2m AGL)" ;
stmp:standard_name = "air_temperature" ;
stmp:units = "K" ;
float spfh(time, ny_grid, nx_grid) ;
spfh:long_name = "Surface Specific Humidity (2m AGL)" ;
spfh:standard_name = "specific_humidity" ;
spfh:units = "1" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

===rad-files (displayed using the ncdump utility)===
<code>
netcdf sflux_rad_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float dlwrf(time, ny_grid, nx_grid) ;
dlwrf:long_name = "Downward Long Wave Radiation Flux" ;
dlwrf:standard_name = "surface_downwelling_longwave_flux_in_air" ;
dlwrf:units = "W/m^2" ;
float dswrf(time, ny_grid, nx_grid) ;
dswrf:long_name = "Downward Short Wave Radiation Flux" ;
dswrf:standard_name = "surface_downwelling_shortwave_flux_in_air" ;
dswrf:units = "W/m^2" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

===prc-files (displayed using the ncdump utility)===
<code>
netcdf sflux_prc_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float prate(time, ny_grid, nx_grid) ;
prate:long_name = "Surface Precipitation Rate" ;
prate:standard_name = "precipitation_flux" ;
prate:units = "kg/m^2/s" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

Atmospheric forcing

2013-01-31T13:27:17Z

Jan: /* rad-files */

==Theory==

[http://ccrm.vims.edu/yinglong/wiki_files/HeatExchange-Theory.pdf Here] is some info about the heat exchange module inside SELFE, based on Dr. Mike Zuluaf's slides.

==File description==
In SELFE, atmoserpic forcings include wind, air pressure and temperature, precipitation, humidity and longwave and shortwave fluxes. These are specified in the netcdf files inside sflux/ dir, and conform to the [http://cf-pcmdi.llnl.gov/documents/cf-conventions/1.0/cf-conventions.html NetCDF Climate and Forecast (CF) Metadata Convention 1.0].

There are 4 types of files in sflux/dir; see [http://www.stccmop.org/CORIE/modeling/selfe/NARR/narr.html this site] for sample files.

<OL>
<LI>sflux_inputs.txt: 
$cat sflux_inputs.txt 
&sflux_inputs ! file name 
start_year = 1998, ! start year 
start_month = 09, ! start month 
start_day = 01, ! start day 
start_hour = 0.0, ! start hour 
utc_start = 0.0, ! offset from UTC in hours, with positive numbers for western hemishpere 
/ 
<LI>sflux_air_1.0*.nc: netcdf files that have time, wind speed, air temperature and pressure and specific humidity;
<LI>sflux_prc_1.0*.nc: netcdf files that have precipitation data;
<LI>sflux_rad_1.0*.nc: netcdf files that have downward long and short (solar) wave radiation fluxes.
</OL>

==Using NARR files for your simulation (North America only)==

First, make sure the [http://ccrm.vims.edu/yinglong/wiki_files/narr.ll NARR grid] covers your hgrid.ll.

In your run directory, mkdir sflux and inside it, create symbolic links to the NARR files. e.g., if you run starts from June 10, 2004 and ends June 20, 2004, then

sflux_air_1.001.nc --> narr_air.2004_06_10.nc 
sflux_air_1.002.nc --> narr_air.2004_06_11.nc 
... 
sflux_air_1.011.nc --> narr_air.2004_06_20.nc 
sflux_air_1.012.nc --> narr_air.2004_06_21.nc (extra day to account for time zone difference) 

Similarly for sflux_rad_*.nc and sflux_prc_*.nc. The number "1" after "air_" denotes first data set used; you can use up to 2 sets in SELFE (which combines them with some given weights set in sflux_subs.F90); here you only use 1 set.

Don't forget to update sflux_inputs.txt also.

==Preparing your own sflux inputs==

After familiarize yourself with the NARR files and their format, you may embark on creating your own nc files. The best way is to modify existing matlab scripts (src/Utility/Sflux_nc/readnc*.m) included in the source code bundle, which have extensively in-line comments to guide you along the way.

'''Conventions for .nc-files in sflux directory and additional files needed for sflux run:'''

'''wind:''' u-component is eastward, v-comp. is northward (normal math convention, not compass convention)

'''windrot_geo2proj.gr3:''' rotates winds in case they do not align with coordinate axes,
i.e. lat/lon

'''watertype.gr3:''' 6 is clear water, 7 is muddiest. Search in selfe_init.F90 for watertype.
===air-files (displayed using the ncdump utility)===
<code>
netcdf sflux_air_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float uwind(time, ny_grid, nx_grid) ;
uwind:long_name = "Surface Eastward Air Velocity (10m AGL)" ;
uwind:standard_name = "eastward_wind" ;
uwind:units = "m/s" ;
float vwind(time, ny_grid, nx_grid) ;
vwind:long_name = "Surface Northward Air Velocity (10m AGL)" ;
vwind:standard_name = "northward_wind" ;
vwind:units = "m/s" ;
float prmsl(time, ny_grid, nx_grid) ;
prmsl:long_name = "Pressure reduced to MSL" ;
prmsl:standard_name = "air_pressure_at_sea_level" ;
prmsl:units = "Pa" ;
float stmp(time, ny_grid, nx_grid) ;
stmp:long_name = "Surface Air Temperature (2m AGL)" ;
stmp:standard_name = "air_temperature" ;
stmp:units = "K" ;
float spfh(time, ny_grid, nx_grid) ;
spfh:long_name = "Surface Specific Humidity (2m AGL)" ;
spfh:standard_name = "specific_humidity" ;
spfh:units = "1" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

===rad-files (displayed using the ncdump utility)===
<code>
netcdf sflux_rad_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float dlwrf(time, ny_grid, nx_grid) ;
dlwrf:long_name = "Downward Long Wave Radiation Flux" ;
dlwrf:standard_name = "surface_downwelling_longwave_flux_in_air" ;
dlwrf:units = "W/m^2" ;
float dswrf(time, ny_grid, nx_grid) ;
dswrf:long_name = "Downward Short Wave Radiation Flux" ;
dswrf:standard_name = "surface_downwelling_shortwave_flux_in_air" ;
dswrf:units = "W/m^2" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

===prc-files===
<code>
netcdf sflux_prc_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float prate(time, ny_grid, nx_grid) ;
prate:long_name = "Surface Precipitation Rate" ;
prate:standard_name = "precipitation_flux" ;
prate:units = "kg/m^2/s" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

Atmospheric forcing

2013-01-31T13:27:09Z

Jan: /* air-files */

==Theory==

[http://ccrm.vims.edu/yinglong/wiki_files/HeatExchange-Theory.pdf Here] is some info about the heat exchange module inside SELFE, based on Dr. Mike Zuluaf's slides.

==File description==
In SELFE, atmoserpic forcings include wind, air pressure and temperature, precipitation, humidity and longwave and shortwave fluxes. These are specified in the netcdf files inside sflux/ dir, and conform to the [http://cf-pcmdi.llnl.gov/documents/cf-conventions/1.0/cf-conventions.html NetCDF Climate and Forecast (CF) Metadata Convention 1.0].

There are 4 types of files in sflux/dir; see [http://www.stccmop.org/CORIE/modeling/selfe/NARR/narr.html this site] for sample files.

<OL>
<LI>sflux_inputs.txt: 
$cat sflux_inputs.txt 
&sflux_inputs ! file name 
start_year = 1998, ! start year 
start_month = 09, ! start month 
start_day = 01, ! start day 
start_hour = 0.0, ! start hour 
utc_start = 0.0, ! offset from UTC in hours, with positive numbers for western hemishpere 
/ 
<LI>sflux_air_1.0*.nc: netcdf files that have time, wind speed, air temperature and pressure and specific humidity;
<LI>sflux_prc_1.0*.nc: netcdf files that have precipitation data;
<LI>sflux_rad_1.0*.nc: netcdf files that have downward long and short (solar) wave radiation fluxes.
</OL>

==Using NARR files for your simulation (North America only)==

First, make sure the [http://ccrm.vims.edu/yinglong/wiki_files/narr.ll NARR grid] covers your hgrid.ll.

In your run directory, mkdir sflux and inside it, create symbolic links to the NARR files. e.g., if you run starts from June 10, 2004 and ends June 20, 2004, then

sflux_air_1.001.nc --> narr_air.2004_06_10.nc 
sflux_air_1.002.nc --> narr_air.2004_06_11.nc 
... 
sflux_air_1.011.nc --> narr_air.2004_06_20.nc 
sflux_air_1.012.nc --> narr_air.2004_06_21.nc (extra day to account for time zone difference) 

Similarly for sflux_rad_*.nc and sflux_prc_*.nc. The number "1" after "air_" denotes first data set used; you can use up to 2 sets in SELFE (which combines them with some given weights set in sflux_subs.F90); here you only use 1 set.

Don't forget to update sflux_inputs.txt also.

==Preparing your own sflux inputs==

After familiarize yourself with the NARR files and their format, you may embark on creating your own nc files. The best way is to modify existing matlab scripts (src/Utility/Sflux_nc/readnc*.m) included in the source code bundle, which have extensively in-line comments to guide you along the way.

'''Conventions for .nc-files in sflux directory and additional files needed for sflux run:'''

'''wind:''' u-component is eastward, v-comp. is northward (normal math convention, not compass convention)

'''windrot_geo2proj.gr3:''' rotates winds in case they do not align with coordinate axes,
i.e. lat/lon

'''watertype.gr3:''' 6 is clear water, 7 is muddiest. Search in selfe_init.F90 for watertype.
===air-files (displayed using the ncdump utility)===
<code>
netcdf sflux_air_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float uwind(time, ny_grid, nx_grid) ;
uwind:long_name = "Surface Eastward Air Velocity (10m AGL)" ;
uwind:standard_name = "eastward_wind" ;
uwind:units = "m/s" ;
float vwind(time, ny_grid, nx_grid) ;
vwind:long_name = "Surface Northward Air Velocity (10m AGL)" ;
vwind:standard_name = "northward_wind" ;
vwind:units = "m/s" ;
float prmsl(time, ny_grid, nx_grid) ;
prmsl:long_name = "Pressure reduced to MSL" ;
prmsl:standard_name = "air_pressure_at_sea_level" ;
prmsl:units = "Pa" ;
float stmp(time, ny_grid, nx_grid) ;
stmp:long_name = "Surface Air Temperature (2m AGL)" ;
stmp:standard_name = "air_temperature" ;
stmp:units = "K" ;
float spfh(time, ny_grid, nx_grid) ;
spfh:long_name = "Surface Specific Humidity (2m AGL)" ;
spfh:standard_name = "specific_humidity" ;
spfh:units = "1" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

===rad-files===
<code>
netcdf sflux_rad_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float dlwrf(time, ny_grid, nx_grid) ;
dlwrf:long_name = "Downward Long Wave Radiation Flux" ;
dlwrf:standard_name = "surface_downwelling_longwave_flux_in_air" ;
dlwrf:units = "W/m^2" ;
float dswrf(time, ny_grid, nx_grid) ;
dswrf:long_name = "Downward Short Wave Radiation Flux" ;
dswrf:standard_name = "surface_downwelling_shortwave_flux_in_air" ;
dswrf:units = "W/m^2" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>
===prc-files===
<code>
netcdf sflux_prc_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float prate(time, ny_grid, nx_grid) ;
prate:long_name = "Surface Precipitation Rate" ;
prate:standard_name = "precipitation_flux" ;
prate:units = "kg/m^2/s" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

Atmospheric forcing

2013-01-31T13:26:28Z

Jan: /* rad-files */

==Theory==

[http://ccrm.vims.edu/yinglong/wiki_files/HeatExchange-Theory.pdf Here] is some info about the heat exchange module inside SELFE, based on Dr. Mike Zuluaf's slides.

==File description==
In SELFE, atmoserpic forcings include wind, air pressure and temperature, precipitation, humidity and longwave and shortwave fluxes. These are specified in the netcdf files inside sflux/ dir, and conform to the [http://cf-pcmdi.llnl.gov/documents/cf-conventions/1.0/cf-conventions.html NetCDF Climate and Forecast (CF) Metadata Convention 1.0].

There are 4 types of files in sflux/dir; see [http://www.stccmop.org/CORIE/modeling/selfe/NARR/narr.html this site] for sample files.

<OL>
<LI>sflux_inputs.txt: 
$cat sflux_inputs.txt 
&sflux_inputs ! file name 
start_year = 1998, ! start year 
start_month = 09, ! start month 
start_day = 01, ! start day 
start_hour = 0.0, ! start hour 
utc_start = 0.0, ! offset from UTC in hours, with positive numbers for western hemishpere 
/ 
<LI>sflux_air_1.0*.nc: netcdf files that have time, wind speed, air temperature and pressure and specific humidity;
<LI>sflux_prc_1.0*.nc: netcdf files that have precipitation data;
<LI>sflux_rad_1.0*.nc: netcdf files that have downward long and short (solar) wave radiation fluxes.
</OL>

==Using NARR files for your simulation (North America only)==

First, make sure the [http://ccrm.vims.edu/yinglong/wiki_files/narr.ll NARR grid] covers your hgrid.ll.

In your run directory, mkdir sflux and inside it, create symbolic links to the NARR files. e.g., if you run starts from June 10, 2004 and ends June 20, 2004, then

sflux_air_1.001.nc --> narr_air.2004_06_10.nc 
sflux_air_1.002.nc --> narr_air.2004_06_11.nc 
... 
sflux_air_1.011.nc --> narr_air.2004_06_20.nc 
sflux_air_1.012.nc --> narr_air.2004_06_21.nc (extra day to account for time zone difference) 

Similarly for sflux_rad_*.nc and sflux_prc_*.nc. The number "1" after "air_" denotes first data set used; you can use up to 2 sets in SELFE (which combines them with some given weights set in sflux_subs.F90); here you only use 1 set.

Don't forget to update sflux_inputs.txt also.

==Preparing your own sflux inputs==

After familiarize yourself with the NARR files and their format, you may embark on creating your own nc files. The best way is to modify existing matlab scripts (src/Utility/Sflux_nc/readnc*.m) included in the source code bundle, which have extensively in-line comments to guide you along the way.

'''Conventions for .nc-files in sflux directory and additional files needed for sflux run:'''

'''wind:''' u-component is eastward, v-comp. is northward (normal math convention, not compass convention)

'''windrot_geo2proj.gr3:''' rotates winds in case they do not align with coordinate axes,
i.e. lat/lon

'''watertype.gr3:''' 6 is clear water, 7 is muddiest. Search in selfe_init.F90 for watertype.
===air-files===
<code>
these lines (displayed using the ncdump utility):
netcdf sflux_air_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float uwind(time, ny_grid, nx_grid) ;
uwind:long_name = "Surface Eastward Air Velocity (10m AGL)" ;
uwind:standard_name = "eastward_wind" ;
uwind:units = "m/s" ;
float vwind(time, ny_grid, nx_grid) ;
vwind:long_name = "Surface Northward Air Velocity (10m AGL)" ;
vwind:standard_name = "northward_wind" ;
vwind:units = "m/s" ;
float prmsl(time, ny_grid, nx_grid) ;
prmsl:long_name = "Pressure reduced to MSL" ;
prmsl:standard_name = "air_pressure_at_sea_level" ;
prmsl:units = "Pa" ;
float stmp(time, ny_grid, nx_grid) ;
stmp:long_name = "Surface Air Temperature (2m AGL)" ;
stmp:standard_name = "air_temperature" ;
stmp:units = "K" ;
float spfh(time, ny_grid, nx_grid) ;
spfh:long_name = "Surface Specific Humidity (2m AGL)" ;
spfh:standard_name = "specific_humidity" ;
spfh:units = "1" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

===rad-files===
<code>
netcdf sflux_rad_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float dlwrf(time, ny_grid, nx_grid) ;
dlwrf:long_name = "Downward Long Wave Radiation Flux" ;
dlwrf:standard_name = "surface_downwelling_longwave_flux_in_air" ;
dlwrf:units = "W/m^2" ;
float dswrf(time, ny_grid, nx_grid) ;
dswrf:long_name = "Downward Short Wave Radiation Flux" ;
dswrf:standard_name = "surface_downwelling_shortwave_flux_in_air" ;
dswrf:units = "W/m^2" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>
===prc-files===
<code>
netcdf sflux_prc_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float prate(time, ny_grid, nx_grid) ;
prate:long_name = "Surface Precipitation Rate" ;
prate:standard_name = "precipitation_flux" ;
prate:units = "kg/m^2/s" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

Atmospheric forcing

2013-01-31T13:26:06Z

Jan: /* Preparing your own sflux inputs */

==Theory==

[http://ccrm.vims.edu/yinglong/wiki_files/HeatExchange-Theory.pdf Here] is some info about the heat exchange module inside SELFE, based on Dr. Mike Zuluaf's slides.

==File description==
In SELFE, atmoserpic forcings include wind, air pressure and temperature, precipitation, humidity and longwave and shortwave fluxes. These are specified in the netcdf files inside sflux/ dir, and conform to the [http://cf-pcmdi.llnl.gov/documents/cf-conventions/1.0/cf-conventions.html NetCDF Climate and Forecast (CF) Metadata Convention 1.0].

There are 4 types of files in sflux/dir; see [http://www.stccmop.org/CORIE/modeling/selfe/NARR/narr.html this site] for sample files.

<OL>
<LI>sflux_inputs.txt: 
$cat sflux_inputs.txt 
&sflux_inputs ! file name 
start_year = 1998, ! start year 
start_month = 09, ! start month 
start_day = 01, ! start day 
start_hour = 0.0, ! start hour 
utc_start = 0.0, ! offset from UTC in hours, with positive numbers for western hemishpere 
/ 
<LI>sflux_air_1.0*.nc: netcdf files that have time, wind speed, air temperature and pressure and specific humidity;
<LI>sflux_prc_1.0*.nc: netcdf files that have precipitation data;
<LI>sflux_rad_1.0*.nc: netcdf files that have downward long and short (solar) wave radiation fluxes.
</OL>

==Using NARR files for your simulation (North America only)==

First, make sure the [http://ccrm.vims.edu/yinglong/wiki_files/narr.ll NARR grid] covers your hgrid.ll.

In your run directory, mkdir sflux and inside it, create symbolic links to the NARR files. e.g., if you run starts from June 10, 2004 and ends June 20, 2004, then

sflux_air_1.001.nc --> narr_air.2004_06_10.nc 
sflux_air_1.002.nc --> narr_air.2004_06_11.nc 
... 
sflux_air_1.011.nc --> narr_air.2004_06_20.nc 
sflux_air_1.012.nc --> narr_air.2004_06_21.nc (extra day to account for time zone difference) 

Similarly for sflux_rad_*.nc and sflux_prc_*.nc. The number "1" after "air_" denotes first data set used; you can use up to 2 sets in SELFE (which combines them with some given weights set in sflux_subs.F90); here you only use 1 set.

Don't forget to update sflux_inputs.txt also.

==Preparing your own sflux inputs==

After familiarize yourself with the NARR files and their format, you may embark on creating your own nc files. The best way is to modify existing matlab scripts (src/Utility/Sflux_nc/readnc*.m) included in the source code bundle, which have extensively in-line comments to guide you along the way.

'''Conventions for .nc-files in sflux directory and additional files needed for sflux run:'''

'''wind:''' u-component is eastward, v-comp. is northward (normal math convention, not compass convention)

'''windrot_geo2proj.gr3:''' rotates winds in case they do not align with coordinate axes,
i.e. lat/lon

'''watertype.gr3:''' 6 is clear water, 7 is muddiest. Search in selfe_init.F90 for watertype.
===air-files===
<code>
these lines (displayed using the ncdump utility):
netcdf sflux_air_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float uwind(time, ny_grid, nx_grid) ;
uwind:long_name = "Surface Eastward Air Velocity (10m AGL)" ;
uwind:standard_name = "eastward_wind" ;
uwind:units = "m/s" ;
float vwind(time, ny_grid, nx_grid) ;
vwind:long_name = "Surface Northward Air Velocity (10m AGL)" ;
vwind:standard_name = "northward_wind" ;
vwind:units = "m/s" ;
float prmsl(time, ny_grid, nx_grid) ;
prmsl:long_name = "Pressure reduced to MSL" ;
prmsl:standard_name = "air_pressure_at_sea_level" ;
prmsl:units = "Pa" ;
float stmp(time, ny_grid, nx_grid) ;
stmp:long_name = "Surface Air Temperature (2m AGL)" ;
stmp:standard_name = "air_temperature" ;
stmp:units = "K" ;
float spfh(time, ny_grid, nx_grid) ;
spfh:long_name = "Surface Specific Humidity (2m AGL)" ;
spfh:standard_name = "specific_humidity" ;
spfh:units = "1" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

===rad-files===
<code>
netcdf sflux_rad_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float dlwrf(time, ny_grid, nx_grid) ;
dlwrf:long_name = "Downward Long Wave Radiation Flux" ;
dlwrf:standard_name = "surface_downwelling_longwave_flux_in_air" ;
dlwrf:units = "W/m^2" ;
float dswrf(time, ny_grid, nx_grid) ;
dswrf:long_name = "Downward Short Wave Radiation Flux" ;
dswrf:standard_name = "surface_downwelling_shortwave_flux_in_air" ;
dswrf:units = "W/m^2" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>
prc-files
<code>
netcdf sflux_prc_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float prate(time, ny_grid, nx_grid) ;
prate:long_name = "Surface Precipitation Rate" ;
prate:standard_name = "precipitation_flux" ;
prate:units = "kg/m^2/s" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

Atmospheric forcing

2013-01-31T13:25:15Z

Jan: /* air-files: */

==Theory==

[http://ccrm.vims.edu/yinglong/wiki_files/HeatExchange-Theory.pdf Here] is some info about the heat exchange module inside SELFE, based on Dr. Mike Zuluaf's slides.

==File description==
In SELFE, atmoserpic forcings include wind, air pressure and temperature, precipitation, humidity and longwave and shortwave fluxes. These are specified in the netcdf files inside sflux/ dir, and conform to the [http://cf-pcmdi.llnl.gov/documents/cf-conventions/1.0/cf-conventions.html NetCDF Climate and Forecast (CF) Metadata Convention 1.0].

There are 4 types of files in sflux/dir; see [http://www.stccmop.org/CORIE/modeling/selfe/NARR/narr.html this site] for sample files.

<OL>
<LI>sflux_inputs.txt: 
$cat sflux_inputs.txt 
&sflux_inputs ! file name 
start_year = 1998, ! start year 
start_month = 09, ! start month 
start_day = 01, ! start day 
start_hour = 0.0, ! start hour 
utc_start = 0.0, ! offset from UTC in hours, with positive numbers for western hemishpere 
/ 
<LI>sflux_air_1.0*.nc: netcdf files that have time, wind speed, air temperature and pressure and specific humidity;
<LI>sflux_prc_1.0*.nc: netcdf files that have precipitation data;
<LI>sflux_rad_1.0*.nc: netcdf files that have downward long and short (solar) wave radiation fluxes.
</OL>

==Using NARR files for your simulation (North America only)==

First, make sure the [http://ccrm.vims.edu/yinglong/wiki_files/narr.ll NARR grid] covers your hgrid.ll.

In your run directory, mkdir sflux and inside it, create symbolic links to the NARR files. e.g., if you run starts from June 10, 2004 and ends June 20, 2004, then

sflux_air_1.001.nc --> narr_air.2004_06_10.nc 
sflux_air_1.002.nc --> narr_air.2004_06_11.nc 
... 
sflux_air_1.011.nc --> narr_air.2004_06_20.nc 
sflux_air_1.012.nc --> narr_air.2004_06_21.nc (extra day to account for time zone difference) 

Similarly for sflux_rad_*.nc and sflux_prc_*.nc. The number "1" after "air_" denotes first data set used; you can use up to 2 sets in SELFE (which combines them with some given weights set in sflux_subs.F90); here you only use 1 set.

Don't forget to update sflux_inputs.txt also.

==Preparing your own sflux inputs==

After familiarize yourself with the NARR files and their format, you may embark on creating your own nc files. The best way is to modify existing matlab scripts (src/Utility/Sflux_nc/readnc*.m) included in the source code bundle, which have extensively in-line comments to guide you along the way.

'''Conventions for .nc-files in sflux directory and additional files needed for sflux run:'''

'''wind:''' u-component is eastward, v-comp. is northward (normal math convention, not compass convention)

'''windrot_geo2proj.gr3:''' rotates winds in case they do not align with coordinate axes,
i.e. lat/lon

'''watertype.gr3:''' 6 is clear water, 7 is muddiest. Search in selfe_init.F90 for watertype.
===air-files===
<code>
these lines (displayed using the ncdump utility):
netcdf sflux_air_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float uwind(time, ny_grid, nx_grid) ;
uwind:long_name = "Surface Eastward Air Velocity (10m AGL)" ;
uwind:standard_name = "eastward_wind" ;
uwind:units = "m/s" ;
float vwind(time, ny_grid, nx_grid) ;
vwind:long_name = "Surface Northward Air Velocity (10m AGL)" ;
vwind:standard_name = "northward_wind" ;
vwind:units = "m/s" ;
float prmsl(time, ny_grid, nx_grid) ;
prmsl:long_name = "Pressure reduced to MSL" ;
prmsl:standard_name = "air_pressure_at_sea_level" ;
prmsl:units = "Pa" ;
float stmp(time, ny_grid, nx_grid) ;
stmp:long_name = "Surface Air Temperature (2m AGL)" ;
stmp:standard_name = "air_temperature" ;
stmp:units = "K" ;
float spfh(time, ny_grid, nx_grid) ;
spfh:long_name = "Surface Specific Humidity (2m AGL)" ;
spfh:standard_name = "specific_humidity" ;
spfh:units = "1" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

===rad-files===
<code>
netcdf sflux_rad_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float dlwrf(time, ny_grid, nx_grid) ;
dlwrf:long_name = "Downward Long Wave Radiation Flux" ;
dlwrf:standard_name = "surface_downwelling_longwave_flux_in_air" ;
dlwrf:units = "W/m^2" ;
float dswrf(time, ny_grid, nx_grid) ;
dswrf:long_name = "Downward Short Wave Radiation Flux" ;
dswrf:standard_name = "surface_downwelling_shortwave_flux_in_air" ;
dswrf:units = "W/m^2" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

Atmospheric forcing

2013-01-31T13:24:55Z

Jan: /* air-files */

==Theory==

[http://ccrm.vims.edu/yinglong/wiki_files/HeatExchange-Theory.pdf Here] is some info about the heat exchange module inside SELFE, based on Dr. Mike Zuluaf's slides.

==File description==
In SELFE, atmoserpic forcings include wind, air pressure and temperature, precipitation, humidity and longwave and shortwave fluxes. These are specified in the netcdf files inside sflux/ dir, and conform to the [http://cf-pcmdi.llnl.gov/documents/cf-conventions/1.0/cf-conventions.html NetCDF Climate and Forecast (CF) Metadata Convention 1.0].

There are 4 types of files in sflux/dir; see [http://www.stccmop.org/CORIE/modeling/selfe/NARR/narr.html this site] for sample files.

<OL>
<LI>sflux_inputs.txt: 
$cat sflux_inputs.txt 
&sflux_inputs ! file name 
start_year = 1998, ! start year 
start_month = 09, ! start month 
start_day = 01, ! start day 
start_hour = 0.0, ! start hour 
utc_start = 0.0, ! offset from UTC in hours, with positive numbers for western hemishpere 
/ 
<LI>sflux_air_1.0*.nc: netcdf files that have time, wind speed, air temperature and pressure and specific humidity;
<LI>sflux_prc_1.0*.nc: netcdf files that have precipitation data;
<LI>sflux_rad_1.0*.nc: netcdf files that have downward long and short (solar) wave radiation fluxes.
</OL>

==Using NARR files for your simulation (North America only)==

First, make sure the [http://ccrm.vims.edu/yinglong/wiki_files/narr.ll NARR grid] covers your hgrid.ll.

In your run directory, mkdir sflux and inside it, create symbolic links to the NARR files. e.g., if you run starts from June 10, 2004 and ends June 20, 2004, then

sflux_air_1.001.nc --> narr_air.2004_06_10.nc 
sflux_air_1.002.nc --> narr_air.2004_06_11.nc 
... 
sflux_air_1.011.nc --> narr_air.2004_06_20.nc 
sflux_air_1.012.nc --> narr_air.2004_06_21.nc (extra day to account for time zone difference) 

Similarly for sflux_rad_*.nc and sflux_prc_*.nc. The number "1" after "air_" denotes first data set used; you can use up to 2 sets in SELFE (which combines them with some given weights set in sflux_subs.F90); here you only use 1 set.

Don't forget to update sflux_inputs.txt also.

==Preparing your own sflux inputs==

After familiarize yourself with the NARR files and their format, you may embark on creating your own nc files. The best way is to modify existing matlab scripts (src/Utility/Sflux_nc/readnc*.m) included in the source code bundle, which have extensively in-line comments to guide you along the way.

'''Conventions for .nc-files in sflux directory and additional files needed for sflux run:'''

'''wind:''' u-component is eastward, v-comp. is northward (normal math convention, not compass convention)

'''windrot_geo2proj.gr3:''' rotates winds in case they do not align with coordinate axes,
i.e. lat/lon

'''watertype.gr3:''' 6 is clear water, 7 is muddiest. Search in selfe_init.F90 for watertype.
===air-files:===
<code>
these lines (displayed using the ncdump utility):
netcdf sflux_air_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float uwind(time, ny_grid, nx_grid) ;
uwind:long_name = "Surface Eastward Air Velocity (10m AGL)" ;
uwind:standard_name = "eastward_wind" ;
uwind:units = "m/s" ;
float vwind(time, ny_grid, nx_grid) ;
vwind:long_name = "Surface Northward Air Velocity (10m AGL)" ;
vwind:standard_name = "northward_wind" ;
vwind:units = "m/s" ;
float prmsl(time, ny_grid, nx_grid) ;
prmsl:long_name = "Pressure reduced to MSL" ;
prmsl:standard_name = "air_pressure_at_sea_level" ;
prmsl:units = "Pa" ;
float stmp(time, ny_grid, nx_grid) ;
stmp:long_name = "Surface Air Temperature (2m AGL)" ;
stmp:standard_name = "air_temperature" ;
stmp:units = "K" ;
float spfh(time, ny_grid, nx_grid) ;
spfh:long_name = "Surface Specific Humidity (2m AGL)" ;
spfh:standard_name = "specific_humidity" ;
spfh:units = "1" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

===rad-files===
<code>
netcdf sflux_rad_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float dlwrf(time, ny_grid, nx_grid) ;
dlwrf:long_name = "Downward Long Wave Radiation Flux" ;
dlwrf:standard_name = "surface_downwelling_longwave_flux_in_air" ;
dlwrf:units = "W/m^2" ;
float dswrf(time, ny_grid, nx_grid) ;
dswrf:long_name = "Downward Short Wave Radiation Flux" ;
dswrf:standard_name = "surface_downwelling_shortwave_flux_in_air" ;
dswrf:units = "W/m^2" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

Atmospheric forcing

2013-01-31T13:24:33Z

Jan: /* Preparing your own sflux inputs */

==Theory==

[http://ccrm.vims.edu/yinglong/wiki_files/HeatExchange-Theory.pdf Here] is some info about the heat exchange module inside SELFE, based on Dr. Mike Zuluaf's slides.

==File description==
In SELFE, atmoserpic forcings include wind, air pressure and temperature, precipitation, humidity and longwave and shortwave fluxes. These are specified in the netcdf files inside sflux/ dir, and conform to the [http://cf-pcmdi.llnl.gov/documents/cf-conventions/1.0/cf-conventions.html NetCDF Climate and Forecast (CF) Metadata Convention 1.0].

There are 4 types of files in sflux/dir; see [http://www.stccmop.org/CORIE/modeling/selfe/NARR/narr.html this site] for sample files.

<OL>
<LI>sflux_inputs.txt: 
$cat sflux_inputs.txt 
&sflux_inputs ! file name 
start_year = 1998, ! start year 
start_month = 09, ! start month 
start_day = 01, ! start day 
start_hour = 0.0, ! start hour 
utc_start = 0.0, ! offset from UTC in hours, with positive numbers for western hemishpere 
/ 
<LI>sflux_air_1.0*.nc: netcdf files that have time, wind speed, air temperature and pressure and specific humidity;
<LI>sflux_prc_1.0*.nc: netcdf files that have precipitation data;
<LI>sflux_rad_1.0*.nc: netcdf files that have downward long and short (solar) wave radiation fluxes.
</OL>

==Using NARR files for your simulation (North America only)==

First, make sure the [http://ccrm.vims.edu/yinglong/wiki_files/narr.ll NARR grid] covers your hgrid.ll.

In your run directory, mkdir sflux and inside it, create symbolic links to the NARR files. e.g., if you run starts from June 10, 2004 and ends June 20, 2004, then

sflux_air_1.001.nc --> narr_air.2004_06_10.nc 
sflux_air_1.002.nc --> narr_air.2004_06_11.nc 
... 
sflux_air_1.011.nc --> narr_air.2004_06_20.nc 
sflux_air_1.012.nc --> narr_air.2004_06_21.nc (extra day to account for time zone difference) 

Similarly for sflux_rad_*.nc and sflux_prc_*.nc. The number "1" after "air_" denotes first data set used; you can use up to 2 sets in SELFE (which combines them with some given weights set in sflux_subs.F90); here you only use 1 set.

Don't forget to update sflux_inputs.txt also.

==Preparing your own sflux inputs==

After familiarize yourself with the NARR files and their format, you may embark on creating your own nc files. The best way is to modify existing matlab scripts (src/Utility/Sflux_nc/readnc*.m) included in the source code bundle, which have extensively in-line comments to guide you along the way.

'''Conventions for .nc-files in sflux directory and additional files needed for sflux run:'''

'''wind:''' u-component is eastward, v-comp. is northward (normal math convention, not compass convention)

'''windrot_geo2proj.gr3:''' rotates winds in case they do not align with coordinate axes,
i.e. lat/lon

'''watertype.gr3:''' 6 is clear water, 7 is muddiest. Search in selfe_init.F90 for watertype.
===air-files===
<code>
these lines (displayed using the ncdump utility):
netcdf sflux_air_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float uwind(time, ny_grid, nx_grid) ;
uwind:long_name = "Surface Eastward Air Velocity (10m AGL)" ;
uwind:standard_name = "eastward_wind" ;
uwind:units = "m/s" ;
float vwind(time, ny_grid, nx_grid) ;
vwind:long_name = "Surface Northward Air Velocity (10m AGL)" ;
vwind:standard_name = "northward_wind" ;
vwind:units = "m/s" ;
float prmsl(time, ny_grid, nx_grid) ;
prmsl:long_name = "Pressure reduced to MSL" ;
prmsl:standard_name = "air_pressure_at_sea_level" ;
prmsl:units = "Pa" ;
float stmp(time, ny_grid, nx_grid) ;
stmp:long_name = "Surface Air Temperature (2m AGL)" ;
stmp:standard_name = "air_temperature" ;
stmp:units = "K" ;
float spfh(time, ny_grid, nx_grid) ;
spfh:long_name = "Surface Specific Humidity (2m AGL)" ;
spfh:standard_name = "specific_humidity" ;
spfh:units = "1" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>
===rad-files===
<code>
netcdf sflux_rad_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float dlwrf(time, ny_grid, nx_grid) ;
dlwrf:long_name = "Downward Long Wave Radiation Flux" ;
dlwrf:standard_name = "surface_downwelling_longwave_flux_in_air" ;
dlwrf:units = "W/m^2" ;
float dswrf(time, ny_grid, nx_grid) ;
dswrf:long_name = "Downward Short Wave Radiation Flux" ;
dswrf:standard_name = "surface_downwelling_shortwave_flux_in_air" ;
dswrf:units = "W/m^2" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

Atmospheric forcing

2013-01-31T13:23:42Z

Jan: /* Preparing your own sflux inputs */

==Theory==

[http://ccrm.vims.edu/yinglong/wiki_files/HeatExchange-Theory.pdf Here] is some info about the heat exchange module inside SELFE, based on Dr. Mike Zuluaf's slides.

==File description==
In SELFE, atmoserpic forcings include wind, air pressure and temperature, precipitation, humidity and longwave and shortwave fluxes. These are specified in the netcdf files inside sflux/ dir, and conform to the [http://cf-pcmdi.llnl.gov/documents/cf-conventions/1.0/cf-conventions.html NetCDF Climate and Forecast (CF) Metadata Convention 1.0].

There are 4 types of files in sflux/dir; see [http://www.stccmop.org/CORIE/modeling/selfe/NARR/narr.html this site] for sample files.

<OL>
<LI>sflux_inputs.txt: 
$cat sflux_inputs.txt 
&sflux_inputs ! file name 
start_year = 1998, ! start year 
start_month = 09, ! start month 
start_day = 01, ! start day 
start_hour = 0.0, ! start hour 
utc_start = 0.0, ! offset from UTC in hours, with positive numbers for western hemishpere 
/ 
<LI>sflux_air_1.0*.nc: netcdf files that have time, wind speed, air temperature and pressure and specific humidity;
<LI>sflux_prc_1.0*.nc: netcdf files that have precipitation data;
<LI>sflux_rad_1.0*.nc: netcdf files that have downward long and short (solar) wave radiation fluxes.
</OL>

==Using NARR files for your simulation (North America only)==

First, make sure the [http://ccrm.vims.edu/yinglong/wiki_files/narr.ll NARR grid] covers your hgrid.ll.

In your run directory, mkdir sflux and inside it, create symbolic links to the NARR files. e.g., if you run starts from June 10, 2004 and ends June 20, 2004, then

sflux_air_1.001.nc --> narr_air.2004_06_10.nc 
sflux_air_1.002.nc --> narr_air.2004_06_11.nc 
... 
sflux_air_1.011.nc --> narr_air.2004_06_20.nc 
sflux_air_1.012.nc --> narr_air.2004_06_21.nc (extra day to account for time zone difference) 

Similarly for sflux_rad_*.nc and sflux_prc_*.nc. The number "1" after "air_" denotes first data set used; you can use up to 2 sets in SELFE (which combines them with some given weights set in sflux_subs.F90); here you only use 1 set.

Don't forget to update sflux_inputs.txt also.

==Preparing your own sflux inputs==

After familiarize yourself with the NARR files and their format, you may embark on creating your own nc files. The best way is to modify existing matlab scripts (src/Utility/Sflux_nc/readnc*.m) included in the source code bundle, which have extensively in-line comments to guide you along the way.

'''Conventions for .nc-files in sflux directory and additional files needed for sflux run:'''

'''wind:''' u-component is eastward, v-comp. is northward (normal math convention, not compass convention)

'''windrot_geo2proj.gr3:''' rotates winds in case they do not align with coordinate axes,
i.e. lat/lon

'''watertype.gr3:''' 6 is clear water, 7 is muddiest. Search in selfe_init.F90 for watertype.
<code>
these lines (displayed using the ncdump utility):
netcdf sflux_air_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float uwind(time, ny_grid, nx_grid) ;
uwind:long_name = "Surface Eastward Air Velocity (10m AGL)" ;
uwind:standard_name = "eastward_wind" ;
uwind:units = "m/s" ;
float vwind(time, ny_grid, nx_grid) ;
vwind:long_name = "Surface Northward Air Velocity (10m AGL)" ;
vwind:standard_name = "northward_wind" ;
vwind:units = "m/s" ;
float prmsl(time, ny_grid, nx_grid) ;
prmsl:long_name = "Pressure reduced to MSL" ;
prmsl:standard_name = "air_pressure_at_sea_level" ;
prmsl:units = "Pa" ;
float stmp(time, ny_grid, nx_grid) ;
stmp:long_name = "Surface Air Temperature (2m AGL)" ;
stmp:standard_name = "air_temperature" ;
stmp:units = "K" ;
float spfh(time, ny_grid, nx_grid) ;
spfh:long_name = "Surface Specific Humidity (2m AGL)" ;
spfh:standard_name = "specific_humidity" ;
spfh:units = "1" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

<code>
netcdf sflux_rad_1.001 {
dimensions:
nx_grid = 349 ;
ny_grid = 277 ;
time = UNLIMITED ; // (8 currently)
variables:
float time(time) ;
time:long_name = "Time" ;
time:standard_name = "time" ;
time:units = "days since 2001-01-01" ;
time:base_date = 2001, 1, 1, 0 ;
float lon(ny_grid, nx_grid) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(ny_grid, nx_grid) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float dlwrf(time, ny_grid, nx_grid) ;
dlwrf:long_name = "Downward Long Wave Radiation Flux" ;
dlwrf:standard_name = "surface_downwelling_longwave_flux_in_air" ;
dlwrf:units = "W/m^2" ;
float dswrf(time, ny_grid, nx_grid) ;
dswrf:long_name = "Downward Short Wave Radiation Flux" ;
dswrf:standard_name = "surface_downwelling_shortwave_flux_in_air" ;
dswrf:units = "W/m^2" ;

// global attributes:
:Conventions = "CF-1.0" ;
}
</code>

Atmospheric forcing

2013-01-30T09:24:19Z

Jan: /* Preparing your own sflux inputs */

Atmospheric forcing

2013-01-30T09:23:00Z

Jan: /* Preparing your own sflux inputs */

Atmospheric forcing

2013-01-30T09:21:48Z

Jan: /* Conventions for .nc-files in sflux directory: */

Atmospheric forcing

2013-01-30T09:15:22Z

Jan: /* Preparing your own sflux inputs */