Differences between revisions 1 and 10 (spanning 9 versions)
⇤ ← Revision 1 as of 2014-12-08 11:50:31 →
Size: 2564
Editor: NicoleThomas
Comment:
← Revision 10 as of 2015-12-11 12:28:14 → ⇥
Size: 9593
Editor: NicoleThomas
Comment:
Deletions are marked like this. Additions are marked like this.
Line 1: Line 1:
= MESSy/CLaMS: Dimensionen und Representationen für Channels =


== Allgemeines ==

 * Channel-Objekte sind immer vom Typ ''REAL(DP)''

 * Verwendete parallele Zerlegungstypen (''dctype'')
  * DC_BC = 0 ! BROADCAST (IO-PE ONLY) -> keine Zerlegung		 = MESSy/CLaMS: dimensions and representations =


 * The type of channel objects is ''REAL(DP)''

 * Used parallel decomposition types (''dctype'')
  * DC_BC = 0 ! BROADCAST (IO-PE ONLY) -> no decomposition
Line 12: Line 10:
== Deklaration der Dimensionen und Representationen ==

In ''bmil/messy_main_channel_bi.f90''

==== Dimensionen ====		 == Deklaration of dimensions and representations ==

In ''bmil/messy_main_channel_bi.f90'':

==== Dimensions ====
Line 23: Line 21:

  INTEGER, SAVE, PUBLIC :: DIMID_THETA = DIMID_UNDEF
Line 24: Line 24:
  INTEGER, SAVE, PUBLIC :: DIMID_MIX_GRID = DIMID_UNDEF
Line 26: Line 27:
}}}

==== Representationen ====		   INTEGER, SAVE, PUBLIC :: DIMID_SEDI_PARTICLE = DIMID_UNDEF

  INTEGER, SAVE, PUBLIC :: DIMID_WAVE = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_DISSOC_LEV = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_DISSOC_LEVC = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_DISSOC_LAT = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_DISSOC_SZA = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_DISSOC_NUMJ = DIMID_UNDEF
}}}

==== Representations ====
Line 34: Line 43:
  INTEGER, SAVE, PUBLIC :: REPR_3DINP_CLAMS = REPR_UNDEF   INTEGER, SAVE, PUBLIC :: REPR_3DINP_CLAMS = REPR_UNDEF 
Line 38: Line 47:
}}}

== Initialisieren der Dimensionen und Representationen ==

In ''bmil/messy_main_channel_clams.inc''

==== Dimensionen ====

Initialisieren mit ''new_dimension''

|| ID || Name || Länge ||		   INTEGER, SAVE, PUBLIC :: REPR_SEDI_PARTICLE = REPR_UNDEF

  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_WAVE = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_LEV = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_LEVC = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_LAT = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_SZA = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_TABS = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_TABS_DAVG = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_TWOD = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_TWODC = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_RATE = REPR_UNDEF
}}}

== Initialize dimensions and representations ==

In ''bmil/messy_main_channel_clams.inc'':

==== Dimensions ====

Initialize with ''new_dimension''

|| ID || Name || Length ||
Line 53: Line 74:
|| DIMID_THETA || THETA || ntheta ||
Line 55: Line 77:
|| DIMID_MIX_GRID || NLEVS || adapt_par%nlevs ||
Line 56: Line 79:
|| DIMID_WAVE || WAVELENGTH || jpschi || || DIMID_WAVE || WAVELENGTH || jpwave ||
Line 62: Line 85:


==== Dimensionsvariablen ====

Anlegen mit ''add_dimension_variable''

||Name der Dimension || Name der Dimensionsvariablen || Werte der Dimensionsvariablen ||
||time|| time || 0.0 ||
||time|| YYYYMMDD || yyyymmdd ||
||time|| dt || delta_time ||
||time|| nstep || REAL(current_time_step, DP) ||
||LON || LON || longrid ||
||LAT || LAT || latgrid ||
||LEV || LEV || levelgrid(1:nz) ||


==== Representationen ====

Anlegen mit ''new_representation''		 || DIMID_SEDI_PARTICLE|| SEDI_PARTICLE || nparticle_max ||


==== Dimension variables ====

Create with ''add_dimension_variable''

