New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
README.md in NEMO/branches/2020/r12377_ticket2386/tests/STATION_ASF – NEMO

source: NEMO/branches/2020/r12377_ticket2386/tests/STATION_ASF/README.md @ 13540

Last change on this file since 13540 was 13540, checked in by andmirek, 3 years ago

Ticket #2386: update to latest trunk

File size: 5.1 KB
Line 
1# *Station Air-Sea Fluxes* demonstration case
2
3Last successful test done with NEMOGCM trunk: `r13263`
4
5Author: Laurent Brodeau, 2020
6
7NOTE: if working with the trunk of NEMO, you are strongly advised to use the same test-case but on the `NEMO-examples` GitHub depo:
8https://github.com/NEMO-ocean/NEMO-examples/tree/master/STATION_ASF
9
10## Objectives
11
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.
13
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).
15
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.
17
18As input ```STATION_ASF``` will require the traditional *bulk* sea surface parameters:
19
20- Bulk sea surface temperature (SST) at _z<sub>SST</sub>_ meters below the surface
21- Surface current vector
22- Sea surface salinity
23
24as well as the usual surface atmospheric state:
25
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
29- Sea level atmospheric pressure (SLP)
30- Downwelling solar radiation
31- Downwelling longwave radiation
32
33### Example of diagnostics from `STATION_ASF`
34
35(Generated with script `./EXPREF/plot_station_asf_simple.py`)
36
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
46## Physical description
47
48### Important namelist parameters specific to STATION_ASF
49
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!
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
62One full year (2018) of processed hourly data from the PAPA station (buoy) is found into the `input_data` directory.
63These three files are everything you need to play with the set of *namelists* provided for this test-case.
64
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
68
69For 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
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.```
79
80
81
82## Playing with STATION_ASF
83
84First compile the test-case as follows (compile with xios-2.5 support → check your ARCH file):
85
86```./makenemo -a STATION_ASF -m <your_arch> -n STATION_ASF2 -j 4```
87
88Then 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:
89
90- ```NEMO_ROOT_DIR``` : NEMO root directory where to fetch compiled STATION_ASF ```nemo.exe``` + setup (such as ```${NEMO_ROOT_DIR}/tests/STATION_ASF```)
91
92- ```PROD_DIR``` :  Directory where to run the simulation
93
94- ```DATA_IN_DIR``` : Directory containing sea-surface + atmospheric forcings (found here in ```input_data/```)
95
96If 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
Note: See TracBrowser for help on using the repository browser.