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 122 (modified by flavoni, 9 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 ......): runoff file

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 example for plots, and example for time series
  3. automatize creation of pdf or html, by shell script example for plots, and example for time series

and all detailed scripts for single : plots and for single : 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 : opa9.card
  2. COMP/opa9.driver : opa9.driver
  3. COMP/lim2.card : lim2.card
  4. PARAM/namelist : namelist
  5. PARAM/namelist_ice : namelist_ice_lim2
  6. PARAM/iodef.xml : iodef.xml
  7. config.card : config.card

NOTE : be careful to opa9.driver & iodef.xml file
in opa9.driver there is a specific sed in iodef.xml to have monthly output only for specifics decades
in iodef.xml file there is :
<group id="1m" output_freq="-1" output_level="10" enabled=".FALSE.">
'

Plots (Climatology)

Job characteristics :

 forcing : COREII 
 2000 years of simulation
 strategy of output: 
         writing frequency of output file : 10 years (one file contains 10 years)
         ocean output frequency: all years annual mean output (T, U, V)
         ocean output frequency: true monthly output for the following decades : 1-10, 41-50, 91-100, 491-500, 991-1000, 1491-1500, 1991-2000
         ice output frequency: true monthly output for all 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 here

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 : 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 change_names.DDMMMYYYY.sh
3.2 change_names.date.sh

NOTE: you can find interpolated files at
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 : opa9.card
  2. COMP/opa9.driver : opa9.driver
  3. COMP/lim2.card : lim2.card
  4. PARAM/namelist : namelist
  5. PARAM/namelist_ice : namelist_ice_lim2
  6. PARAM/iodef.xml : iodef.xml
  7. config.card : config.card

Plots (Interannual)

Job characteristics :

 forcing : COREII 
 60 years of simulation (repeted 3 times)
 strategy of output: 
             writing frequency of output file : 10 years (one file contains 10 years)
             ocean output frequency: all years annual mean output (T, U, V)
             ocean output frequency: true monthly output for the following decades : 1948-1957 and 1998-2007
             ice output frequency: true monthly output for all 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 here

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
sowaflup wfo
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)