The start.nc and startfi.nc input files - Revision history

Aurelien.falco: controle is read from starfi.nc

2024-10-18T17:02:14Z

controle is read from starfi.nc

Jleconte at 09:58, 16 October 2024

2024-10-16T09:58:41Z

Aurelien.falco: /* controle */ more on controle

2024-02-06T10:15:55Z

‎controle: more on controle

Aurelien.falco: /* start.nc file Contents */ explain some of the variables of "controle"

2024-02-06T09:59:28Z

‎start.nc file Contents: explain some of the variables of "controle"

Emillour at 08:16, 26 May 2023

2023-05-26T08:16:05Z

Martin.turbet at 16:13, 20 March 2023

2023-03-20T16:13:10Z

Martin.turbet at 10:25, 9 March 2023

2023-03-09T10:25:00Z

Emillour: Created page with "== GCM initial condition files == To run, the GCM requires initial conditions (i.e. values of state variables at a given instant). These are provided as two NetCDF files, '''..."

2022-05-04T15:01:35Z

Created page with "== GCM initial condition files == To run, the GCM requires initial conditions (i.e. values of state variables at a given instant). These are provided as two NetCDF files, '''..."

New page

== GCM initial condition files ==

To run, the GCM requires initial conditions (i.e. values of state variables at a given instant). These are provided as two NetCDF files, '''start.nc''' (which contains state variables for the dynamics) and '''startfi.nc''' (state variables for the physics). As these files contain information about fields at given grid points, they must be ''at the same resolution'' as the GCM; in other words, the GCM does not do any interpolation of initial conditions which must thus be provided on the exact same grid that the model uses.

== start.nc file Contents ==
A dump of the header of a typical '''start.nc''' (obtained via '''ncdump -h start.nc''') will show contents of the likes of:
<pre>
netcdf start {
dimensions:
index = 100 ;
rlonu = 33 ;
latitude = 33 ;
longitude = 33 ;
rlatv = 32 ;
altitude = 15 ;
interlayer = 16 ;
Time = UNLIMITED ; // (1 currently)
variables:
float controle(index) ;
controle:title = "Parametres de controle" ;
float rlonu(rlonu) ;
rlonu:title = "Longitudes des points U" ;
float rlatu(latitude) ;
rlatu:title = "Latitudes des points U" ;
float rlonv(longitude) ;
rlonv:title = "Longitudes des points V" ;
float rlatv(rlatv) ;
rlatv:title = "Latitudes des points V" ;
float ap(interlayer) ;
ap:title = "Coef A: hybrid pressure levels" ;
float bp(interlayer) ;
bp:title = "Coef B: hybrid sigma levels" ;
float aps(altitude) ;
aps:title = "Coef AS: hybrid pressure at midlayers" ;
float bps(altitude) ;
bps:title = "Coef BS: hybrid sigma at midlayers" ;
float presnivs(altitude) ;
float latitude(latitude) ;
latitude:units = "degrees_north" ;
latitude:long_name = "North latitude" ;
float longitude(longitude) ;
longitude:long_name = "East longitude" ;
longitude:units = "degrees_east" ;
float altitude(altitude) ;
altitude:long_name = "pseudo-alt" ;
altitude:units = "km" ;
altitude:positive = "up" ;
float cu(latitude, rlonu) ;
cu:title = "Coefficient de passage pour U" ;
float cv(rlatv, longitude) ;
cv:title = "Coefficient de passage pour V" ;
float aire(latitude, longitude) ;
aire:title = "Aires de chaque maille" ;
float phisinit(latitude, longitude) ;
phisinit:title = "Geopotentiel au sol" ;
float Time(Time) ;
Time:title = "Temps de simulation" ;
Time:units = "days since 1-01-01 00:00:00" ;
float ucov(Time, altitude, latitude, rlonu) ;
ucov:title = "Vitesse U" ;
float vcov(Time, altitude, rlatv, longitude) ;
vcov:title = "Vitesse V" ;
float teta(Time, altitude, latitude, longitude) ;
teta:title = "Temperature" ;
float co2_ice(Time, altitude, latitude, longitude) ;
co2_ice:title = "Traceur co2_ice" ;
float h2o_ice(Time, altitude, latitude, longitude) ;
h2o_ice:title = "Traceur h2o_ice" ;
float h2o_vap(Time, altitude, latitude, longitude) ;
h2o_vap:title = "Traceur h2o_vap" ;
float masse(Time, altitude, latitude, longitude) ;
masse:title = "C est quoi ?" ;
float ps(Time, latitude, longitude) ;
ps:title = "Pression au sol" ;

// global attributes:
:title = "Dynamic start file" ;
}
</pre>

