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	Name of dimension variable	values of dimension variable
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)

Representations

Create with new_representation

SCALAR

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

Variablen mit dieser Representation:

In messy_clams_si.f90
- grid_switch
- year_uvfirst, month_uvfirst, day_uvfirst, hour_uvfirst
In messy_clamstraj_si.f90
- 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---'                                             &
     )

Variable mit dieser Representation in messy_clams_si.f90:
- dnparts
Eindimensionales Feld der Länge NTASKS
Index-verteilt auf die Tasks (dctype=DC_IX), so dass jede Task dnparts (Anz. der von der Task zu bearbeitenden Trajektorien) als Skalar nutzen kann

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_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_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