Changeset 13682

Timestamp:

2020-10-27T17:12:39+01:00 (3 years ago)

Author:

laurent

Message:

Still improving the oce+ice capable STATION_ASF...

Location:

NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF

Files:

• EXPREF/analyze_output.py (modified) (2 diffs)
• EXPREF/context_nemo.xml (deleted)
• EXPREF/domain_def_nemo.xml (deleted)
• EXPREF/field_def_nemo-oce.xml (deleted)
• EXPREF/grid_def_nemo.xml (deleted)
• EXPREF/launch_sasf.sh (modified) (8 diffs)
• EXPREF/namelist_cfg (deleted)
• EXPREF/namelist_ice_cfg (added)
• EXPREF/namelist_ref (deleted)
• EXPREF/sbcblk_sanity_check.sh (modified) (4 diffs)
• README.md (modified) (1 diff)

NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/analyze_output.py

@@ -17 +17 @@
 import math
 import numpy as nmp
-from netCDF4 import Dataset,num2date
-import matplotlib as mpl
-mpl.use('Agg')
-import matplotlib.pyplot as plt
-import matplotlib.dates as mdates
+from netCDF4 import Dataset
 
 l_t_shift = False  ; # because time interp. is set to FALSE into "&namsbc_blk" of NEMO...
@@ -59 +55 @@
 if narg==4:
     l_more = ( sys.argv[3] in ['m','M'] )
+    import matplotlib as mpl
+    mpl.use('Agg')
+    import matplotlib.pyplot as plt
+
 
 

NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/launch_sasf.sh

@@ -15 +15 @@
 
 # Root directory NEMOGCM reference depository where to fetch compiled STATION_ASF nemo.exe + default namelists:
-NEMO_REF_DIR="/home/laurent/DEV/NEMO/NEMOGCM_dev_r13648_ASINTER-04_laurent_bulk_ice"
+NEMO_REF_DIR="`dirname ${PWD} | sed -e 's|/tests/STATION_ASF||g'`" ; # that should normally do the trick!
 
 # NEMOGCM root directory where to fetch compiled STATION_ASF nemo.exe:
@@ -32 +32 @@
 i_sea_ice=0 ; # default does not compute fluxes over sea-ice
 #
-#NFORC="ERA5_arctic" ; i_sea_ice=1 ; FORC="ERA5_arctic_surface_81N_36p75E_1h" ; # string sufficient to copy relevant files as in "*${FORC}*.nc"
-NFORC="PAPA" ; FORC="Station_PAPA_50N-145W" ; # string sufficient to copy relevant files as in "*${FORC}*.nc"
+NFORC="ERA5_arctic" ; i_sea_ice=1 ; FORC="ERA5_arctic_surface_81N_36p75E_1h" ; # string sufficient to copy relevant files as in "*${FORC}*.nc"
+#NFORC="PAPA" ; FORC="Station_PAPA_50N-145W" ; # string sufficient to copy relevant files as in "*${FORC}*.nc"
 
 
@@ -57 +57 @@
 i_si3=0
 FCPP="${SASF_WRK_DIR}/cpp_${CONFIG_BLD}.fcm"
-if [ ! -f ${FCPP} ]; then echo " Mhhh, we did not find 'cpp_STATION_ASF.fcm' into `dirname ${FCPP}` !"; exit; fi
+if [ ! -f ${FCPP} ]; then echo " Mhhh, we could not find 'cpp_STATION_ASF.fcm' into `dirname ${FCPP}` !"; exit; fi
 ca=`cat ${FCPP} | grep 'key_si3'`
 
@@ -78 +78 @@
     LIST_OA_ALGOS="ECMWF" ;          # list of air-sea algorithms to test
     LIST_IA_ALGOS="LG15 LU12 CSTC" ; # list of air-ice algorithms to test
-    DIR_NL=${NFORC}/ice ; # where to fetch the namelists from...
+    DIR_NL=${NFORC}/oce+ice ; # where to fetch the namelists from...
 else
     # Only OPEN-OCEAN/AIR flux computation