TODO: EXPLAIN SOME MAIN VARIABLES/DIMENSIONS

== startfi.nc File Contents ==
A dump of the header of a typical '''startfi.nc''' (obtained via '''ncdump -h startfi.nc''') will show contents of the likes of:
<pre>
netcdf startfi {
dimensions:
index = 100 ;
physical_points = 994 ;
subsurface_layers = 18 ;
nlayer_plus_1 = 16 ;
number_of_advected_fields = 3 ;
nlayer = 15 ;
variables:
float controle(index) ;
controle:title = "Control parameters" ;
float soildepth(subsurface_layers) ;
soildepth:title = "Soil mid-layer depth" ;
float longitude(physical_points) ;
longitude:title = "Longitudes of physics grid" ;
float latitude(physical_points) ;
latitude:title = "Latitudes of physics grid" ;
float area(physical_points) ;
area:title = "Mesh area" ;
float phisfi(physical_points) ;
phisfi:title = "Geopotential at the surface" ;
float albedodat(physical_points) ;
albedodat:title = "Albedo of bare ground" ;
float ZMEA(physical_points) ;
ZMEA:title = "Relief: mean relief" ;
float ZSTD(physical_points) ;
ZSTD:title = "Relief: standard deviation" ;
float ZSIG(physical_points) ;
ZSIG:title = "Relief: sigma parameter" ;
float ZGAM(physical_points) ;
ZGAM:title = "Relief: gamma parameter" ;
float ZTHE(physical_points) ;
ZTHE:title = "Relief: theta parameter" ;
float inertiedat(subsurface_layers, physical_points) ;
inertiedat:title = "Soil thermal inertia" ;
float tsurf(physical_points) ;
tsurf:title = "Surface temperature" ;
float tsoil(subsurface_layers, physical_points) ;
tsoil:title = "Soil temperature" ;
float emis(physical_points) ;
emis:title = "Surface emissivity" ;
float q2(nlayer_plus_1, physical_points) ;
q2:title = "pbl wind variance" ;
float cloudfrac(nlayer, physical_points) ;
cloudfrac:title = "Cloud fraction" ;
float totcloudfrac(physical_points) ;
totcloudfrac:title = "Total fraction" ;
float hice(physical_points) ;
hice:title = "Height of oceanic ice" ;
float co2_ice(physical_points) ;
co2_ice:title = "tracer on surface" ;
float h2o_ice(physical_points) ;
h2o_ice:title = "tracer on surface" ;
float h2o_vap(physical_points) ;
h2o_vap:title = "tracer on surface" ;

// global attributes:
:title = "Physics start file" ;
}
</pre>

TODO: EXPLAIN SOME MAIN VARIABLES/DIMENSIONS/GRID STRUCTURE

[[Category:Inputs]]
[[Category:WhatIs]]

@@ Line 79: / Line 79: @@
 </pre>
-=== controle ===
+'''controle''' is mainly useless, as the model reads the '''controle''' variable from startfi.nc (see below).
-Then there is either for the generic/mars/pluto model:
+TODO: EXPLAIN SOME MAIN VARIABLES/DIMENSIONS
-Then idecal is set to 4 in the generic/mars/pluto model and 5 for earth/titan.
+* cu:
 == startfi.nc File Contents ==
@@ Line 180: / Line 153: @@
+}
 </pre>
 TODO: EXPLAIN SOME MAIN VARIABLES/DIMENSIONS/GRID STRUCTURE

@@ Line 99: / Line 99: @@
 * idecal+3   : g         (surface gravity)
 * idecal+4   : cpp       (heat capacity at constant pressure)
