-
Notifications
You must be signed in to change notification settings - Fork 0
Running Coupled Ocean
Steps to run AOGCM
As of this writing the stable tag is 'yuri-coupled-Fortuna-2_5_obio' and module is 'Fortuna'. Check with Yury Vikhliaev for the latest stable tag and module.
$ cvs co -r yuri-coupled-Fortuna-2_5_obio Fortuna
$ cd GEOSagcm/src
$ cvs up -r Ganymed-2_0_UNSTABLE g5_modules
$ source g5_modules
Compile time specification of resolution/layout
The seaice component has the resolution and layout specified in the GNUmakefile in:
$ESMADIR/src/GEOSgcs_GridComp/GEOSgcm_GridComp/GEOSagcm_GridComp/
GEOSphysics_GridComp/GEOSsurface_GridComp/GEOSCICEThermo_GridComp/cice/GNUmakefile
Inside this are lines:
NXGLOB = 360
NYGLOB = 200
BLCKX = 45
BLCKY = 20
MXBLCKS = 1
where (NXGLOB, NYGLOB) define the ocean resolution. It is the user's responsibility to ensure that BLCKX (BLCKY) is a divisor of NXGLOB (NYGLOB) and the quotient defines the number of processors OGCM_NX (OGCM_NY) in the X (Y) direction, i.e.
OGCM_NX = NXGLOB/BLCKX (= 8 in our example)
OGCM_NY = NYGLOB/BLCKY (= 10 in our example)
[I had tried with NXGLOB = 360, BLCKX = 40 => OGCM_NX = 9. It didn't work. Apparently, odd numbered processes don't. Update: I could get it to work on 1x1 layout.]
Since the resolution/layout are specified at compile time (really!), the model needs to be re-compiled for a new resolution/layout combination.
Finally, compile:
$ make install
$ cd Applications/GEOSgcm_App
Set correct values of OGCM_NX and OGCM_NY in gcm_setup (if they are not already set):
set OGCM_NX = 8
set OGCM_NY = 10
Then run gcm_setup:
$ ./gcm_setup
Experiment ID: Choose
Experiement Description: Choose
Atmospheric Horizontal Resolution: 144 91
Atmospheric Vertical Resolution: Default
COUPLED Ocean/Sea-Ice Model: YES
Ocean Lat/Lon Horizontal Resolution: Default
Ocean Model Vertical Resolution: Default
GOCART: Default
HISTORY template: HISTORY_COUPLED.rc.tmpl
HOME directory: Choose
EXPERIMENT directory: Same as HOME
BUILD directory: Default
GROUP ID: Default
$ cd /path/to/exp/dir
To run a 1 day experiment etc. Also timers and memutils can be enabled here.
Need to ensure:
OGCM_NX = NX = 8
OGCM_NY = NY = 10
CICE_NPROCS = 80
layout = 8,10,
Comment out the line
'RUNOFF' , 'SURFACE' ,
in geosgcm_surf.fields (model crashes with this option).
#PBS -l select=7:ncpus=12:mpiprocs=12
(need 80 cpus, hence 7 nodes) Also commenting out the 2 qsub lines prevents gcm_run.j from running post processing and re-submitting jobs.
Example restarts can be found in ~yvikhlia/geos5/test180/expdir (for atmosphere resolution of 144x91 and ocean resolution of 360x200). Several new restarts are needed but these can be bootstrapped. The new ones are orad_import, seaice_import, and seaice_internal. [The required restarts are fvcore_internal_rst, lake_internal_rst, landice_internal_rst, catch_internal_rst, moist_internal_rst, saltwater_internal_rst, seaice_internal_rst, seaice_import_rst, orad_import_rst. If using GOCART, then gocart_internal_rst, carma_internal_rst, stratchem_internal_rst, gmichem_internal_rst as well.]
$ cp ~yvikhlia/geos5/test180/expdir/*_rst .
$ cp ~yvikhlia/geos5/test180/expdir/cap_restart .
cap_restart may need to edited as well to specify a particular year (no changes in month, day or time).
When the experiment directory is created (by gcm_setup) a new directory is now there called RESTART. Inside are restart files [only ocean_temp_salt.res.nc is required, rest can be bootstrapped].
$ cp ~yvikhlia/geos5/test180/expdir/RESTART/* RESTART/
We need ocean grid_spec, ice grid, KPAR data and river routing file during a run. They are all in directory pointed to by $GRIDDIR in the runscript (gcm_run.j).
$ qsub gcm_run.j
Category:SI Team Category:Running the Model