[[TOC(heading=Integrating COREII forcing in nemo-libigcm)]] = '''Integrating COREII forcing in nemo tags v3.2 / libigcm tags v1.4''' = [[BR]] == Purposes: == [[BR]] * Document the using of COREII forcing for NEMO (tag v32, revision : 1803), in configuration ORCA2_LIM, with libIGCM (tags v1.4, revision : 200). * [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 === 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/core2_cmip/runoff_core_monthly.nc runoff_core_monthly.nc]] NOTE: you can find interpolated files at [http://dodsp.idris.fr/reee605/CORE2_clima/ IDRIS] [[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/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]] 6. 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 all single : [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/plots.tar.gz plots ]] and for all single : [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/ts.tar.gz time series ]] 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/iodef.xml : [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/iodef.xml iodef.xml ]] 6. 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) }}} [http://dodsp.idris.fr/reee605/plots_CORE2/CLIMATO/plots_pdf/core2000.pdf plots.pdf ] [[BR]] [http://dodsp.idris.fr/reee605/ts_CORE2/CLIMATO/clima_ts.pdf timeseries.pdf ] '''Supplementary plots''' :[[BR]] PLOTS:[[BR]] Plots after 10 years : [http://dodsp.idris.fr/reee605/plots_CORE2/CLIMATO/plots_pdf/core2000_10y.pdf core2000_10y.pdf] [[BR]] Plots after 100 years : [http://dodsp.idris.fr/reee605/plots_CORE2/CLIMATO/plots_pdf/core2000_100y.pdf core2000_100y.pdf] [[BR]] Plots after 500 years : [http://dodsp.idris.fr/reee605/plots_CORE2/CLIMATO/plots_pdf/core2000_500y.pdf core2000_500y.pdf] [[BR]] Plots after 1000 years : [http://dodsp.idris.fr/reee605/plots_CORE2/CLIMATO/plots_pdf/core2000_1000y.pdf core2000_1000y.pdf] [[BR]] [[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/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/core2_cmip/change_names.DDMMMYYYY.sh change_names.DDMMMYYYY.sh]] [[BR]] 3.2 [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/change_names.date.sh change_names.date.sh]] [[BR]] NOTE: you can find interpolated files at [[BR]] [[http://dodsp.idris.fr/reee605/CORE2_interan/ IDRIS]] [[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 (same as climatological case) 8. 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/iodef.xml : [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/iodef.xml iodef.xml ]] 6. 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) }}} Plots : [[BR]] [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/plots_interan_run3.pdf plots_interan_run3.pdf ]] [[BR]] Time Series : [[BR]] For the third run of 60 years ( and differences between run3 and run2) [[BR]] [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/timeseries_interan_run3.pdf timeseries_interan_run3.pdf ]] '''Supplementary plots & time series''' :[[BR]] PLOTS:[[BR]] Differences between run3 and run1 :[[BR]] [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/plots_interan_run31.pdf plots_interan_run31.pdf ]] [[BR]] Differences between run3 and run2 :[[BR]] [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/plots_interan_run32.pdf plots_interan_run32.pdf ]] [[BR]] Differences between run2 and run1 :[[BR]] [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/plots_interan_run21.pdf plots_interan_run21.pdf ]] [[BR]] TIME SERIES: [[BR]] For the first run of 60 years : run1 [[BR]] [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/timeseries_interan_run1.pdf timeseries_interan_run1.pdf ]] Differences between run2 and run1 [[BR]] [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/timeseries_interan_run2.pdf timeseries_interan_run2.pdf ]] Differences between run3 and run1 [[BR]] [[https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/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]]