||Name of dimension variable|| Name of dimension || values of dimension variable ||
|| time || time || days since starttime ||
|| YYYYMMDD || time || yyyymmdd ||
|| dt || time || delta_time ||
|| nstep || time || REAL(current_time_step, DP) ||
|| JULIAN_SECONDS|| time || time in julian seconds ||
|| LON || LON || longrid ||
|| LAT || LAT || latgrid ||
|| LEV || LEV || levelgrid(1:nz) ||
|| THETA || THETA|| thetagrid(1:ntheta) ||


==== Representations ====

Create with ''new_representation''

 * '''SCALAR'''
 {{{
CALL new_representation(status, SCALAR, 'SCALAR' &
     , rank = 0, link = '----', dctype = 0 &
     , dimension_ids = scdim &
     , ldimlen = scdim &
     , axis = '----' &
     )
}}}
 Variables with representation SCALAR:
  * In ''messy_clams_si.f90''
   * '''grid_switch'''
   * '''pre_year''', '''pre_month''', '''pre_day''', '''pre_sec'''
   * '''JULTIME'''

 * '''REPR_NTASKS'''
 {{{
CALL new_representation(status, REPR_NTASKS, 'REPR_NTASKS' &
     , rank = 1, link = 'x---', dctype = DC_IX &
     , dimension_ids = (/ DIMID_NTASKS /) &
     , ldimlen = (/ 1 /) & INTEGER, SAVE, PUBLIC :: REPR_SEDI_TBOX = REPR_UNDEF
 
     , axis = 'N---' &
     )
}}}
  * Variables with representation REPR_NTASKS in ''messy_clams_si.f90'':
    * '''dnparts'''
    * '''dnparts_max'''
  * one-dimensional array of length NTASKS
  * index distributed on tasks (dctype=DC_IX), every task can use ''dnparts'' and ''dnparts_max'' as a scalar

 * '''REPR_3DINP_CLAMS'''
 {{{
CALL new_representation(status, REPR_3DINP_CLAMS, 'REPR_3DINP_CLAMS' &
     , rank = 3, link = 'xxx-', dctype = 0 &
     , dimension_ids = (/ DIMID_LON, DIMID_LAT, DIMID_LEV /) &
     , ldimlen = (/ AUTO, AUTO, AUTO /) &
     , axis = 'XYZ-' &
     )
}}}
  * Variablen mit dieser Representation in ''messy_clamstraj_si.f90'':
   * '''UDT, VDT, WDT, LEVELDT'''
   * '''UFUT, VFUT, WFUT, LEVELFUT'''
  * 3-dimensionale Felder
  * Dimensionierung (nlon,nlat,nlev) wie sie innerhalb des TRAJ genutzt wird
  * Die Windfelder müssen auf allen Tasks vollständig vorhanden sein => keine parallele Zerlegung der Felder (dctype=0)

 * '''REPR_3DINP_CLAMS2'''
 {{{
CALL new_representation(status, REPR_3DINP_CLAMS2, 'REPR_3DINP_CLAMS2' &
     , rank = 3, link = 'xxx-', dctype = 0 &
     , dimension_ids = (/ DIMID_LEV, DIMID_LAT, DIMID_LON /) &
     , ldimlen = (/ AUTO, AUTO, AUTO /) &
     , axis = 'ZYX-' &
     )
}}}
  * Variablen mit dieser Representation in ''messy_clamstraj_si.f90'':
    * '''PREDATA_param, FUTDATA_param''' (für TEMP,PRESS etc)
  * 3-dimensionale Felder
  * Dimensionierung wie auf Netcdf-Winddateien (nlev,nlat,nlon)
  * Die Felder müssen auf allen Tasks vollständig vorhanden sein => keine parallele Zerlegung der Felder (dctype=0)


 * '''REPR_LG_CLAMS'''
 {{{
CALL new_representation(status, REPR_LG_CLAMS, 'REPR_LG_CLAMS' &
     , rank = 1, link = 'x---', dctype = DC_IX &
     , dimension_ids = (/ DIMID_TRAJ /) &
     , ldimlen = (/ dnparts_max /) &
     , axis = 'N---' &
     )
}}}
  * Variablen mit dieser Representation in ''messy_clams_si.f90'':
    * '''SPECARR''' für alle Spezies (1-nspec)
  * Variablen mit dieser Representation in ''messy_clamstraj_si.f90'':
    * '''JULSEC'''
    * '''LAT, LAT_OLD, LAT_OLD_MIX'''
    * '''LON, LON_OLD, LON_OLD_MIX'''
    * '''LEV, LEV_OLD'''
    * '''param, param_OLD ''' (PRESS,TEMP etc)
  * Variablen mit dieser Representation in ''messy_dissoc_si.f90'':
    * '''jresult'''
  * Eindimensionale Felder der Länge '''nparts_max'''
  * Index-Zerlegung (dctype=DC_IX) => jede Task erhält ein eindimensionales Feld der Länge '''dnparts_max'''

 * '''REPR_3DINP_CLAMSTHETA'''
 {{{
CALL new_representation(status, REPR_3DINP_CLAMSTHETA, 'REPR_3DINP_CLAMSTHETA' &
     , rank = 3, link = 'xxx-', dctype = 0 &
     , dimension_ids = (/ DIMID_THETA, DIMID_LAT, DIMID_LON /) &
     , ldimlen = (/ AUTO, AUTO, AUTO /) &
     , axis = 'ZYX-' &
     )
}}}

 * '''REPR_CLAMS_SHUFFLE'''
 {{{
CALL new_representation(status, REPR_CLAMS_SHUFFLE, 'REPR_CLAMS_SHUFFLE' &
     , rank = 1, link = 'x---', dctype = DC_IX &
     , dimension_ids = (/ DIMID_SHUFFLE /) &
     , ldimlen = (/ dnparts_max_shuffle /) &
     , axis = 'N---' &
     )
}}}
  * Variablen mit dieser Representation in ''messy_clamsmix_si.f90'':
    * '''ZETA_SHUFFLED'''
    * '''LAT_SHUFFLED'''
    * '''LON_SHUFFLED'''
    * '''LAT_OLD_MIX_SHUFFLED'''
    * '''LON_OLD_MIX_SHUFFLED'''
    * '''param_SHUFFLED '''(für alle Spezies)
  * Eindimensionale Felder der Länge '''nparts_max_shuffle'''
  * Index-Zerlegung (dctype=DX_IX) => jede Task erhält ein eindimensionales Feld der Länge '''dnparts_max_shuffle'''


 * '''REPR_MIX_GRID'''
 {{{
CALL new_representation(status, REPR_MIX_GRID, 'REPR_MIX_GRID' &
     , rank = 1, link = 'x---', dctype = 0 &
     , dimension_ids = (/ DIMID_MIX_GRID /) &
     , ldimlen = (/ adapt_par%nlevs /) &
     , axis = 'N---' &
     )
}}}


 * REPR_DISSOC_WAVE
 * REPR_DISSOC_LEV
 * REPR_DISSOC_LEVC
 * REPR_DISSOC_SZA
 * REPR_DISSOC_LAT
 * REPR_DISSOC_TABS
 * REPR_DISSOC_TABS_DAVG
 * REPR_DISSOC_TWOD
 * REPR_DISSOC_TWODC
 * REPR_DISSOC_RATE
 * REPR_SEDI_PARTICLE

