Generate Forcing Files

1. Grab the entire region’s stations (https://coastwatch.glerl.noaa.gov/marobs/stations.lmd) or for a specific lake “stations-superior.lmd”. And save it to the target path.

2. In the directory (target path) where you want to create the forcing, edit the file “metedit.in” and update the .lmd file and the glsea-temps.dat files on lines 2, 3 and 4. Also edit the year in lines 5 and 6, and then change line 8 to the appropriate # of hours (8760 in normal years, 8784 in leap years).
   metedit.in shows like below:

                       # leave first line blank
stations_superior.lmd
2008.lmd
glsea-temps2008_1024.dat
e200800100
2008 001 00    IYR,IJD,IHR
1.             DTWIND
8784           NHRS
2              SELECTMOD

• The *lmd file for line 3 can be downloaded here: https://coastwatch.glerl.noaa.gov/AnnualLMD/
• The glsea-temps files for line 4 can be downloaded here: https://coastwatch.glerl.noaa.gov/ftp/glsea/avgtemps/1995/glsea-temps1995_1024.dat (Note: remove all the lines of headers)

wget https://coastwatch.glerl.noaa.gov/ftp/glsea/avgtemps/1995/glsea-temps1995_1024.dat

3. Compile metedit-crmod.f

module load intel/2016.1
ifort -o metedit-crmod.exe metedit-crmod.f

./metedit-crmod.exe <metedit.in >metedit.out

• This program will create both an output file (metedit.out) and an adjusted meteorology file (metobs.dat)

5. edit “interp.in” to update the year (Line 5), and the # of hours on Line 7 (8761 for normal years, 8785 for leap years)

superior_grd.dat
metobs.dat

e200800100
2008 001 00    IYR,IJD,IHR
1.             DTWIND
8785           NHRS

6. Compile interpun.f

module load intel/2016.1
icc -c stack.c stackpair.c volume.c
ifort -O -o interpun-new.exe interpun.f stack.o stackpair.o volume.o

7. qsub interpun.sh
8. rename *_uvw.dat to *_wnd.dat
9. Edit your “*_forcing.nml” file (superior_forcing.nml in example) with appropriate dates, file locations, and file names (note: this assumes wind speed is being used instead of stress)

&NML_SURFACE_FORCE
 INPUT_DIR       = './',
 OUTPUT_DIR      = './',
 INPUT_FILE             = 'e200800100',
 INPUT_FILE_TYPE        = 'binary speed',
 OUTPUT_FILE            = 'e200800100_forcing_marobs.nc',
 OUTPUT_FILE_TYPE       = 'unstructured',
 GRID_FILE              = 'superior_grd.dat',
 GRID_FILE_UNITS        = 'degrees',
 PROJECTION_REFERENCE   = 'none',
 DATE_FORMAT            = 'YMD',
 START_DATE             = '2008-01-01 00:00:00',
 TIMEZONE               = 'UTC',
 END_DATE               = '2009-01-01 00:00:00',
 OUT_INTERVAL           = 'seconds=3600.0',
 ZERO_TIME              = '2008-01-01 00:00:00',
 DATE_REFERENCE = '1858-11-17 00:00:00',
/

10. Make sure there is “xsurfaceforce” executable file in the present directory.
    To generate “xsurfaceforce”, first compile the fvcom code (/pexue5/chuyan/generate_forcing/forcing/code_xsurfaceforce), then compile (/pexue5/chuyan/generate_forcing/forcing/code_xsurfaceforce/input/surface_forcing).
11. qsub xsurfaceforce.sh
12. The new generated nc file is the target forcing input file.