Changes between Version 4 and Version 5 of PMIP4Page

Timestamp:

04/30/19 17:51:22 (5 years ago)

Author:

flhardy

Comment:

--

PMIP4Page

@@ -6 +6 @@
 The high resolution topography etopo1 is used to generate a new bathymetry and land-sea mask for CLIO in PI conditions. [[BR]]
 
-__Pre-processing :__ [[BR]]
-1. etopo1 --> ''connectivity program '' --> clio_rotated.txt and clio_regular.txt --> ''land-sea mask program '' --> bath.txt [[BR]]
+__Semi-automated pre-processing steps :__ [[BR]]
 
-The connectivity program writes the chosen variables (mean bathymetry, hypsometry...) into a text file, either on the rotated or regular clio grid. [[BR]] The two grids are put together in the land-sea mask program, which writes a text file containing the bathymetry with the land-sea mask (0 in every land grid cells). [[BR]] The mask is generated using the hypsometry, a chosen sea-level (-0.5 m for the PI), and a chosen threshold (% of surface of a grid cell above which the cell is defined as ocean - here 40 %). The mask of a few ocean grid cells is manually forced so that all the critical traits stay open. Small seas are closed. [[BR]] 
-The bathymetry is converted into vertical levels using the bath.om file, which contains the old depth as a vertical level (from 1 to 20, 20 being the deepest cell) and the definition of these irregular vertical levels (center and height). The old vertical level is forced for a few problematic grid cells in order to get realistic salinity values in the Mediterranean Sea and Hudson Bay. The vertical level 1 is avoided (either forced to 0 or 2), because the model cannot deal with the very shallow grid cells. [[BR]] 
+[[Image(preprocessing_bathychange.png)]]
 
-2. bath.txt --> ''pre-poids '' (modified defgrid.f) --> clio.msk and clio.srf --> ''poids '' --> fractoc.dat and mozaic.w [[BR]]
+1. The connectivity program writes the chosen variables (mean bathymetry, hypsometry...) into a text file, either on the rotated or regular CLIO grid. It can also produce the connections between ocean basins thanks to the computation of subgrid sills. [[BR]] 
+2. The two grids are put together in the bathymetry program, which writes a text file containing the bathymetry with the land-sea mask (0 in every land grid cells). [[BR]] The mask is generated using the hypsometry, a chosen sea-level (-0.5 m for the PI), and a chosen threshold (% of surface of a grid cell above which the cell is defined as ocean - here 40 %). The mask of a few ocean grid cells is manually forced so that all the critical traits stay open. Small isolated seas are closed. [[BR]] The bathymetry is converted into vertical levels using the bath.om file, which contains the old depth as a vertical level (from 1 to 20, 20 being the deepest cell) and the definition of these irregular vertical levels (center and height). The old vertical level is forced for a few problematic grid cells in order to get realistic salinity values in the Mediterranean Sea and Hudson Bay. The vertical level 1 is avoided (either forced to 0 or 2), because the model cannot deal with these very shallow grid cells. As the model also cannot deal with 'isolated' oceanic grid cells for which the last vertical level is isolated (e.g. deep grid cells with shallower neighbours), a procedure similar to a smoothing filter is applied. [[BR]] 
+3. The pre-poids program loads the new bathymetry (bath.txt file) and generate two files necessary to the next step. It contains the defgrid.f routine, slightly modified to deal the newly generated 'isolated' cells at the edges (longitudes 2 and 121). [[BR]]
+4. The poids program uses the clio.msk and clio.srf files to compute the interpolation points between the CLIO and the T21 grid (mozaic.w), as well as the fraction of ocean seen by each T21 grid cells (fractoc.dat, aka fort.8). When running this program, kamax = 14 and komax = 17 should be printed. [[BR]]
 
-''Note '' : kamax = 14 and komax = 17 when running the poids program. [[BR]]
+__How to run a PI simulation with the new bathymetry and land-sea mask :__ [[BR]]
 
-__How to run :__ [[BR]]
+Run with -F -I (-P), using a scenario (pi-etopo1 - commited but outdated) containing :
+- inputdata : fractoc.dat, clio/bath.txt, clio/mozaic.w, clio/bath.om (definition of vertical levels still needed)
+- parameters : run.param, with a BERING_VALUE tuned to get a realistic Bering transport (around 0.8 Sv)
+- sources : defgrid.f (with the right bath.txt name and modification for isolated cells at the edges)
 
-To run a PI simulation with the new bathymetry and land-sea mask, a directory scenario/ has to be created. The fractoc.dat and bath.txt (but also the old bath.om) have to be in inputdata/, whereas mozaic.w has to be in inputdata/clio. Don't forget the modified defgrid.f in sources/. Run with -I -F. [[BR]]
+__Potential sources of error when running :__ [[BR]]
+- "no ocean basin defined in x,y" : change the ocbas.dat file in ecbilt_clio/parameters or include a modified one in scenario/inputdata.
 
-__3 potential sources of error when running :__ [[BR]]
-- "no ocean basin defined in x,y" : change the ocbas.dat file in ecbilt_clio/parameters/ [[BR]]
+- "error in runoff budget" : set bricole_Caspienne=.false. in ludus-code/lbm/sources/routageEAU.f90.
 
-- "error in runoff budget" : set bricole_Caspienne=.false. in ludus-code/lbm/sources/routageEAU.f90 [[BR]]
+- "Albedo out of range" : conflict with alb.dat, comment the LGMSWITCH.EQ.1 condition in landmodel0.f for now.
 
-- "rain larger than maximum set by rainmax" and "Albedo of land out of range" : variables are diverging due to a non-zero sea ice fraction in a few tropical grid cells (which were formerly land but are now ocean). The initialization of all the restart variables in new ocean grid cells in redrunb.f has been modified.
+- "rain larger than maximum set by rainmax" : you're in trouble, variables are diverging (e.g. due to a non-zero sea ice fraction in a few tropical grid cells - which were formerly land but are now ocean ?). Check if the number of 'isolated' grid cells has fallen to 0 in step 2.
+
+- NaN : whatever, try a restart.
 
 == '''2. LGM simulations''' ==
@@ -33 +39 @@
 __Pre-processing :__ [[BR]]
 
-The same pre-processing is applied with the following LGM topography, computed with a regridding : [[BR]]
+The same pre-processing is applied with the following LGM topography, computed with a regridding on the etopo1 grid : [[BR]]
 {{{
 LGM topography = PI (etopo1) + LGM Peltier (ICE-6G-C 21kyr) - PI Peltier (ICE-6G-C 0kyr)
@@ -67 +73 @@
 
 == '''3. Transient simulations''' ==
+
+To accomodate restarts with a modified bathymetry (and in the long run an interactive bathymetry for transient simulations), the model has to generate realistic values of the tracers content of new oceanic grid cells. A new module has been commited (clio/sources/update_clio_bathy_tools.f90). The initialization of all the restart variables in new ocean grid cells (in redrunb.f) is now done by either setting 0 or averaging the values in neighbouring oceanic grid cells.