MESSy/CLaMS: dimensions and representations

  • The type of channel objects is REAL(DP)

  • Used parallel decomposition types (dctype)

    • DC_BC = 0 ! BROADCAST (IO-PE ONLY) -> no decomposition

    • DC_IX = 4 ! INDEX

Deklaration of dimensions and representations

In bmil/messy_main_channel_bi.f90:

Dimensions

  •   INTEGER, SAVE, PUBLIC :: DIMID_TIME = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_LEV  = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_LON  = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_LAT  = DIMID_UNDEF

  INTEGER, SAVE, PUBLIC :: DIMID_THETA   = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_TRAJ    = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_MIX_GRID = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_NTASKS  = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_SHUFFLE = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_SEDI_PARTICLE = DIMID_UNDEF

  INTEGER, SAVE, PUBLIC :: DIMID_WAVE        = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_DISSOC_LEV  = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_DISSOC_LEVC = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_DISSOC_LAT  = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_DISSOC_SZA  = DIMID_UNDEF
  INTEGER, SAVE, PUBLIC :: DIMID_DISSOC_NUMJ = DIMID_UNDEF

Representations

  •   INTEGER, SAVE, PUBLIC :: REPR_LG_CLAMS     = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_CLAMS_SHUFFLE= REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_LEVARR_CLAMS = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_3DINP_CLAMS  = REPR_UNDEF 
  INTEGER, SAVE, PUBLIC :: REPR_3DINP_CLAMS2 = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_2D_CLAMSTEST = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_NTASKS = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_SEDI_PARTICLE = REPR_UNDEF

  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_WAVE  = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_LEV   = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_LEVC  = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_LAT   = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_SZA   = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_TABS  = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_TABS_DAVG = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_TWOD  = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_TWODC = REPR_UNDEF
  INTEGER, SAVE, PUBLIC :: REPR_DISSOC_RATE  = REPR_UNDEF

