Running CSIM 5.0 at UD

Transcription

1 Running CSIM 5.0 at UD The following instructions are for running NCAR s Community Sea Ice Model, version 5.0, uncoupled on a SunFire 6800 running Solaris 9 with Studio 9 compiler suite. Environmental variables have been set by other scripts that are user specific ( Download the source code Very roughly, the source code is about 12 Mb while the input data and scripts require about 75 Mb once the executable has been built without atmospheric forcing data Atmospheric forcing data will require an additional 1 Gb. Daily output for a standard uncoupled, CSIM run will use about 4 Mb a day. Download the source code and datasets from the CSIM 5.0 public release homepage ( Source Code: o csim5_code.tar.gz (0.17 Mb compressed) (Sept 2004) Input Data: o csim5_inputdata_gx3v5.tar.gz (5.2 Mb compressed) (6/23/2004) o csim5_inputdata_gx1v3.tar.gz (39 Mb compressed) (6/23/2004) *Atmospheric Forcing Data: o csim5_atmforcing_gx3v5.tar.gz (253 Mb compressed) (6/23/2004) * This forcing is year 1997 of NCEP state variables, MXA precipitation and ISCCPM radiative fluxes. These datasets are meant for testing the model and are not the best observational data for research. CSIM 5.0 does not require atmospheric forcing data to setup and run (NCAR). Uncompress Files Uncompress and untar the files: gunzip filename.tar.gz tar xvf filename.tar

2 Modify run script (csim_run) 1.) Change all directory paths to your directories. CSIMDIR home directory of untarred CSIM CSIMDATA location of input data (created when the inputdata tarball is untarred) CASE name of the current experiment EXEROOT working directory and location of output 2.) Set the run configuration. CASESTR description of run OCEANMIXED_ICE whether or not to use the ocean mixed layer model (pg 21 CSIM users guide) RUNTYPE start up or continue (pg 22 CSIM users guide) NCAT # of ice categories (pg 23 CSIM users guide; pg 13 CSIM scientific description) NX & NY # of processors - set to one for a single processor run (only option at UD) - BINTYPE single for single processor run 3.) Delete NCAR specific stuff. There are a lot of references to NCAR specific computers that can be removed if it makes it easier to look through the code. This is not necessary and not recommended until you are certain you know what lines of code you need. 4.) Change run command. if ( $OS = = AIX ) csim if ( $OS = = SunOS ) csim Modify setup script (csim.setup.csh) 1.) Switch the order of the first set of if statements located six lines down from Position the input datasets. If this isn t changed it is impossible to enter the loop. Originally: If ($RESTART = =.true.) then If ($RUNTYPE = = startup) then Change to: If ($RUNTYPE = = startup) then If ($RESTART = =.true.) then

3 2.) Set the namelist variables for each simulation. Namelist variables determine time management, output frequency, filenames, model physics and options (pg 12-21, CSIM user s guide). Comment out year_int, fyear_init, ycycle, and atm_data_dir unless you are using the atmospheric forcing data. If you are not using the atmospheric forcing information you must also comment out call init_getflux and call getflux in ice.f If you are using the forcing information leave the variables active and change the directory path to the data to it s location in your space. 3.) Change the directory paths under Filepath to the appropriate path to the location of the source directory in your space. There should only be one directory in the CSIM directory tree that contains source files. Therefore there should only be one path listed under the filepath command. This command creates a Srcfile that is used by the Makefile to create the object files. This should work but if it doesn t, create a srcfile in the $CSIMDIR/source directory and add a line after the filepath command that copies the srcfile you created to the OBJDIR \cat >! Filepath << EOF $CSIMDIR/source EOF cp $CSIMDIR/source/Srcfile Scrfile 4.) Set DEBUG to true if it is a test run of false if it is a long run. 5.) Add an alias for the make command: alias gmake / path / to /make Set variables used in the makefile after the gmake command. gmake d VPFILE=Filepath MODEL=csim EXEC=$EXEROOT/csim NX=$NX NY=$NY \ f $CBLD/Makefile MACFILE=$CBLD/Macros.$OS SRCFILE=Srcfile DEBUG=$DEBUG

