[[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). * [http://data1.gfdl.noaa.gov/nomads/forms/core/COREv2.html Documentation of CORE :] == 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]] 1. 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] 1. 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 file] NOTE: you can find interpolated files at [http://prodn.idris.fr/thredds/catalog/ipsl_public/romr009/CORE2_clima/catalog.html IDRIS][[BR]] 4. removed part "if ORCA2" in sbcrnf.F90 4. 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] 4. 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] 1. COMP/opa9.driver : [https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/opa9.driver opa9.driver] 1. COMP/lim2.card : [https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/lim2.card lim2.card] 1. PARAM/namelist : [https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/namelist namelist] 1. PARAM/namelist_ice : [https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/namelist_ice_lim2 namelist_ice_lim2] 1. PARAM/iodef.xml : [https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/iodef.xml iodef.xml] 1. config.card : [https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/config.card config.card] ''' NOTE : be careful to opa9.driver & iodef.xml file''' [[BR]] ''' in opa9.driver there is a specific sed in iodef.xml to have monthly output only for specifics decades''' [[BR]] ''' in iodef.xml file there is : [[BR]] ''''' ''' ''''' ' === 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 '''[http://prodn.idris.fr/thredds/catalog/ipsl_public/romr009/REFERENCE_CLIMATO/catalog.html web page] Reference is named "tmx_2000" [[BR]] [http://prodn.idris.fr/thredds/fileServer/ipsl_public/romr009/plots_CORE2/CLIMATO/plots_pdf/tmx2000_2000y.pdf plots.pdf] [[BR]] [http://prodn.idris.fr/thredds/fileServer/ipsl_public/romr009/ts_CORE2/CLIMATO/ts_tmx2000.pdf timeseries.pdf] '''Note: experience named "core2000" is like reference but without key key_zdftmx '''[[BR]] '''Supplementary plots''' :[[BR]][http://dodsp.idris.fr/romr009/plots_CORE2/CLIMATO/plots_pdf plots] :[[BR]] Plots after 10 years : [http://dodsp.idris.fr/romr009/plots_CORE2/CLIMATO/plots_pdf/tmx2000_10y.pdf tmx2000_10y.pdf] [[BR]] Plots after 100 years : [http://dodsp.idris.fr/romr009/plots_CORE2/CLIMATO/plots_pdf/tmx2000_100y.pdf tmx2000_100y.pdf] [[BR]] Plots after 500 years : [http://dodsp.idris.fr/romr009/plots_CORE2/CLIMATO/plots_pdf/tmx2000_500y.pdf tmx2000_500y.pdf] [[BR]] Plots after 1000 years : [http://dodsp.idris.fr/romr009/plots_CORE2/CLIMATO/plots_pdf/tmx2000_1000y.pdf tmx2000_1000y.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]] 1. 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] 1. 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/romr009/CORE2_interan/ IDRIS] [[BR]] 4. removed part "if ORCA2" in sbcrnf.F90 4. in fldread.F90 changed: "IF( llprev .AND. sdjf%num == 0 ) " in "IF( llprev .AND. sdjf%num <= 0 )" 4. 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) 4. 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] 1. COMP/opa9.driver : [https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/opa9_inter.driver opa9.driver] 1. COMP/lim2.card : [https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/lim2.card lim2.card] 1. PARAM/namelist : [https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/namelist_inter namelist] 1. PARAM/namelist_ice : [https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/namelist_ice_lim2 namelist_ice_lim2] 1. PARAM/iodef.xml : [https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/core2_cmip/iodef.xml iodef.xml] 1. 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 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 '''[http://dodsp.idris.fr/romr009/REFERENCE_INTERAN/ DODS] Reference is named "tmxint3" [http://dodsp.idris.fr/reee605/plots_CORE2/INTERAN/plots_pdf/tmx_int3.pdf plots.pdf] [[BR]] [http://dodsp.idris.fr/reee605/ts_CORE2/INTERAN/ts_tmxint3_tmxint3.pdf time_series.pdf] is time series of the third run of 60 years.[[BR]] ( and time series of difference between run3 and run2): [http://dodsp.idris.fr/reee605/ts_CORE2/INTERAN/ts_tmxint3_tmxint2.pdf ts_tmxint3_tmxint2.pdf] [[BR]] '''Note: experience named "core360i" is like reference but without key key_zdftmx'''[[BR]] '''Supplementary plots ''' :[[BR]] Differences between run3 and run1 : [http://dodsp.idris.fr/romr009/plots_CORE2/INTERAN/plots_pdf/tmxint3_tmxint1.pdf tmxint3_tmxint1.pdf] [[BR]] Differences between run3 and run2 : [http://dodsp.idris.fr/romr009/plots_CORE2/INTERAN/plots_pdf/tmxint3_tmxint2.pdf tmxint3_tmxint2.pdf] [[BR]] Differences between run2 and run1 : [http://dodsp.idris.fr/romr009/plots_CORE2/INTERAN/plots_pdf/tmxint2_tmxint1.pdf tmxint2_tmxint1.pdf] [[BR]] ''' & time series ''': [[BR]] For the first run of 60 years : run1 [http://dodsp.idris.fr/romr009/ts_CORE2/INTERAN/ts_tmxint1_tmxint1.pdf ts_tmxint1_tmxint1.pdf] Differences between run2 and run1 : [http://dodsp.idris.fr/romr009/ts_CORE2/INTERAN/ts_tmxint2_tmxint1.pdf ts_tmxint2_tmxint1.pdf] Differences between run3 and run1 : [http://dodsp.idris.fr/romr009/ts_CORE2/INTERAN/ts_tmxint3_tmxint1.pdf 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]]