[13540] | 1 | # *Station Air-Sea Fluxes* demonstration case |
---|
[11850] | 2 | |
---|
[13540] | 3 | Last successful test done with NEMOGCM trunk: `r13263` |
---|
[11850] | 4 | |
---|
[13540] | 5 | Author: Laurent Brodeau, 2020 |
---|
[11850] | 6 | |
---|
[13540] | 7 | NOTE: if working with the trunk of NEMO, you are strongly advised to use the same test-case but on the `NEMO-examples` GitHub depo: |
---|
| 8 | https://github.com/NEMO-ocean/NEMO-examples/tree/master/STATION_ASF |
---|
| 9 | |
---|
[11850] | 10 | ## Objectives |
---|
| 11 | |
---|
[13540] | 12 | ```STATION_ASF``` is a demonstration test-case that mimics a (static) in-situ station (buoy, platform) dedicated to the estimation of surface air-sea fluxes by means of *widely-measured* (bulk) meteorological surface parameters. |
---|
[11850] | 13 | |
---|
[13540] | 14 | ```STATION_ASF``` has been constructed by merging the *single column* and the *standalone surface module* configurations of NEMO. In short, it can be defined as "SAS meets C1D". As such, the spatial domain of ```STATION_ASF``` is punctual (1D, well actually 3 x 3 as in C1D). |
---|
[11850] | 15 | |
---|
[13540] | 16 | ```STATION_ASF``` is therefore a versatile tool, and extremely lightweight in terms of computing requirements, to test the different bulk algorithms and cool-skin/warm-layer parameterization options included in NEMO. |
---|
[11850] | 17 | |
---|
[11892] | 18 | As input ```STATION_ASF``` will require the traditional *bulk* sea surface parameters: |
---|
[11850] | 19 | |
---|
[13540] | 20 | - Bulk sea surface temperature (SST) at _z<sub>SST</sub>_ meters below the surface |
---|
[11850] | 21 | - Surface current vector |
---|
| 22 | - Sea surface salinity |
---|
| 23 | |
---|
| 24 | as well as the usual surface atmospheric state: |
---|
| 25 | |
---|
[13540] | 26 | - air temperature at _z<sub>t</sub>_ meters above the surface |
---|
| 27 | - air humidity at _z<sub>t</sub>_ meters above the surface (specific humidity or relative humidity or dew-point temperature) |
---|
| 28 | - wind speed vector at _z<sub>u</sub>_ meters above the surface |
---|
[11850] | 29 | - Sea level atmospheric pressure (SLP) |
---|
| 30 | - Downwelling solar radiation |
---|
| 31 | - Downwelling longwave radiation |
---|
| 32 | |
---|
[13540] | 33 | ### Example of diagnostics from `STATION_ASF` |
---|
[11850] | 34 | |
---|
[13540] | 35 | (Generated with script `./EXPREF/plot_station_asf_simple.py`) |
---|
[11850] | 36 | |
---|
[13540] | 37 | ![plot](https://github.com/NEMO-ocean/NEMO-examples/blob/master/STATION_ASF/figs/01_temperatures_ECMWF.svg) |
---|
| 38 | |
---|
| 39 | ![plot](https://github.com/NEMO-ocean/NEMO-examples/blob/master/STATION_ASF/figs/Cd.svg) |
---|
| 40 | |
---|
| 41 | ![plot](https://github.com/NEMO-ocean/NEMO-examples/blob/master/STATION_ASF/figs/dT_skin.svg) |
---|
| 42 | |
---|
| 43 | ![plot](https://github.com/NEMO-ocean/NEMO-examples/blob/master/STATION_ASF/figs/Qlat.svg) |
---|
| 44 | |
---|
| 45 | |
---|
[11850] | 46 | ## Physical description |
---|
| 47 | |
---|
[13540] | 48 | ### Important namelist parameters specific to STATION_ASF |
---|
[11850] | 49 | |
---|
[13540] | 50 | * ```rn_dept1@namusr_def:``` depth (m) at which the prescribed SST is taken (*i.e.* depth of first T-point); important due to impact on warm-layer estimate, the deeper, the more pronounced! |
---|
[11850] | 51 | |
---|
| 52 | * ```rn_lat1d,rn_lon1d@namc1d:``` fixed coordinates of the location of the station (buoy, platform, etc). |
---|
| 53 | |
---|
| 54 | * ```namsbc_blk:``` to be filled carefully, just as for "C1D", the prescribed surface ATMOSPHERIC state (files) are time series of shape 3x3 in space |
---|
| 55 | |
---|
| 56 | * ```namsbc_sas:``` to be filled carefully, just as for "C1D", the prescribed surface OCEAN state (files) are time series of shape 3x3 in space |
---|
| 57 | |
---|
| 58 | |
---|
| 59 | |
---|
| 60 | ## Input files to test STATION ASF |
---|
| 61 | |
---|
[13540] | 62 | One full year (2018) of processed hourly data from the PAPA station (buoy) is found into the `input_data` directory. |
---|
| 63 | These three files are everything you need to play with the set of *namelists* provided for this test-case. |
---|
[11850] | 64 | |
---|
[13540] | 65 | - ```Station_PAPA_50N-145W_atm_hourly_y2018.nc``` → contains hourly surface atmospheric state |
---|
| 66 | - ```Station_PAPA_50N-145W_precip_daily_y2018.nc``` → contains daily precipitation |
---|
| 67 | - ```Station_PAPA_50N-145W_oce_hourly_y2018.nc``` → contains hourly sea surface state |
---|
[11850] | 68 | |
---|
[11892] | 69 | For station PAPA (50.1 N, 144.9 W), air temperature and humidity are measured at 2.5 m, the wind speed at 4 m, and the SST at 1 m below the surface, hence the following namelist parameters are given: |
---|
| 70 | |
---|
[13540] | 71 | - `&namusr_def` |
---|
| 72 | - ```rn_dept1 = 1. ``` |
---|
| 73 | - `&namc1d` |
---|
| 74 | - ```rn_lat1d = 50.1 ``` |
---|
| 75 | - ```rn_lon1d = 215.1``` |
---|
| 76 | - `&namsbc_blk` |
---|
| 77 | - ```rn_zqt = 2.5``` |
---|
| 78 | - ```rn_zu = 4.``` |
---|
[11892] | 79 | |
---|
| 80 | |
---|
| 81 | |
---|
[11930] | 82 | ## Playing with STATION_ASF |
---|
[11892] | 83 | |
---|
[11930] | 84 | First compile the test-case as follows (compile with xios-2.5 support → check your ARCH file): |
---|
[11892] | 85 | |
---|
[13540] | 86 | ```./makenemo -a STATION_ASF -m <your_arch> -n STATION_ASF2 -j 4``` |
---|
[11892] | 87 | |
---|
[11930] | 88 | Then you can use the script ``launch_sasf.sh`` found in ```EXPREF/``` to launch 3 simulations (one for each bulk parameterization available). You need to adapt the following variable to your environment in the script: |
---|
[11892] | 89 | |
---|
[13540] | 90 | - ```NEMO_ROOT_DIR``` : NEMO root directory where to fetch compiled STATION_ASF ```nemo.exe``` + setup (such as ```${NEMO_ROOT_DIR}/tests/STATION_ASF```) |
---|
[11892] | 91 | |
---|
[13540] | 92 | - ```PROD_DIR``` : Directory where to run the simulation |
---|
[11892] | 93 | |
---|
[13540] | 94 | - ```DATA_IN_DIR``` : Directory containing sea-surface + atmospheric forcings (found here in ```input_data/```) |
---|
[11930] | 95 | |
---|
[13540] | 96 | If everything goes according to plan, ``launch_sasf.sh`` should have generated the 3 following sets of output files into `${PROD_DIR}/output`: |
---|
| 97 | |
---|
| 98 | STATION_ASF-COARE3p6_1h_20180101_20181231_gridT.nc |
---|
| 99 | STATION_ASF-COARE3p6_1h_20180101_20181231_gridU.nc |
---|
| 100 | STATION_ASF-COARE3p6_1h_20180101_20181231_gridV.nc |
---|
| 101 | STATION_ASF-ECMWF_1h_20180101_20181231_gridT.nc |
---|
| 102 | STATION_ASF-ECMWF_1h_20180101_20181231_gridU.nc |
---|
| 103 | STATION_ASF-ECMWF_1h_20180101_20181231_gridV.nc |
---|
| 104 | STATION_ASF-NCAR_1h_20180101_20181231_gridT.nc |
---|
| 105 | STATION_ASF-NCAR_1h_20180101_20181231_gridU.nc |
---|
| 106 | STATION_ASF-NCAR_1h_20180101_20181231_gridV.nc |
---|
| 107 | |
---|
| 108 | --- |
---|
| 109 | |
---|
| 110 | */Laurent, July 2020.* |
---|
| 111 | |
---|