Differences between revisions 5 and 16 (spanning 11 versions)
Revision 5 as of 2015-11-24 12:40:39
Size: 3365
Editor: NicoleThomas
Comment:
Revision 16 as of 2023-03-20 11:55:32
Size: 4244
Editor: NicoleThomas
Comment:
Deletions are marked like this. Additions are marked like this.
Line 1: Line 1:
= Randersetzung in MESSy-CLaMS = = Update of the boundary conditions in MESSy-CLaMS =
Line 3: Line 3:
Alle Randersetzungen werden im Unterprogramm '''clamsbmix_global_end''' (messy_clamsbmix_si.f90) aufgerufen.

Die Ersetzungen im CLaMS-Paket BMIX wurden in die Messy-Routine '''bmix''' übernommen. Die Randersetzungen erfolgen in der in MIX/BMIX festgelegten untersten bzw. obersten Schicht.
 
Die Randersetzungen aus ''replace_lbound'', ''pos_chem'' und ''pos_dyn'' werden in der Routine '''replace_boundaries''' vorgenommen (in messy_clamsbmix_replace_bounds). Diese Routine wird in ''clamsbmix_global_end'' nach BMIX aufgerufen. Die Ersetzungen erfolgen in einem vorgegebenen Bereich (unabhängig von den in MIX/BMIX verwendeten Schichten).
The routines for updating the boundary conditions are called in ''clamsbmix_global_end'' (''messy_clamsbmix_si.f90'').
Line 10: Line 6:
Angaben in '''clamsbmix.nml''':
{{{
There are two routines for updating the boundary conditions:

 * '''bmix''': replace lowest, highest and/or vertical layers (layers used for mixing)
 * '''replace_boundaries''': set boundary conditions below or above a specified zeta-level and within a specified period


=== Specifications in clamsbmix.nml ===
 {{{
&CTRL
timestep_bmix=24, ! timestep of bmix in hours
lev_in_down=50., ! lower boundary
lev_in_up=2500., ! upper boundary
lat_in_down=-90., ! south vertical boundary
lat_in_up=90., ! north vertical boundary
delta_lev=0.,
interpol_from_init=0,
Line 14: Line 24:
bmix_boundlist='/usr/nfs/local/messy_config/bounds_bmix_clim.inp',
clams_boundlist='/usr/nfs/local/messy_config/bounds_clams_clim.inp'
bmix_boundlist='/usr/nfs/software/icg1/messy_config/bounds_bmix_clim.inp',
clams_boundlist='/usr/nfs/software/icg1/messy_config/bounds_clams_clim.inp'
loutput_bmix=.false.,
/
Line 18: Line 30:
 * ''file_bounds'': Boundfile/Init-File für BMIX
 * ''dir_boundfiles'': Verzeichnis mit Boundfiles (mit Ausnahme des Init-Files)
 * ''bmix_boundlist'': Datei, in der die in BMIX zu ersetzenden Spezies aufgeführt sind
 * ''clams_boundlist'': Datei, in der alle weiteren Randersetzungen (in replace_boundaries) aufgelistet sind
 * ''file_bounds'': boundfile/init file used in ''bmix''
 * ''dir_boundfiles'': directory with all other boundfiles used in ''bmix'' or ''replace_boundaries''
 * ''bmix_boundlist'': file with specifications for all species replaced in ''bmix''
 * ''clams_boundlist'': file with specifications for all species replaced in ''replace_boundaries''
Line 25: Line 37:

Spezifikation der Randersetzungen für '''bmix''' in '''bmix_boundlist''' (z.B. ''bounds_bmix_clim.inp''):
{{{
=== Specifications for bmix in bmix_boundlist (e.g. bounds_bmix_clim.inp) ===
 {{{
ALL 1 1 chem_bound_0050K.nc
ALL 2 1 chem_bound_2500K.nc
Line 30: Line 43:
CO2 1 1 chem_bound_noaa_cmdl.nc CO2 1 1 chem_bound_noaa_cmdl.nc   1970 2000
Line 35: Line 48:
}}}

Angaben pro Zeile:

 * Spezies
 * Rand:
  * 1 = Unterrand
  * 2 = Oberrand
  * 3 = Unter- und Oberrand
 * Aktion:
  * 1 = interpoliere + ersetze
  * 2 = interpoliere + addiere
 * Boundfile: Datei, aus der die Ränder ersetzt werden. Die Datei muss im angegebenen Verzeichnis (dir_boundfiles) liegen.
 }}}
Line 50: Line 51:
Spezifikation der Randersetzungen für '''replace_boundaries''' in '''clams_boundlist''' (z.B. ''bounds_clams_clim.inp'') :
{{{
specifications per line:

 * name of species; if '''''ALL''''' is specified, all species from boundfile are read in
 * boundary, which should be replaced:
     * 1 = lower boundary
     * 2 = upper boundary
     * 3 = lower and upper boundary
 * action:
     * 1 = interpolate + replace
     * 2 = interpolate + add
 * boundfile: file with boundary conditions (in directory ''dir_boundfiles'')
 * optional: startyear and endyear: replacement of boundaries only between start- and endyear
 

=== Specifications for replace_boundaries in clams_boundlist (e.g. bounds_clams_clim.inp) ===
 {{{
Line 59: Line 74:
}}}  }}}
Line 61: Line 76:
Angaben pro Zeile: specifications per line:
Line 63: Line 78:
   * Speziesname in CLaMS
   * Speziesname im Boundfile
   * Unterrand in ZETA
   * Aktion am Unterrand
     * 0: keine Randersetzung
     * 1: interpolieren + ersetzen
     * 9: setze auf 0.
   * Oberrand in ZETA
   * Aktion am Oberrand
     * 0: keine Randersetzung
     * 1: interpolieren + ersetzen
     * 9: setze auf 0.
   * Startjahr und Endjahr: nur im vorgegebenen Zeitraum erfolgt eine Ersetzung der Ränder
   * Boundfile: Die Datei muss im angegebenen Verzeichnis (dir_boundfiles) liegen
   * speciesname in CLaMS
   * speciesname in boundfile
   * lower boundary [ZETA]
   * action at lower boundary
     * 0: no replacement
     * 1: interpolate + replace
     * 2: interpolate * add
     * 9: set to 0.
   * upper boundary [ZETA]
   * action at upper boundary
     * 0: no replacement
     * 1: interpolate + replace
     * 2: interpolate + add
     * 9: set to 0.
   * startyear and endyear: replacement of boundaries only between start- and endyear
   * boundfile: file with boundary conditions (in directory ''dir_boundfiles'')
     * '''AIRS/mopitt''' files: "_yyyy010112_yyyy123112.nc" will be attached to filename
     * '''METDATA''':
          * meteorological data is used
          * filename: ''met_prefix''_yymmddhh.nc (''met_prefix'' is specified in clams.nml)
          * file must be located in ''met_dir'' (specified in clams.nml)
     * possible dimensioning of species in boundfile:
        * species (lat,time)
        * species (lat,lon,time)
        * species (lat/eqlat,zeta/theta,month)

Update of the boundary conditions in MESSy-CLaMS

The routines for updating the boundary conditions are called in clamsbmix_global_end (messy_clamsbmix_si.f90).

There are two routines for updating the boundary conditions:

  • bmix: replace lowest, highest and/or vertical layers (layers used for mixing)

  • replace_boundaries: set boundary conditions below or above a specified zeta-level and within a specified period

Specifications in clamsbmix.nml

  • &CTRL
    timestep_bmix=24,     ! timestep of bmix in hours
    lev_in_down=50.,      ! lower boundary 
    lev_in_up=2500.,      ! upper boundary 
    lat_in_down=-90.,     ! south vertical boundary
    lat_in_up=90.,        ! north vertical boundary
    delta_lev=0.,
    interpol_from_init=0, 
    file_bounds='/usr/nfs/data/clams/clim/init/init_KY2_79010112.nc',
    dir_boundfiles='/usr/nfs/data/clams/clim/config/',
    bmix_boundlist='/usr/nfs/software/icg1/messy_config/bounds_bmix_clim.inp',
    clams_boundlist='/usr/nfs/software/icg1/messy_config/bounds_clams_clim.inp'
    loutput_bmix=.false.,
    /
  • file_bounds: boundfile/init file used in bmix

  • dir_boundfiles: directory with all other boundfiles used in bmix or replace_boundaries

  • bmix_boundlist: file with specifications for all species replaced in bmix

  • clams_boundlist: file with specifications for all species replaced in replace_boundaries

Specifications for bmix in bmix_boundlist (e.g. bounds_bmix_clim.inp)

  • ALL   1  1  chem_bound_0050K.nc
    ALL   2  1  chem_bound_2500K.nc
    CH4   1  1  chem_bound_noaa_cmdl.nc
    CO    1  1  chem_bound_noaa_cmdl.nc
    CO2   1  1  chem_bound_noaa_cmdl.nc   1970 2000
    N2O   1  1  chem_bound_noaa_cmdl.nc
    F11   1  1  chem_bound_noaa_cmdl.nc 
    F12   1  1  chem_bound_noaa_cmdl.nc 
    HCN   1  2  hcn_boundary_2_0_deg_50p_2000010112_2007123112.nc

specifications per line:

  • name of species; if ALL is specified, all species from boundfile are read in

  • boundary, which should be replaced:
    • 1 = lower boundary
    • 2 = upper boundary
    • 3 = lower and upper boundary
  • action:
    • 1 = interpolate + replace
    • 2 = interpolate + add
  • boundfile: file with boundary conditions (in directory dir_boundfiles)

  • optional: startyear and endyear: replacement of boundaries only between start- and endyear

Specifications for replace_boundaries in clams_boundlist (e.g. bounds_clams_clim.inp)

  • CO       CO    200    1  0      0  2001  2011  mopitt_chem_bound_2_0_0200K     
    CO       CO    200    1  0      0  2012  2015  AIRS_chem_bound_CO_1_5_0200K   
    CH4      CH4   200    1  0      0  2012  2015  AIRS_chem_bound_CH4_0200K      
    HCl      HCl   0.     0  500.   1  1970  2020  haloe_eqlat_climat_theta_update.nc      
    O3       O3    0.     0  500.   1  1970  2020  haloe_eqlat_climat_theta_update.nc      
    TRACER   O3    0.     0  500.   1  1970  2020  haloe_eqlat_climat_theta_update.nc      
    H2O      SH    250.   1  0      0  1970  2020  METDATA       

specifications per line:

  • speciesname in CLaMS
  • speciesname in boundfile
  • lower boundary [ZETA]
  • action at lower boundary
    • 0: no replacement
    • 1: interpolate + replace
    • 2: interpolate * add
    • 9: set to 0.
  • upper boundary [ZETA]
  • action at upper boundary
    • 0: no replacement
    • 1: interpolate + replace
    • 2: interpolate + add
    • 9: set to 0.
  • startyear and endyear: replacement of boundaries only between start- and endyear
  • boundfile: file with boundary conditions (in directory dir_boundfiles)

    • AIRS/mopitt files: "_yyyy010112_yyyy123112.nc" will be attached to filename

    • METDATA:

      • meteorological data is used
      • filename: met_prefix_yymmddhh.nc (met_prefix is specified in clams.nml)

      • file must be located in met_dir (specified in clams.nml)

    • possible dimensioning of species in boundfile:
      • species (lat,time)
      • species (lat,lon,time)
      • species (lat/eqlat,zeta/theta,month)

messy/bmix (last edited 2023-03-20 11:55:56 by NicoleThomas)