[[TOC(heading=Integrating COREII forcing in nemo-libigcm)]] = '''Integrating COREII forcing in nemo v3.2 / libigcm''' = [[BR]] == Purposes: == [[BR]] * Document the using of COREII forcing for NEMO (tag v32, revision : 1803), in configuration ORCA2_LIM, with libIGCM. * [http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6VPS-4TGHNK7-1&_user=1082852&_coverDate=12%2F31%2F2009&_rdoc=1&_fmt=high&_orig=browse&_srch=doc-info%28%23toc%236214%232009%23999739998%23739137%23FLA%23display%23Volume%29&_cdi=6214&_sort=d&_docanchor=&_ct=6&_acct=C000051401&_version=1&_urlVersion=0&_userid=1082852&md5=89e823a638c6e3c91a61dac716fa3eea Documentation of CORE :] [[BR]] [[BR]] == Step by step: == Preparing COREII forcing files on ORCA2 grid '''TO DO : '''[[BR]] '''understand how to sort rain without key_coupled (can we change iodef.xml file?)''' === climatological files === 1. Download COREII forcing files : [http://data1.gfdl.noaa.gov/nomads/forms/mom4/COREv2.html CYNF] [[BR]] 2. Interpolation of forcing files (with cdo (version 1.4.0.1) & nco (version 3.9.5) ) : [http://www.nemo-ocean.eu/Using-NEMO/Pre-and-post-processing-packages/Pre-Processing/cdo-interpolation 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/private/core2_cmip/runoff_core_monthly.nc runoff_core_monthly.nc]] NOTE: for the moment you can find interpolated files at [/u/rech/eee/reee605/CORE2_clima/OCE/ORCA2_LIM2 IDRIS] and [/dmnfs13/cont003/flavoni/CORE2_clima/OCE/ORCA2_LIM2 CCRT] [[BR]] 4. removed part "if ORCA2" in sbcrnf.F90 5. preparing script ".pro" for post processing with idl [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/std_plots.pro.clima example for plots]], and [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/std_ts.pro.clima example for time series]] 6. automatize creation of pdf or html, by shell script [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/std_plots.sh.clima example for plots ]], and [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/time_series.sh.clima example for time series ]] Configuration of libIGCM files: 1. COMP/opa9.card : [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/opa9.card opa9.card]] 2. COMP/opa9.driver : [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/opa9.driver opa9.driver ]] 3. PARAM/namelist : [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/namelist namelist ]] 4. PARAM/iodef.xml : [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/iodef.xml iodef.xml ]] 5. config.card : [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/config.card config.card ]] === Plots (Climatology) === Job characteristics : {{{ forcing : COREII 2000 years of simulation output: annual mean (T,U,V) monthly mean (ice) monthly mean (T,U,V) for year 50, 100, 500, 1000, 1500, 2000 (used ipcc's names, see table above) }}} [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/all_plots.pdf plots.pdf ]] [[BR]] [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/all_ts.pdf timeseries.pdf ]] [[BR]] === 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 : [http://data1.gfdl.noaa.gov/nomads/forms/mom4/COREv2/CIAF_v2.html CIAF] [[BR]] 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/private/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) ) : [http://www.nemo-ocean.eu/Using-NEMO/Pre-and-post-processing-packages/Pre-Processing/cdo-interpolation example of script] [[BR]] 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'''[[BR]] There are 2 cases: [[BR]] 3.1 If this script of example is used you remove part with date (04MAR2010) [[BR]] 3.2 Or you can leave part with date of file's creation [[BR]] Here you have corresponding examples of script to move "_year" in filename :[[BR]] 3.1 [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/change_names.DDMMMYYYY.sh change_names.DDMMMYYYY.sh]] [[BR]] 3.2 [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/change_names.date.sh change_names.date.sh]] [[BR]] NOTE: '''for the moment''' you can find interpolated files at [[BR]] [[/u/rech/eee/reee605/CORE2_interan_lastvers/OCE/ORCA2_LIM2 IDRIS]] (version 04MARS2010) [[BR]] [[/dmnfs13/cont003/flavoni/CORE2_interan_lastvers/OCE/ORCA2_LIM2 CCRT]] (version 04MARS2010) [[BR]] ( note: it is also available previous version : [[BR]] [[/u/rech/eee/reee605/CORE2_interan/OCE/ORCA2_LIM2 IDRIS: version 26JAN2010]] [[BR]] [[/dmnfs13/cont003/flavoni/CORE2_interan/OCE/ORCA2_LIM2 CCRT: version 26JAN2010]] [[BR]] 4. removed part "if ORCA2" in sbcrnf.F90 5. in fldread.F90 changed: "IF( llprev .AND. sdjf%num == 0 ) " in "IF( llprev .AND. sdjf%num <= 0 )" 6. added in subbasins_orca21.nc mask for indopacific ocean and indian ocean without red sea. 7. preparing script ".pro" for post processing with idl (still not decide all plots) 8. automatize creation of pdf or html, by shell script Configuration of libIGCM files: 1. COMP/opa9.card : [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/opa9_inter.card opa9.card]] 2. COMP/opa9.driver : [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/opa9_inter.driver opa9.driver ]] 3. PARAM/namelist : [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/namelist_inter namelist ]] 4. PARAM/iodef.xml : [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/iodef.xml iodef.xml ]] 5. config.card : [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/config_inter.card config.card ]] === Plots (Interannual) === Job characteristics : {{{ output: annual mean (T,U,V) monthly mean (ice) (used ipcc's names, see table above) }}} For the first run of 60 years (core2_run1) [[BR]] [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/timeseries_interan_run1.pdf timeseries_interan_run1.pdf ]] For the second run of 60 years (core2_run2) [[BR]] [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/timeseries_interan_run2.pdf timeseries_interan_run2.pdf ]] For the third run of 60 years (core2_run3) [[BR]] [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/timeseries_interan_run3.pdf timeseries_interan_run3.pdf ]] [[BR]] [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/plots_interan_run3.pdf plots_interan_run3.pdf ]] (RUN3 - RUN1) [[BR]] [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/private/core2_cmip/timeseries_interan_run3_1.pdf timeseries_interan_run3_1.pdf ]] [[BR]] [[BR]] === 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 || [[BR]] [[BR]] Simona Flavoni last update: [[Timestamp]]