4 Create a macros file for the Sun operating system (Macros.SunOS) CSIM does not come with a macros file for a Sun operating system so one must be created from scratch. This may not the most efficient macros file but it works. Paths to directories are specific to the system at UD and should be changed to your directory paths. Note, the Makefile doesn t need to be changed. #********************************************************************** #Macros.SUN --> Mary D. Lemcke (5/2/2005) #********************************************************************** # Set up special characters null := space := $(null) $(null) comma := $(null),$(null) # Determine distribution root directory ifeq ($(ROOTDIR),$(null)) ROOTDIR := $(CSIMDIR) #NETCDF Definitions, $NETCDF_ROOT is /hsm/users/hanson/netcdf ifeq ($(MOD_NETCDF),$(null)) MOD_NETCDF := $(NETCDF_LIB) # Build the ESMF library (if there isn t an ESMF directory in your root # directory you can copy it over by adding a copy statement in the # setup file, which really doesn t seem to matter) ifeq ($(ESMF_ROOT),$(null)) ESMF_ROOT := ${ROOTDIR}/esmf ifeq ($(ESMF_BLD),$(null)) ESMF_BLD := $(shell (mkdir./esmf 2>/dev/null;cd./esmf;pwd)) ESMF_BOPT := O ifeq ($(DEBUG),TRUE) ESMF_BOPT := g ESMF_MOD := $(ESMF_BLD)/mod/mod$(ESMF_BOPT) ESMF_LIB := $(ESMF_BLD)/lib/lib$(ESMF_BOPT) # Check if SPMD is defined in "misc.h" # Ensure that it is defined and not just "undef SPMD" set in file #ifeq ($(SPMD),$(null)) # SPMDSET := $(shell grep SPMD misc.h) # ifneq (,$(findstring define,$(spmdset))) # SPMD := TRUE # else SPMD := FALSE # #

5 # Determine platform --> SunOS ifeq ($(UNAMES),$(null)) UNAMES := $(shell uname -s) #SunOS ifeq ($(UNAMES),SunOS) ESMF_ARCH = solaris ##NOTE: Include the paths in green as they are INCLDIR := -I. -I/usr/local/include -I/usr/include \ -I$(CSIMDIR)/source -I$(NETCDF_INC) SLIBS := -L $(NETCDF_LIB) -lnetcdf -L/usr/bin #ULIBS := -L $(ESMF_LIB)/$(ESMF_ARCH) FC := /opt/studio9/sunwspro/bin/f90 FC77 := /opt/studio9/sunwspro/bin/f77 CC := /opt/studio9/sunwspro/bin/cc CPP := /usr/ccs/lib/cpp CPPFLAGS := -P CPPDEFS := -DSUNOS -Duncoupled -DNPROC_X=$(NX) -DNPROC_Y=$(NY) CFLAGS := -c -DFORTRANUNDERSCORE FIXEDFLAGS := FREEFLAGS := # Don't use OpenMP on Solaris as it currently causes problems FFLAGS := -c -M$(CSIMDIR)/source -e -r8const - xtypemap=real:64,integer:32 MOD_SUFFIX := mod LD := $(FC) LDFLAGS := -L/opt/studio9/SUNWspro/lib -lf77compat -openmp -fast debug: $(OBJS) echo "FFLAGS: $(FFLAGS)" echo "LDFLAGS: $(LDFLAGS)" echo "OBJS: $(OBJS)" test_fc: test_fc.o $(FC) -o $@ test_fc.o $(LDFLAGS) test_nc: test_nc.o $(FC) -o $@ test_nc.o -L$(NETCDF_LIB) -lnetcdf $(LDFLAGS) ifeq ($(DEBUG),TRUE) FFLAGS += -g -dalign else # Check for override of default Fortran compiler optimizations ifeq ($(F_OPTIMIZATION_OVERRIDE),$(null)) # Inline code when not debugging -inline # Use symbol table (-g) as make execution stable and can # find out where problems are. # Turn off aggressive optimization options (fsimple=1, fns=no)

6 -O4 # so the climate will be similar to climate on the IBM FORTRAN_OPTIMIZATION := -inline=%auto -fast -g -fsimple=1 -fns=no FFLAGS += $(FORTRAN_OPTIMIZATION) ifeq ($(BINTYPE), MPI) CPPDEFS += -D_MPI else CPPDEFS += -DHIDE_MPI ifeq ($(SPMD),TRUE) FFLAGS += -I$(INC_MPI) LDFLAGS += -L$(LIB_MPI) -lmpich -lnsl -lsocket else FFLAGS += -DHIDE_MPI # Targets/rules that depend on architecture specific variables. # The ESMF library is not made to be built in parallel, so specify only # one job is to run. $(ESMF_LIB)/$(ESMF_ARCH)/libesmf.a: cd $(ESMF_ROOT); \ echo "Build the ESMF library."; \ echo "ESMF is NOT supported by the CCSM project, but by the ESMF core team in NCAR/SCD"; \ echo "See \ $(MAKE) -j 1 BOPT=$(ESMF_BOPT) ESMF_BUILD=$(ESMF_BLD) ESMF_DIR=$(ESMF_ROOT) ESMF_ARCH=$(ESMF_ARCH); time_manager.o : $(ESMF_LIB)/$(ESMF_ARCH)/libesmf.a Edit the Real Time Clock command The the subroutine timers(t1) located at the end of the file ice_timers.f uses the unix real time clock command rtc but this function isn t in the UD system. At UD there is an equivalent command called time that can be used instead of rtc. Add real (kind=dbl_kind) :: time after real (kind=dbl_kind) :: rtc Change t1=rtc( ) to t1=time( ) Mary Lemcke, Department of Geography, University of Delaware (6/17/2005)