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.
user/flavoni/libIGCM/change_BC – NEMO
wiki:user/flavoni/libIGCM/change_BC

Version 13 (modified by flavoni, 15 years ago) (diff)

--

Change Surface Boundary Conditions


How to change Surface Boundary forcing files

By default from nemo web site :

http://www.nemo-ocean.eu/index.php//Using-NEMO/Configurations/ORCA2_LIM

you have forcing files for bulk CLIO formulae

If you want to use other atmospheric forcing methods you can chose :

analytical
flux
bulk ( CLIO formulae )
bulk ( CORE formulae )
coupled

Example : using CORE bulk formulae

You have to change

  • in modipsl/config/ORCA2_LIM/IGCM00/COMP/opa9.card :
  1. use ${year} for loop
  2. create directory to put new forcing files ( or use common account : https://forge.ipsl.jussieu.fr/nemo/wiki/private/libIGCM/common_account )
  3. Add new forcing file to list of files for initial conditions

ListNonDel?= (${R_BC}/OCE/${config_UserChoices_TagName}/radlw_DFS4-ORCA2_${yearm1}.nc, qlw_core_y${yearm1}.nc), \

(${R_BC}/OCE/${config_UserChoices_TagName}/radlw_DFS4-ORCA2_${year}.nc, qlw_core_y${year}.nc), \
(${R_BC}/OCE/${config_UserChoices_TagName}/radlw_DFS4-ORCA2_${yearp1}.nc, qlw_core_y${yearp1}.nc), \
(${R_BC}/OCE/${config_UserChoices_TagName}/radsw_DFS4-ORCA2_${yearm1}.nc, qsw_core_y${yearm1}.nc), \
(${R_BC}/OCE/${config_UserChoices_TagName}/radsw_DFS4-ORCA2_${year}.nc, qsw_core_y${year}.nc), \
(${R_BC}/OCE/${config_UserChoices_TagName}/radsw_DFS4-ORCA2_${yearp1}.nc, qsw_core_y${yearp1}.nc), \
(${R_BC}/OCE/${config_UserChoices_TagName}/q2_DFS4-ORCA2_${yearm1}.nc, q2_core_y${yearm1}.nc), \
(${R_BC}/OCE/${config_UserChoices_TagName}/q2_DFS4-ORCA2_${year}.nc, q2_core_y${year}.nc), \
(${R_BC}/OCE/${config_UserChoices_TagName}/q2_DFS4-ORCA2_${yearp1}.nc, q2_core_y${yearp1}.nc), \
(${R_BC}/OCE/${config_UserChoices_TagName}/t2_DFS4-ORCA2_${yearm1}.nc, t2_core_y${yearm1}.nc), \
(${R_BC}/OCE/${config_UserChoices_TagName}/t2_DFS4-ORCA2_${year}.nc, t2_core_y${year}.nc), \
(${R_BC}/OCE/${config_UserChoices_TagName}/t2_DFS4-ORCA2_${yearp1}.nc, t2_core_y${yearp1}.nc), \
(${R_BC}/OCE/${config_UserChoices_TagName}/u10_DFS4-ORCA2_${yearm1}.nc, u10_core_y${yearm1}.nc),\
(${R_BC}/OCE/${config_UserChoices_TagName}/u10_DFS4-ORCA2_${year}.nc, u10_core_y${year}.nc),\
(${R_BC}/OCE/${config_UserChoices_TagName}/u10_DFS4-ORCA2_${yearp1}.nc, u10_core_y${yearp1}.nc),\
(${R_BC}/OCE/${config_UserChoices_TagName}/v10_DFS4-ORCA2_${yearm1}.nc, v10_core_y${yearm1}.nc),\
(${R_BC}/OCE/${config_UserChoices_TagName}/v10_DFS4-ORCA2_${year}.nc, v10_core_y${year}.nc),\
(${R_BC}/OCE/${config_UserChoices_TagName}/v10_DFS4-ORCA2_${yearp1}.nc, v10_core_y${yearp1}.nc),\
(${R_BC}/OCE/${config_UserChoices_TagName}/precip_DFS4-ORCA2_${yearm1}_12.nc, precip_core_y${yearm1}.nc),\
(${R_BC}/OCE/${config_UserChoices_TagName}/precip_DFS4-ORCA2_${year}_12.nc, precip_core_y${year}.nc),\
(${R_BC}/OCE/${config_UserChoices_TagName}/precip_DFS4-ORCA2_${yearp1}_12.nc, precip_core_y${yearp1}.nc),\
(${R_BC}/OCE/${config_UserChoices_TagName}/snow_DFS4-ORCA2_${yearm1}_12.nc, snow_core_y${yearm1}.nc),\
(${R_BC}/OCE/${config_UserChoices_TagName}/snow_DFS4-ORCA2_${year}_12.nc, snow_core_y${year}.nc),\
(${R_BC}/OCE/${config_UserChoices_TagName}/snow_DFS4-ORCA2_${yearp1}_12.nc, snow_core_y${yearp1}.nc)

NOTA BENE : It is very important to control variable name in netcdf files. In namelist you have to put same variable name that you have in netcdf files!!! So you will change also output name file in namelist ( like qw_lcore, etc... )

  • in modipsl/config/ORCA2_LIM/IGCM00/PARAM/namelist :
  1. activation of bulk core formulae in Surface Boundary Condition function ( namsbc )

ln_blk_core = .true.

  1. set correct frequency, name and interpolation in CORE bulk formulae ( namsbc_core )
file name frequency (hours) (if <0 months) variable name time interpol. (logical) clim (T/F) 'yearly' or monthly weights filename rotation pairing
sn_wndi 'u10_core' 6. 'u10' .true. .false. 'yearly' ' ' ' '
sn_wndj 'v10_core' 6. 'v10' .true. .false. 'yearly' ' ' ' '
sn_qsr 'qsw_core' 6. 'radsw' .true. .false. 'yearly' ' ' ' '
sn_qlw 'qlw_core' 6. 'radlw' .true. .false. 'yearly' ' ' ' '
sn_tair 't2_core' 6. ' t2' .true. .false. 'yearly' ' ' ' '
sn_humi 'q2_core' 6. 'q2' .true. .false. 'yearly' ' ' ' '
sn_prec 'precip_core' -1. 'precip' .true. .false. 'monthly' ' ' ' '
sn_snow 'snow_core' -1. 'snow' .true. .false. 'monthly' ' ' ' '

!

cn_dir = './' ! root directory for the location of the bulk files
ln_2m = .true. ! air temperature and humidity referenced at 2m (T) instead 10m (F)
alpha_precip= 1. ! multiplicative factor for precipitation (total & snow)