Initialize dimensions and representations

In bmil/messy_main_channel_clams.inc:

Dimensions

Initialize with new_dimension

ID

Name

Length

DIMID_TIME

time

1

DIMID_LON

LON

nx

DIMID_LAT

LAT

ny

DIMID_LEV

LEV

nz

DIMID_THETA

THETA

ntheta

DIMID_TRAJ

NPARTS

nparts_max

DIMID_SHUFFLE

NPARTS_SHUFFLE

nparts_max_shuffle

DIMID_MIX_GRID

NLEVS

adapt_par%nlevs

DIMID_NTASKS

ntasks

ntasks

DIMID_WAVE

WAVELENGTH

jpwave

DIMID_DISSOC_LEV

DISSOC_LEV

jpslevall

DIMID_DISSOC_LEVC

DISSOC_LEVC

jpslev

DIMID_DISSOC_LAT

DISSOC_LAT

jpschi

DIMID_DISSOC_SZA

DISSOC_SZA

jpschi

DIMID_DISSOC_NUMJ

DISSOC_NUMJ

numj_max

DIMID_SEDI_PARTICLE

SEDI_PARTICLE

nparticle_max

Dimension variables

Create with add_dimension_variable

Name of dimension variable

Name of dimension

values of dimension variable

time

time

days since starttime

YYYYMMDD

time

yyyymmdd

dt

time

delta_time

nstep

time

REAL(current_time_step, DP)

JULIAN_SECONDS

time

time in julian seconds

LON

LON

longrid

LAT

LAT

latgrid

LEV

LEV

levelgrid(1:nz)

THETA

THETA

thetagrid(1:ntheta)

Representations