-* idecal+5   : kappa     (potential temperature)
+* idecal+5   : rcp       (R/cp of the dry gas, used to compute R and mu from cpp)
 * idecal+6   : daysec    (seconds in a day)
 * idecal+7   : dtvr      (seconds in a dynamic step?)
@@ Line 107: / Line 107: @@
 * idecal+11  : stot0     (entropie?)
 * idecal+12  : ang0
-* idecal+13  : pa
+* idecal+13  : pa        (transition pressure between sigma and p coordinates in hybrid coordinates. Used in computation of ap and bp)
 * idecal+14  : preff     (used in ap(l)+bp(l)*preff)
 * idecal+15  : clon

@@ Line 82: / Line 82: @@
 The variable "controle" lists a set of parameters.
 The first three parameters are always:
-* 1          : im
+* 1          : im         (physical points)
-* 2          : jm
+* 2          : jm         (nb of layers)
 * 3          : lllm
@@ Line 95: / Line 95: @@
 The list follows as such:
-* idecal+1   : rad
+* idecal+1   : rad       (planet radius)
-* idecal+2   : omeg
+* idecal+2   : omeg      (angular speed)
-* idecal+3   : g
+* idecal+3   : g         (surface gravity)
-* idecal+4   : cpp
+* idecal+4   : cpp       (heat capacity at constant pressure)
-* idecal+5   : kappa
+* idecal+5   : kappa     (potential temperature)
 * idecal+6   : daysec    (seconds in a day)
 * idecal+7   : dtvr      (seconds in a dynamic step?)
-* idecal+8   : etot0
+* idecal+8   : etot0     (energie)
 * idecal+9   : ptot0
-* idecal+10  : ztot0     (periastron)
+* idecal+10  : ztot0     (enstrophie?)
-* idecal+11  : stot0     (apoastron)
+* idecal+11  : stot0     (entropie?)
 * idecal+12  : ang0
 * idecal+13  : pa
@@ Line 112: / Line 112: @@
 * idecal+16  : clat
 * idecal+17  : grossismx
 * idecal+18  : grossismy
 === Other parameters ===

← Older revision		Revision as of 09:59, 6 February 2024
Line 79:		Line 79:
	</pre>		</pre>

		+	=== controle ===
		+	The variable "controle" lists a set of parameters.
		+	The first three parameters are always:
		+	* 1 : im
		+	* 2 : jm
		+	* 3 : lllm
		+
		+	Then there is either for the generic/mars/pluto model:
		+	* 4 : day_ref
		+	or for earth/titan:
		+	* 4 : day_ini
		+	* 5 : annee_ref
		+
		+	Then idecal is set to 4 in the generic/mars/pluto model and 5 for earth/titan.
		+
		+	The list follows as such:
		+	* idecal+1 : rad
		+	* idecal+2 : omeg
		+	* idecal+3 : g
		+	* idecal+4 : cpp
		+	* idecal+5 : kappa
		+	* idecal+6 : daysec (seconds in a day)
		+	* idecal+7 : dtvr (seconds in a dynamic step?)
		+	* idecal+8 : etot0
		+	* idecal+9 : ptot0
		+	* idecal+10 : ztot0 (periastron)
		+	* idecal+11 : stot0 (apoastron)
		+	* idecal+12 : ang0
		+	* idecal+13 : pa
		+	* idecal+14 : preff (used in ap(l)+bp(l)*preff)
		+	* idecal+15 : clon
		+	* idecal+16 : clat
		+	* idecal+17 : grossismx
		+	* idecal+18 : grossismy
		+
		+	=== Other parameters ===
	TODO: EXPLAIN SOME MAIN VARIABLES/DIMENSIONS		TODO: EXPLAIN SOME MAIN VARIABLES/DIMENSIONS

@@ Line 149: / Line 149: @@
 [[Category:Inputs]]
 [[Category:WhatIs]]
-[[Category:Generic-LMDz]]
+[[Category:Generic-LMDZ]]

← Older revision		Revision as of 10:25, 9 March 2023
Line 149:		Line 149:
	[[Category:Inputs]]		[[Category:Inputs]]
	[[Category:WhatIs]]		[[Category:WhatIs]]
		+	[[Category:Generic-Model]]