Changes between Version 5 and Version 6 of RunningDynamico

Timestamp:

09/08/14 09:19:59 (10 years ago)

Author:

millour

Comment:

--

RunningDynamico

@@ -17 +17 @@
 By default DYNAMICO solves the traditional, shallow-atmosphere, hydrostatic equations. An example is provided in [source:codes/icosagcm/trunk/param_sets/dcmip2012/ dcmip2012/test4/test4.1/test4.1-0  ]. This example runs a dry baroclinic instability test case (Jablonowski & Williamson, 2006). Sample configuration files for the climate-like Held and Suarez (1994) benchmark can be found in [source:codes/icosagcm/trunk/param_sets/climate/Held_Suarez climate/Held_Suarez ]
 
+== Running DYNAMICO with Saturn (LMDZ-GENERIC) physics ==
+This is for those playing with the SATURN_DYNAMICO branch.
+A quick tutorial on setting up the Saturn test case on Ada:
+* Download the SATURN_DYNAMICO branch:
+{{{
+svn co svn+ssh://yourlogin@forge.ipsl.jussieu.fr/ipsl/forge/projets/dynamico/svn/codes/icosagcm/branches/SATURN_DYNAMICO
+}}}
+* compile XIOS using script in directory XIOS:
+{{{
+./compile_ada
+}}}
+* Compile the model using script (current settings are for MPI compilation) in directory ICOSAGCM:
+{{{
+./compile_dynlmdz_ada
+}}}
+* Set up the test case in some directory:
+ 1. Copy over icosa_gcm.exe from ICOSAGCM/bin
+ 2. Copy over apbp.txt (vertical coordinates specifications) and temp_profile.txt (input initial temperature profile) from TEST directory
+ 3. Copy over all the *.def files from TEST directory
+ 4. Copy over all the *.xml files from TEST directory (these control the XIOS ouputs)
+ 5. Modify the path "datadir" in callphys.def:
+{{{
+datadir = /ccc/scratch/cont003/dsm/p86yann/SATURNE_128x96x64/DATAGENERIC/
+}}}
+to point to the TEST/DATAGENERIC directory
+ 6. Set run parameters in run_icosa.def (e.g. nqtot=0 since there are no tracers, run_length=..., etc.)
+ 7. Run the model using a job (see sample script launch.ada in TEST
+ 8. With XIOS outputs, the output file xios_diagfi.nc is on native icosahedral grid, so it is usually better to reinterpolate onto a lon-lat grid, e.g.:
+{{{
+% cat mygrid
+gridtype = lonlat
+xsize    = 90
+ysize    = 45
+xfirst   = -180
+xinc     = 4
+yfirst   = -90
+yinc     = 4
+% cdo remapdis,mygrid xios_diagfi.nc xios_latlon.nc
+}}} 
+
 == Horizontal resolution ==
 
-Horinzontal resolution is controlled by the parameter nbp defined in run.def. The total number of hexagonal cells is about 10 x nbp x nbp, corresponding to subdividing each main triangle of the icosahedron in nbp x nbp sub-triangles (there are about twice as many triangles as there are hexagons).
+Horizontal resolution is controlled by the parameter nbp defined in run.def. The total number of hexagonal cells is about 10 x nbp x nbp, corresponding to subdividing each main triangle of the icosahedron in nbp x nbp sub-triangles (there are about twice as many triangles as there are hexagons).
 Notice that, everything else being equal, the time step (dt) should be inversely proportional to nbp for numerical stability.