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.
IS-ENES/core2_cmip – NEMO
wiki:IS-ENES/core2_cmip

Version 104 (modified by flavoni, 14 years ago) (diff)

--

TOC(heading=Integrating COREII forcing in nemo-libigcm)?

Integrating COREII forcing in nemo tags v3.2 / libigcm tags v1.4


Purposes:


  • Document the using of COREII forcing for NEMO (tag v32, revision : 1803), in configuration ORCA2_LIM, with libIGCM (tags v1.4, revision : 200).



Step by step:

Preparing COREII forcing files on ORCA2 grid

climatological files

  1. Download COREII forcing files : CYNF
  2. Interpolation of forcing files (with cdo (version 1.4.0.1) & nco (version 3.9.5) ) : example of script
  3. create coefr weights for runoffs ( NOTE: in our case we created weights = 0, because max of runoff in COREII is < 50 Sv ) File runoff with 2 weight's coefficient (coefr0 is 0 averywhere, and coefr ......): https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/runoff_core_monthly.nc runoff_core_monthly.nc

NOTE: you can find interpolated files at IDRIS

  1. removed part "if ORCA2" in sbcrnf.F90
  2. preparing script ".pro" for post processing with idl https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/std_plots.pro.clima example for plots, and https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/std_ts.pro.clima example for time series
  3. automatize creation of pdf or html, by shell script https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/std_plots.sh.clima example for plots , and https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/time_series.sh.clima example for time series

and all detailed scripts for single : https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/plots.tar.gz plots and for single : https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/ts.tar.gz time series

NOTE : it is better to do not give name with "_" because to have timeseries the script looks for experience name and dates in filenames, and to split filename it uses "_" like separator.

Configuration of libIGCM files:

  1. COMP/opa9.card : https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/opa9.card opa9.card
  2. COMP/opa9.driver : https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/opa9.driver opa9.driver
  3. COMP/lim2.card : https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/lim2.card lim2.card
  4. PARAM/namelist : https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/namelist namelist
  5. PARAM/namelist_ice : https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/namelist_ice_lim2 namelist_ice_lim2
  6. PARAM/iodef.xml : https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/iodef.xml iodef.xml
  7. config.card : https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/config.card config.card

Plots (Climatology)

Job characteristics :

 forcing : COREII 
 2000 years of simulation
 strategy of output: 
         writing frequency of output : 10 years
         ocean output : all years annual mean output (T, U, V)
         ocean output : true monthly output for the following decades : 1-10, 41-50, 91-100, 491-500, 991-1000, 1491-1500, 1991-2000
         oce restart : 1 year every ten years 
         ice output : 1 month for all years 
         ice restart : 1 year every ten years 
         restart of the system saved each ten years

 keys list: key_trabbl_dif key_vectopt_loop key_vectopt_memory key_orca_r2 key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv
            key_dynldf_c3d key_dtatem key_dtasal key_tradmp key_trabbc key_zdftke key_zdfddm key_diahth key_iomput key_zdftmx
 (used ipcc's names, see table above)

Outputs & Restart on DODS

Reference is named "tmx_2000"

plots.pdf
timeseries.pdf

Note: experience named "core2000" is like reference but without key key_zdftmx

Supplementary plots :
plots :
Plots after 10 years : tmx2000_10y.pdf

Plots after 100 years : tmx2000_100y.pdf

Plots after 500 years : tmx2000_500y.pdf

Plots after 1000 years : tmx2000_1000y.pdf


Inter Annual files

NOTE: runoff for interannual is the same of climatological so be careful to change (if necessary) the name in opa9.card!!!!

  1. Download COREII forcing files : CIAF
  2. For Inter Annual case you can download "Combined years data sets" from 1948 to 2007 (60 years). Here there is a script to extract every year from global file : https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/script_extract.sh example of extract script
  3. Interpolation of forcing files (with cdo (version 1.4.0.1) & nco (version 3.9.5) ) : example of script

In this case I needed also a script to move "_year" at the end of file to work with Inter Annual forcing files, because in namelist you have to put your file name and NEMO add the year at the end of the file name

There are 2 cases:
3.1 If this script of example is used you remove part with date (04MAR2010)
3.2 Or you can leave part with date of file's creation

Here you have corresponding examples of script to move "_year" in filename :
3.1 https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/change_names.DDMMMYYYY.sh change_names.DDMMMYYYY.sh
3.2 https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/change_names.date.sh change_names.date.sh

NOTE: you can find interpolated files at
http://dodsp.idris.fr/reee605/CORE2_interan/ IDRIS

  1. removed part "if ORCA2" in sbcrnf.F90
  2. in fldread.F90 changed: "IF( llprev .AND. sdjf%num == 0 ) " in "IF( llprev .AND. sdjf%num <= 0 )"
  3. added in subbasins_orca21.nc mask for indopacific ocean and indian ocean without red sea.
  4. preparing script ".pro" for post processing with idl (same as climatological case)
  5. automatize creation of pdf or html, by shell script (same as climatological case)

Configuration of libIGCM files:

  1. COMP/opa9.card : https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/opa9_inter.card opa9.card
  2. COMP/opa9.driver : https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/opa9_inter.driver opa9.driver
  3. COMP/lim2.card : https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/lim2.card lim2.card
  4. PARAM/namelist : https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/namelist_inter namelist
  5. PARAM/namelist_ice : https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/namelist_ice_lim2 namelist_ice_lim2
  6. PARAM/iodef.xml : https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/iodef.xml iodef.xml
  7. config.card : https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/config_inter.card config.card

Plots (Interannual)

Job characteristics :

 forcing : COREII 
 60 years of simulation (repeted 3 times)
 strategy of output: 
             writing frequency of output : 10 years
             ocean output : all years annual mean output
             ocean output : true monthly mean output for the following decades : 1948-1957 and 1998-2007
             oce restart : 1 year every ten years 
             ice output : 1 month for all years 
             ice restart : 1 year every ten years restart of the system saved each ten years

 keys list: key_trabbl_dif key_vectopt_loop key_vectopt_memory key_orca_r2 key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv
            key_dynldf_c3d key_dtatem key_dtasal key_tradmp key_trabbc key_zdftke key_zdfddm key_diahth key_iomput key_zdftmx
 (used ipcc's names, see table above)

Outputs & Restart on DODS

Reference is named "tmxint3"

plots.pdf
time_series.pdf is time series of the third run of 60 years.

( and time series of difference between run3 and run2): ts_tmxint3_tmxint2.pdf

Note: experience named "core360i" is like reference but without key key_zdftmx

Supplementary plots :

Differences between run3 and run1 : tmxint3_tmxint1.pdf

Differences between run3 and run2 : tmxint3_tmxint2.pdf

Differences between run2 and run1 : tmxint2_tmxint1.pdf

& time series :
For the first run of 60 years : run1 ts_tmxint1_tmxint1.pdf

Differences between run2 and run1 : ts_tmxint2_tmxint1.pdf

Differences between run3 and run1 : ts_tmxint3_tmxint1.pdf



table of conversion of variable's names

Here is a table of conversion from old names to ipcc's names for variables used in these examples with COREII forcing files:

old name ipcc's name
votemper thetao
vosaline so
vozocrtx uoce_eff
vomecrty voce_eff
sosstsst tos
sohefldo qns+qsr
sowafldp wfcorr
somxl010 mldr10_1
sosaline sos
sossheig zos
soshfldo rsntds
sozotaux tauuo
sozotauy tauvo
iicethic sit
ice_cover sic
isnowthi snd



Simona Flavoni last update: Timestamp?

Attachments (39)