@@ -96 +96 @@
 
 
-if [ ! -d ${DATA_IN_DIR} ]; then echo "PROBLEM!!! => did not find directory 'input_data' with input forcing..."; exit; fi
+CFGS_SHARED="${NEMO_REF_DIR}/cfgs/SHARED"
+if [ ! -d ${CFGS_SHARED} ]; then echo "PROBLEM!!! => could not find directory ${CFGS_SHARED} !"; exit; fi
+
+if [ ! -d ${DATA_IN_DIR} ]; then echo "PROBLEM!!! => could not find directory 'input_data' with input forcing..."; exit; fi
 
 cdt_cmpl="`\ls -l ${NEMO_EXE} | cut -d' ' -f 6,7,8`"
@@ -125 +128 @@
 ################
 
-list_xml="iodef.xml file_def_nemo-oce.xml field_def_nemo-oce.xml"
+list_xml_ref="field_def_nemo-oce.xml domain_def_nemo.xml grid_def_nemo.xml"
+list_xml_cfg="iodef.xml file_def_nemo-oce.xml"
 fcntxt="context_nemo_OCE.xml"
 if [ ${i_sea_ice} -eq 1 ]; then
-    list_xml+=" file_def_nemo-ice.xml field_def_nemo-ice.xml"
+    list_xml_ref+=" field_def_nemo-ice.xml"
+    list_xml_cfg+=" file_def_nemo-ice.xml"
     fcntxt="context_nemo_OCE+ICE.xml"
 fi
@@ -137 +142 @@
 
 # All remaining "*.xml" files:
-for ff in ${list_xml} ; do
+for ff in ${list_xml_cfg} ; do
     if [ ! -f ${SASF_EXPREF}/${ff} ]; then echo " Mhhh, ${ff} not found into ${SASF_EXPREF} !"; exit; fi
     rsync -avPL ${SASF_EXPREF}/${ff} ${PROD_DIR}/
 done
-
-
-# Getting reference/defaults files from reference NEMO distro:
-rdir="${NEMO_REF_DIR}/cfgs/SHARED"
-for ff in "domain_def_nemo.xml" "grid_def_nemo.xml"; do
-    if [ ! -f ${rdir}/${ff} ]; then echo " Mhhh, ${ff} not found into ${rdir} !"; exit; fi
-    ln -sf ${rdir}/${ff} ${PROD_DIR}/.
+for ff in ${list_xml_ref} ; do
+    if [ ! -f ${CFGS_SHARED}/${ff} ]; then echo " Mhhh, ${ff} not found into ${CFGS_SHARED} !"; exit; fi
+    rsync -avPL ${CFGS_SHARED}/${ff} ${PROD_DIR}/
 done
 
@@ -186 +187 @@
         # The namelists:
         rsync -avPL ${SASF_EXPREF}/${DIR_NL}/namelist_${scase}_cfg ${PROD_DIR}/namelist_cfg
-        rsync -avPL ${SASF_EXPREF}/namelist_ref                   ${PROD_DIR}/namelist_ref
+        rsync -avPL ${CFGS_SHARED}/namelist_ref                    ${PROD_DIR}/namelist_ref
         if [ ${i_sea_ice} -eq 1 ]; then
-            rsync -avPL ${SASF_EXPREF}/namelist_ice_ref ${PROD_DIR}/namelist_ice_ref
             rsync -avPL ${SASF_EXPREF}/namelist_ice_cfg ${PROD_DIR}/namelist_ice_cfg
+            rsync -avPL ${CFGS_SHARED}/namelist_ice_ref ${PROD_DIR}/namelist_ice_ref
         fi
         

NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/sbcblk_sanity_check.sh

@@ -25 +25 @@
 
 # Root directory NEMOGCM reference depository where to fetch compiled STATION_ASF nemo.exe + default namelists:
-NEMO_REF_DIR="/home/laurent/DEV/NEMO/NEMOGCM_dev_r13648_ASINTER-04_laurent_bulk_ice"
+NEMO_REF_DIR="`dirname ${PWD} | sed -e 's|/tests/STATION_ASF||g'`" ; # that should normally do the trick!
 
 # NEMOGCM root directory where to fetch compiled STATION_ASF nemo.exe:
@@ -41 +41 @@
 #############################################################################################
 i_sea_ice=0 ; # default does not compute fluxes over sea-ice
-NFORC="IDEALIZED" ; i_sea_ice=0 ; FORC="input_output_VALIDATION_IDEALIZED.nc"
+NFORC="IDEALIZED" ; i_sea_ice=0 ; FORC="input_output_VALIDATION_IDEALIZED.nc" ; # for now the sanity test does not test ice-atmo bulk parameterizations
 
 
@@ -47 +47 @@
 
 #================================================================================
+
+# Should the `analyze_output.py` script provide more output ???
+cmore=""
+if [ "$1" != "" ]; then
+    if [ "$1" = "more" ]; then
+        cmore="m"
+    else
+        echo "Only 'more' is accepted as argument for script `basename $0` !"
+        exit
+    fi
+fi
+
+
 
 rm -f SBCBLK.success SBCBLK.fail
@@ -191 +204 @@
 
 cd ${HERE}
-python3 ./analyze_output.py ${DATA_IN_DIR}/${FORC} ${PROD_DIR}/output
-
-
+python3 ./analyze_output.py ${DATA_IN_DIR}/${FORC} ${PROD_DIR}/output ${cmore}
+
+

NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/README.md

@@ -80 +80 @@
 
 
+## Testing the sanity of the SBCBLK (bulk atmospheric forcing) interface of NEMO
+
+`STATION_ASF` can be used to perform a sanity test of the SBCBLK interface of
+NEMO.  It will test all the bulk-parameterization algorithms using an idealized
+forcing that includes a wide range of *SSX / surface atmospheric state*
+conditions to detect potential error / inconsitencies.  Both a short report and
+boolean output: *passed* or *failed* is provided as an output.
+
+First compile the `STATION_ASF` test-case as follows (compile with xios-2.5 support → check your ARCH file):
+    
+```./makenemo -a STATION_ASF -m <your_arch> -n STATION_ASF2 -j 4```
+
+(successfull compilation generates ```tests/STATION_ASF2/BLD/bin/nemo.exe```
+
+Move to: ```tests/STATION_ASF/EXPREF/```
+
+There, in script ``sbcblk_sanity_check.sh``, double check that the variables given a value in the first lines are consistent with your setup (it should).
+
+Launch the script:
+
+    $ ./sbcblk_sanity_check.sh
+
+
+The report is both interactively shown in the standard output and spwaned in the report file:
+
+- `SBCBLK.success` the test passed :)
+- `SBCBLK.fail` the test failed :(
+
+In case of failure, read the rapport to know which algorithm/computed variables did not pass the test... You should also re-run the script ``sbcblk_sanity_check.sh`` with the "more" argument: 
+
+    $ ./sbcblk_sanity_check.sh more
+
+ This will generate more output such as figure of time series for the tested variables.
+
+### Dependencies
+
+In order for the python script `analyze_output.py` to work, you need `Python 3` with the following modules:
+- numpy
+- netCDF4
+- (matplotlib)
+
+
+
+
+
+
 ## Playing with STATION_ASF
 
 First compile the test-case as follows (compile with xios-2.5 support → check your ARCH file):
-
+    
 ```./makenemo -a STATION_ASF -m <your_arch> -n STATION_ASF2 -j 4```
 
-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:
+Move to: ```tests/STATION_ASF/EXPREF/```
+
+
+
+
+to launch 3 simulations (one for each bulk parameterization available). You need to adapt the following variable to your environment in the script:
 
 - ```NEMO_ROOT_DIR``` : NEMO root directory where to fetch compiled STATION_ASF ```nemo.exe``` + setup (such as ```${NEMO_ROOT_DIR}/tests/STATION_ASF```)