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 37 (modified by flavoni, 14 years ago) (diff)

--

Change Surface Boundary Conditions


How to change Surface Boundary forcing files

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

When you download
./model NEMO_DEV
by default you have a configuration with climatological forcing, bulk clio formulae.

Change from climatolgical (bulk_clio)

to interannual (bulk_core) (with DRAKKAR forcing files) ===

If you want to change forcing method to have an Interannaul run with bulk_core formulae
you have to :

  • in modipsl/config/ORCA2_LIM/IGCM00/COMP/opa9.card:
  1. Choose interannual run

Interannual_Run='''y''' [[BR]]

  1. List Surface Boundary Condition original files name ( get by the job )
List_jobsbc=(radlw_DFS4-ORCA2, radsw_DFS4-ORCA2, q2_DFS4-ORCA2, t2_DFS4-ORCA2, u10_DFS4-ORCA2, v10_DFS4-ORCA2, snow_DFS4-ORCA2, precip_DFS4-ORCA2)

NOTA BENE : It is very important to control variable name in netcdf files. In namelist you have to put exactly the same variable name that you have in netcdf files.

  1. List Surface Boundary Condition files name expected by NEMO ( same as the ones in namelist )
List_runsbc=(qlw_core, qsw_core, q2_core, t2_core, u10_core, v10_core, snow_core, precip_core)

NOTA BENE : If you leave List_runsbc=() this means that List_runsbc=List_jobsbc=(radlw_DFS4-ORCA2, radsw_DFS4-ORCA2, q2_DFS4-ORCA2, t2_DFS4-ORCA2, u10_DFS4-ORCA2, v10_DFS4-ORCA2, snow_DFS4-ORCA2, precip_DFS4-ORCA2) and in this case it CAN'T run because these are not same names that in the namelist!

  • in modipsl/config/ORCA2_LIM/IGCM00/PARAM/namelist:
  1. Activation of bulk core formulae in Surface Boundary Condition function ( namsbc )
ln_blk_core = .true.  [[BR]]
(and ln_blk_clio = .false.)
  • in modipsl/config/ORCA2_LIM/IGCM00/config.card:
  1. Change CalendarType? and PerdioLenght?
CalendarType=noleap
PeriodLength=1M
R_BC=/path/directory/with/DFS4_files_and_sbc

(note: DFS4_files_and_sbc directory has to have DFS4 forcing files, and boundary conditions files. [BR]] In my directory I changed precip and snow file's names, because I removed suffix _12, to have same names for all variables)

  1. Rebuild job
  cd modipsl
  ./util/ins_job

Change from interannaul, bulk_core to climatolgical, bulk_clio

  1. Choose climatological run
  • in modipsl/config/ORCA2_LIM/IGCM00/COMP/opa9.card you have to :

Interannual_Run='''n''' 
  1. List Surface Boundary Condition original files name ( get by the job )
List_jobsbc'''=(flx, taux_1m, tauy_1m)

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.

  1. List Surface Boundary Condition files name expected by NEMO ( same as the ones in namelist )
List_runsbc'''=()

NOTA BENE : List_runsbc=() this means that List_runsbc=List_jobsbc=(flx, taux_1m, tauy_1m) and in this case it is ok. If you prefer you can put also List_runsbc=(flx, taux_1m, tauy_1m); is exactly the same thing.

  • in modipsl/config/ORCA2_LIM/IGCM00/PARAM/namelist :
  1. Activation of bulk clio formulae in Surface Boundary Condition function ( namsbc )
ln_blk_clio = .true. 
(and ln_blk_core = .false.)
  • in modipsl/config/ORCA2_LIM/IGCM00/config.card:
  1. Change CalendarType? and PerdioLenght?
CalendarType=noleap
PeriodLength=1Y
  1. Rebuild job
  cd modipsl
  ./util/ins_job

Change from climatolgical, bulk_clio to climatological, bulk_core (with CORE2 forcing files)

If you want to change forcing method to have a Climatological run with bulk_core formulae
you have to :

NOTE: if forcing files are not on ORCA2 grid you have to interpolate them. (see: http://www.nemo-ocean.eu/Using-NEMO/Pre-and-post-processing-packages/Pre-Processing/cdo-interpolation )

  1. List Surface Boundary Condition original files name ( get by the job )
List_jobsbc=(u_10.15JUNE2009_orca2, v_10.15JUNE2009_orca2, ncar_rad.15JUNE2009_orca2, t_10.15JUNE2009_orca2, q_10.15JUNE2009_orca2, ncar_precip.15JUNE2009_orca2)

NOTA BENE : It is very important to control variable name in netcdf files. In namelist you have to put exactly the same variable name that you have in netcdf files.

  1. List Surface Boundary Condition files name expected by NEMO ( same as the ones in namelist )
List_runsbc=()

  • in modipsl/config/ORCA2_LIM/IGCM00/PARAM/namelist :
  1. Activation of bulk core formulae in Surface Boundary Condition function ( namsbc )
ln_blk_core = .true.  [[BR]]
(and ln_blk_clio = .false.)
  1. Example of namelist :
&namsbc_core namsbc_core CORE bulk formulea
file name frequency (hours) variable time interpol. clim 'yearly'/ weights rotation
(if <0 months) name (logical) (T/F) 'monthly' filename pairing
sn_wndi ='u_10.15JUNE2009_orca2' 6 'U_10' .false. .true. 'yearly'
sn_wndj ='v_10.15JUNE2009_orca2' 6 'V_10' .false. .true. 'yearly'
sn_qsr ='ncar_rad.15JUNE2009_orca2' 24 'SWDN' .false. .true. 'yearly'
sn_qlw ='ncar_rad.15JUNE2009_orca2' 24 'LWDN' .false. .true. 'yearly'
sn_tair ='t_10.15JUNE2009_orca2' 6 'T_10' .false. .true. 'yearly'
sn_humi ='q_10.15JUNE2009_orca2' 6 'Q_10' .false. .true. 'yearly'
sn_prec ='ncar_precip.15JUNE2009_orca2' -1 'RAIN' .false. .true. 'yearly'
sn_snow ='ncar_precip.15JUNE2009_orca2' -1 'SNOW' .false. .true. 'yearly'
sn_tdif ='taudif_core' 24 'taudif' .false. .true. 'yearly'
cn_dir ='./' root directory for the location of the bulk files
ln_2m = .false. air temperature and humidity referenced at 2m (T) instead 10m (F)
ln_taudif = .false. HF tau contribution: use "mean of stress module - module of the mean stress" data ?
rn_pfac = 1. multiplicative factor for precipitation (total & snow)
  • in modipsl/config/ORCA2_LIM/IGCM00/config.card:
  1. Change CalendarType? and PerdioLenght?
CalendarType=noleap
PeriodLength=1M
R_BC=/path/directory/with/CORE2_climatological_interpolates_files_and_sbc

(note: CORE2_climatological_interpolates_files_and_sbc directory has to have CORE2 climatological, interpolated forcing files, and boundary conditions files.)

  1. Rebuild job
  cd modipsl
  ./util/ins_job