Skip to content
Snippets Groups Projects

Compare revisions

Changes are shown as if the source revision was being merged into the target revision. Learn more about comparing revisions.

Source

Select target project
No results found

Target

Select target project
  • nemo/nemo
  • sparonuz/nemo
  • hatfield/nemo
  • extdevs/nemo
4 results
Show changes
Hide whitespace changes
Inline Side-by-side
Showing
with 186 additions and 149 deletions
cfgs/C1D/EXP_SASICE/grid_def_nemo.xml 0 → 120000
View file @ 43f40c58
../../SHARED/grid_def_nemo.xml
\ No newline at end of file
tests/STATION_ASF/EXPREF/ERA5/oce+ice/namelist_ecmwf-cstc_cfg → cfgs/C1D/EXP_SASICE/namelist_cfg
View file @ 43f40c58
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! NEMO/OCE Configuration namelist : overwrite default values defined in SHARED/namelist_ref
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! STATION_ASF configuration with ICE, <FORCING> forcing
!! C1D - SAS - ICE configuration !!
!!======================================================================
!! *** Domain & Run management namelists *** !!
!! !!
......@@ -10,22 +10,19 @@
!! namcfg parameters of the configuration (default: user defined GYRE)
!! namwad Wetting and drying (default: OFF)
!! namtsd data: temperature & salinity (default: OFF)
!! namcrs coarsened grid (for outputs and/or TOP) (ln_crs =T)
!! namc1d 1D configuration options (ln_c1d =T)
!! namc1d_dyndmp 1D newtonian damping applied on currents (ln_c1d =T)
!! namc1d_uvd 1D data (currents) (ln_c1d =T)
!!======================================================================
!
!-----------------------------------------------------------------------
&namusr_def ! STATION_ASF user defined namelist
&namusr_def ! C1D user defined namelist
!-----------------------------------------------------------------------
rn_dept1 = 1.0 ! depth (m) at which the prescribed SST is taken == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!)
! ! [ERA5_NorthGreenland station = 1 m]
rn_bathy = 1. ! depth in meters
! rn_dept1 = 1.0 ! depth (m) at which the prescribed SST is taken == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!)
/
!-----------------------------------------------------------------------
&namrun ! parameters of the run
!-----------------------------------------------------------------------
cn_exp = 'STATION_ASF-ECMWF-CSTC_<FORCING>' ! experience name
cn_exp = 'C1D_SASICE' ! experience name
nn_it000 = 1 ! first time step
nn_itend = 8760 ! last time step => 3 years (including 1 leap!) at dt=3600s
nn_date0 = 20180101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
......@@ -33,37 +30,30 @@
!
nn_stock = 8760 ! 1 year at dt=3600s / frequency of creation of a restart file (modulo referenced to 1)
nn_write = 8760 ! 1 year at dt=3600s / frequency of write in the output file (modulo referenced to nn_it000)
ln_top = .false. ! Consider (T) or bypass (F) the TOP routines when the key_top is activated
/
!-----------------------------------------------------------------------
&namdom ! time and space domain / STATION_ASF
&namdom ! time and space domain
!-----------------------------------------------------------------------
ln_linssh = .true. ! =T linear free surface ==>> model level are fixed in time
rn_Dt = 3600. ! time step for the dynamics and tracer => frequency of atmospheric forcing
ln_meshmask = .false.
!
ln_c1d = .true. ! Single column domain (1x1pt) (T => fill namc1d)
/
!-----------------------------------------------------------------------
&namcfg ! parameters of the configuration (default: use namusr_def in namelist_cfg)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namc1d ! 1D configuration options / STATION_ASF
&namc1d ! 1D configuration options (ln_c1d =T default: PAPA station)
!-----------------------------------------------------------------------
rn_lat1d = 84. ! Column latitude [ERA5_NorthGreenland station]
rn_lon1d = 324. ! Column longitude [ERA5_NorthGreenland station]
/
!-----------------------------------------------------------------------
&namc1d_dyndmp ! U & V newtonian damping (ln_c1d =T default: OFF)
!-----------------------------------------------------------------------
ln_dyndmp = .false. ! add a damping term (T) or not (F)
/
!-----------------------------------------------------------------------
&namc1d_uvd ! data: U & V currents (ln_c1d =T default: OFF)
!-----------------------------------------------------------------------
! ! =T read U-V fields for:
ln_uvd_init = .false. ! ocean initialisation
ln_uvd_dyndmp = .false. ! U-V restoring
cn_dir = './' ! root directory for the U-V data location
cn_dir = './SAS/' ! root directory for the U-V data location
!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
......@@ -74,8 +64,6 @@
!!======================================================================
!! *** Surface Boundary Condition namelists *** !!
!! !!
!! STATION_ASF test-case !!
!! !!
!! namsbc surface boundary condition manager (default: NO selection)
!! namsbc_flx flux formulation (ln_flx =T)
!! namsbc_blk Bulk formulae formulation (ln_blk =T)
......@@ -88,34 +76,34 @@
!! namsbc_apr Atmospheric Pressure (ln_apr_dyn =T)
!! namsbc_wave external fields from wave model (ln_wave =T)
!! namberg iceberg floats (ln_icebergs=T)
!! namsbc_fwb freshwater-budget adjustment (nn_fwb > 0)
!!======================================================================
!
!-----------------------------------------------------------------------
&namsbc ! Surface Boundary Condition manager (default: NO selection)
!-----------------------------------------------------------------------
nn_fsbc = 1 !!LOLO! frequency of SBC module call
nn_fsbc = 1 ! frequency of SBC module call
ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk )
!
nn_ice = 2 !LOLO ! =0 no ice boundary condition
nn_ice = 2 ! =0 no ice boundary condition
! ! =1 use observed ice-cover ( => fill namsbc_iif )
! ! =2 or 3 for SI3 and CICE, respectively
!
ln_dm2dc = .false. !#lolo daily mean to diurnal cycle on short wave
/
!-----------------------------------------------------------------------
&namsbc_blk ! namsbc_blk generic Bulk formula / STATION_ASF (ln_blk =T)
&namsbc_blk ! namsbc_blk generic Bulk formula (ln_blk =T)
!-----------------------------------------------------------------------
! ! bulk algorithm :
ln_NCAR = .false. ! "NCAR" algorithm (Large and Yeager 2008)
ln_COARE_3p0 = .false. ! "COARE 3.0" algorithm (Fairall et al. 2003)
ln_COARE_3p6 = .false. ! "COARE 3.6" algorithm (Edson et al. 2013)
ln_ECMWF = .true. ! "ECMWF" algorithm (IFS cycle 45r1)
ln_MFS = .false. ! "MFS" algorithm (MFS/BS Copernicus, Petenuzzo et al 2010)
ln_ANDREAS = .false. ! "ANDREAS" algorithm (Andreas et al. 2015)
rn_zqt = 2. ! Air temperature & humidity reference height (m)
rn_zu = 10. ! Wind vector reference height (m)
nn_iter_algo = 20 ! Number of iterations in bulk param. algo ("stable ABL + weak wind" requires more)
ln_skin_cs = .true. ! use the cool-skin parameterization
ln_skin_wl = .true. ! use the warm-layer parameterization
ln_skin_cs = .false. ! use the cool-skin parameterization => use at least nn_iter_algo > 10
ln_skin_wl = .false. ! use the warm-layer parameterization => use at least nn_iter_algo > 10
!
ln_humi_sph = .false. ! humidity "sn_humi" is specific humidity [kg/kg]
ln_humi_dpt = .true. ! humidity "sn_humi" is dew-point temperature [K]
......@@ -124,26 +112,26 @@
!!
!! Bulk transfer coefficients over sea-ice: (relevant IF: nn_ice >=1 )
ln_Cx_ice_cst = .true. ! use constant ice-air bulk transfer coefficients (value given below)
rn_Cd_i = 1.4e-3 ! sea-ice drag coefficient
rn_Ce_i = 1.4e-3 ! " sublimation coefficient
rn_Ch_i = 1.4e-3 ! " sensible heat flux coefficient
rn_Cd_ia = 1.4e-3 ! sea-ice drag coefficient
rn_Ce_ia = 1.4e-3 ! " sublimation coefficient
rn_Ch_ia = 1.4e-3 ! " sensible heat flux coefficient
ln_Cx_ice_AN05 = .false. ! (Andreas et al. 2005)
ln_Cx_ice_LU12 = .false. ! (Lupkes et al. 2012)
ln_Cx_ice_LG15 = .false. ! (Lupkes & Gryanik 2015)
!
cn_dir = './' ! root directory for the bulk data location
cn_dir = './SAS/' ! root directory for the bulk data location
!___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
sn_wndi = '<SFORC>', 1. , 'u10' , .true. , .false. , 'yearly' , '' ,'' , ''
sn_wndj = '<SFORC>', 1. , 'v10' , .true. , .false. , 'yearly' , '' ,'' , ''
sn_qsr = '<SFORC>', 1. , 'ssrd' , .false. , .false. , 'yearly' , '' , '' , ''
sn_qlw = '<SFORC>', 1. , 'strd' , .false. , .false. , 'yearly' , '' , '' , ''
sn_tair = '<SFORC>', 1. , 't2m' , .true. , .false. , 'yearly' , '' , '' , ''
sn_humi = '<SFORC>', 1. , 'd2m' , .true. , .false. , 'yearly' , '' , '' , ''
sn_prec = '<SFORC>', 1. , 'tp' , .false. , .false. , 'yearly' , '' , '' , ''
sn_snow = '<SFORC>', 1. , 'sf' , .false. , .false. , 'yearly' , '' , '' , ''
sn_slp = '<SFORC>', 1. , 'msl' , .true. , .false. , 'yearly' , '' , '' , ''
sn_wndi = 'ERA5_NorthGreenland_surface_84N_-36E_1h', 1. , 'u10' , .true. , .false. , 'yearly' , '' ,'' , ''
sn_wndj = 'ERA5_NorthGreenland_surface_84N_-36E_1h', 1. , 'v10' , .true. , .false. , 'yearly' , '' ,'' , ''
sn_qsr = 'ERA5_NorthGreenland_surface_84N_-36E_1h', 1. , 'ssrd' , .false. , .false. , 'yearly' , '' , '' , ''
sn_qlw = 'ERA5_NorthGreenland_surface_84N_-36E_1h', 1. , 'strd' , .false. , .false. , 'yearly' , '' , '' , ''
sn_tair = 'ERA5_NorthGreenland_surface_84N_-36E_1h', 1. , 't2m' , .true. , .false. , 'yearly' , '' , '' , ''
sn_humi = 'ERA5_NorthGreenland_surface_84N_-36E_1h', 1. , 'd2m' , .true. , .false. , 'yearly' , '' , '' , ''
sn_prec = 'ERA5_NorthGreenland_surface_84N_-36E_1h', 1. , 'tp' , .false. , .false. , 'yearly' , '' , '' , ''
sn_snow = 'ERA5_NorthGreenland_surface_84N_-36E_1h', 1. , 'sf' , .false. , .false. , 'yearly' , '' , '' , ''
sn_slp = 'ERA5_NorthGreenland_surface_84N_-36E_1h', 1. , 'msl' , .true. , .false. , 'yearly' , '' , '' , ''
/
!-----------------------------------------------------------------------
&namsbc_abl ! Atmospheric Boundary Layer formulation (ln_abl = T)
......@@ -154,24 +142,24 @@
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namsbc_sas ! Stand-Alone Surface boundary condition / STATION_ASF
&namsbc_sas ! Stand-Alone Surface module: ocean data (SAS_SRC only)
!-----------------------------------------------------------------------
l_sasread = .true. ! =T Read in file ; =F set all to 0. (see sbcssm)
!
cn_dir = './' ! root directory for the ocean data location
cn_dir = './SAS/' ! root directory for the ocean data location
!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
sn_usp = '<SFORC>', 1. , 'ssu' , .true. , .false. , 'yearly' , '' , '' , ''
sn_vsp = '<SFORC>', 1. , 'ssv' , .true. , .false. , 'yearly' , '' , '' , ''
sn_tem = '<SFORC>', 1. , 'sst' , .true. , .false. , 'yearly' , '' , '' , ''
sn_sal = '<SFORC>', 1. , 'sss', .true. , .false. , 'yearly' , '' , '' , ''
sn_ssh = '<SFORC>', 1. , 'ssh', .true. , .false. , 'yearly' , '' , '' , ''
sn_usp = 'ERA5_NorthGreenland_surface_84N_-36E_1h', 1. , 'ssu' , .true. , .false. , 'yearly' , '' , '' , ''
sn_vsp = 'ERA5_NorthGreenland_surface_84N_-36E_1h', 1. , 'ssv' , .true. , .false. , 'yearly' , '' , '' , ''
sn_tem = 'ERA5_NorthGreenland_surface_84N_-36E_1h', 1. , 'sst' , .true. , .false. , 'yearly' , '' , '' , ''
sn_sal = 'ERA5_NorthGreenland_surface_84N_-36E_1h', 1. , 'sss', .true. , .false. , 'yearly' , '' , '' , ''
sn_ssh = 'ERA5_NorthGreenland_surface_84N_-36E_1h', 1. , 'ssh', .true. , .false. , 'yearly' , '' , '' , ''
!!
!! Following only needed with STATION_ASF compiled with "sea-ice" support: "key_si3" (ice fraction, ice surface temperature and sea-ice albedo:
sn_ifr = '<SFORC>', 1. , 'siconc' , .true. , .false. , 'yearly' , '' , '' , ''
sn_tic = '<SFORC>', 1. , 'istl1' , .true. , .false. , 'yearly' , '' , '' , ''
sn_ial = '<SFORC>', 1. , 'ialb' , .true. , .false. , 'yearly' , '' , '' , ''
!! Following only needed with SAS compiled with sea-ice model (ice fraction, ice surface temperature and sea-ice albedo)
sn_ifr = 'ERA5_NorthGreenland_surface_84N_-36E_1h', 1. , 'siconc' , .true. , .false. , 'yearly' , '' , '' , ''
sn_tic = 'ERA5_NorthGreenland_surface_84N_-36E_1h', 1. , 'istl1' , .true. , .false. , 'yearly' , '' , '' , ''
sn_ial = 'ERA5_NorthGreenland_surface_84N_-36E_1h', 1. , 'ialb' , .true. , .false. , 'yearly' , '' , '' , ''
/
!!======================================================================
!! *** Lateral boundary condition *** !!
......@@ -193,8 +181,8 @@
!! *** Top/Bottom boundary condition *** !!
!! !!
!! namdrg top/bottom drag coefficient (default: NO selection)
!! namdrg_top top friction (ln_OFF=F & ln_isfcav=T)
!! namdrg_bot bottom friction (ln_OFF=F)
!! namdrg_top top friction (ln_drg_OFF=F & ln_isfcav=T)
!! namdrg_bot bottom friction (ln_drg_OFF=F)
!! nambbc bottom temperature boundary condition (default: OFF)
!! nambbl bottom boundary layer scheme (default: OFF)
!!======================================================================
......@@ -246,7 +234,6 @@
!! namhsb Heat and salt budgets (default: OFF)
!! namdiu Cool skin and warm layer models (default: OFF)
!! namdiu Cool skin and warm layer models (default: OFF)
!! namflo float parameters (default: OFF)
!! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)
!! nam_diadct transports through some sections (default: OFF)
!! nam_diatmb Top Middle Bottom Output (default: OFF)
......
cfgs/C1D/EXP_SASICE/namelist_ice_ref 0 → 120000
View file @ 43f40c58
../../SHARED/namelist_ice_ref
\ No newline at end of file
cfgs/C1D/EXP_SASICE/namelist_ref 0 → 120000
View file @ 43f40c58
../../SHARED/namelist_ref
\ No newline at end of file
cfgs/C1D/EXP_SASOCE/axis_def_nemo.xml 0 → 120000
View file @ 43f40c58
../../SHARED/axis_def_nemo.xml
\ No newline at end of file
tests/STATION_ASF/EXPREF/context_nemo_OCE.xml → cfgs/C1D/EXP_SASOCE/context_nemo.xml
View file @ 43f40c58
......@@ -4,8 +4,6 @@
==============================================================================================
-->
<context id="nemo">
<variable_definition>
<!-- Year/Month/Day of time origin for NetCDF files; defaults to 1800-01-01 -->
......@@ -22,7 +20,6 @@
<!-- Files definition -->
<file_definition src="./file_def_nemo-oce.xml"/> <!-- NEMO ocean dynamics -->
<!--
============================================================================================================
= grid definition = = DO NOT CHANGE =
......@@ -34,7 +31,6 @@
<axis id="depthu" long_name="Vertical U levels" unit="m" positive="down" />
<axis id="depthv" long_name="Vertical V levels" unit="m" positive="down" />
<axis id="depthw" long_name="Vertical W levels" unit="m" positive="down" />
<axis id="nfloat" long_name="Float number" unit="-" />
<axis id="icbcla" long_name="Iceberg class" unit="1" />
<axis id="ncatice" long_name="Ice category" unit="1" />
<axis id="iax_20C" long_name="20 degC isotherm" unit="degC" />
......
cfgs/C1D/EXP_SASOCE/domain_def_nemo.xml 0 → 120000
View file @ 43f40c58
../../SHARED/domain_def_nemo.xml
\ No newline at end of file
cfgs/C1D/EXP_SASOCE/field_def_nemo-oce.xml 0 → 120000
View file @ 43f40c58
../../SHARED/field_def_nemo-oce.xml
\ No newline at end of file
cfgs/C1D/EXP_SASOCE/file_def_nemo-oce.xml 0 → 100644
View file @ 43f40c58
<?xml version="1.0"?>
<file_definition type="one_file" name="@expname@_@freq@_@startdate@_@enddate@" sync_freq="1mo" min_digits="4">
<file_group id="oce_1h" output_freq="1h" output_level="10" enabled=".TRUE."> <!-- 1h files -->
<file id="file1" name_suffix="_gridT" description="ocean T grid variables" compression_level="3" >
<field field_ref="t_skin" name="t_skin" />
<field field_ref="sst_m" name="sst" />
<field field_ref="sss_m" name="sss" />
<field field_ref="ssu_m" name="ssu" />
<field field_ref="ssv_m" name="ssv" />
<!-- -->
<field field_ref="rho_air" name="rho_air" />
<field field_ref="dt_skin" name="dt_skin" />
<!-- -->
<field field_ref="qt_oce" name="qt_oce" />
<field field_ref="qsr_oce" name="qsr_oce" />
<field field_ref="qns_oce" name="qns_oce" />
<field field_ref="qla_oce" name="qla_oce" />
<field field_ref="qsb_oce" name="qsb_oce" />
<field field_ref="qlw_oce" name="qlw_oce" />
<field field_ref="qemp_oce" name="qemp_oce" />
<!-- -->
<field field_ref="evap_oce" name="evap_oce" />
<field field_ref="precip" name="precip" />
<!-- <field field_ref="empmr" name="empmr" /> -->
<!-- -->
<field field_ref="taum_oce" name="taum" />
<field field_ref="utau_oce" name="utau" />
<field field_ref="vtau_oce" name="vtau" />
<field field_ref="wspd" name="windsp" />
<!-- -->
<field field_ref="Cd_oce" name="Cd_oce" />
<field field_ref="Ce_oce" name="Ce_oce" />
<field field_ref="Ch_oce" name="Ch_oce" />
<!-- -->
<field field_ref="theta_zt" name="theta_zt" />
<field field_ref="q_zt" name="q_zt" />
<field field_ref="theta_zu" name="theta_zu" />
<field field_ref="q_zu" name="q_zu" />
<field field_ref="ssq" name="ssq" />
<field field_ref="wspd_blk" name="wspd_blk" />
</file>
</file_group>
<file_group id="oce_1ts" output_freq="1ts" output_level="10" enabled=".TRUE."/> <!-- 1 time step files -->
<file_group id="oce_2h" output_freq="2h" output_level="10" enabled=".TRUE."/> <!-- 2h files -->
<file_group id="oce_3h" output_freq="3h" output_level="10" enabled=".TRUE."/> <!-- 3h files -->
<file_group id="oce_4h" output_freq="4h" output_level="10" enabled=".TRUE."/> <!-- 4h files -->
<file_group id="oce_6h" output_freq="6h" output_level="10" enabled=".TRUE."/> <!-- 6h files -->
<file_group id="oce_1d" output_freq="1d" output_level="10" enabled=".TRUE."/> <!-- 1d files -->
<file_group id="oce_3d" output_freq="3d" output_level="10" enabled=".TRUE."/> <!-- 3d files -->
<file_group id="oce_5d" output_freq="5d" output_level="10" enabled=".TRUE."/> <!-- 5d files -->
<file_group id="oce_1m" output_freq="1mo" output_level="10" enabled=".TRUE."/> <!-- real monthly files -->
<file_group id="oce_2y" output_freq="2y" output_level="10" enabled=".TRUE."/> <!-- real 2y files -->
<file_group id="oce_5y" output_freq="5y" output_level="10" enabled=".TRUE."/> <!-- real 5y files -->
<file_group id="oce_10y" output_freq="10y" output_level="10" enabled=".TRUE."/> <!-- real 10y files -->
</file_definition>
cfgs/C1D/EXP_SASOCE/grid_def_nemo.xml 0 → 120000
View file @ 43f40c58
../../SHARED/grid_def_nemo.xml
\ No newline at end of file
tests/ISOMIP/EXPREF/iodef.xml → cfgs/C1D/EXP_SASOCE/iodef.xml
View file @ 43f40c58
......@@ -12,7 +12,7 @@
<variable id="info_level" type="int">10</variable>
<variable id="using_server" type="bool">false</variable>
<variable id="using_oasis" type="bool">false</variable>
<variable id="oasis_codes_id" type="string" >oceanx</variable>
<variable id="oasis_codes_id" type="string">oceanx</variable>
</variable_definition>
</context>
......
tests/STATION_ASF/EXPREF/PAPA/oce/namelist_ecmwf_cfg → cfgs/C1D/EXP_SASOCE/namelist_cfg
View file @ 43f40c58
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! NEMO/OCE Configuration namelist : overwrite default values defined in SHARED/namelist_ref
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! STATION_ASF configuration, PAPA forcing
!! C1D - SAS - OCE configuration !!
!!======================================================================
!! *** Domain & Run management namelists *** !!
!! !!
......@@ -10,22 +10,19 @@
!! namcfg parameters of the configuration (default: user defined GYRE)
!! namwad Wetting and drying (default: OFF)
!! namtsd data: temperature & salinity (default: OFF)
!! namcrs coarsened grid (for outputs and/or TOP) (ln_crs =T)
!! namc1d 1D configuration options (ln_c1d =T)
!! namc1d_dyndmp 1D newtonian damping applied on currents (ln_c1d =T)
!! namc1d_uvd 1D data (currents) (ln_c1d =T)
!!======================================================================
!
!-----------------------------------------------------------------------
&namusr_def ! STATION_ASF user defined namelist
&namusr_def ! C1D user defined namelist
!-----------------------------------------------------------------------
rn_dept1 = 1.0 ! depth (m) at which the prescribed SST is taken == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!)
! ! [PAPA station = 1 m]
rn_bathy = 1. ! depth in meters
! rn_dept1 = 1.0 ! depth (m) at which the prescribed SST is taken == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!)
/
!-----------------------------------------------------------------------
&namrun ! parameters of the run
!-----------------------------------------------------------------------
cn_exp = 'STATION_ASF-ECMWF_PAPA' ! experience name
cn_exp = 'C1D_SASOCE' ! experience name
nn_it000 = 1 ! first time step
nn_itend = 8760 ! last time step => 3 years (including 1 leap!) at dt=3600s
nn_date0 = 20180101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
......@@ -33,37 +30,30 @@
!
nn_stock = 8760 ! 1 year at dt=3600s / frequency of creation of a restart file (modulo referenced to 1)
nn_write = 8760 ! 1 year at dt=3600s / frequency of write in the output file (modulo referenced to nn_it000)
ln_top = .false. ! Consider (T) or bypass (F) the TOP routines when the key_top is activated
/
!-----------------------------------------------------------------------
&namdom ! time and space domain / STATION_ASF
&namdom ! time and space domain
!-----------------------------------------------------------------------
ln_linssh = .true. ! =T linear free surface ==>> model level are fixed in time
rn_Dt = 3600. ! time step for the dynamics and tracer => frequency of atmospheric forcing
ln_meshmask = .false.
!
ln_c1d = .true. ! Single column domain (1x1pt) (T => fill namc1d)
/
!-----------------------------------------------------------------------
&namcfg ! parameters of the configuration (default: use namusr_def in namelist_cfg)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namc1d ! 1D configuration options / STATION_ASF (ln_c1d =T default: PAPA station)
&namc1d ! 1D configuration options (ln_c1d =T default: PAPA station)
!-----------------------------------------------------------------------
rn_lat1d = 50.1 ! Column latitude [PAPA station]
rn_lon1d = 215.1 ! Column longitude [PAPA station]
/
!-----------------------------------------------------------------------
&namc1d_dyndmp ! U & V newtonian damping (ln_c1d =T default: OFF)
!-----------------------------------------------------------------------
ln_dyndmp = .false. ! add a damping term (T) or not (F)
/
!-----------------------------------------------------------------------
&namc1d_uvd ! data: U & V currents (ln_c1d =T default: OFF)
!-----------------------------------------------------------------------
! ! =T read U-V fields for:
ln_uvd_init = .false. ! ocean initialisation
ln_uvd_dyndmp = .false. ! U-V restoring
cn_dir = './' ! root directory for the U-V data location
cn_dir = './SAS/' ! root directory for the U-V data location
!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
......@@ -74,8 +64,6 @@
!!======================================================================
!! *** Surface Boundary Condition namelists *** !!
!! !!
!! STATION_ASF test-case !!
!! !!
!! namsbc surface boundary condition manager (default: NO selection)
!! namsbc_flx flux formulation (ln_flx =T)
!! namsbc_blk Bulk formulae formulation (ln_blk =T)
......@@ -88,22 +76,17 @@
!! namsbc_apr Atmospheric Pressure (ln_apr_dyn =T)
!! namsbc_wave external fields from wave model (ln_wave =T)
!! namberg iceberg floats (ln_icebergs=T)
!! namsbc_fwb freshwater-budget adjustment (nn_fwb > 0)
!!======================================================================
!
!-----------------------------------------------------------------------
&namsbc ! Surface Boundary Condition manager (default: NO selection)
!-----------------------------------------------------------------------
nn_fsbc = 1 !!LOLO! frequency of SBC module call
nn_fsbc = 1 ! frequency of SBC module call
ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk )
!
nn_ice = 0 !LOLO ! =0 no ice boundary condition
! ! =1 use observed ice-cover ( => fill namsbc_iif )
! ! =2 or 3 for SI3 and CICE, respectively
!
ln_dm2dc = .false. !#lolo daily mean to diurnal cycle on short wave
/
!-----------------------------------------------------------------------
&namsbc_blk ! namsbc_blk generic Bulk formula / STATION_ASF (ln_blk =T)
&namsbc_blk ! namsbc_blk generic Bulk formula (ln_blk =T)
!-----------------------------------------------------------------------
! ! bulk algorithm :
ln_NCAR = .false. ! "NCAR" algorithm (Large and Yeager 2008)
......@@ -118,11 +101,9 @@
ln_skin_wl = .true. ! use the warm-layer parameterization
!
ln_humi_sph = .false. ! humidity "sn_humi" is specific humidity [kg/kg]
ln_humi_dpt = .false. ! humidity "sn_humi" is dew-point temperature [K]
ln_humi_rlh = .true. ! humidity "sn_humi" is relative humidity [%]
ln_tair_pot = .false. ! air temperature read in "sn_tair" is already POTENTIAL TEMPERATURE, NOT ABSOLUTE (ECMWF => ln_tair_pot=.false.)
!
cn_dir = './' ! root directory for the bulk data location
cn_dir = './SAS/' ! root directory for the bulk data location
!___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
......@@ -145,11 +126,11 @@
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namsbc_sas ! Stand-Alone Surface boundary condition / STATION_ASF
&namsbc_sas ! Stand-Alone Surface module: ocean data (SAS_SRC only)
!-----------------------------------------------------------------------
l_sasread = .true. ! =T Read in file ; =F set all to 0. (see sbcssm)
!
cn_dir = './' ! root directory for the ocean data location
cn_dir = './SAS/' ! root directory for the ocean data location
!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
......@@ -179,8 +160,8 @@
!! *** Top/Bottom boundary condition *** !!
!! !!
!! namdrg top/bottom drag coefficient (default: NO selection)
!! namdrg_top top friction (ln_OFF=F & ln_isfcav=T)
!! namdrg_bot bottom friction (ln_OFF=F)
!! namdrg_top top friction (ln_drg_OFF=F & ln_isfcav=T)
!! namdrg_bot bottom friction (ln_drg_OFF=F)
!! nambbc bottom temperature boundary condition (default: OFF)
!! nambbl bottom boundary layer scheme (default: OFF)
!!======================================================================
......@@ -232,7 +213,6 @@
!! namhsb Heat and salt budgets (default: OFF)
!! namdiu Cool skin and warm layer models (default: OFF)
!! namdiu Cool skin and warm layer models (default: OFF)
!! namflo float parameters (default: OFF)
!! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)
!! nam_diadct transports through some sections (default: OFF)
!! nam_diatmb Top Middle Bottom Output (default: OFF)
......
cfgs/C1D/EXP_SASOCE/namelist_ref 0 → 120000
View file @ 43f40c58
../../SHARED/namelist_ref
\ No newline at end of file
tests/C1D_ASICS/MY_SRC/usrdef_hgr.F90 → cfgs/C1D/MY_SRC/usrdef_hgr.F90
View file @ 43f40c58
......@@ -2,7 +2,7 @@ MODULE usrdef_hgr
!!======================================================================
!! *** MODULE usrdef_hgr ***
!!
!! === C1D_ASICS configuration ===
!! === C1D configuration ===
!!
!! User defined : mesh and Coriolis parameter of a user configuration
!!======================================================================
......@@ -27,8 +27,7 @@ MODULE usrdef_hgr
PUBLIC usr_def_hgr ! called in domhgr.F90
!!----------------------------------------------------------------------
!! NEMO/OCE 4.0 , NEMO Consortium (2018)
!! $Id: usrdef_hgr.F90 10072 2018-08-28 15:21:50Z nicolasmartin $
!! NEMO/OCE 5.0, NEMO Consortium (2024)
!! Software governed by the CeCILL license (see ./LICENSE)
!!----------------------------------------------------------------------
CONTAINS
......@@ -46,9 +45,8 @@ CONTAINS
!!
!! ** Method : set all intent(out) argument to a proper value
!!
!! Here C1D_ASICS configuration :
!! Rectangular 3x3 domain
!! - Located at 10E-42N
!! Here C1D configuration :
!! Rectangular 1x1 domain
!! - a constant horizontal resolution
!!
!! ** Action : - define longitude & latitude of t-, u-, v- and f-points (in degrees)
......@@ -101,6 +99,10 @@ CONTAINS
!
! !== Coriolis parameter ==!
kff = 0 ! indicate to compute Coriolis parameter afterward
IF( .FALSE.) THEN ! just to avoid compilation warning...
pff_f(1,1) = -1._wp
pff_t(1,1) = -1._wp
ENDIF
!
END SUBROUTINE usr_def_hgr
......
tests/C1D_ASICS/MY_SRC/usrdef_nam.F90 → cfgs/C1D/MY_SRC/usrdef_nam.F90
View file @ 43f40c58
......@@ -2,7 +2,7 @@ MODULE usrdef_nam
!!======================================================================
!! *** MODULE usrdef_nam ***
!!
!! === C1D_ASICS configuration ===
!! === C1D configuration ===
!!
!! User defined : set the domain characteristics of a user configuration
!!======================================================================
......@@ -13,8 +13,8 @@ MODULE usrdef_nam
!! usr_def_nam : read user defined namelist and set global domain size
!! usr_def_hgr : initialize the horizontal mesh
!!----------------------------------------------------------------------
USE oce , ONLY: l_SAS
USE dom_oce , ONLY: nimpp, njmpp ! ocean space and time domain
USE dom_oce , ONLY: ln_zco, ln_zps, ln_sco ! flag of type of coordinate
USE par_oce ! ocean space and time domain
USE phycst ! physical constants
!
......@@ -31,9 +31,10 @@ MODULE usrdef_nam
INTEGER, PUBLIC :: nn_GYRE ! 1/nn_GYRE = the resolution chosen in degrees and thus defining the horizontal domain size
REAL(wp), PUBLIC:: rn_bathy ! Depth in meters for 1D configuration
!! * Substitutions
# include "read_nml_substitute.h90"
!!----------------------------------------------------------------------
!! NEMO/OCE 4.0 , NEMO Consortium (2018)
!! $Id: usrdef_nam.F90 12377 2020-02-12 14:39:06Z acc $
!! NEMO/OCE 5.0, NEMO Consortium (2024)
!! Software governed by the CeCILL license (see ./LICENSE)
!!----------------------------------------------------------------------
CONTAINS
......@@ -62,8 +63,7 @@ CONTAINS
NAMELIST/namusr_def/ rn_bathy
!!----------------------------------------------------------------------
!
READ ( numnam_cfg, namusr_def, IOSTAT = ios, ERR = 902 )
902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namusr_def in configuration namelist' )
READ_NML_(numnam_cfg,cfg,namusr_def,.TRUE.)
!
IF(lwm) WRITE( numond, namusr_def )
!
......@@ -73,7 +73,11 @@ CONTAINS
! Global Domain size: C1D domain is 1 x 1 grid-points x 75 or vertical levels
kpi = 1
kpj = 1
kpk = 75
IF( .NOT.l_sas ) THEN
kpk = 75
ELSE
kpk = 2
ENDIF
! ! Set the lateral boundary condition of the global domain
ldIperio = .TRUE. ; ldJperio = .TRUE. ! C1D configuration : 1x1 basin with cyclic Est-West and Norht-South condition
ldNFold = .FALSE. ; cdNFtype = '-'
......@@ -81,16 +85,12 @@ CONTAINS
! ! control print
IF(lwp) THEN
WRITE(numout,*) ' '
WRITE(numout,*) 'usr_def_nam : read the user defined namelist (namusr_def) in namelist_cfg'
WRITE(numout,*) 'usr_def_nam : read the user defined namelist (namusr_def) in namelist_cfg'
WRITE(numout,*) '~~~~~~~~~~~ '
WRITE(numout,*) ' Namelist namusr_def : C1 case'
WRITE(numout,*) ' type of vertical coordinate : '
WRITE(numout,*) ' z-coordinate flag ln_zco = ', ln_zco
WRITE(numout,*) ' z-partial-step coordinate flag ln_zps = ', ln_zps
WRITE(numout,*) ' s-coordinate flag ln_sco = ', ln_sco
WRITE(numout,*) ' C1D domain = 1 x 1 x 75 grid-points '
WRITE(numout,*) ' resulting global domain size : jpiglo = ', kpi
WRITE(numout,*) ' jpjglo = ', kpj
WRITE(numout,*) ' Namelist namusr_def : C1D case'
WRITE(numout,*) ' C1D domain = ',kpi,' x ',kpj,' x ',kpk,' grid-points'
WRITE(numout,*) ' resulting global domain size : Ni0glo = ', kpi
WRITE(numout,*) ' Nj0glo = ', kpj
WRITE(numout,*) ' jpkglo = ', kpk
WRITE(numout,*) ' '
ENDIF
......
cfgs/C1D_PAPA/MY_SRC/usrdef_zgr.F90 → cfgs/C1D/MY_SRC/usrdef_zgr.F90
View file @ 43f40c58
......@@ -2,11 +2,12 @@ MODULE usrdef_zgr
!!======================================================================
!! *** MODULE usrdef_zgr ***
!!
!! === C1D_PAPA configuration ===
!! === C1D configuration ===
!!
!! User defined : vertical coordinate system of a user configuration
!!======================================================================
!! History : 4.0 ! 2016-06 (R. Bourdalle-Badie) Original code
!! 4.3 ! 2023-01 (S. Techene,G. Madec) New zps : depth and w-level scale factors are horizontally uniform
!!----------------------------------------------------------------------
!!----------------------------------------------------------------------
......@@ -32,32 +33,31 @@ MODULE usrdef_zgr
!! * Substitutions
# include "do_loop_substitute.h90"
!!----------------------------------------------------------------------
!! NEMO/OCE 4.0 , NEMO Consortium (2018)
!! $Id: usrdef_zgr.F90 13295 2020-07-10 18:24:21Z acc $
!! NEMO/OCE 5.0, NEMO Consortium (2024)
!! Software governed by the CeCILL license (see ./LICENSE)
!!----------------------------------------------------------------------
CONTAINS
CONTAINS
SUBROUTINE usr_def_zgr( ld_zco , ld_zps , ld_sco , ld_isfcav, & ! type of vertical coordinate
& k_top , k_bot , & ! top & bottom ocean level
& pdept_1d, pdepw_1d, pe3t_1d , pe3w_1d , & ! 1D reference vertical coordinate
& pe3t , pe3u , pe3v , pe3f , & ! 3D t-level vertical scale factors
& pdept , pdepw , & ! 3D t & w-points depth
& pe3t , pe3u , pe3v , pe3f , & ! vertical scale factors
& pe3w , pe3uw , pe3vw , & ! - - -
& k_top , k_bot ) ! top & bottom ocean level
& pe3w , pe3uw , pe3vw ) ! 3D w-level vertical scale factors
!!---------------------------------------------------------------------
!! *** ROUTINE usr_def_zgr ***
!!
!! ** Purpose : User defined the vertical coordinates
!!
!!----------------------------------------------------------------------
LOGICAL , INTENT(out) :: ld_zco, ld_zps, ld_sco ! vertical coordinate flags
LOGICAL , INTENT(out) :: ld_isfcav ! under iceshelf cavity flag
REAL(wp), DIMENSION(:) , INTENT(out) :: pdept_1d, pdepw_1d ! 1D grid-point depth [m]
REAL(wp), DIMENSION(:) , INTENT(out) :: pe3t_1d , pe3w_1d ! 1D grid-point depth [m]
REAL(wp), DIMENSION(:,:,:), INTENT(out) :: pdept, pdepw ! grid-point depth [m]
REAL(wp), DIMENSION(:,:,:), INTENT(out) :: pe3t , pe3u , pe3v , pe3f ! vertical scale factors [m]
REAL(wp), DIMENSION(:,:,:), INTENT(out) :: pe3w , pe3uw, pe3vw ! i-scale factors
INTEGER , DIMENSION(:,:) , INTENT(out) :: k_top, k_bot ! first & last ocean level
LOGICAL , INTENT(out) :: ld_zco, ld_zps, ld_sco ! vertical coordinate flags
LOGICAL , INTENT(out) :: ld_isfcav ! under iceshelf cavity flag
INTEGER , DIMENSION(:,:) , INTENT(out) :: k_top, k_bot ! first & last ocean level
REAL(wp), DIMENSION(:) , INTENT(out) :: pdept_1d, pdepw_1d ! 1D grid-point depth [m]
REAL(wp), DIMENSION(:) , INTENT(out) :: pe3t_1d , pe3w_1d ! 1D grid-point depth [m]
REAL(wp), DIMENSION(:,:,:), OPTIONAL, INTENT(out) :: pe3t , pe3u , pe3v , pe3f ! t-level vertical scale factors [m]
REAL(wp), DIMENSION(:,:,:), OPTIONAL, INTENT(out) :: pdept, pdepw ! grid-point depth [m]
REAL(wp), DIMENSION(:,:,:), OPTIONAL, INTENT(out) :: pe3w , pe3uw, pe3vw ! w-level vertical scale factors [m]
!
INTEGER :: ji, jj, jk ! dummy indices
INTEGER :: ik ! local integers
......@@ -142,40 +142,43 @@ CONTAINS
! !* bottom ocean compute from the depth of grid-points
k_bot(:,:) = jpkm1
DO jk = jpkm1, 1, -1
ze3min = 0.1_wp * pe3t_1d (jk)
WHERE( zht(:,:) < pdepw_1d(jk) + ze3min ) k_bot(:,:) = jk-1
ze3min = 0.1_wp * pe3t_1d (jk)
WHERE( zht(:,:) < pdepw_1d(jk) + ze3min ) k_bot(:,:) = MAX( 1 , jk-1)
END DO
!
! !* vertical coordinate system
DO jk = 1, jpk ! initialization to the reference z-coordinate
pdept(:,:,jk) = pdept_1d(jk)
pdepw(:,:,jk) = pdepw_1d(jk)
pe3t (:,:,jk) = pe3t_1d (jk)
pe3u (:,:,jk) = pe3t_1d (jk)
pe3v (:,:,jk) = pe3t_1d (jk)
pe3f (:,:,jk) = pe3t_1d (jk)
!!st : new zps framework toward penalization only fluxed at t-level are affected by scale factor change
#if defined key_vco_3d
pdept(:,:,jk) = pdept_1d(jk)
pdepw(:,:,jk) = pdepw_1d(jk)
pe3w (:,:,jk) = pe3w_1d (jk)
pe3uw(:,:,jk) = pe3w_1d (jk)
pe3vw(:,:,jk) = pe3w_1d (jk)
#endif
END DO
! bottom scale factors and depth at T- and W-points
DO_2D( 1, 1, 1, 1 )
ik = k_bot(ji,jj)
pdepw(ji,jj,ik+1) = MIN( zht(ji,jj) , pdepw_1d(ik+1) )
pe3t (ji,jj,ik ) = pdepw(ji,jj,ik+1) - pdepw(ji,jj,ik)
!!st pdepw(ji,jj,ik+1) = MIN( zht(ji,jj) , pdepw_1d(ik+1) )
pe3t (ji,jj,ik ) = MIN( zht(ji,jj) , pdepw_1d(ik+1) )- pdepw_1d(ik)
pe3t (ji,jj,ik+1) = pe3t (ji,jj,ik )
!
pdept(ji,jj,ik ) = pdepw(ji,jj,ik ) + pe3t (ji,jj,ik ) * 0.5_wp
pdept(ji,jj,ik+1) = pdepw(ji,jj,ik+1) + pe3t (ji,jj,ik+1) * 0.5_wp
pe3w (ji,jj,ik+1) = pdept(ji,jj,ik+1) - pdept(ji,jj,ik) ! = pe3t (ji,jj,ik )
!!st pdept(ji,jj,ik ) = pdepw(ji,jj,ik ) + pe3t (ji,jj,ik ) * 0.5_wp
!!st pdept(ji,jj,ik+1) = pdepw(ji,jj,ik+1) + pe3t (ji,jj,ik+1) * 0.5_wp
!!st pe3w (ji,jj,ik+1) = pdept(ji,jj,ik+1) - pdept(ji,jj,ik) ! = pe3t (ji,jj,ik )
END_2D
! ! bottom scale factors and depth at U-, V-, UW and VW-points
! ! usually Computed as the minimum of neighbooring scale factors
pe3u (:,:,:) = pe3t(:,:,:) ! HERE C1D configuration :
pe3v (:,:,:) = pe3t(:,:,:) ! e3 increases with k-index
pe3f (:,:,:) = pe3t(:,:,:) ! so e3 minimum of (i,i+1) points is (i) point
pe3uw(:,:,:) = pe3w(:,:,:) ! in j-direction e3v=e3t and e3f=e3v
pe3vw(:,:,:) = pe3w(:,:,:) ! ==>> no need of lbc_lnk calls
!!st pe3uw(:,:,:) = pe3w(:,:,:) ! in j-direction e3v=e3t and e3f=e3v
!!st pe3vw(:,:,:) = pe3w(:,:,:) ! ==>> no need of lbc_lnk calls
!
!
END SUBROUTINE usr_def_zgr
......