Create with new_representation

  • SCALAR 

    CALL new_representation(status, SCALAR, 'SCALAR'                 &
     , rank = 0, link = '----', dctype = 0                       &
     , dimension_ids = scdim &
     , ldimlen       = scdim &
     , axis = '----'                                             &
     )
    Variables with representation SCALAR:

    • In messy_clams_si.f90

      • grid_switch

      • pre_year, pre_month, pre_day, pre_sec

      • JULTIME

  • REPR_NTASKS 

    CALL new_representation(status, REPR_NTASKS, 'REPR_NTASKS'       &
     , rank = 1, link = 'x---', dctype = DC_IX                   &
     , dimension_ids = (/ DIMID_NTASKS /) &
     , ldimlen       = (/ 1 /) &  INTEGER, SAVE, PUBLIC :: REPR_SEDI_TBOX = REPR_UNDEF
 
     , axis = 'N---'                                             &
     )

    • Variables with representation REPR_NTASKS in messy_clams_si.f90:

      • dnparts

      • dnparts_max

    • one-dimensional array of length NTASKS

    • index distributed on tasks (dctype=DC_IX), every task can use dnparts and dnparts_max as a scalar

  • REPR_3DINP_CLAMS 

    CALL new_representation(status, REPR_3DINP_CLAMS, 'REPR_3DINP_CLAMS' &
     , rank = 3, link = 'xxx-', dctype = 0                           &
     , dimension_ids = (/ DIMID_LON, DIMID_LAT, DIMID_LEV /)         &
     , ldimlen = (/ AUTO, AUTO, AUTO /)                              &
     , axis = 'XYZ-'                                                 &
     )

    • Variablen mit dieser Representation in messy_clamstraj_si.f90:

      • UDT, VDT, WDT, LEVELDT

      • UFUT, VFUT, WFUT, LEVELFUT

    • 3-dimensionale Felder
    • Dimensionierung (nlon,nlat,nlev) wie sie innerhalb des TRAJ genutzt wird

    • Die Windfelder müssen auf allen Tasks vollständig vorhanden sein => keine parallele Zerlegung der Felder (dctype=0)

  • REPR_3DINP_CLAMS2 

    CALL new_representation(status, REPR_3DINP_CLAMS2, 'REPR_3DINP_CLAMS2' &
     , rank = 3, link = 'xxx-', dctype = 0                           &
     , dimension_ids = (/ DIMID_LEV, DIMID_LAT, DIMID_LON /)         &
     , ldimlen = (/ AUTO, AUTO, AUTO /)                              &
     , axis = 'ZYX-'                                                 &
     )

    • Variablen mit dieser Representation in messy_clamstraj_si.f90:

      • PREDATA_param, FUTDATA_param (für TEMP,PRESS etc)

    • 3-dimensionale Felder
    • Dimensionierung wie auf Netcdf-Winddateien (nlev,nlat,nlon)

    • Die Felder müssen auf allen Tasks vollständig vorhanden sein => keine parallele Zerlegung der Felder (dctype=0)

  • REPR_LG_CLAMS 

    CALL new_representation(status, REPR_LG_CLAMS, 'REPR_LG_CLAMS'         &
     , rank = 1, link = 'x---', dctype = DC_IX                 &
     , dimension_ids = (/ DIMID_TRAJ /) &
     , ldimlen       = (/ dnparts_max /) &
     , axis = 'N---'                                           &
     )

    • Variablen mit dieser Representation in messy_clams_si.f90:

      • SPECARR für alle Spezies (1-nspec)

    • Variablen mit dieser Representation in messy_clamstraj_si.f90:

      • JULSEC

      • LAT, LAT_OLD, LAT_OLD_MIX

      • LON, LON_OLD, LON_OLD_MIX

      • LEV, LEV_OLD

      • param, param_OLD (PRESS,TEMP etc)

    • Variablen mit dieser Representation in messy_dissoc_si.f90:

      • jresult

    • Eindimensionale Felder der Länge nparts_max

    • Index-Zerlegung (dctype=DC_IX) => jede Task erhält ein eindimensionales Feld der Länge dnparts_max

  • REPR_3DINP_CLAMSTHETA 

    CALL new_representation(status, REPR_3DINP_CLAMSTHETA, 'REPR_3DINP_CLAMSTHETA' &
     , rank = 3, link = 'xxx-', dctype = 0                           &
     , dimension_ids = (/ DIMID_THETA, DIMID_LAT, DIMID_LON /)         &
     , ldimlen = (/ AUTO, AUTO, AUTO /)                              &
     , axis = 'ZYX-'                                                 &
     )

  • REPR_CLAMS_SHUFFLE 

    CALL new_representation(status, REPR_CLAMS_SHUFFLE, 'REPR_CLAMS_SHUFFLE'  &
     , rank = 1, link = 'x---', dctype = DC_IX                       &
     , dimension_ids = (/ DIMID_SHUFFLE /)                              &
     , ldimlen       = (/ dnparts_max_shuffle /)                              &
     , axis = 'N---'                                                 &
     )

    • Variablen mit dieser Representation in messy_clamsmix_si.f90:

      • ZETA_SHUFFLED

      • LAT_SHUFFLED

      • LON_SHUFFLED

      • LAT_OLD_MIX_SHUFFLED

      • LON_OLD_MIX_SHUFFLED

      • param_SHUFFLED (für alle Spezies)

    • Eindimensionale Felder der Länge nparts_max_shuffle

    • Index-Zerlegung (dctype=DX_IX) => jede Task erhält ein eindimensionales Feld der Länge dnparts_max_shuffle

  • REPR_MIX_GRID 

    CALL new_representation(status, REPR_MIX_GRID, 'REPR_MIX_GRID'         &
     , rank = 1, link = 'x---', dctype = 0                     &
     , dimension_ids = (/ DIMID_MIX_GRID /) &
     , ldimlen       = (/ adapt_par%nlevs /) &
     , axis = 'N---'                                           &
     )
  • REPR_DISSOC_WAVE
  • REPR_DISSOC_LEV
  • REPR_DISSOC_LEVC
  • REPR_DISSOC_SZA
  • REPR_DISSOC_LAT
  • REPR_DISSOC_TABS
  • REPR_DISSOC_TABS_DAVG
  • REPR_DISSOC_TWOD
  • REPR_DISSOC_TWODC
  • REPR_DISSOC_RATE
  • REPR_SEDI_PARTICLE

messy/ClamsRepr (last edited 2018-04-24 11:04:12 by NicoleThomas)