Index: L/script_install_new =================================================================== --- /BOL/script_install_new (revision 6163) +++ (revision ) @@ -1,1539 +1,0 @@ -#!/bin/bash -set -eu # error on command failure, and on unset variables -export LC_ALL=C # standardize awk format - -########################################################################### -# Author : Laurent Fairhead et Frédéric Hourdin -# Usage : install_lmdz.sh -help -# -# bash installation script of the LMDZ model on different computer types : -# Linux PC, "mesocentre" (IPSL-UPMC, IPSL-X), super-computer (IDRIS) -# -# The model is downloaded in the following directory tree -# $MODEL/modipsl/modeles/... -# using the "modipsl" infrastructure created by the "IPSL" -# for coupled (atmosphere/ocean/vegetation/chemistry) climate modeling -# activities. -# Here we only download atmospheric (LMDZ) and vegetation (ORCHIDEE) -# components. -# -# The sources of the models can be found in the "modeles" directory. -# In the present case, LMDZ, ORCHIDEE, and IOIPSL or XIOS (handling of -# input-outputs using the NetCDF library). -# -# The script downloads various source files (including a version of NetCDF) -# and utilities, compiles the model, and runs a test simulation in a -# minimal configuration. -# -# Prerequisites : pgf90/gfortran, bash or ksh, wget , gunzip, tar, ... -# -# Modif 18/11/2011 -# changes for option real 8. -# We compile with -r8 (or equivalent) and -DNC_DOUBLE for the GCM -# but with -r4 for netcdf. Variable real must be set to -# r4 or r8 at the beginning of the script below. -# -########################################################################### - -echo install_lmdz.sh DEBUT "$(date)" - -################################################################ -# Choice of installation options -################################################################ - -################################################################ -# A function to fetch files either locally or on the internet -################################################################ -function myget { #1st and only argument should be file name - # Path on local computer where to look for the datafile - if [[ -f /u/lmdz/WWW/LMDZ/pub/$1 ]]; then - \cp -f -p "/u/lmdz/WWW/LMDZ/pub/$1" . - elif [[ -f ~/LMDZ/pub/$1 ]]; then - \cp -f -p ~"/LMDZ/pub/$1" . - else - wget --no-check-certificate -nv "http://lmdz.lmd.jussieu.fr/pub/$1" - save_pub_locally=0 - if [[ $save_pub_locally = 1 ]]; then # saving wget files on ~/LMDZ/pub - dir=~"/LMDZ/pub/$(dirname "$1")"; mkdir -p "$dir" - cp -r "$(basename "$1")" "$dir" - fi - fi -} - -# 04_2021 : tester si r4 marche encore ! -real="r8" # ou "r4" - - -######################################################################### -# Valeur par défaut des parametres -######################################################################### -svn="" -version="20231022.trunk" - -getlmdzor=1 -netcdf=1 # 1: for automatic installation; -# or 0: do not install NetCDF and look for it in standard locations; -# or absolute path: look for NetCDF there -check_linux=1 -ioipsl=1 -bench=1 -pclinux=1 -pcmac=0 # default: not on a Mac -compiler="gfortran" -if [[ $(gfortran -dumpversion | cut -d. -f1) -ge 10 ]]; then allow_arg_mismatch="-fallow-argument-mismatch"; fi -SCM=0 -# surface/vegetation scheme treatment -# controlled by the single variable veget which can have the following values -# - NONE: bucket scheme (default) -# - CMIP6: orchidee version used in CMIP exercise, rev 5661 -# - number: orchidee version number -veget="NONE" -# choose the resolution for the bench runs -# grid_resolution= 32x24x11 or 48x36x19 for tests (test without ORCHIDEE) -# 96x71x19 standard configuration -grid_resolution="144x142x79" -grid_resolution="96x95x39" -grid_resolution="48x36x19" -grid_resolution="32x32x39" -# choose the physiq version you want to test -#physiq=NPv6.0.14splith -physiq="" - -## parallel can take the values none/mpi/omp/mpi_omp -parallel="mpi_omp" -parallel="none" -idris_acct="lmd" -trusting="testing" -OPT_GPROF="" -OPT_MAKELMDZ="" -MODEL="" - -## also compile XIOS? (and more recent NetCDF/HDF5 libraries) Default=no -with_xios="n" -opt_makelmdz_xios="" - -## compile with oldrad/rrtm/ecrad radiatif code (Default=rrtm) -rad="rrtm" - -## compile_with_fcm=1 : use makelmdz_fcm (1) or makelmdz (0) -compile_with_fcm=1 - -#Compilation with Cosp (cosp=NONE/v1/v2; default=NONE) -cosp="NONE" -opt_cosp="" - -# Check if on a Mac -if [[ $(uname) = "Darwin" ]]; then - pcmac=1 - export MAKE=make -fi - -env_file="" -optim="" - -######################################################################### -# Options interactives -######################################################################### -while (($# > 0)) -do - case $1 in - "-h") cat <<........fin - $0 [ -v version ] [ -r svn_release ] - [ -parallel PARA ] [ -d GRID_RESOLUTION ] [ -bench 0/1 ] - [-name LOCAL_MODEL_NAME] [-gprof] [-opt_makelmdz] [-rad RADIATIF] - - -v "version" like 20150828.trunk - see http://www.lmd.jussieu.fr/~lmdz/Distrib/LISMOI.trunk - - -r "svn_release" : either the svn release number or "last" - - -compiler gfortran|ifort|pgf90 (default: gfortran) - - -parallel PARA : can be mpi_omp (mpi with openMP) or none (for sequential) - - -d GRID_RESOLUTION should be among the available benchs if -bench 1 - among which : 48x36x19, 48x36x39 - if wanting to run a bench simulation in addition to compilation - default : 48x36x19 - - -bench activating the bench or not (0/1). Default 1 - - -testing/unstable - - -name LOCAL_MODEL_NAME : default = LMDZversion.release - - -netcdf 0, 1 or PATH - 0: do not download NetCDF, look for it in standard locations - 1: download and compile NetCDF - PATH: full path to an existing installed NetCDF library - - -xios also download and compile the XIOS library - (requires the NetCDF4-HDF5 library, also installed by default) - (requires to also have -parallel mpi_omp) - - -gprof to compile with -pg to enable profiling with gprof - - -cosp to run without our with cospv1 or cospv2 [none/v1/v2] - - -rad RADIATIF can be oldrad, rrtm or ecrad radiatif code - - -nofcm to compile without fcm - - -SCM install 1D version automatically - - -debug compile everything in debug mode - - -opt_makelmdz to call makelmdz or makelmdz_fcm with additional options - - -physiq to choose which physics package to use - - -env_file specify an arch.env file to overwrite the existing one - - -veget surface model to run [NONE/CMIP6/xxxx] - -........fin - exit 1;; - "-v") version=$2; shift; shift;; - "-r") svn=$2; shift; shift;; - "-compiler") compiler=$2 - case $compiler in - "gfortran"|"ifort"|"pgf90") compiler=$2; shift - shift;; - *) echo "Only gfortran , ifort or pgf90 for the " \ - "compiler option" - exit 1 - esac;; - "-d") grid_resolution=$2; shift; shift;; - "-gprof") OPT_GPROF="-pg"; shift;; - "-unstable"|"-testing") trusting=$(echo "$1" | cut -c2-) ; shift;; - "-cosp") cosp=$2 - case $cosp in - "none"|"v1"|"v2") cosp=$2; shift; shift;; - *) echo Only none v1 v2 for cosp option; exit 1 - esac;; - "-nofcm") compile_with_fcm=0; echo This option will be reactivated soon '(promesse du 8dec2022)'; exit 1; shift;; - "-SCM") SCM=1; shift;; - "-opt_makelmdz") OPT_MAKELMDZ="$2"; shift; shift;; - "-rrtm") if [[ $2 = "false" ]]; then - rad="oldrad" - else - rad="rrtm" - fi - shift; shift;; - "-rad") rad=$2 - case $rad in - "oldrad"|"rrtm"|"ecrad") rad=$2; shift; shift;; - *) echo "Only oldrad rrtm ecrad for rad option"; exit 1 - esac;; - "-parallel") parallel=$2 - case $parallel in - "none"|"mpi"|"omp"|"mpi_omp") parallel=$2; shift - shift;; - *) echo "Only none mpi omp mpi_omp for the parallel option" - exit 1 - esac;; - "-bench") bench=$2; shift; shift;; - "-debug") optim=-debug; shift;; - "-name") MODEL=$2; shift; shift;; - "-netcdf") netcdf=$2; shift; shift;; - "-physiq") physiq=$2; shift; shift;; - "-xios") with_xios="y"; shift;; - "-env_file") env_file=$2; shift; shift;; - "-veget") veget=$2; shift; shift;; - *) bash install_lmdz.sh -h; exit 1 - esac -done - -# Option de compilation du rayonnement : depend de $mysvn ><= r4185, -# sera donc definie plus bas - -#opt_rad="" -#case $rad in -# rrtm) opt_rad="-rad rrtm";; -# ecrad) opt_rad="-rad ecrad";; -#esac - - -# Option de compilation pour Cosp -opt_cosp="" -case $cosp in - v1) opt_cosp="-cosp true";; - v2) opt_cosp="-cospv2 true";; -esac - -# Check on veget version -if [[ $veget != 'NONE' && $veget != "CMIP6" ]]; then - re='^[0-9]+$' - if ! [[ $veget =~ $re ]]; then - echo 'Valeur de l option veget non valable' - exit 1 - fi -fi - -#Define veget-related suffix for gcm name -if [[ $veget = 'NONE' ]]; then - suff_orc='' - #For use with tutorial, orchidee_rev is also defined (will be - #written in surface_env at the end of the script) - orchidee_rev='' -else - suff_orc='_orch' -fi - - -if [[ $parallel = none ]]; then - sequential=1; suff_exe='_seq' -else - sequential=0; suff_exe='_para_mem' -fi - -#Chemin pour placer le modele -if [[ $MODEL = "" ]]; then MODEL=./LMDZ$version$svn$optim; fi - - -arch=local - - -if [[ $compiler = g95 ]]; then echo g95 is not supported anymore; exit 1; fi - -################################################################ -# Specificite des machines -################################################################ - -hostname=$(hostname) -if [[ $pclinux = 1 ]]; then o_ins_make="-t g95"; else o_ins_make=""; fi - -case ${hostname:0:5} in - - jean-) compiler="mpiifort"; - par_comp="mpiifort"; - o_ins_make="-t jeanzay"; - make=gmake; - module purge - module load intel-compilers/19.0.4; - #module load intel-mpi/19.0.4; - #module load intel-mkl/19.0.4; - module load hdf5/1.10.5-mpi; - module load netcdf/4.7.2-mpi; - module load netcdf-fortran/4.5.2-mpi; - module load subversion/1.9.7; - module load cmake - export NETCDF_LIBDIR=./ - export NETCDFFORTRAN_INCDIR=./ - export NETCDFFORTRAN_LIBDIR=./ - arch=X64_JEANZAY;; - - cicla|camel) compiler="gfortran"; - module purge - module load gnu/10.2.0 - module load openmpi/4.0.5 - module load hdf5/1.10.7-mpi - module load netcdf-c/4.7.4-mpi - module load netcdf-fortran/4.5.3-mpi - netcdf=/net/nfs/tools/PrgEnv/linux-scientific6-x86_64/gcc-10.2.0/netcdf-fortran-4.5.3-k3drgfqok3lip62hnm3tsyof4cjen5sk - module load svn/1.14.0 - - if [[ $parallel != none ]]; then - root_mpi=/net/nfs/tools/meso-sl6/openmpi/4.0.5-gcc-10.2.0 - path_mpi=$root_mpi/bin; - par_comp=${path_mpi}/mpif90; - mpirun=${path_mpi}/mpirun; - fi; - arch="local" ; - make="make"; - o_ins_make="-t g95";; - - *) if [[ $parallel = none || -f /usr/bin/mpif90 ]]; then - path_mpi=$(which mpif90 | sed -e s:/mpif90::); - if [[ -d /usr/lib64/openmpi ]]; then - root_mpi="/usr/lib64/openmpi" - else - root_mpi="/usr" - fi - else - echo "Cannot find mpif90"; - if [[ $parallel = none ]]; then exit 1; fi; - fi; - if [[ $parallel != none ]]; then - root_mpi=$(which mpif90 | sed -e s:/bin/mpif90::) - path_mpi=$(which mpif90 | sed -e s:/mpif90::) - export LD_LIBRARY_PATH=${root_mpi}/lib:$LD_LIBRARY_PATH - fi - par_comp=${path_mpi}/mpif90; - mpirun=${path_mpi}/mpirun; - arch="local" ; - make="make"; - o_ins_make="-t g95" -esac - -# Flags for parallelism: -if [[ $parallel != none ]]; then - # MPI_LD are the flags needed for linking with MPI - MPI_LD="-L${root_mpi}/lib -lmpi" - if [[ "$compiler" = "gfortran" ]]; then - # MPI_FLAGS are the flags needed for compilation with MPI - MPI_FLAGS="-fcray-pointer" - # OMP_FLAGS are the flags needed for compilation with OpenMP - OMP_FLAGS="-fopenmp -fcray-pointer" - # OMP_LD are the flags needed for linking with OpenMP - OMP_LD="-fopenmp" - elif [[ "$compiler" = "ifort" ]]; then - MPI_FLAGS="" - OMP_FLAGS="-openmp" - OMP_LD="-openmp" - else # pgf90 - MPI_FLAGS="" - OMP_FLAGS="-mp" - OMP_LD="-mp" - fi -fi - -##################################################################### -# Test for old gfortran compilers -# If the compiler is too old (older than 4.3.x) we test if the -# temporary gfortran44 patch is available on the computer in which -# case the compiler is changed from gfortran to gfortran44 -# Must be aware than parallelism can not be activated in this case -##################################################################### - -if [[ "$compiler" = "gfortran" ]]; then - gfortran=gfortran - gfortranv=$(gfortran --version | \ - head -1 | awk ' { print $NF } ' | awk -F. ' { print $1 * 10 + $2 } ') - if [[ $gfortranv -le 43 ]]; then - \cat < tt.f90 -print*,'coucou' -end -eod - $compiler tt.f90 || a.out - ./a.out >| tt - if [[ $(< tt sed -e 's/ //g' ) != "coucou" ]]; then - echo "problem installing with compiler $compiler"; exit 1; fi - \rm tt a.out tt.f90 - fi -fi - -########################################################################### -if [[ $getlmdzor = 1 && ! -d $MODEL/modipsl ]]; then -########################################################################### - echo "##########################################################" - echo "Download a slightly modified version of LMDZ" - echo "##########################################################" - cd "$MODEL" - getlog=$(pwd)/get.log - echo logfile : "$getlog" - myget "src_archives/$trusting/modipsl.$version.tar.gz" >> get.log 2>&1 - echo install_lmdz.sh wget_OK "$(date)" - gunzip "modipsl.$version.tar.gz" >> get.log 2>&1 - tar xf "modipsl.$version.tar" >> get.log 2>&1 - \rm "modipsl.$version.tar" -fi - -########################################################################### -echo Installing Netcdf -########################################################################### - -if [[ $netcdf = 0 ]]; then - ncdfdir=/usr - -else - cd "$MODEL" - - case $compiler in - gfortran) opt1="-compiler gnu"; opt2="-CC gcc -FC gfortran -CXX g++";; - ifort) opt1="-compiler intel"; opt2="-CC icc -FC ifort -CXX icpc";; - pgf90) opt1="-compiler pgf90"; opt2="-CC pgcc -FC pgf90 -CXX pgCC";; - *) echo "unexpected compiler $compiler" for netcdf; exit 1 - esac - - case $with_xios in - n) script_install_netcdf=install_netcdf4_hdf5_seq.bash - ncdfdir=netcdf4_hdf5_seq - opt_=$opt1;; - y) script_install_netcdf=install_netcdf4_hdf5.bash - ncdfdir=netcdf4_hdf5 - opt_="$opt2 -MPI $root_mpi";; - *) echo with_xios=$with_xios, should be n or y; exit 1 - esac - if [[ $netcdf = 1 ]]; then - ncdfdir=$MODEL/$ncdfdir - else - mkdir -p $netcdf; ncdfdir=$netcdf/$ncdfdir - fi - - echo "Repertoire netcdf $ncdfdir" - if [[ ! -d $ncdfdir ]]; then - netcdflog=$(pwd)/netcdf.log - echo "----------------------------------------------------------" - echo "Compiling the Netcdf library" - echo "----------------------------------------------------------" - echo "log file : $netcdflog" - myget script_install/$script_install_netcdf >> "$netcdflog" 2>&1 - chmod u=rwx $script_install_netcdf - # shellcheck disable=SC2086 - ./$script_install_netcdf -prefix "$ncdfdir" $opt_ >> "$netcdflog" 2>&1 - fi - - #---------------------------------------------------------------------------- - # LF rajout d'une verrue, pour une raison non encore expliquee, - # la librairie est parfois rang�e dans lib64 et non dans lib - # par certains compilateurs - if [[ ! -e lib && -d lib64 ]]; then ln -s lib64 lib; fi - #---------------------------------------------------------------------------- - - echo install_lmdz.sh netcdf_OK "$(date)" - -fi - -cat >test_netcdf90.f90 < tmp - sed -e "s:F_L = g95:F_L = $compiler:" \ - -e "s:F_C = g95 -c -cpp:F_C = $F_C": \ - -e 's/g95.*.w_w.*.(F_D)/g95 w_w = '"$OPTIMGCM"'/' \ - -e 's:g95.*.NCDF_INC.*.$:g95 NCDF_INC= '"$ncdfdir"'/include:' \ - -e 's:g95.*.NCDF_LIB.*.$:g95 NCDF_LIB= -L'"$ncdfdir"'/lib -lnetcdff -lnetcdf:' \ - -e 's:g95 L_O =:g95 L_O = -Wl,-rpath='"$ncdfdir"'/lib:' \ - -e "s:-fmod=:-$fmod:" -e 's/-fno-second-underscore//' \ - -e 's:#-Q- g95 M_K = gmake:#-Q- g95 M_K = make:' \ - tmp >| AA_make.gdef - - if [[ $pcmac == 1 ]] - then - cp AA_make.gdef tmp - sed -e 's/rpath=/rpath,/g' tmp > AA_make.gdef - fi - - - # We use lines for g95 even for the other compilers to run ins_make - if [[ "$use_shell" = "ksh" ]]; then - ./ins_make "$o_ins_make" - else # bash - sed -e s:/bin/ksh:/bin/bash:g ins_make > ins_make.bash - if [[ $(grep jeanzay AA_make.gdef) = "" ]]; then - # Bidouille pour compiler sur ada des vieux modipsl.tar - echo 'Warning jean-zay not in AA_make.gdef' - echo 'Think about updating' - exit 1 - fi - - chmod u=rwx ins_make.bash - ./ins_make.bash "$o_ins_make" - fi # of if [[ "$use_shell" = "ksh" ]] - - echo "install_lmdz.sh MODIPSL_OK $(date)" - - cd "$MODEL/modipsl/modeles/IOIPSL/src" - ioipsllog=$(pwd)/ioipsl.log - echo '##########################################################' - echo 'Compiling IOIPSL, the interface library with Netcdf' - echo '##########################################################' - echo "log file : $ioipsllog $(date)" - - if [[ "$use_shell" = "bash" ]]; then - cp Makefile Makefile.ksh - sed -e s:/bin/ksh:/bin/bash:g Makefile.ksh > Makefile - fi - ### if [[ "$pclinux" = 1 ]]; then - # Build IOIPSL modules and library - $make clean - $make > "$ioipsllog" 2>&1 - if [[ $compiler = "$gfortran" || $compiler = "mpif90" ]]; then - # copy module files to lib - cp -f -- *.mod ../../../lib - fi - # Build IOIPSL tools (ie: "rebuild", if present) - # For IOIPSLv_2_2_2, "rebuild" files are in IOIPSL/tools - rebuild_dir="" - if [[ -f $MODEL/modipsl/modeles/IOIPSL/tools/rebuild ]]; then - rebuild_dir=$MODEL/modipsl/modeles/IOIPSL/tools - elif [[ -d $MODEL/modipsl/modeles/IOIPSL/rebuild ]]; then - rebuild_dir=$MODEL/modipsl/modeles/IOIPSL/rebuild - fi - if [[ $rebuild_dir != "" ]]; then - cd "$rebuild_dir" - # adapt Makefile & rebuild script if in bash - if [[ "$use_shell" = "bash" ]]; then - cp Makefile Makefile.ksh - sed -e s:/bin/ksh:/bin/bash:g Makefile.ksh > Makefile - cp rebuild rebuild.ksh - sed -e 's:/bin/ksh:/bin/bash:g' \ - -e 's:print -u2:echo:g' \ - -e 's:print:echo:g' rebuild.ksh > rebuild - fi - $make clean - $make > "$ioipsllog" 2>&1 - fi - ### fi # of if [[ "$pclinux" = 1 ]] which is commented out - - else # of if [[ $ioipsl = 1 ]] - if [[ ${hostname:0:5} = jean- ]]; then - cd "$MODEL/modipsl" - cd util - if [[ $(grep jeanzay AA_make.gdef) = "" ]]; then - echo 'Warning jean-zay not in AA_make.gdef' - echo 'Think about updating' - exit 1 - fi - ./ins_make "$o_ins_make" - # Compile IOIPSL on jean-zay - cd "$MODEL/modipsl/modeles/IOIPSL/src" - gmake > ioipsl.log - # For IOIPSLv_2_2_2, "rebuild" files are in IOIPSL/tools, so "gmake" in IOIPSL/tools is enough - # For IOIPSLv_2_2_5, "rebuild" files are in a separate IOIPSL/rebuild folder , while "tools" only contains "FCM" - if [[ -f $MODEL/modipsl/modeles/IOIPSL/tools/Makefile ]]; then - cd "$MODEL/modipsl/modeles/IOIPSL/tools" - gmake > ioipsl.log - fi - if [[ -d $MODEL/modipsl/modeles/IOIPSL/rebuild ]]; then - cd "$MODEL/modipsl/modeles/IOIPSL/rebuild" - gmake > ioipsl.log - fi - - fi - echo "install_lmdz.sh ioipsl_OK $(date)" - fi # of if [[ $ioipsl = 1 ]] -fi # of if [[ ! -f $MODEL/modipsl/lib/libioipsl.a ]] - -#========================================================================= -if [[ "$with_xios" = "y" ]]; then - echo '##########################################################' - echo 'Compiling XIOS' - echo '##########################################################' - cd "$MODEL/modipsl/modeles" - xioslog="$(pwd)/xios.log" - #wget http://www.lmd.jussieu.fr/~lmdz/Distrib/install_xios.bash - myget script_install/install_xios.bash - chmod u=rwx install_xios.bash -# following will be recalculated later on once LMDZ is updated -# mysvn=`svnversion LMDZ | egrep || "[0-9]+" 2>/dev/null` - mysvn=$(grep 'Revision: [0-9]' "$MODEL"/Read*.md | awk ' { print $2 } ' 2>/dev/null) - if [[ $svn != "" ]]; then mysvn=$svn; fi - echo "mysvn $mysvn" - - if [[ ${hostname:0:5} = jean- ]]; then - if [[ $mysvn -ge 4619 ]]; then - svn co http://forge.ipsl.jussieu.fr/ioserver/svn/XIOS2/branches/xios-2.6 \ - XIOS - else - svn co http://forge.ipsl.jussieu.fr/ioserver/svn/XIOS2/branches/xios-2.5 \ - XIOS - fi - cd XIOS/arch - svn update - cd .. - echo "Compiling XIOS, start $(date) (it takes about 20 min on Jean-Zay)" - echo "log file: $xioslog" - ./make_xios --prod --arch $arch --job 4 > xios.log 2>&1 - else - ./install_xios.bash -prefix "$MODEL/modipsl/modeles" \ - -netcdf "$ncdfdir" -hdf5 "$ncdfdir" \ - -MPI "$root_mpi" -arch $arch > xios.log 2>&1 - fi # of case Jean-Zay - if [[ -f XIOS/lib/libxios.a ]]; then - echo "XIOS library successfully generated" - echo install_lmdz.sh XIOS_OK "$(date)" - fi -fi - -#============================================================================ -veget_version=false -if [[ "$veget" != 'NONE' ]]; then - cd "$MODEL/modipsl/modeles/ORCHIDEE" - set +e; svn upgrade; set -e - if [[ "$veget" = "CMIP6" ]]; then - veget_version=orchidee2.0 - orchidee_rev=6592 - else # specific orchidee revision newer than CMIP6, on 2_1 or 2_2 branches - veget_version=orchidee2.1 - orchidee_rev=$veget - if [[ $veget -lt 4465 ]]; then - echo 'Stopping, ORCHIDEE version too old, script needs work on ' \ - 'the CPP flags to pass to makelmdz' - exit 1 - fi - set +e - # which branch is my version on? - orcbranch=$(svn log -v -q svn://forge.ipsl.jussieu.fr/orchidee/ -r "$veget" |grep ORCHIDEE |head -1| sed -e 's:ORCHIDEE/.*$:ORCHIDEE:' | awk '{print $2}') - # switch to that branch - echo IF YOU INSTALL ORCHIDEE THE VERY FIRST TIME, ASK for PASSWORD at \ - orchidee-help@listes.ipsl.fr - svn switch -r "$veget" --accept theirs-full "svn://forge.ipsl.jussieu.fr/orchidee/$orcbranch" - if ! (svn log -r "$veget" | grep "$veget"); then - echo 'Cannot update ORCHIDEE as not on the right branch for ORCHIDEE' - exit 1 - fi - set -e - set +e; svn update -r "$veget"; set -e - fi - # Correctif suite debug Jean-Zay - sed -i -e 's/9010 FORMAT(A52,F17.14)/9010 FORMAT(A52,F20.14)/' \ - src_stomate/stomate.f90 - opt_orc="-prod"; if [[ "$optim" = "-debug" ]]; then opt_orc="-debug"; fi - - orchideelog=$(pwd)/orchidee.log - echo '########################################################' - echo 'Compiling ORCHIDEE, the continental surface model ' - echo '########################################################' - echo "Start of the first compilation of orchidee, in sequential mode: $(date)" - echo "log file : $orchideelog" - - ORCHPATH=$(pwd) - export ORCHPATH=$ORCHPATH - xios_orchid="-noxios" - if [[ "$with_xios" = "y" ]]; then - xios_orchid="-xios" - fi - if [[ -d tools ]]; then - ################################################################### - # Pour les experts qui voudraient changer de version d'orchidee. - # Attention : necessite d'avoir le password pour orchidee - - # Correctif suite debug Jean-Zay - if [[ -f src_global/time.f90 ]]; then - sed -i -e 's/CALL tlen2itau/\!CALL tlen2itau/' src_global/time.f90 - fi - ################################################################### - if [[ "$veget_version" == "false" ]]; then - veget_version=orchidee2.0 - fi - cd arch - sed -e s:"%COMPILER .*.$:%COMPILER $compiler":1 \ - -e s:"%LINK .*.$:%LINK $compiler":1 \ - -e s:"%FPP_FLAGS .*.$:%FPP_FLAGS $fpp_flags":1 \ - -e s:"%PROD_FFLAGS .*.$:%PROD_FFLAGS $OPTIM":1 \ - -e s:"%DEV_FFLAGS .*.$:%DEV_FFLAGS $OPTDEV":1 \ - -e s:"%DEBUG_FFLAGS .*.$:%DEBUG_FFLAGS $OPTDEB":1 \ - -e s:"%BASE_FFLAGS .*.$:%BASE_FFLAGS $OPTPREC":1 \ - -e s:"%BASE_LD .*.$:%BASE_LD $BASE_LD":1 \ - -e s:"%ARFLAGS .*.$:%ARFLAGS $ARFLAGS":1 \ - arch-gfortran.fcm > arch-local.fcm - cat < arch-local.path -NETCDF_LIBDIR="-L$ncdfdir/lib -lnetcdff -lnetcdf" -NETCDF_INCDIR="$ncdfdir/include" -IOIPSL_INCDIR="$ORCHPATH/../../lib" -IOIPSL_LIBDIR="$ORCHPATH/../../lib" -XIOS_INCDIR="\$ORCHDIR/../XIOS/inc" -XIOS_LIBDIR="\$ORCHDIR/../XIOS/lib -lxios" -EOF - cd ../ - - echo ./makeorchidee_fcm -j $xios_orchid $opt_orc -parallel none \ - -arch $arch - ./makeorchidee_fcm -j 8 $xios_orchid $opt_orc -parallel none \ - -arch $arch > "$orchideelog" 2>&1 - pwd - else # of "if [[ -d tools ]]" - if [[ -d src_parallel ]]; then - liste_src="parallel parameters global stomate sechiba driver" - if [[ "$veget_version" == "false" ]]; then - veget_version=orchidee2.0 - fi - fi - for d in $liste_src; do - src_d=src_$d - echo src_d "$src_d" - echo ls; ls - if [[ ! -d $src_d ]]; then - echo "Problem orchidee : no $src_d"; exit 1 - fi - cd "$src_d"; \rm -f -- *.mod make; $make clean - $make > "$orchideelog" 2>&1 - if [[ "$compiler" = "$gfortran" || "$compiler" = "mpif90" ]]; then - cp -f -- *.mod ../../../lib - fi - cd .. - done - fi # of "if [[ -d tools ]]" - echo "install_lmdz.sh orchidee_compil_seq_OK $(date)" -fi # of if [[ "$veget" != 'NONE' ]] - - -#============================================================================ -# Ehouarn: the directory name LMDZ* depends on version/tar file... -if [[ -d $MODEL/modipsl/modeles/LMD* ]]; then - echo '###############################################################' - echo 'Preparing LMDZ compilation : arch file, svn switch if needed...' - echo '###############################################################' - cd "$MODEL/modipsl/modeles/"LMD* - LMDZPATH=$(pwd) -else - echo "ERROR: No LMD* directory !!!" - exit 1 -fi - -########################################################### -# For those who want to use fcm to compile via : -# makelmdz_fcm -arch local ..... -############################################################ - - - -cd "$MODEL/modipsl/modeles"/LMDZ* -lmdzlog=$(pwd)/lmdz.log - -################################################################## -# Possibly update LMDZ if a specific svn release is requested -################################################################## - -set +e; svn upgrade; set -e - -if [[ $svn = "last" ]]; then svnopt=""; else svnopt="-r $svn"; fi -if [[ $svn != "" ]]; then - if set +e; svn info | grep -q 'https:'; then - svn switch --relocate https://svn.lmd.jussieu.fr/LMDZ \ - http://svn.lmd.jussieu.fr/LMDZ - fi - svn update "$svnopt" - set -e -fi - -#--------------------------------------------------------------------- -# Retrieve the final svn release number, and adjust compilation -# options accordingly -# If svn not available, will use the svn writen in $MODEL/Readm*md -# For old version it assumes that it is before 4185 (the version -# for which the test was introduced -#--------------------------------------------------------------------- - -set +e; mysvn=$(svnversion . | grep -E || "[0-9]+" 2>/dev/null); set -e -if [[ "$mysvn" = "" ]]; then mysvn=$(grep 'Revision: [0-9]' "$MODEL"/Read*.md | awk ' { print $2 } ' 2>/dev/null); fi -if [[ "$mysvn" = "" ]]; then mysvn=4190; fi - -if [[ "$pclinux" = "1" && ! -f arch/arch-local.path ]]; then - - # create local 'arch' files (if on Linux PC): - cd arch - # arch-local.path file - # test for version as arch.pth file changed format with rev 4426 - if [[ "$mysvn" -gt 4425 ]]; then - cat < arch-local.path -NETCDF_LIBDIR="-L$ncdfdir/lib" -NETCDF_LIB="-lnetcdff -lnetcdf\" -NETCDF_INCDIR=-I$ncdfdir/include -NETCDF95_INCDIR=-I\$LMDGCM/../../include -NETCDF95_LIBDIR=-L\$LMDGCM/../../lib -NETCDF95_LIB=-lnetcdf95 -IOIPSL_INCDIR=-I\$LMDGCM/../../lib -IOIPSL_LIBDIR=-L\$LMDGCM/../../lib -IOIPSL_LIB=-lioipsl -XIOS_INCDIR=-I\$LMDGCM/../XIOS/inc -XIOS_LIBDIR=-L\$LMDGCM/../XIOS/lib -XIOS_LIB="-lxios -lstdc++" -ORCH_INCDIR=-I\$LMDGCM/../../lib -ORCH_LIBDIR=-L\$LMDGCM/../../lib -EOF - else - cat < arch-local.path -NETCDF_LIBDIR="-L$ncdfdir/lib -lnetcdff -lnetcdf" -NETCDF_INCDIR=-I$ncdfdir/include -NETCDF95_INCDIR=\$LMDGCM/../../include -NETCDF95_LIBDIR=\$LMDGCM/../../lib -IOIPSL_INCDIR=\$LMDGCM/../../lib -IOIPSL_LIBDIR=\$LMDGCM/../../lib -XIOS_INCDIR=\$LMDGCM/../XIOS/inc -XIOS_LIBDIR=\$LMDGCM/../XIOS/lib -ORCH_INCDIR=\$LMDGCM/../../lib -ORCH_LIBDIR=\$LMDGCM/../../lib -EOF - fi - - if [[ $pcmac == 1 ]]; then - BASE_LD="$BASE_LD -Wl,-rpath,${ncdfdir}/lib" - else - BASE_LD="$BASE_LD -Wl,-rpath=${ncdfdir}/lib" - fi - # Arch-local.fcm file (adapted from arch-linux-32bit.fcm) - - if [[ $real = r8 ]]; then FPP_DEF=NC_DOUBLE; else FPP_DEF=""; fi - sed -e s:"%COMPILER .*.$:%COMPILER $compiler":1 \ - -e s:"%LINK .*.$:%LINK $compiler":1 \ - -e s:"%PROD_FFLAGS .*.$:%PROD_FFLAGS $OPTIM":1 \ - -e s:"%DEV_FFLAGS .*.$:%DEV_FFLAGS $OPTDEV":1 \ - -e s:"%DEBUG_FFLAGS .*.$:%DEBUG_FFLAGS $OPTDEB":1 \ - -e s:"%BASE_FFLAGS .*.$:%BASE_FFLAGS $OPTPREC":1 \ - -e s:"%FPP_DEF .*.$:%FPP_DEF $FPP_DEF":1 \ - -e s:"%BASE_LD .*.$:%BASE_LD $BASE_LD":1 \ - -e s:"%ARFLAGS .*.$:%ARFLAGS $ARFLAGS":1 \ - arch-linux-32bit.fcm > arch-local.fcm - - cd .. - ### Adapt "bld.cfg" (add the shell): - #whereisthatshell=$(which ${use_shell}) - #echo "bld::tool::SHELL $whereisthatshell" >> bld.cfg - -fi # of if [[ "$pclinux" = 1 ]] -#--------------------------------------------------------------------- -# Option de compilation du rayonnement : depend de $mysvn ><= r4185 -#--------------------------------------------------------------------- -opt_rad="" - -case $rad in - oldrad) iflag_rrtm=0; NSW=2; opt_rad="";; - rrtm) iflag_rrtm=1; NSW=6 - if [[ $mysvn -le 4185 ]]; then - opt_rad="-rrtm true" - else - opt_rad="-rad rrtm" - fi;; - ecrad) iflag_rrtm=2; NSW=6; opt_rad="-rad ecrad";; - *) echo "Only oldrad rrtm ecrad for rad option"; exit 1 -esac - -if [[ $mysvn -le 4185 && $rad = "ecrad" ]]; then - echo "ecrad only available for LMDZ rev starting with 4186 "; exit 1 -fi - -################################################################## - - -if [[ ! -f libf/misc/netcdf95.F90 && ! -d $MODEL/NetCDF95-0.3 ]]; then - cd "$MODEL" - myget src_archives/netcdf/NetCDF95-0.3.tar.gz - tar -xf NetCDF95-0.3.tar.gz - rm NetCDF95-0.3.tar.gz - cd NetCDF95-0.3 - mkdir build - cd build - cmake .. -DCMAKE_BUILD_TYPE=Release -DCMAKE_PREFIX_PATH="$ncdfdir" \ - -DCMAKE_INSTALL_PREFIX="$MODEL/modipsl" - make install - cd "$MODEL/modipsl/modeles/"LMDZ* -fi - -echo '##################################################################' -echo "Preparing script compile.sh for LMDZ compilation" -echo "It will only be run automatically if bench=1/tuto" -echo "Here bench=$bench" -echo '##################################################################' - -if [[ $env_file != "" ]]; then - mv arch/arch-${arch}.env arch/arch-${arch}.orig - \cp -f "$env_file" arch/arch-${arch}.env -fi - -if [[ $compile_with_fcm = 1 ]]; then - makelmdz="makelmdz_fcm $optim -arch $arch -j 8 " -else - makelmdz="makelmdz $optim -arch $arch" -fi - -# sequential compilation -if [[ "$sequential" = 1 ]]; then - echo Sequential compilation command, saved in compile.sh: - echo "./$makelmdz $optim $OPT_MAKELMDZ $optim $opt_rad $opt_cosp " \ - "-d ${grid_resolution} -v $veget_version gcm " - echo "./$makelmdz $optim $OPT_MAKELMDZ $optim $opt_rad $opt_cosp " \ - "-d ${grid_resolution} -v $veget_version gcm " > compile.sh - chmod +x ./compile.sh - if [[ $bench = 1 ]]; then - echo "install_lmdz.sh start_lmdz_seq_compilation $(date)" - echo "log file: $lmdzlog" - ./compile.sh > "$lmdzlog" 2>&1 - echo "install_lmdz.sh end_lmdz_seq_compilation $(date)" - fi -fi # fin sequential - -# compiling in parallel mode -if [[ $parallel != "none" ]]; then - echo '##########################################################' - echo ' Parallel compile ' - echo '##########################################################' - echo "(after saving the sequential libs and binaries)" - cd "$MODEL/modipsl" - tar cf sequential.tar bin/ lib/ - # - # Orchidee - # - cd "$ORCHPATH" - if [[ -d src_parallel && $veget != 'NONE' ]]; then - cd arch - sed \ - -e s:"%COMPILER.*.$:%COMPILER $par_comp":1 \ - -e s:"%LINK.*.$:%LINK $par_comp":1 \ - -e s:"%MPI_FFLAG.*.$:%MPI_FFLAGS $MPI_FLAGS":1 \ - -e s:"%OMP_FFLAG.*.$:%OMP_FFLAGS $OMP_FLAGS":1 \ - -e s:"%MPI_LD.*.$:%MPI_LD $MPI_LD":1 \ - -e s:"%OMP_LD.*.$:%OMP_LD $OMP_LD":1 \ - arch-local.fcm > tmp.fcm - - mv tmp.fcm arch-local.fcm - cd ../ - echo "Compiling ORCHIDEE in parallel mode $(date)" - echo "logfile $orchideelog" - echo "NOTE : to recompile it when necessary, use ./compile_orc.sh " \ - "in modipsl/modeles/ORCHIDEE" - echo ./makeorchidee_fcm -j 8 -clean $xios_orchid $opt_orc \ - -parallel "$parallel" -arch $arch > compile_orc.sh - echo ./makeorchidee_fcm -j 8 $xios_orchid $opt_orc \ - -parallel "$parallel" -arch $arch >> compile_orc.sh - echo echo Now you must also recompile LMDZ, by running ./compile.sh \ - in modeles/LMDZ >> compile_orc.sh - chmod u+x compile_orc.sh - ./makeorchidee_fcm -j 8 -clean $xios_orchid $opt_orc \ - -parallel "$parallel" -arch $arch > "$orchideelog" 2>&1 - ./makeorchidee_fcm -j 8 $xios_orchid $opt_orc -parallel "$parallel" \ - -arch $arch >> "$orchideelog" 2>&1 - echo "End of ORCHIDEE compilation in parallel mode $(date)" - elif [[ $veget != 'NONE' ]]; then - echo '##########################################################' - echo ' Orchidee version too old ' - echo ' Please update to new version ' - echo '##########################################################' - exit 1 - fi #�of [ -d src_parallel && $veget != 'NONE' ] - - # LMDZ - cd "$LMDZPATH" - if [[ $arch = local ]]; then - cd arch - sed -e s:"%COMPILER.*.$:%COMPILER $par_comp":1 \ - -e s:"%LINK.*.$:%LINK $par_comp":1 \ - -e s:"%MPI_FFLAG.*.$:%MPI_FFLAGS $MPI_FLAGS":1 \ - -e s:"%OMP_FFLAG.*.$:%OMP_FFLAGS $OMP_FLAGS":1 \ - -e s:"%ARFLAGS.*.$:%ARFLAGS $ARFLAGS":1 \ - -e "s§%BASE_LD.*.$§%BASE_LD -Wl,-rpath=${root_mpi}/lib:${ncdfdir}/lib§1" \ - -e s:"%MPI_LD.*.$:%MPI_LD $MPI_LD":1 \ - -e s:"%OMP_LD.*.$:%OMP_LD $OMP_LD":1 \ - arch-local.fcm > tmp.fcm - mv tmp.fcm arch-local.fcm - cd ../ - fi - rm -f compile.sh - echo resol="$grid_resolution" >> compile.sh - if [[ ${hostname:0:5} = jean- && "$cosp" = "v2" ]]; then - - echo "LMDZ compilation command in parallel mode, saved in compile.sh, is :" - echo "(ATTENTION le probleme de cospv2 sur jean-zay en mode prod " \ - "n est pas corrige ! )" - # ATTENTION le probleme de cospv2 sur jean-zay en mode prod n - # est pas corrige - echo ./"$makelmdz" -dev $optim "$OPT_MAKELMDZ" "$opt_rad" "$opt_cosp" \ - "$opt_makelmdz_xios" -d \$resol -v $veget_version -mem \ - -parallel "$parallel" gcm >> compile.sh - echo ./"$makelmdz" -dev $optim "$OPT_MAKELMDZ" "$opt_rad" "$opt_cosp" \ - "$opt_makelmdz_xios" -d \$resol -v $veget_version -mem \ - -parallel "$parallel" gcm - else - echo ./"$makelmdz" $optim "$OPT_MAKELMDZ" "$opt_rad" "$opt_cosp" \ - "$opt_makelmdz_xios" -d \$resol -v $veget_version -mem \ - -parallel "$parallel" gcm >> compile.sh - echo ./"$makelmdz" $optim "$OPT_MAKELMDZ" "$opt_rad" "$opt_cosp" \ - "$opt_makelmdz_xios" -d \$resol -v $veget_version -mem \ - -parallel "$parallel" gcm - fi - chmod +x ./compile.sh - - if [[ $bench = 1 ]]; then - echo "Compiling LMDZ in parallel mode $(date), LMDZ log file: $lmdzlog" - ./compile.sh > "$lmdzlog" 2>&1 - fi - -fi # of if [[ $parallel != "none" ]] - - -################################################################## -# Verification du succes de la compilation -################################################################## - -# Recherche de l'executable dont le nom a change au fil du temps ... -# suffix contains radiative option starting with revision 4186 -if [[ $mysvn -ge 4186 ]]; then suff_exe=_${rad}${suff_exe}; fi -gcm="" - -for exe in gcm.e bin/gcm_${grid_resolution}_phylmd${suff_exe}${suff_orc}.e -do - if [[ -f $exe ]]; then gcm=$exe; fi -done - -if [[ "$gcm" = "" ]]; then - if [[ $bench = 1 ]]; then - echo 'Compilation failed !! Cannot run the benchmark;' - exit 1 - else - echo 'Compilation not done (only done when bench=1)' - fi -else - echo '##########################################################' - echo "Compilation successfull !! $(date)" - echo '##########################################################' - echo "The executable is $gcm" -fi - -################################################################## -# Below, we run a benchmark if bench=1 or tuto -################################################################## - -if [[ $bench = tuto ]]; then - myget Training/tutorial.tar; tar xf tutorial.tar; cd TUTORIAL - ./init.sh - -elif [[ $bench = 1 && ! -d BENCH${grid_resolution} ]]; then - # TOUTE CETTE SECTION DEVRAIT DISPARAITRE POUR UNE COMMANDE - # OU DES BENCHS PAR MOTS CLES COMME tuto - - echo '##########################################################' - echo ' Running a test run ' - echo '##########################################################' - - \rm -rf "BENCH${grid_resolution}" - bench=bench_lmdz_${grid_resolution} - echo "install_lmdz.sh before bench download $(date)" - #wget http://www.lmd.jussieu.fr/~lmdz/Distrib/$bench.tar.gz - myget "3DBenchs/$bench.tar.gz" - echo "install_lmdz.sh after bench download $(date)" - tar xf "$bench.tar.gz" - - if [[ "$cosp" = "v1" || "$cosp" = "v2" ]]; then - cd "BENCH${grid_resolution}" - # copier les fichiers namelist input et output our COSP - cp ../DefLists/cosp*_input_nl.txt . - cp ../DefLists/cosp*_output_nl.txt . - # Activer la cles ok_cosp pour tourner avec COSP - sed -e 's@ok_cosp=n@ok_cosp=y@' config.def > tmp - \mv -f tmp config.def - cd .. - fi - - if [[ -n "$physiq" ]]; then - cd "BENCH${grid_resolution}" - if [[ -f "physiq.def_${physiq}" ]]; then - cp "physiq.def_${physiq}" physiq.def - echo using "physiq.def_${physiq}" - else - echo using standard physiq.def - fi - cd .. - else - echo using standard physiq.def - fi - - if [[ "$with_xios" = "y" ]]; then - cd "BENCH${grid_resolution}" - cp ../DefLists/iodef.xml . - cp ../DefLists/context_lmdz.xml . - cp ../DefLists/field_def_lmdz.xml . - # A raffiner par la suite - echo "A FAIRE : Copier les *xml en fonction de l option cosp" - cp ../DefLists/field_def_cosp*.xml . - cp ../DefLists/file_def_hist*xml . - # adapt iodef.xml to use attached mode - sed -e 's@"using_server" type="bool">true@"using_server" type="bool">false@' \ - iodef.xml > tmp - \mv -f tmp iodef.xml - - # and convert all the enabled="_AUTO_" (for libIGCM) to enabled=.FALSE. - # except for histday - for histfile in file_def_hist*xml - do - if [[ "$histfile" = "file_def_histday_lmdz.xml" ]]; then - sed -e 's@enabled="_AUTO_"@type="one_file" enabled=".TRUE."@' \ - "$histfile" > tmp - \mv -f tmp "$histfile" - sed -e 's@output_level="_AUTO_"@output_level="5"@' "$histfile" \ - > tmp - \mv -f tmp "$histfile" - sed -e 's@compression_level="2"@compression_level="0"@' \ - "$histfile" > tmp - \mv -f tmp "$histfile" - else - sed -e 's@enabled="_AUTO_"@type="one_file" enabled=".FALSE."@' \ - "$histfile" > tmp - \mv -f tmp "$histfile" - fi - done - # and add option "ok_all_xml=y" in config.def - echo "### XIOS outputs" >> config.def - echo 'ok_all_xml=.true.' >> config.def - - #activer les sorties pour Cosp - if [[ "$cosp" = "v1" ]]; then - sed -i'' -e 's@enabled=".FALSE."@enabled=".TRUE."@' \ - -e 's@output_level="_AUTO_"@output_level="5"@' \ - -e 's@compression_level="2"@compression_level="0"@' \ - file_def_histdayCOSP_lmdz.xml - fi - if [[ "$cosp" = "v2" ]]; then - sed -e 's@compression_level="2"@compression_level="0"@' file_def_histdayCOSPv2_lmdz.xml - for type_ in hf day mth; do - file=file_def_hist${type_}COSP - sed -i'' -e 's@src="./'${file}'_lmdz.xml"@src="./'${file}'v2_lmdz.xml"@' context_lmdz.xml - done - sed -i '' -e 's@field_def_cosp1.xml@field_def_cospv2.xml@' field_def_lmdz.xml - fi - - cd .. - fi - - # Cas Bensh avec ecrad - if [[ "$rad" = "ecrad" ]]; then - cd "BENCH${grid_resolution}" - cp ../DefLists/namelist_ecrad . - cp -r ../libf/phylmd/ecrad/data . - cd .. - fi - - # Adjusting bench physiq.def to radiative code chosen - cd "BENCH${grid_resolution}" - sed -e 's/iflag_rrtm=.*.$/iflag_rrtm='$iflag_rrtm'/' \ - -e 's/NSW=.*.$/NSW='$NSW'/' physiq.def > tmpdef - \mv tmpdef physiq.def - cd .. - - cp "$gcm" "BENCH${grid_resolution}/gcm.e" - - cd "BENCH${grid_resolution}" - # On cree le fichier bench.sh au besoin - # Dans le cas 48x36x39 le bench.sh existe deja en parallele - - if [[ "$grid_resolution" = "48x36x39" ]]; then - echo On ne touche pas au bench.sh - # But we have to adapt "run_local.sh" for $mpirun - sed -e "s@mpirun@$mpirun@g" run_local.sh > tmp - mv -f tmp run_local.sh - chmod u=rwx run_local.sh - elif [[ "${parallel:0:3}" = "mpi" ]]; then - # Lancement avec deux procs mpi et 2 openMP - echo "export OMP_STACKSIZE=800M" > bench.sh - if [[ "${parallel:4:3}" = "omp" ]]; then - echo "export OMP_NUM_THREADS=2" >> bench.sh - fi - if [[ "$cosp" = "v1" || "$cosp" = "v2" ]]; then - if [[ ${hostname:0:5} = jean- ]]; then - chmod +x ../arch.env - ../arch.env - echo "ulimit -s 2000000" >> bench.sh - else - echo "ulimit -s 200000" >> bench.sh - fi - else - echo "ulimit -s unlimited" >> bench.sh - fi - if [[ ${hostname:0:5} = jean- ]]; then - . ../arch/arch-${arch}.env - echo "srun -n 2 && $idris_acct@cpu gcm.e > listing 2>&1" \ - >> bench.sh - else - echo "$mpirun -np 2 gcm.e > listing 2>&1" >> bench.sh - fi - # Add rebuild, using reb.sh if it is there - cat <> bench.sh -if [[ -f reb.sh ]]; then - ./reb.sh histday; ./reb.sh histmth; ./reb.sh histhf; - ./reb.sh histins; ./reb.sh stomate_history; - ./reb.sh sechiba_history; ./reb.sh sechiba_out_2 -fi -EOF - else - echo "./gcm.e > listing 2>&1" > bench.sh - fi - # Getting orchidee stuff - if [[ $veget == 'CMIP6' ]]; then - #echo 'myget 3DBenchs/BENCHorch11.tar.gz' - #myget 3DBenchs/BENCHorch11.tar.gz - #tar xvzf BENCHorch11.tar.gz - echo 'myget 3DBenchs/BENCHCMIP6.tar.gz' - myget 3DBenchs/BENCHCMIP6.tar.gz - tar xvzf BENCHCMIP6.tar.gz - sed -e "s:VEGET=n:VEGET=y:" config.def > tmp - mv -f tmp config.def - if [[ "$with_xios" = "y" ]]; then - cp ../../ORCHIDEE/src_xml/context_orchidee.xml . - echo '' \ - > add.tmp - cp ../../ORCHIDEE/src_xml/field_def_orchidee.xml . - cp ../../ORCHIDEE/src_xml/file_def_orchidee.xml . - cp ../../ORCHIDEE/src_xml/file_def_input_orchidee.xml . - if [[ -f ../../ORCHIDEE/src_xml/context_input_orchidee.xml ]]; then - cp ../../ORCHIDEE/src_xml/context_input_orchidee.xml . - echo '' >> add.tmp - fi - sed -e '/id="LMDZ"/r add.tmp' iodef.xml > tmp - mv tmp iodef.xml - sed -e'{/sechiba1/ s/enabled="_AUTO_"/type="one_file" enabled=".TRUE."/}' \ - file_def_orchidee.xml > tmp - \mv -f tmp file_def_orchidee.xml - sed -e 's@enabled="_AUTO_"@type="one_file" enabled=".FALSE."@' \ - file_def_orchidee.xml > tmp - \mv -f tmp file_def_orchidee.xml - sed -e 's@output_level="_AUTO_"@output_level="1"@' \ - file_def_orchidee.xml > tmp - \mv -f tmp file_def_orchidee.xml - sed -e 's@output_freq="_AUTO_"@output_freq="1d"@' \ - file_def_orchidee.xml > tmp - \mv -f tmp file_def_orchidee.xml - sed -e 's@compression_level="4"@compression_level="0"@' \ - file_def_orchidee.xml > tmp - \mv -f tmp file_def_orchidee.xml - sed -e 's@XIOS_ORCHIDEE_OK = n@XIOS_ORCHIDEE_OK = y@' \ - orchidee.def > tmp - \mv -f tmp orchidee.def - fi - fi - - if [[ -f ../arch.env ]] - then - source ../arch.env - fi - - echo EXECUTION DU BENCH - set +e - date; ./bench.sh > out.bench 2>&1; date - set -e - tail listing - - - echo "##########################################################" - echo "Simulation finished in $(pwd)" - echo "You have compiled with:" - cat ../compile.sh - if [[ $parallel = "none" ]]; then - echo "You may re-run it with : cd $(pwd) ; gcm.e" - echo "or ./bench.sh" - else - echo "You may re-run it with : " - echo "cd $(pwd) ; ./bench.sh" - # echo 'ulimit -s unlimited' - # echo 'export OMP_NUM_THREADS=2' - # echo 'export OMP_STACKSIZE=800M' - # echo "$mpirun -np 2 gcm.e " - fi - echo '##########################################################' - -fi # bench - - -################################################################# -# Installation eventuelle du 1D -################################################################# - -if [[ $SCM = 1 ]]; then - cd "$MODEL" - myget 1D/1D.tar.gz - tar xf 1D.tar.gz - cd 1D - if [[ $rad = oldrad ]]; then - sed -i'' -e 's/^rad=.*$/rad=oldrad/' run.sh - sed -i'' -e 's/^rad=.*$/rad=oldrad/' bin/compile - fi - echo "Running 1D/run.sh, log in $(pwd)/run1d.log" - ./run.sh > "$(pwd)/run1d.log" 2>&1 -fi - - -################################################################# -# sauvegarde des options veget pour utilisation eventuelle tutorial_prod -################################################################# -cd "$MODEL/modipsl/modeles" -#echo surface_env file created in $MODEL -cat <surface_env -veget=$veget -opt_veget="-v $veget_version" -orchidee_rev=$orchidee_rev -suforch=$suff_orc -EOF Index: DZ.3.3/trunk/README =================================================================== --- /LMDZ.3.3/trunk/README (revision 6163) +++ (revision ) @@ -1,78 +1,0 @@ -13/04 -Installation de LMDZ.3.3 sur Etoile avec M-A Foujols -Installation de ioipsl sur Etoile (M-A Foujols) -Modification de create_gcm, makegcm en consequence -Pour tenir compte de la taille des fichiers *.o et -des librairies en F90 on compile seulement les routines -utilisant ioipsl en F90, le reste en F77. <- en cours - -Creation d'un directory phylmd.old (sans traceurs) -dans phylmd suppression des *new.F pour eviter erreurs. <- en cours -modif dans dyn3d pour compatibilite 2 phylmd. <- OK - -14/04 -Dans phylmd.old mise a jour physiq.F / ioipsl F90 -dans phylmd suppression des *new.F <- OK -Dans clmain il faut supprimer cltrac car appels differents -dans phylmd et plylmd.old - -Modification de makegcm et create_make_gcm (M-A Foujols) -pour tenir compte de la compil F90 de ioipsl. <-OK - -15/04 -Creation dans RUN de job_init et job_gcm pour lancer en batch -les creations de fichiers et le modele. - -16/04 -correction de cltracrn.F (rajout d'une variable d'appel). -correction de makegcm / option debug - -19/04 -correction de clafis.F initialisation de dq. -appel de fluxstoke.F dans gcm.F - -20/04 -suppression de idadvtrac -rajout de iapp_phy et iapp_dyn les pdt de sorties dans la dyn et la phy. - -21/04 -correction de fluxstoke.F rajout save phic,tetac -correction fisrtilp_tr.F manque init de zoliq si rneb=0 - -22/04 -iapp_dyn rempalce par istdyn et iapp_phy par istphy - -23/04 -Correction dans physiq dans sorties ioipsl = supression de d_q_ajs -qui n'est plus def a ce niveau. -corrections alex dans clmain, fisrtilpr_tr, physiq - -26/04 -correction cltracrn / def de hsoltr, tautr, vdept, scavtr, masktr, fshtr -?? a quoi set lessi* dans phytrac -phytrac rajout du calcul de trs (+SAVE) - -rajout de assign pour option CRAY pour ecriture de fluxmass et physique. - -Ioipsl dans la partie traceurs (histrac.nc) => phytrac en F90 - -27/04 -plante quand option convection activee = dans nflxtr -zmfxd mal def -dx calcule seulement pour klev - -03/05 -correction de nflxtr -zmfdx et dx - -A FAIRE -Verifier que run on-line off-line donne meme resultat / traceurs - -Verifier que run sans traceurs donne meme resultat meteo / traceurs - -Phaser / ioispl de L. Fairhead. - -Ioipsl dans la partie traceurs quand OK pour entrees-sorties -- dans phylmd = phystoke (Kz et flux convect) -- dans dyn3d = fluxstoke (flux masse, cumul) -(garder poss. anciennes sorties.) Index: DZ.3.3/trunk/create_make_gcm =================================================================== --- /LMDZ.3.3/trunk/create_make_gcm (revision 6163) +++ (revision ) @@ -1,227 +1,0 @@ -#!/bin/sh -#set -xv -machine=`hostname` -os=`uname` -gcm=`pwd` -libf=$gcm/libf -libo=$gcm/libo -CRAY=0 -if [ "$machine" = "atlas" -o "$machine" = "etoile" -o "$machine" = "axis" ] ; then - CRAY=1 -fi -XNEC=0 -if [ "$machine" = "rhodes" ] ; then - XNEC=1 -fi -VPP=0 -if [ "$machine" = "nymphea0" ] ; then - VPP=1 -fi -# -echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" -echo "# Definitions de Macros pour Make" -echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" -echo -echo "# Repertoires :" -echo -echo "GCM = "$gcm -if [ "$CRAY" = '0' ] ; then - echo "MACHINE = sun4" -fi -echo 'LIBF = $(GCM)/libf' -if [ "$CRAY" = '0' ] ; then -# echo 'LIBO = $(GCM)/libo/$(MACHINE)' - echo 'LIBO = $(LIBOGCM)/$(MACHINE)' -else - echo 'LIBO = $(GCM)/libo' -fi -#echo 'LOCAL_DIR=$(GCM)' -#echo $localdir -echo "LOCAL_DIR=`echo $localdir`" -echo 'BIBIO = $(LIBF)/bibio' -echo "FILTRE = filtre" -echo "PHYS = " -echo "DYN = dyn " -echo 'LIBPHY = $(LIBO)/libphy$(PHYS).a' -echo 'DIRMAIN=dyn$(DIM)d' -echo 'RM=rm' -echo -echo "OPLINK = " -echo -echo '# Les differentes librairies pour l"edition des liens:' -echo -if [ "$XNEC" = '1' ] ; then - echo 'dyn3d = $(LIBO)/libsxdyn3d.a $(LIBO)/libsx$(FILTRE).a' - echo 'dyn2d = $(LIBO)/libsxdyn2d.a' - echo 'dyn1d = $(LIBO)/libsxdyn1d.a' - echo 'L_DYN = -lsxdyn$(DIM)d' - echo 'L_FILTRE = -lsx$(FILTRE)' - echo 'L_PHY = -lsxphy$(PHYS) ' - echo 'L_BIBIO = -lsxbibio' - echo 'L_ADJNT =' -else - echo 'dyn3d = $(LIBO)/libdyn3d.a $(LIBO)/lib$(FILTRE).a' - echo 'dyn2d = $(LIBO)/libdyn2d.a' - echo 'dyn1d = $(LIBO)/libdyn1d.a' - echo 'L_DYN = -ldyn$(DIM)d' - echo 'L_FILTRE = -l$(FILTRE)' - echo 'L_PHY = -lphy$(PHYS) ' - echo 'L_BIBIO = -lbibio' - echo 'L_ADJNT =' -fi -echo -echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" -echo "# Option de compilation FORTRAN" -echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" -echo - echo 'COMPILE = $(F77) $(OPTIM) $(INCLUDE) -c' - echo 'COMPILE90 = $(F90) $(OPTIM90) $(INCLUDE) -c' - echo 'COMPTRU90 = $(F90) $(OPTIMTRU90) $(INCLUDE) -c' - echo "LINK = $LINK" - echo "AR = $AR" -echo -echo -echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" -echo "# Creation des differents executables" -echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" -echo -echo "# Executables:" -echo "# ------------" -echo -echo "PROG = code" -echo -echo 'main : $(DYN) bibio phys $(OPTION_DEP) ' -echo ' cd $(LIBO) ; $(RANLIB) lib*.a ; cd $(GCM) ;\' -echo ' cd $(LOCAL_DIR); \' -echo ' $(COMPILE90) $(LIBF)/$(DIRMAIN)/$(PROG).F -o $(PROG).o ; \' -if [ "$CRAY" = '0' ] ; then -echo ' $(LINK) $(PROG).o -L$(LIBO) $(L_DYN) $(L_ADJNT) $(L_FILTRE) $(L_PHY) $(L_BIBIO) $(L_DYN) $(OPLINK) $(OPTION_LINK) -o $(LOCAL_DIR)/$(PROG).e ; $(RM) $(PROG).o ' -else -echo ' $(LINK) $(PROG).o -L$(LIBO) $(L_DYN) $(L_ADJNT) $(L_FILTRE) $(L_PHY) $(L_BIBIO) $(OPLINK) $(OPTION_LINK) -o $(LOCAL_DIR)/$(PROG).e ; $(RM) $(PROG).o ' -fi -echo -echo 'dyn : $(LIBO)/libdyn$(DIM)d.a $(FILTRE)$(DIM)d' -echo -echo 'phys : $(LIBPHY)' -echo -echo 'bibio : $(LIBO)/libbibio.a' -echo -echo 'adjnt : $(LIBO)/libadjnt.a' -echo -echo '$(FILTRE)3d : $(LIBO)/lib$(FILTRE).a' -echo -echo '$(FILTRE)2d :' -echo -echo '$(FILTRE)1d :' -echo -echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" -echo "# Contenu des differentes bibliotheques" -echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" -echo -echo -cd $libf >/dev/null 2>&1 -for diri in ` ls ` -do - if [ -d $diri ] ; then - if [ "`ls $diri/*.F`" != "" ] || [ "`ls $diri/*.F90`" != "" ] ; then - cd $diri >/dev/null 2>&1 - echo - listlib="" - for i in `ls *.F` - do - fili=`basename $i .F` - test=` ( head $i | grep ' PROGRAM' ) ` - if [ "$test" = "" ] ; then - listlib=$listlib" "$fili - fi - done - for i in `ls *.F90` - do - fili=`basename $i .F90` - test=` ( head $i | grep ' PROGRAM' ) ` - if [ "$test" = "" ] ; then - listlib=$listlib" "$fili - fi - done -# - echo - echo - echo "#=======================================================================" - echo "# Contenu de la bibliotheque correspondant au Directory "$diri - echo "#=======================================================================" - echo - for fili in $listlib - do - echo '$(LIBO)/lib'$diri".a : " '$(LIBO)/lib'$diri".a("$fili".o)" - echo - done - echo '.PRECIOUS : $(LIBO)/lib'$diri'.a' - echo - echo - echo "# Compilation des membres de la bibliotheque lib"$diri".a" - echo - for fili in $listlib - do - if [ -f $fili.F90 ] ; then - trufile=$fili.F90 - else - trufile=$fili.F - fi - F90=0 ; egrep -i '^ *use ' $trufile > /dev/null 2>&1 && F90=1 - egrep -i '^ *module ' $trufile > /dev/null 2>&1 && F90=1 - str1='$(LIBO)/lib'$diri'.a('$fili'.o) : $(LIBF)/'$diri/$trufile - [ "$fili" = "chem.subs" ] && str1=$str1' $(LIBF)/'$diri/chem.mods.F - for stri in ` ( sed -n "/\#include/s/\#include//p" $trufile | sed 's/\"//g' ; egrep -i '^ *use ' $trufile | awk ' { print $2 } ' ) ` - do - stri=`echo $stri | tr [A-Z] [a-z]` - if [ -f $stri ] ; then - echo $str1 \\ - str1='$(LIBF)/'$diri'/'$stri - else - if [ -f $stri.F ] || [ -f $stri.F90 ] ; then - echo $str1 \\ - str1='$(LIBO)/lib'$diri'.a('$stri'.o)' - else - for dirinc in dyn3d grid bibio - do - if [ -f ../$dirinc/$stri ] ; then - echo $str1 \\ - str1='$(LIBF)/'`cd .. ; ls */$stri | head -1` - fi - done - fi - fi - done - echo $str1 - if [ "$F90" -eq '0' ] ; then - echo ' cd $(LOCAL_DIR); \' - echo ' $(COMPILE) $(LIBF)/'$diri'/'$trufile' ; \' - else - echo ' cd $(LOCAL_DIR); \' - if [ -f $fili.F90 ] ; then - echo ' $(COMPTRU90) $(LIBF)/'$diri'/'$trufile' ; \' - else - echo ' $(COMPILE90) $(LIBF)/'$diri'/'$trufile' ; \' - fi - MODU=0; egrep -i ' *module ' $trufile> /dev/null 2>&1 && MODU=1 - if [ "$MODU" -eq '1' -a "$CRAY" != '1' ] ; then - if [ "$os" = 'UNIX_System_V' ] ; then - echo ' cp $(MOD_LOC_DIR)/*.$(MOD_SUFFIX) $(LIBO)/ ; \' - else - echo ' mv $(MOD_LOC_DIR)/'$fili'.$(MOD_SUFFIX) $(LIBO)/'$fili'.$(MOD_SUFFIX) ; \' - fi - fi - fi - if [ "$XNEC" -eq '1' ] ; then - echo ' sxar r $(LIBO)/libsx'$diri'.a '$fili'.o ; \' - fi - echo ' $(AR) r $(LIBO)/lib'$diri'.a '$fili'.o ; $(RM) '$fili'.o ; \' - echo ' cd $(GCM)' - echo - done -# - echo - cd $libf - fi - fi -done Index: DZ.3.3/trunk/gcm.def =================================================================== --- /LMDZ.3.3/trunk/gcm.def (revision 6163) +++ (revision ) @@ -1,81 +1,0 @@ -# -## $Header$ -# -## Jour de l'etat initial ( = 350 si 20 Decembre ,par expl. ,comme ici ) -dayref=1 -## Annee de l'etat initial ( avec 4 chiffres ) -anneeref=1979 -## Nombre de jours d'integration -nday=5 -## nombre de pas par jour (multiple de iperiod) ( ici pour dt = 1 min ) -day_step=480 -## periode pour le pas Matsuno (en pas) -iperiod=5 -## periode de sortie des variables de controle (en pas) -iconser=240 -## periode d'ecriture du fichier histoire (en jour) -iecri=1 -## periode de stockage fichier histmoy (en jour) -periodav=1. -## periode de la dissipation (en pas) -idissip=5 -## choix de l'operateur de dissipation (star ou non star ) -lstardis=y -## nombre d'iterations de l'operateur de dissipation gradiv -nitergdiv=1 -## nombre d'iterations de l'operateur de dissipation nxgradrot -nitergrot=2 -## nombre d'iterations de l'operateur de dissipation divgrad -niterh=2 -## temps de dissipation des plus petites long.d ondes pour u,v (gradiv) -tetagdiv=3600. -## temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot) -tetagrot=3600. -## temps de dissipation des plus petites long.d ondes pour h ( divgrad) -tetatemp=3600. -## coefficient pour gamdissip -coefdis=0. -## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog) -purmats=n -## avec ou sans physique -physic=y -## periode de la physique (en pas) -iphysiq=10 -## frequence (en jours ) de l'ecriture du fichier histphy -ecritphy=30 -## Cycle diurne ou non -cycle_diurne=y -## Soil Model ou non -soil_model=y -## Choix ou non de New oliq -new_oliq=y -## Orodr ou non pour l orographie -ok_orodr=y -## Orolf ou non pour l orographie -ok_orolf=y -## Si = .T. , lecture du fichier limit avec la bonne annee -ok_limitvrai=n -## Nombre d'appels des routines de rayonnements ( par jour) -nbapp_rad=12 -## Flag pour la convection ( 1 pour LMD, 2 pour Tiedtke, 3 CCM ) -iflag_con=2 -## longitude en degres du centre du zoom -clon=0. -## latitude en degres du centre du zoom -clat=0. -## facteur de grossissement du zoom,selon longitude -grossismx=1.0 -## facteur de grossissement du zoom ,selon latitude -grossismy=1.0 -## Fonction f(y) hyperbolique si = .true. , sinon sinusoidale -fxyhypb=y -## extension en longitude de la zone du zoom ( fraction de la zone totale) -dzoomx=0.0 -## extension en latitude de la zone du zoom ( fraction de la zone totale) -dzoomy=0.0 -##raideur du zoom en X -taux=3. -##raideur du zoom en Y -tauy=3. -## Fonction f(y) avec y = Sin(latit.) si = .true. , sinon y = latit. -ysinus=y Index: DZ.3.3/trunk/initialisation =================================================================== --- /LMDZ.3.3/trunk/initialisation (revision 6163) +++ (revision ) @@ -1,52 +1,0 @@ -#!/bin/csh -vx - -# en local - -set dir=~lmdz/DATA/Netcdf -set cp="ln -s" -setenv LIBOGCM /tmp2/fairhead - -#=================================================================== -# Fichiers necessaires a la creations des etats initiaux -#=================================================================== -$cp $dir/ECDYN.nc . # Etat initial dynamique -$cp $dir/ECPHY.nc . # Etat initial physique -$cp $dir/Relief.nc . # Fichier de relief -$cp $dir/Rugos.nc . # Fichier rugosite sur un an -$cp $dir/AMIP.nc . # Fichier SST/Seaice sur un an -$cp $dir/Albedo.nc . # Fichier albedo sur un an -#=================================================================== - - -if ( ! -f run.def ) then -echo Il faut avoir sur le repertoire courant un fichier run.def -exit -endif - -# Compilation etat0_netcdf.e -# ========================== -$LMDGCM/makegcm etat0 - -# Creation des etats initiaux -# =========================== -echo "execution de etat0.e, sortie dans etat0.l" -./etat0.e >! etat0.l - -# Compilation limit_netcdf.e -# ========================== -$LMDGCM/makegcm limit_netcdf - -# Creation des conditions aux limites -# =================================== -echo "execution de limit_lmd.e, sortie dans limit_lmd.l" -./limit_netcdf.e >! limit_netcdf.l - -# Conditions aux limites pour les especes traces -# ============================================== -#if ( -f limit_tr.e ) then -#$cp $dir/so2_1x1 so2_1x1 -#limit_tr.e > limit_tr.l -#else -#echo Vous n avez pas cree le fichier limit_tr utilise par les versions -#echo du modele avec traceurs -#endif Index: DZ.3.3/trunk/makegcm =================================================================== --- /LMDZ.3.3/trunk/makegcm (revision 6163) +++ (revision ) @@ -1,814 +1,0 @@ -#!/bin/csh -# $Header$ -#set verbose echo -######################################################################## -# options par defaut pour la commande make -######################################################################## -set dim="72x45x19" -set physique=lmd -set phys="PHYS=$physique" -set include='-I$(LIBF)/grid -I$(LIBF)/bibio -I. ' -set ntrac = 4 -set filtre=filtrez -set grille=reg -set couple=false -set veget=true -######################################################################## -# path a changer contenant les sources et les objets du modele -######################################################################## - -###### VERSION LMDZ.3.3 -# set LMDGCM=$HOME/LMDZ.3.3 -# - -setenv localdir "`pwd`" -set MODIPSL=0 -echo $localdir | grep modipsl >& /dev/null -if ( ! $status ) then - set MODIPSL=1 - setenv LMDGCM $localdir - cd ../.. - setenv LIBOGCM "`pwd`/lib" - setenv IOIPSLDIR $LIBOGCM - setenv MODIPSLDIR $LIBOGCM - cd $localdir - if ( `hostname` == rhodes ) then - set NCDFINC=`grep sxnec ../../util/AA_make.gdef| grep NCDF_INC|sed -e "s/^.* =//"` - set NCDFLIB=`grep sxnec ../../util/AA_make.gdef| grep NCDF_LIB|sed -e 's/^.* =//'` - else - if ( `hostname` == nymphea0 ) then - set NCDFINC=`grep fjvpp ../../util/AA_make.gdef| grep NCDF_INC|sed -e "s/^.* =//"` - set NCDFLIB=`grep fjvpp ../../util/AA_make.gdef| grep NCDF_LIB|sed -e 's/^.* =//'` - else - echo 'Probleme de definition des variables NCDFINC et NCDFLIB' - endif - endif -else - if ( ! $?LMDGCM ) then - echo You must initialize the variable LMDGCM in your environnement - echo for instance: "setenv LMDGCM /usr/myself/supergcm" in .cshrc - exit - endif - if ( ! $?LIBOGCM ) then - set LIBOGCM=$LMDGCM/libo - endif - if ( ! $?IOIPSLDIR ) then - echo You must initialize the variable IOIPSLDIR in your environnement - echo for instance: "setenv IOIPSLDIR /usr/myself/ioipsl" in .cshrc - exit - else - setenv MODIPSLDIR $IOIPSLDIR - endif - if ( ! $?NCDFLIB ) then - echo You must initialize the variable NCDFLIB in your environnement - echo for instance: "setenv NCDFLIB /usr/myself/netcdf" in .cshrc - exit - endif - if ( ! $?NCDFINC ) then - echo You must initialize the variable NCDFINC in your environnement - echo for instance: "setenv NCDFINC /usr/myself/netcdf" in .cshrc - exit - endif -endif -set model=$LMDGCM -set libo=$LIBOGCM - -######################################################################## -# Les differentes platformes reconnues -######################################################################## - -set HP=0 -set IBM=0 -set SUN=0 -set VPP=0 -set CRAY=0 -set DEC=0 -set LINUX=0 -set NEC=0 -set XNEC=0 -if ( `uname` == HP-UX ) then - set machine=HP - set HP=1 -else if (`uname` == UNIX_System_V ) then - set machine=VPP - set VPP=1 -else if (`uname` == SunOS ) then - set machine=SUN - set SUN=1 -else if ( `uname` == AIX ) then - set machine=IBM - set IBM=1 -else if ( `uname` == OSF1 ) then - set machine=ALPHA - set DEC=1 -else if ( `uname` == Linux ) then - set machine=LINUX - set LINUX=1 -else if ( `hostname` == atlas || `hostname` == axis || `hostname` == etoile ) then - set machine=CRAY - set CRAY=1 -else if ( `uname` == SUPER-UX ) then - set machine=NEC - set NEC=1 -else if ( `hostname` == rhodes) then - set machine=XNEC - set XNEC=1 -else - echo Vous travaillez sur une machine non prevue par le reglement - exit -endif - -if ( ! -d $libo ) then - mkdir $libo -endif - - -if $VPP then -set netcdf=netcdf_v -else -set netcdf=netcdf -endif -######################################################################## -# Quelques initialisations de variables du shell. -######################################################################## - -set dyn= -set opt_link="" -set adjnt="" -set opt_dep="" - -set optim="" -set optimbis="" -set optim90="" -set oplink="" - -######################################################################## -# Optimisations par defaut suivant les machines -######################################################################## - -echo "Optimisations par defaut suivant les machines" -set libf=$model/libf -#setenv localdir "LOCAL_DIR=`pwd`" -#setenv localdir "`pwd`" -cd $model -if $CRAY then -# set optim="-Wf'-ei' -dp -Wf'-a static'" - set optimbis=" -DCRAY " - set optim90="-Wp'-P' -DCRAY -p$IOIPSLDIR "'-p$(LIBO) -eiv ' - set optim="$optim90" - set oplink="-Wl'-DSTACK=128 -f indef' -L$IOIPSLDIR -lioipsl -L$NCDFLIB -lnetcdf " - set mod_loc_dir=" " - set mod_suffix=" " -else if $SUN then - set optim=" -fast " - set optimbis=" " - set optim90=" -fast -fixed " - set optimtru90=" -fast -free " - set opt_link="-lf77compat -L$MODIPSLDIR -lioipsl -L$NCDFLIB -lnetcdf " - set mod_loc_dir=$localdir - set mod_suffix=mod -else if $HP then - set optim=" +U77 -O +E1 " - set optimbis=" " -else if $IBM then - set optim=" -O3 -qtune=pwr2 -qarch=pwr2" - set optimbis=" " -else if $VPP then -# set optim="-Dasuxm -On, -g -Ad -Potilax -Eciplume -Si" -# set optimbis=" -Wv,-m3 -Wp,-DVPP -Z $LMDGCM/listage" - set optimbis=" -Wp,-DNC_DOUBLE -Ad -Z $LMDGCM/listage" - set optim90="$optim $optimbis -X9 -w" - set mod_loc_dir=$MODIPSLDIR - set mod_suffix=mod -else if $DEC then - set optim=" " - set optimbis=" " -else if $LINUX then - set optim="-fast " - set optim90=" -fast " - set optimtru90=" -fast -c -Mfree -module $MODIPSLDIR " - set opt_link=" -Mfree -L /usr/local/pgi/linux86/lib -lpgf90 -lpgftnrtl -lpghpf -lpghpf2 -L$MODIPSLDIR -L$NCDFLIB -lnetcdf -lioipsl -Wl,-Bstatic -L/usr/lib/gcc-lib/i386-linux/2.95.2/" - set mod_loc_dir=$MODIPSLDIR - set mod_suffix=mod -else if $NEC then - set optim90=' -clear -C hopt -float0 -ew -P static -Wf,"-pvctl fullmsg noassume "' - set optimtru90=' -clear -f4 -C hopt -float0 -ew -P static -Wf,"-pvctl fullmsg noassume "' - set optim="$optim90" - set optimbis=" " - set opt_link=" -C hopt -float0 -ew -P static -L$MODIPSLDIR -lioipsl -L$NCDFLIB -lnetcdf_i8r8_v " - set mod_loc_dir="." - set mod_suffix="mod" -else if $XNEC then - set optim90=" -clear -float0 -ew -R5 " - set optimtru90=" -clear -f4 -float0 -ew -R2 -R3 -R4 -R5 " - set optim="$optim90" - set optimbis=" " - set mod_suffix="mod" -else - set optim="" - set optimbis=" " -endif - -set nomlib=$machine - -######################################################################## -# lecture des options de mymake -######################################################################## - -top: -if ($#argv > 0) then - switch ($1:q) - - case -h: - -######################################################################## -# Manuel en ligne -######################################################################## -more <> tmp -end -set suf= -foreach i ( `sort tmp | uniq ` ) - set suf=$suf$i -end -if ( ! $IBM ) then - set nomlib="$nomlib$suf" -endif -if ( $DEC ) then - set nomlib=DEC -endif -if ( $IBM ) then - set dim=`echo $dim | sed -en 's/[^0-9]/ /g'` - set dim_=`echo $dim | sed -en 's/[^0-9]/_/g'` -else if ( $SUN ) then - set dim=`echo $dim | sed -e 's/[^0-9]/ /g'` - set dim_=`echo $dim | sed -e 's/[^0-9]/_/g'` -else - set dim=`echo $dim | sed -n -e 's/[^0-9]/ /gp'` - set dim_=`echo $dim | sed -n -e 's/[^0-9]/_/gp'` -endif -set nomlib=${nomlib}_${dim_}_t${ntrac}_$grille -## M-A-F nomlib trop long sur CRAY pour ar -if ( $CRAY ) then - set nomlib=F90_${dim_}_t${ntrac} -endif -if ( $NEC || $XNEC ) then - set nomlib=F90_${dim_}_t${ntrac} -endif -echo calcul de la dimension -set dimc=`echo $dim | wc -w` - -if ( "$dimc" == "2" ) then -set include="$include "'-I$(LIBF)/dyn2d ' -set dimh=$dim -else -set include="$include "'-I$(LIBF)/dyn3d ' -set dimh=`echo $dim | awk ' { print $1 "." $2 } '` -endif -echo $dimc - -######################################################################## -# path pour les #include -######################################################################## - -set include="$include -I$NCDFINC -I$IOIPSLDIR" -echo $include - -######################################################################## -# Gestion des dimensions du modele. -# on cree ou remplace le fichier des dimensions/nombre de traceur -######################################################################## - -cd $libf/grid -if ( -f dimensions.h ) then -echo 'ATTENTION: vous etes sans doute en train de compiler le modele par ailleurs' -echo "Attendez que la premiere compilation soit terminee pour relancer la suivante." -echo "Si vous etes sur que vous ne compilez pas le modele par ailleurs," -echo vous pouvez continuer en repondant oui. -echo "Voulez-vous vraiment continuer?" -if ( $< == "oui" ) then -\rm -f $libf/grid/dimensions.h -else -exit -endif -endif - -cd dimension -./makdim $ntrac $dim -cat $libf/grid/dimensions.h - -cd $LMDGCM -set libo=$libo/$nomlib -if ( ! -d $libo ) then - mkdir $libo - cd $model -endif - -######################################################################## -# Differentes dynamiques (3d, 2d, 1d) -######################################################################## - -set dimension=`echo $dim | wc -w` -echo dimension $dimension -if ( $dimension == 1 ) then -echo pas de dynamique -set dyn="L_DYN= DYN= L_FILTRE= DIRMAIN=phy$physique " -endif -endif -cd $model -if ( $dimension == 3 ) then -cd libf/grid -\rm fxyprim.h -cp -p fxy_${grille}.h fxyprim.h -endif - -###################################################################### -# Traitement special pour le nouveau rayonnement de Laurent Li. -###################################################################### - -if ( -f $libf/phy$physique/raddim.h ) then - if ( -f $libf/phy$physique/raddim.$dimh.h ) then - \rm -f $libf/phy$physique/raddim.h - cp -p $libf/phy$physique/raddim.$dimh.h $libf/phy$physique/raddim.h - echo $libf/phy$physique/raddim.$dimh.h - cat $libf/phy$physique/raddim.$dimh.h - cat $libf/phy$physique/raddim.h - else - echo On peut diminuer la taille de l executable en creant - echo le fichier $libf/phy$physique/raddim.$dimh.h - \cp -p $libf/phy$physique/raddim.defaut.h $libf/phy$physique/raddim.h - endif -endif - -###################################################################### -# Gestion du filtre qui n'existe qu'en 3d. -###################################################################### - -if ( `expr $dimc \> 2` == 1 ) then - set filtre="FILTRE=$filtre" -else - set filtre="FILTRE= L_FILTRE= " -endif -echo MACRO FILTRE $filtre - -echo $dimc - -######################################################################## -# utilisation des vraies routines de couplage si on est en couple -######################################################################## -if ( $couple == 'true' ) then - banner couple - pwd - if ( `diff $libf/phylmd/oasis.F $libf/phylmd/oasis.true | wc -w` ) then - \cp $libf/phylmd/oasis.F $libf/phylmd/oasis.dummy - \cp $libf/phylmd/oasis.true $libf/phylmd/oasis.F - endif -else - if ( `diff $libf/phylmd/oasis.F $libf/phylmd/oasis.dummy | wc -w` ) then - \cp $libf/phylmd/oasis.F $libf/phylmd/oasis.true - \cp $libf/phylmd/oasis.dummy $libf/phylmd/oasis.F - endif -endif -######################################################################## -# Avant de lancer le make, on recree le makefile si necessaire -######################################################################## -# c'est a dire dans 3 cas: -# 1. si la liste des fichiers .F et .h a ete modifiee depuis la -# derniere creation du makefile -# 2. si le fichier contenant cette liste "liste_des_sources" -# n'existe pas. -# 3. Si le makefile n'existe pas. -######################################################################## - -cd $model -find libf -name '*.[Fh]' -print >! tmp77 -find libf -name '*.[Fh]' -exec egrep -i " *use *ioipsl" {} \; -print >! tmp90 -find libf -name '*.[Fh]90' -print >> tmp90 - -if ( `diff tmp77 liste_des_sources_f77 | wc -w` \ - || `diff tmp90 liste_des_sources_f90 | wc -w` \ - || ! -f makefile \ - || ! -f liste_des_sources_f90 \ - || ! -f liste_des_sources_f77 ) then - echo les fichiers suivants ont ete crees ou detruits - echo ou les fichiers suivants sont passes ou ne sont plus en Fortran 90 - diff liste_des_sources_f77 tmp77 - diff liste_des_sources_f90 tmp90 - \cp tmp77 liste_des_sources_f77 - \cp tmp90 liste_des_sources_f90 - echo On recree le makefile - ./create_make_gcm >! tmp - \mv tmp makefile - echo Nouveau makefile cree. -endif - -######################################################################## -# Execution de la comande make -######################################################################## - -echo PHYSIQUE $phys -echo dynamique $dyn $dimension -echo OPTIM="$optim" $filtre LIBO=$libo $dyn PHYS=$phys DIM=$dimc PROG=$code -echo PATH pour les fichiers INCLUDE $include -echo OPLINK="$oplink" - -if $HP then - set f77='fort77 +OP' - set f90='jensaisrien' - set opt_link="$opt_link -lm" -else if $VPP then - set f77=frt - set f90=$f77 - if ($couple == true) then - set opt_link="-Wg,-c $MODIPSLDIR/liboasis2.4_mpi2.a /usr/lang/mpi2/lib64/libmpi.a /usr/lang/mpi2/lib64/libmp.a /usr/local/lib/lib64/libnetcdf_cc_31.a -L$MODIPSLDIR -lioipsl" - set oplink="-Wl,-t,-P,-dy " - else - set opt_link="-Wg,-c /usr/local/lib/lib64/libnetcdf_cc_31.a -L$MODIPSLDIR -lioipsl" - set oplink="-Wl,-t,-dy " - endif - if ($veget == true) then - set opt_link="$opt_link -lsechiba -lparameters -lstomate -lioipsl /usr/local/lib/lib64/libnetcdf_cc_31.a" - endif -else if $CRAY then - set f77=f90 - set f90=f90 -else if $LINUX then - set f77=pgf90 - set f90=pgf90 -else if $SUN then - set f77=f90 - set f90=f90 -else if $NEC then - set f77=f90 - set f90=f90 - set opt_link="-L$MODIPSLDIR" - if ($veget == true) then - set opt_link="$opt_link -lsechiba -lparameters -lstomate" - endif - if ($couple == true) then - set opt_link="$opt_link -lioipsl -loasis2.4_mpi2 -float0 -ew -P static $NCDFLIB " - else - set opt_link="$opt_link -L$MODIPSLDIR -lioipsl -float0 -ew -P static $NCDFLIB " - endif - set mod_loc_dir="./" -else if $XNEC then - set f77=sxmpif90 - set f90=sxmpif90 - if $MODIPSL then - set opt_link="-L$MODIPSLDIR" - if ($veget == true) then - set opt_link="$opt_link -lsxsechiba -lsxparameters -lsxstomate" - endif - if ($couple == true) then - set opt_link="$opt_link -lsxioipsl -loasis2.4_mpi2 -float0 -ew -P static $NCDFLIB " - else - set opt_link="$opt_link -lsxioipsl -float0 -ew -P static $NCDFLIB " - endif - else - set opt_link=" -C hopt -float0 -ew -P static -L$MODIPSLDIR -lsxioipsl $NCDFLIB -I$NCDFINC " - endif - set mod_loc_dir="./" -else - set f77=f77 - set f90=f90 -endif - -cd $model - -if $VPP then -set make="gmake RANLIB=ls" -else if $CRAY then -set make="make RANLIB=ls" -else if $NEC then -set make="make RANLIB=ls" -else if $LINUX then -set make="make -k RANLIB=ranlib" -else if $XNEC then -set make="gmake RANLIB=ls" -else -set make="make RANLIB=ranlib" -endif - - - -# -# les deux test suivants sont "temporaires" pour pallier des "faiblesses" du -# compilateur fortran Sun: f90: SC4.0 11 Sep 1995 FORTRAN 90 1.1 -# -#if ($code == 'create_limit' && $SUN) then -# set link=f77 -# set opt_link="-L$NCDFLIB -lnetcdf" -#endif - -#if ($code == 'create_etat0' && $SUN) then -# if ( ! -f $libo/libdyn3d.a ) then -# echo "Priere de compiler gcm en premier pour des raisons d'optimisation" -# \rm $libf/grid/dimensions.h -# exit -# endif -# set optim=" -dalign " -# set optim90=" -dalign -fixed " -# set opt_link="-L$IOIPSLDIR -lioipsl -L$NCDFLIB -lnetcdf" -# set link="$f90 $optim90" -# touch $LMDGCM/libf/dyn3d/startvar.F -# touch $LMDGCM/libf/dyn3d/etat0_netcdf.F -#endif - -# -# etat0_netcdf a besoin d'info de la physique -# A revoir -set include="$include"' -I$(LIBF)/phylmd' -# -# le programme principal create_limit a besoin de l'info du module -# startvar: on met donc libo dans les include pour Nec -set include="$include"' -I$(LIBO)' - - -################################################################# -# Execution de la comande make... ENFIN! -################################################################# - -if $VPP then - set optim90=" $optim90 -Am -M$libo" - set optimtru90="$optim90" - \cp $IOIPSLDIR/*.mod $libo -else if $SUN then - set optim90=" $optim90 -M$libo -M$MODIPSLDIR " - set optimtru90=" $optimtru90 -M$libo -M$MODIPSLDIR " - set optim="$optim90" - \cp $IOIPSLDIR/*.mod $libo -else if $NEC then - set optim90=" $optim90 -I$libo " -else if $XNEC then - set optim90=" $optim90 -I$libo " - set optimtru90=" $optimtru90 -I$libo " -else if $LINUX then - set optim90=" $optim90 -module $libo " - set optim="$optim90" - set mod_loc_dir=$libo - \cp $IOIPSLDIR/*.mod $libo -endif - -set link="$f90 $optim90" - -set ar=ar - -if $XNEC then - set link="sxld $opt_link" - set link="$f90 " -# set ar=sxar -endif - - -cd $localdir - -echo $make -f $LMDGCM/makefile \ -OPTION_DEP="$opt_dep" OPTION_LINK="$opt_link" \ -OPTIM90="$optim90" \ -OPTIMTRU90="$optimtru90" \ -OPTIM="$optim$optimbis" \ -INCLUDE="$include" \ -$filtre \ -LIBO=$libo \ -$dyn \ -$phys \ -DIM=$dimc \ -L_ADJNT="$adjnt" \ -LOCAL_DIR="$localdir" \ -F77="$f77" \ -F90="$f90" \ -OPLINK="$oplink" \ -LINK="$link" \ -GCM="$LMDGCM" \ -MOD_LOC_DIR=$mod_loc_dir \ -MOD_SUFFIX=$mod_suffix \ -AR=$ar \ -PROG=$code - -$make -f $LMDGCM/makefile \ -OPTION_DEP="$opt_dep" OPTION_LINK="$opt_link" \ -OPTIM90="$optim90" \ -OPTIMTRU90="$optimtru90" \ -OPTIM="$optim$optimbis" \ -INCLUDE="$include" \ -$filtre \ -LIBO=$libo \ -$dyn \ -$phys \ -DIM=$dimc \ -L_ADJNT="$adjnt" \ -LOCAL_DIR="$localdir" \ -F77="$f77" \ -F90="$f90" \ -OPLINK="$oplink" \ -LINK="$link" \ -GCM="$LMDGCM" \ -MOD_LOC_DIR=$mod_loc_dir \ -MOD_SUFFIX=$mod_suffix \ -AR=$ar \ -PROG=$code - -\rm -f $libf/grid/dimensions.h Index: DZ.3.3/trunk/offline.def =================================================================== --- /LMDZ.3.3/trunk/offline.def (revision 6163) +++ (revision ) @@ -1,9 +1,0 @@ -T -4 -T --2. -48.1 -1 -T -6 -2 Index: DZ.3.3/trunk/orchidee.def =================================================================== --- /LMDZ.3.3/trunk/orchidee.def (revision 6163) +++ (revision ) @@ -1,65 +1,0 @@ -# -# $Header$ -# - -DEBUG_INFO=n -# WEATHER_GENERATOR is not set -#FORCING_FILE=/home/polcher/ISLSCP/islscp_Euro_87.nc -HEIGHT_LEV1=2.0 -SPLIT_DT=12 -RESTART_FILEIN=default -RESTART_FILEOUT=driver_rest.nc -TIME_SKIP=0 -TIME_LENGTH=30d -# IMPLICIT is not set -NO_INTER=y -SPRED_PREC=1 -NETRAD_CONS=n -# -# IOIPSL -# - -# -# PARAMETERS -# - -# -# SECHIBA -# -# STOMATE_OK_CO2 is not set -# STOMATE_OK_STOMATE is not set -# STOMATE_OK_DGVM is not set -# STOMATE_WATCHOUT is not set -SECHIBA_restart_in=start_sech.nc -SECHIBA_rest_out=sechiba_rest.nc -SECHIBA_reset_time=y -# SECHIBA_reset_time is not set -OUTPUT_FILE=sechiba_out.nc -WRITE_STEP=86400 -SECHIBA_HISTLEVEL=10 -STOMATE_OUTPUT_FILE=stomate_history.nc -STOMATE_HIST_DT=10. -STOMATE_HISTLEVEL=0 -SECHIBA_DAY=0.0 -SECHIBA_ZCANOP=0.5 -DT_SLOW=86400. -# IMPOSE_VEG is not set -VEGETATION_FILE=carteveg5km.nc -DIFFUCO_LEAFCI=233. -CONDVEG_SNOWA=default -# IMPOSE_AZE is not set -SOILALB_FILE=soils_param.nc -SOILTYPE_FILE=soils_param.nc -ENERBIL_TSURF=280. -HYDROL_SNOW=0.0 -HYDROL_SNOWAGE=0.0 -HYDROL_SNOWICE=0.0 -HYDROL_SNOWICEAGE=0.0 -HYDROL_HDRY=1.0 -HYDROL_HUMR=1.0 -HYDROL_BQSB=default -HYDROL_GQSB=0.0 -HYDROL_DSG=0.0 -HYDROL_DSP=default -HYDROL_QSV=0.0 -THERMOSOIL_TPRO=280. Index: DZ.3.3/trunk/physiq.def =================================================================== --- /LMDZ.3.3/trunk/physiq.def (revision 6163) +++ (revision ) @@ -1,13 +1,0 @@ -# -## $Header$ -# -# -# Automatically generated make config: don't edit -# - -OCEAN=force -VEGET=y -OK_journe=y -OK_mensuel=y -OK_instan=y -SCATTER=n Index: DZ.3.3/trunk/proc =================================================================== --- /LMDZ.3.3/trunk/proc (revision 6163) +++ (revision ) @@ -1,16 +1,0 @@ -#! /bin/sh -touch difdyn -for i in /users/fairhead/modele/source/LMDZ.3.2/libf/dyn3d/*.F - do - ls $i - echo "======== $i =================">>difdyn - sdiff -s $i libf/dyn3d/. >>difdyn - echo "=============================">>difdyn -done -for i in /users/fairhead/modele/source/LMDZ.3.2/libf/dyn3d/*.h - do - ls $i - echo "======== $i =================">>difdyn - sdiff -s $i libf/dyn3d/. >>difdyn - echo "=============================">>difdyn -done Index: DZ.3.3/trunk/run.def =================================================================== --- /LMDZ.3.3/trunk/run.def (revision 6163) +++ (revision ) @@ -1,161 +1,0 @@ - ------------------------------------------------------------------------ -Parametres de controle du run: ------------------------------- - -- Jour de l'etat initial ( = 350 si 20 Decembre ,par expl. ,comme ici ) - dayref - 110 - -- Annee de l'etat initial ( avec 4 chiffres ) - anneeref - 1988 - -- Nombre de jours d'integration - nday - 1 - -- nombre de pas par jour (multiple de iperiod) ( ici , dt= 10 secondes ) - day_step - 8640 - -- periode pour le pas Matsuno (en pas) - iperiod - 5 - -- periode de sortie des variables de controle (en pas) - iconser - 20 - -- periode d'ecriture du fichier histoire (en jour) - iecri - 100 - -- periode de stockage fichier histmoy (en jour) - periodav - 30 - -- periode de la dissipation (en pas) - idissip - 5 - -- choix de l'operateur de dissipation (star ou non star ) - lstardis - T - -- nombre d'iterations de l'operateur de dissipation gradiv - nitergdiv - 1 - -- nombre d'iterations de l'operateur de dissipation nxgradrot - nitergrot - 2 - -- nombre d'iterations de l'operateur de dissipation divgrad - niterh - 2 - -- temps de dissipation des plus petites long.d ondes pour u,v (gradiv) - tetagdiv - 1700. - -- temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot) - tetagrot - 1700. - -- temps de dissipation des plus petites long.d ondes pour h ( divgrad) - tetatemp - 1700. - -- coefficient pour gamdissip - coefdis - 0. - -- choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog) - purmats - F - -- avec ou sans physique - physic - T - -- periode de la physique (en pas) ( toutes les 30 minutes ) - iphysiq - 180 - -- frequence (en jours ) de l'ecriture du fichier histphy - ecritphy - 1. - -- Cycle diurne ou non - cycle_diurne - T - -- Soil Model ou non - soil_model - F - -- Choix ou non de New oliq - new_oliq - T - -- Orodr ou non pour l orographie - ok_orodr - F - -- Orolf ou non pour l orographie - ok_orolf - F - -- Si = .T. , lecture du fichier limit avec la bonne annee - ok_limitvrai - T - -- Nombre d'appels des routines de rayonnements ( par jour) - nbapp_rad - 12 - -- Flag pour la convection ( 1 pour LMD, 2 pour Tiedtke, 3 CCM ) - iflag_con - 2 - -- longitude en degres du centre du zoom - clon - 75. - - latitude en degres du centre du zoom - clat - 15. - -- facteur de grossissement du zoom,selon longitude - grossismx - 1.0 - -- facteur de grossissement du zoom ,selon latitude - grossismy - 1.0 - -- Fonction f(y) hyperbolique ( si .true. ) ou sinusoidale - fxyhypb - T - -- extension en longitude de la zone du zoom ( en degres ) ou fraction( <1 ) - dzoomx - 0.0 - -- extension en latitude de la zone du zoom ( en degres ) ou fraction( <1 ) - dzoomy - 0.0 - -- raideur du zoom en X - taux - 3. - -- raideur du zoom en Y - tauy - 3. - -- Fonction f(y) avec y = Sin(latit.) si = .TRUE. , Sinon y = latit. - ysinus - T - Index: DZ.3.3/trunk/run.def.10jours =================================================================== --- /LMDZ.3.3/trunk/run.def.10jours (revision 6163) +++ (revision ) @@ -1,187 +1,0 @@ - ------------------------------------------------------------------------ -Parametres de controle du run: ------------------------------- - -- Jour de l'etat initial ( = 350 si 20 Decembre ,par expl. ,comme ici ) - dayref - 181 - -- Annee de l'etat initial ( avec 4 chiffres ) - anneeref - 1998 - -- Nombre de jours d'integration - nday - 10 - -- nombre de pas par jour (multiple de iperiod) ( ici pour dt = 1 min ) - day_step - 240 - -- periode pour le pas Matsuno (en pas) - iperiod - 5 - -- periode de sortie des variables de controle (en pas) - iconser - 240 - -- periode d'ecriture du fichier histoire (en jour) - iecri - 1 - -- periode de stockage fichier histmoy (en jour) - periodav - 1. - -- periode de la dissipation (en pas) - idissip - 10 - -- choix de l'operateur de dissipation (star ou non star ) - lstardis - T - -- nombre d'iterations de l'operateur de dissipation gradiv -nitergdiv - 1 - -- nombre d'iterations de l'operateur de dissipation nxgradrot -nitergrot - 2 - -- nombre d'iterations de l'operateur de dissipation divgrad - niterh - 2 - -- temps de dissipation des plus petites long.d ondes pour u,v (gradiv) - tetagdiv - 7200. - -- temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot) - tetagrot - 7200. - -- temps de dissipation des plus petites long.d ondes pour h ( divgrad) - tetatemp - 7200. - -- coefficient pour gamdissip - coefdis - 0. - -- choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog) - purmats - F - -- avec ou sans physique - physic - T - -- periode de la physique (en pas) - iphysiq - 5 - -- frequence (en jours ) de l'ecriture du fichier histphy - ecritphy - 1 - -- Cycle diurne ou non - cycle_diurne - T - -- Soil Model ou non - soil_model - T - -- Choix ou non de New oliq - new_oliq - T - -- Orodr ou non pour l orographie - ok_orodr - T - -- Orolf ou non pour l orographie - ok_orolf - T - -- Si = .T. , lecture du fichier limit avec la bonne annee - ok_limitvrai - F - -- Nombre d'appels des routines de rayonnements ( par jour) - nbapp_rad - 12 - -- Flag pour la convection ( 1 pour LMD, 2 pour Tiedtke, 3 CCM ) - iflag_con - 2 - -- longitude en degres du centre du zoom - clon - 0. - -- latitude en degres du centre du zoom - clat - 0. - -- facteur de grossissement du zoom,selon longitude - grossismx - 1.0 - -- facteur de grossissement du zoom ,selon latitude - grossismy - 1.0 - -- Fonction f(y) hyperbolique si = .true. , sinon sinusoidale - fxyhypb - T - -- extension en longitude de la zone du zoom ( fraction de la zone totale) - dzoomx - 0.0 - -- extension en latitude de la zone du zoom ( fraction de la zone totale) - dzoomy - 0.0 - -- Fonction f(y) avec y = Sin(latit.) si = .true. , sinon y = latit. - ysinus - T - -c -cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc -cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc -cc -c Autres parametres ou cles , non places dans run.def mais -c dans les routines . -cc -ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc -cc -cc A ) iadv(iq ) , iq =1, nqmx , avec iadv = 1 , shema humid. -cc specifique LMD -cc -cc avec iadv = 2 , shema amont -cc -cc ( avec a l'interieur de amont , la cle amont_qs ) -cc -cc avec iadv = 3 , shema Van_Leer -cc -cc dans *** gcm.F **** -cc -cc B ) ok_oasis , pour coupleur Oasis , dans *** physiq.F *** -cc -cc C ) ok_ocean dans *** physiq.F *** -cc -cc D ) ok_journe , ok_instan, ok_region pour sorties analyses -c -cc dans *** physiq.F *** -cc -cc E ) nvm : nbre de vegetations -cc dans *** physiq.F *** -CC -cc F ) deltay : deplacement ( en degres ) en y de la zone du zoom -cc dans *** inigeom.F *** - Index: DZ4/trunk/000-README =================================================================== --- /LMDZ4/trunk/000-README (revision 6163) +++ (revision ) @@ -1,20 +1,0 @@ - -Logiciel LMDZ -------------- - -La documentation relative à LMDZ est accessible sur : -http://lmdz.lmd.jussieu.fr/documentation - -Les quides d'installation et utilisation de LMDZ sont accessibles sur : -http://lmdz.lmd.jussieu.fr/documentation/guides - -========================================================================== - -LMDZ software -------------- - -Documentation about the LMDZ software is available on the web at this address: -http://lmdz.lmd.jussieu.fr/documentation - -Practical installation and user guides are available here: -http://lmdz.lmd.jussieu.fr/documentation/guides Index: DZ4/trunk/bld.cfg =================================================================== --- /LMDZ4/trunk/bld.cfg (revision 6163) +++ (revision ) @@ -1,105 +1,0 @@ -# ----------------------- FCM extract configuration file ----------------------- -cfg::type bld -cfg::version 1.0 - - -# ------------------------------------------------------------------------------ -# Build information -# ------------------------------------------------------------------------------ - -#Default value of FPP fortran preprocessor -%FPP cpp - -inc arch.fcm -inc config.fcm - -%CONFIG_NAME %{ARCH}%SUFF_NAME -%BASE_CONFIG_PATH %LIBO/%CONFIG_NAME -%CONFIG_PATH %BASE_CONFIG_PATH/.config -%SRC_PATH %LIBF - -%FFLAGS %BASE_FFLAGS %COMPIL_FFLAGS %PARA_FFLAGS -%LD_FLAGS %BASE_LD %PARA_LD - -src::dyn %SRC_PATH/%DYN -src::phys %SRC_PATH/%PHYS -src::grid %SRC_PATH/grid -src::filtrez %SRC_PATH/filtrez -src::bibio %SRC_PATH/bibio -src::cosp %COSP - -bld::lib::dyn %DYN -bld::lib::phys %PHYS -bld::lib::grid grid -bld::lib::filtrez filtrez -bld::lib::bibio bibio -bld::lib::cosp cosp - - -bld::outfile_ext::exe %SUFF_NAME.e -bld::target lib%{DYN}.a lib%{PHYS}.a libgrid.a libfiltrez.a libbibio.a libcosp.a -bld::target %EXEC%SUFF_NAME.e -bld::exe_dep %{DYN} %{PHYS} grid filtrez bibio cosp - - -dir::root %CONFIG_PATH -dir::lib %BASE_CONFIG_PATH -dir::bin %ROOT_PATH/bin - -#search_src 1 - -bld::tool::fpp %FPP -bld::tool::fc %COMPILER -bld::tool::ld %LINK -bld::tool::ar %AR -bld::tool::make %MAKE -bld::tool::fflags %FFLAGS %INCDIR -bld::tool::ldflags %LD_FLAGS %LIB - -bld::tool::cppflags %FPP_FLAGS %INCDIR -bld::tool::fppflags %FPP_FLAGS %INCDIR -bld::tool::fppkeys %CPP_KEY %FPP_DEF - - -#bld::tool::fflags::phys::readaerosol %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt -pi auto -#bld::tool::fflags::phys::aeropt_2bands %BASE_FFLAGS %PROD_FFLAGS %INCDIR -#bld::tool::fflags::phys::radiation_AR4 %BASE_FFLAGS %PROD_FFLAGS1 %INCDIR -C hopt -Wf,-O,extendreorder -#bld::tool::fflags::phys::radiation_AR4_param %BASE_FFLAGS %PROD_FFLAGS1 %INCDIR -C hopt -f3 -#bld::tool::fflags::phys::fisrtilp %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt -#bld::tool::fflags::phys::cv30_routines %BASE_FFLAGS %PROD_FFLAGS %INCDIR -Wf,-O,extendreorder -#bld::tool::fflags::phys::cvltr %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt -#bld::tool::fflags::phys::clouds_gno %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt -#bld::tool::fflags::dyn::vlsplt_p %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt -#bld::tool::fflags::dyn::groupeun_p %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt - - -inc arch.opt - -# Pre-process code before analysing dependencies -bld::pp 1 - - -# Ignore the following dependencies -bld::excl_dep inc::netcdf.inc -bld::excl_dep use::netcdf -bld::excl_dep use::typesizes -bld::excl_dep h::netcdf.inc -bld::excl_dep h::mpif.h -bld::excl_dep inc::mpif.h -bld::excl_dep use::ioipsl -bld::excl_dep use::intersurf -bld::excl_dep use::mod_prism_proto -bld::excl_dep use::mod_prism_def_partition_proto -bld::excl_dep use::mod_prism_get_proto -bld::excl_dep use::mod_prism_put_proto -bld::excl_dep use::mkl_dfti - -# Don't generate interface files -bld::tool::geninterface none - -# Allow ".inc" as an extension for CPP include files -bld::infile_ext::inc CPP::INCLUDE - -# extension for module output -bld::outfile_ext::mod .mod -bld::tool::SHELL /bin/bash Index: DZ4/trunk/build_gcm =================================================================== --- /LMDZ4/trunk/build_gcm (revision 6163) +++ (revision ) @@ -1,24 +1,0 @@ -#!/bin/csh - -if ( -f '.lock' ) then - echo 'ATTENTION: vous etes sans doute en train de compiler le modele par ailleurs' - echo "Attendez que la premiere compilation soit terminee pour relancer la suivante." - echo "Si vous etes sur que vous ne compilez pas le modele par ailleurs," - echo vous pouvez continuer en repondant oui. - echo "Voulez-vous vraiment continuer?" - - if ( $< == "oui" ) then - - else - exit - endif -else - echo "compilation en cours..." > '.lock' -endif - -#set arch=$1 - - -fcm build - -\rm -f '.lock' Index: DZ4/trunk/cosp_input_nl.txt =================================================================== --- /LMDZ4/trunk/cosp_input_nl.txt (revision 6163) +++ (revision ) @@ -1,104 +1,0 @@ -! (c) British Crown Copyright 2008, the Met Office. -! All rights reserved. -! -! Redistribution and use in source and binary forms, with or without modification, are permitted -! provided that the following conditions are met: -! -! * Redistributions of source code must retain the above copyright notice, this list -! of conditions and the following disclaimer. -! * Redistributions in binary form must reproduce the above copyright notice, this list -! of conditions and the following disclaimer in the documentation and/or other materials -! provided with the distribution. -! * Neither the name of the Met Office nor the names of its contributors may be used -! to endorse or promote products derived from this software without specific prior written -! permission. -! -! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR -! IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND -! FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR -! CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL -! DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, -! DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER -! IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT -! OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. - - -! Namelist that sets up the main COSP options -&COSP_INPUT - CMOR_NL='./cmor/cosp_cmor_nl.txt', ! CMOR namelist - NPOINTS=9026,! Number of gridpoints (klon dans LMDZi : ici correspond a klon de 96x95) - NPOINTS_IT=10000,! Max number of gridpoints to be processed in one iteration - NCOLUMNS=20, ! Number of subcolumns - NLEVELS=39, ! Number of model levels - USE_VGRID=.true., ! Use fixed vertical grid for outputs? (if .true. then you need to define number of levels with Nlr) - NLR=40, ! Number of levels in statistical outputs (only used if USE_VGRID=.true.) - CSAT_VGRID=.true., ! CloudSat vertical grid? (if .true. then the CloudSat standard grid is used for the outputs. - ! USE_VGRID needs also be .true.) - FINPUT='histday.nc', ! NetCDF file with 1D inputs -! FINPUT='cosp_input_um_2d.nc', ! NetCDF file with 2D inputs - !---------------------------------------------------------------------------------- - !--------------- Inputs related to radar simulations - !---------------------------------------------------------------------------------- - RADAR_FREQ=94.0, ! CloudSat radar frequency (GHz) - SURFACE_RADAR=0, ! surface=1, spaceborne=0 - use_mie_tables=0,! use a precomputed lookup table? yes=1,no=0 - use_gas_abs=1, ! include gaseous absorption? yes=1,no=0 - do_ray=0, ! calculate/output Rayleigh refl=1, not=0 - melt_lay=0, ! melting layer model off=0, on=1 - k2=-1, ! |K|^2, -1=use frequency dependent default - use_reff=.true., ! True if you want effective radius to be used by radar simulator (always used by lidar) - use_precipitation_fluxes=.true., ! True if precipitation fluxes are input to the algorithm - !---------------------------------------------------------------------------------- - !---------------- Inputs related to lidar simulations - !---------------------------------------------------------------------------------- - Nprmts_max_hydro=12, ! Max number of parameters for hydrometeor size distributions - Naero=1, ! Number of aerosol species (Not used) - Nprmts_max_aero=1, ! Max number of parameters for aerosol size distributions (Not used) - lidar_ice_type=0, ! Ice particle shape in lidar calculations (0=ice-spheres ; 1=ice-non-spherical) - OVERLAP=3, ! overlap type: 1=max, 2=rand, 3=max/rand - !---------------------------------------------------------------------------------- - !---------------- Inputs related to ISCCP simulator - !---------------------------------------------------------------------------------- - ISCCP_TOPHEIGHT=1, ! 1 = adjust top height using both a computed - ! infrared brightness temperature and the visible - ! optical depth to adjust cloud top pressure. Note - ! that this calculation is most appropriate to compare - ! to ISCCP data during sunlit hours. - ! 2 = do not adjust top height, that is cloud top - ! pressure is the actual cloud top pressure - ! in the model - ! 3 = adjust top height using only the computed - ! infrared brightness temperature. Note that this - ! calculation is most appropriate to compare to ISCCP - ! IR only algortihm (i.e. you can compare to nighttime - ! ISCCP data with this option) - ISCCP_TOPHEIGHT_DIRECTION=1, ! direction for finding atmosphere pressure level - ! with interpolated temperature equal to the radiance - ! determined cloud-top temperature - ! 1 = find the *lowest* altitude (highest pressure) level - ! with interpolated temperature equal to the radiance - ! determined cloud-top temperature - ! 2 = find the *highest* altitude (lowest pressure) level - ! with interpolated temperature equal to the radiance - ! determined cloud-top temperature - ! ONLY APPLICABLE IF top_height EQUALS 1 or 3 - ! 1 = default setting, and matches all versions of - ! ISCCP simulator with versions numbers 3.5.1 and lower - ! 2 = experimental setting - !---------------------------------------------------------------------------------- - !-------------- RTTOV inputs - !---------------------------------------------------------------------------------- - Platform=1, ! satellite platform - Satellite=15, ! satellite - Instrument=0, ! instrument - Nchannels=8, ! Number of channels to be computed - Channels=1,3,5,6,8,10,11,13, ! Channel numbers (please be sure that you supply Nchannels) - Surfem=0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0, ! Surface emissivity (please be sure that you supply Nchannels) - ZenAng=50.0, ! Satellite Zenith Angle - CO2=5.241e-04, ! Mixing ratios of trace gases - CH4=9.139e-07, - N2O=4.665e-07, - CO=2.098e-07 -/ - - Index: DZ4/trunk/cosp_output_nl.txt =================================================================== --- /LMDZ4/trunk/cosp_output_nl.txt (revision 6163) +++ (revision ) @@ -1,63 +1,0 @@ -! (c) British Crown Copyright 2008, the Met Office. -! All rights reserved. -! -! Redistribution and use in source and binary forms, with or without modification, are permitted -! provided that the following conditions are met: -! -! * Redistributions of source code must retain the above copyright notice, this list -! of conditions and the following disclaimer. -! * Redistributions in binary form must reproduce the above copyright notice, this list -! of conditions and the following disclaimer in the documentation and/or other materials -! provided with the distribution. -! * Neither the name of the Met Office nor the names of its contributors may be used -! to endorse or promote products derived from this software without specific prior written -! permission. -! -! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR -! IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND -! FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR -! CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL -! DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, -! DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER -! IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT -! OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. - -! Namelist that sets up output-related variables. It controls -! the instrument simulators to run and the list of variables -! to be written to file -&COSP_OUTPUT - ! Simulator flags - Lradar_sim=.false., - Llidar_sim=.true., - Lisccp_sim=.true., - Lmisr_sim=.false., - ! Output variables - Lalbisccp=.true., - Latb532=.true., - Lboxptopisccp=.true., - Lboxtauisccp=.true., - Lcfad_dbze94=.false., - Lcfad_lidarsr532=.true., - Lclcalipso=.true., - Lclhcalipso=.true., - Lclisccp2=.true., - Lcllcalipso=.true., - Lclmcalipso=.true., - Lcltcalipso=.true., - Lctpisccp=.true., - Ldbze94=.false., - Ltauisccp=.true., - Ltclisccp=.true., - Llongitude=.false., - Llatitude=.false., - Lparasol_refl=.true., - LclMISR=.false., - Lmeantbisccp=.true., - Lmeantbclrisccp=.true., - ! Use lidar and radar - Lclcalipso2=.false., - Lcltlidarradar=.false., - ! These are provided for debugging or special purposes - Lfrac_out=.false., - Lbeta_mol532=.true., -/ Index: DZ4/trunk/create_make_gcm =================================================================== --- /LMDZ4/trunk/create_make_gcm (revision 6163) +++ (revision ) @@ -1,253 +1,0 @@ -#!/bin/sh -# -# $Header$ -# -#set -xv -machine=`hostname` -os=`uname` -gcm=`pwd` -libf=$gcm/libf -libo=$gcm/libo -CRAY=0 -if [ "$machine" = "atlas" -o "$machine" = "etoile" -o "$machine" = "axis" ] ; then - CRAY=1 -fi -XNEC=0 -if [ "$machine" = "rhodes" ] ; then - XNEC=1 -fi -X6NEC=0 -if [ "$machine" = "mercure" ] ; then - X6NEC=1 -fi -X8BRODIE=0 -if [ "$machine" = "brodie" ] ; then - X8BRODIE=1 -fi -VPP=0 -if [ "$machine" = "nymphea0" ] ; then - VPP=1 -fi -# -echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" -echo "# Definitions de Macros pour Make" -echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" -echo -echo "# Repertoires :" -echo -echo "GCM = "$gcm -echo 'LIBF = $(GCM)/libf' -if [ "$CRAY" = '0' ] ; then -# echo 'LIBO = $(GCM)/libo/$(MACHINE)' - echo 'LIBO = $(LIBOGCM)/$(MACHINE)' -else - echo 'LIBO = $(GCM)/libo' -fi -#echo 'LOCAL_DIR=$(GCM)' -#echo $localdir -echo "LOCAL_DIR=`echo $localdir`" -echo 'BIBIO = $(LIBF)/bibio' -echo "FILTRE = filtre" -echo "PHYS = " -echo "DYN = dyn " -echo 'LIBPHY = $(LIBO)/libphy$(PHYS).a' -echo 'DIRMAIN=dyn$(DIM)d$(FLAG_PARA)' -echo 'RM=rm' -echo -echo "OPLINK = " -echo -echo '# Les differentes librairies pour l"edition des liens:' -echo -if ( [ "$XNEC" = '1' ] || [ "$X6NEC" = '1' ] || [ "$X8BRODIE" = '1' ] ) ; then - echo 'dyn3d = $(LIBO)/libsxdyn3d.a $(LIBO)/libsx$(FILTRE).a' - echo 'dyn3dpar = $(LIBO)/libsxdyn3dpar.a $(LIBO)/libsx$(FILTRE).a' - echo 'dyn2d = $(LIBO)/libsxdyn2d.a' - echo 'dyn1d = $(LIBO)/libsxdyn1d.a' - echo 'L_DYN = -lsxdyn$(DIM)d$(FLAG_PARA)' - echo 'L_FILTRE = -lsx$(FILTRE)' - echo 'L_PHY = -lsxphy$(PHYS) ' - echo 'L_BIBIO = -lsxbibio' - echo 'L_ADJNT =' - echo 'L_COSP = -lsxcosp' -else - echo 'dyn3d = $(LIBO)/libdyn3d.a $(LIBO)/lib$(FILTRE).a' - echo 'dyn3dpar = $(LIBO)/libdyn3dpar.a $(LIBO)/lib$(FILTRE).a' - echo 'dyn2d = $(LIBO)/libdyn2d.a' - echo 'dyn1d = $(LIBO)/libdyn1d.a' - echo 'L_DYN = -ldyn$(DIM)d$(FLAG_PARA)' - echo 'L_FILTRE = -l$(FILTRE)' - echo 'L_PHY = -lphy$(PHYS) ' - echo 'L_BIBIO = -lbibio' - echo 'L_ADJNT =' - echo 'L_COSP = -lcosp' -fi - -echo -echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" -echo "# Option de compilation FORTRAN" -echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" -echo - echo 'COMPILE = $(F77) $(OPTIM) $(INCLUDE) -c' - echo 'COMPILE90 = $(F90) $(OPTIM90) $(INCLUDE) -c' - echo 'COMPTRU90 = $(F90) $(OPTIMTRU90) $(INCLUDE) -c' - echo "LINK = $LINK" - echo "AR = $AR" -echo -echo -echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" -echo "# Creation des differents executables" -echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" -echo -echo "# Executables:" -echo "# ------------" -echo -echo "PROG = code" -echo -#echo 'main : chimie $(DYN) bibio phys $(OPTION_DEP) ' -echo 'main : $(DYN) bibio phys $(OPTION_DEP) ' -echo ' cd $(LIBO) ; $(RANLIB) lib*.a ; cd $(GCM) ;\' -echo ' cd $(LOCAL_DIR); \' -echo ' $(COMPILE90) $(LIBF)/$(DIRMAIN)/$(SOURCE) -o $(PROG).o ; \' -echo ' $(LINK) $(PROG).o -L$(LIBO) $(L_DYN) $(L_ADJNT) $(L_COSP) $(L_FILTRE) $(L_PHY) $(L_DYN) $(L_BIBIO) $(L_DYN) $(OPLINK) $(OPTION_LINK) -o $(LOCAL_DIR)/$(PROG).e ; $(RM) $(PROG).o ' -echo -echo 'dyn : $(LIBO)/libdyn$(DIM)d$(FLAG_PARA).a $(FILTRE)$(DIM)d' -echo -echo 'phys : $(LIBPHY)' -echo -#echo 'chimie : $(LIBO)/libchimie.a' -echo -echo 'bibio : $(LIBO)/libbibio.a' -echo -echo 'adjnt : $(LIBO)/libadjnt.a' -echo -echo 'cosp : $(LIBO)/libcosp.a' -echo -echo '$(FILTRE)3d : $(LIBO)/lib$(FILTRE).a' -echo -echo '$(FILTRE)2d :' -echo -echo '$(FILTRE)1d :' -echo -echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" -echo "# Contenu des differentes bibliotheques" -echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" -echo -echo -cd $libf >/dev/null 2>&1 -for diri in ` ls ` -do - if [ -d $diri ] ; then - if [ "`ls $diri/*.F`" != "" ] || [ "`ls $diri/*.F90`" != "" ] ; then - cd $diri >/dev/null 2>&1 - echo - listlib="" - for i in `ls *.F` - do - fili=`basename $i .F` - test=` ( head $i | grep ' PROGRAM' ) ` - if [ "$test" = "" ] ; then - listlib=$listlib" "$fili - fi - done - for i in `ls *.F90` - do - fili=`basename $i .F90` - test=` ( head $i | grep ' PROGRAM' ) ` - if [ "$test" = "" ] ; then - listlib=$listlib" "$fili - fi - done -# - echo - echo - echo "#=======================================================================" - echo "# Contenu de la bibliotheque correspondant au Directory "$diri - echo "#=======================================================================" - echo - for fili in $listlib - do - echo '$(LIBO)/lib'$diri".a : " '$(LIBO)/lib'$diri".a("$fili".o)" - echo - done - echo '.PRECIOUS : $(LIBO)/lib'$diri'.a' - echo - echo - echo "# Compilation des membres de la bibliotheque lib"$diri".a" - echo - for fili in $listlib - do - if [ -f $fili.F90 ] ; then - trufile=$fili.F90 - else - trufile=$fili.F - fi - F90=0 ; egrep -i '^ *use ' $trufile > /dev/null 2>&1 && F90=1 - egrep -i '^ *module ' $trufile > /dev/null 2>&1 && F90=1 - egrep -i '#include*.inc ' $trufile > /dev/null 2>&1 && F90=1 - str1='$(LIBO)/lib'$diri'.a('$fili'.o) : $(LIBF)/'$diri/$trufile - [ "$fili" = "chem.subs" ] && str1=$str1' $(LIBF)/'$diri/chem.mods.F - for stri in ` ( sed -n "/\#include/s/\#include//p" $trufile | sed 's/\"//g' ; egrep -i '^ *use ' $trufile | sed -e 's/,/ /' | awk ' { print $2 } ' ) ` - do - - -# Differents cas de dependance correspondant a des include ou des -# use module. -# soit dans le repertoire local soit dans un autre. - - stri=`echo $stri | tr [A-Z] [a-z]` - if [ -f $stri ] ; then - echo $str1 \\ - str1='$(LIBF)/'$diri'/'$stri - else - if [ -f $stri.F ] || [ -f $stri.F90 ] ; then - echo $str1 \\ - str1='$(LIBO)/lib'$diri'.a('$stri'.o)' - else - for dirinc in dyn3d grid bibio filtrez - do - if [ -f ../$dirinc/$stri ] ; then - echo $str1 \\ - str1='$(LIBF)/'`cd .. ; ls */$stri | head -1` - fi - if [ -f ../$dirinc/$stri.F90 ] ; then - echo $str1 \\ - str1='$(LIBO)/lib'$dirinc'.a('$stri'.o)' - fi - done - fi - fi - done - echo $str1 - if [ "$F90" -eq '0' ] ; then - echo ' cd $(LOCAL_DIR); \' - echo ' $(COMPILE) $(LIBF)/'$diri'/'$trufile' ; \' - else - echo ' cd $(LOCAL_DIR); \' - if [ -f $fili.F90 ] ; then - echo ' $(COMPTRU90) $(LIBF)/'$diri'/'$trufile' ; \' - else - echo ' $(COMPILE90) $(LIBF)/'$diri'/'$trufile' ; \' - fi - MODU=0; egrep -i '^ *module ' $trufile> /dev/null 2>&1 && MODU=1 - if [ "$MODU" -eq '1' -a "$CRAY" != '1' ] ; then - if [ "$os" = 'UNIX_System_V' ] ; then - echo ' cp $(MOD_LOC_DIR)/*.$(MOD_SUFFIX) $(LIBO)/ ; \' - else - echo ' mv $(MOD_LOC_DIR)/'$fili'.$(MOD_SUFFIX) $(LIBO)/'$fili'.$(MOD_SUFFIX) ; \' -# echo ' if [ "$(MOD_LOC_DIR)" ne "$(LIBO)" ] ; then mv $(MOD_LOC_DIR)/'*'.$(MOD_SUFFIX) $(LIBO) ; fi ; \' - fi - fi - fi - if ( [ "$XNEC" -eq '1' ] || [ "$X6NEC" = '1' ] || [ "$X8BRODIE" = '1' ] ) ; then - echo ' sxar r $(LIBO)/libsx'$diri'.a '$fili'.o ; \' - fi - echo ' $(AR) r $(LIBO)/lib'$diri'.a '$fili'.o ; $(RM) '$fili'.o ; \' - echo ' cd $(GCM)' - echo - done -# - echo - cd $libf - fi - fi -done Index: DZ4/trunk/gcm.def =================================================================== --- /LMDZ4/trunk/gcm.def (revision 6163) +++ (revision ) @@ -1,62 +1,0 @@ -## $Id$ -## nombre de pas par jour (multiple de iperiod) ( ici pour dt = 1 min ) -day_step=480 -## periode pour le pas Matsuno (en pas) -iperiod=5 -## periode de la dissipation (en pas) -idissip=5 -## choix de l'operateur de dissipation (star ou non star ) -lstardis=y -## nombre d'iterations de l'operateur de dissipation gradiv -nitergdiv=1 -## nombre d'iterations de l'operateur de dissipation nxgradrot -nitergrot=2 -## nombre d'iterations de l'operateur de dissipation divgrad -niterh=2 -## temps de dissipation des plus petites long.d ondes pour u,v (gradiv) -tetagdiv=5400. -## temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot) -tetagrot=5400. -## temps de dissipation des plus petites long.d ondes pour h ( divgrad) -tetatemp=5400. -## coefficient pour gamdissip -coefdis=0. -## choix du shema d'integration temporelle (Matsuno:y ou Matsuno-leapfrog:n) -purmats=n -## avec ou sans physique -## 0: pas de physique (e.g. en mode Shallow Water) -## 1: avec physique (e.g. physique phylmd) -## 2: avec rappel newtonien dans la dynamique -iflag_phys=1 -## avec ou sans fichiers de demarrage (start.nc, startphy.nc) ? -## (sans fichiers de demarrage, initialisation des champs par iniacademic -## dans la dynamique) -read_start=y -## periode de la physique (en pas dynamiques, n'a de sens que si iflag_phys=1) -iphysiq=10 -## Avec ou sans strato -ok_strato=n -# Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche -iflag_top_bound=2 -# Coefficient pour la couche eponge (valeur derniere couche) -tau_top_bound=5.e-5 -## longitude en degres du centre du zoom -clon=0. -## latitude en degres du centre du zoom -clat=0. -## facteur de grossissement du zoom,selon longitude -grossismx=1.0 -## facteur de grossissement du zoom ,selon latitude -grossismy=1.0 -## Fonction f(y) hyperbolique si = .true. , sinon sinusoidale -fxyhypb=y -## extension en longitude de la zone du zoom ( fraction de la zone totale) -dzoomx=0.0 -## extension en latitude de la zone du zoom ( fraction de la zone totale) -dzoomy=0.0 -##raideur du zoom en X -taux=3. -##raideur du zoom en Y -tauy=3. -## Fonction f(y) avec y = Sin(latit.) si = .true. , sinon y = latit. -ysinus=y Index: DZ4/trunk/guide.def =================================================================== --- /LMDZ4/trunk/guide.def (revision 6163) +++ (revision ) @@ -1,37 +1,0 @@ -ok_guide=y -# guidage sur niveaux modèle (y) ou standard -guide_modele=y -# inversion de l'ordre des niveaux verticaux -ok_invertp=y -ncep=y - ###################################### - #### guidage de u ##### - guide_u=y - ###################################### - #### guidage de v ##### - guide_v=y - ###################################### - #### guidage de T ##### - guide_T=y - ###################################### - #### guidage de P ##### - guide_P=n - ###################################### - #### guidage de Q (hr=y:hum.rel,n:hum.spec) ##### - guide_Q=n - guide_hr=n - ###################################### - ## guidage dans la couche limite - guide_BL=n - ###################################### -ini_anal=n -tau_min_u=0.04166667 -tau_max_u=0.125 -tau_min_v=0.04166667 -tau_max_v=0.125 -tau_min_T=0.04166667 -tau_max_T=10. -tau_min_Q=0.2 -tau_max_Q=10. -# gamma limité -gamma4=n Index: DZ4/trunk/makegcm =================================================================== --- /LMDZ4/trunk/makegcm (revision 6163) +++ (revision ) @@ -1,1102 +1,0 @@ -#!/bin/csh -# -# $Id$ -# -#set verbose echo -######################################################################## -# options par defaut pour la commande make -######################################################################## -set dim="96x71x19" -set physique=lmd -set phys="PHYS=$physique" -set include='-I$(LIBF)/grid -I$(LIBF)/bibio -I$(LIBF)/filtrez -I. ' -set filtre=filtrez -set grille=reg -set couple=false -set veget=false -set chimie=false -set psmile=true -set parallel=false -set vampir=false -set OPT_STACK='-Wf,-init stack=nan' -set OPT_STACK=' ' -set OPTIMI='-C debug -eC ' -set OPTIMI=' -ftrace ' -set OPT_LINUX='-O3' -set OPT_LINUX="-i4 -r8 -O3" -set io=ioipsl -set cosp=false - -set FC_LINUX=g95 -set FC_LINUX=gfortran -#set FC_LINUX=pgf90 -if ( $FC_LINUX == g95 ) then - set OPT_LINUX="-O3" - set OPT_LINUX="-O3" -else if ( $FC_LINUX == gfortran ) then - set OPT_LINUX="-fdefault-real-8 -O3" -# set OPT_LINUX="-O3 -fno-second-underscore" - set OPT_LINUX="-O3 " -else - # pgf90 options - set OPT_LINUX="-i4 -r8 -O2 -Munroll -Mnoframe -Mautoinline -Mcache_align" -endif - -######################################################################## -# path a changer contenant les sources et les objets du modele -######################################################################## - -###### VERSION LMDZ.4 -set INCALIB=../INCA3/config/lib -set LMDGCM="`pwd`" -setenv LIBOGCM $LMDGCM/libo -# -# -setenv IOIPSLDIR /d4/fairhead/gfortran/ioipsl_v2_1_9 -setenv MODIPSLDIR /d4/fairhead/gfortran/ioipsl_v2_1_9 -setenv NCDFINC /d4/fairhead/gfortran_4.4/netcdf-4.1.1/include -setenv NCDFLIB /d4/fairhead/gfortran_4.4/netcdf-4.1.1/lib - - - - - -setenv localdir "`pwd`" -set MODIPSL=0 -echo $localdir | grep modipsl >& /dev/null -if ( ! $status ) then - set MODIPSL=1 - setenv LMDGCM $localdir - cd ../.. - setenv LIBOGCM "`pwd`/lib" - setenv IOIPSLDIR $LIBOGCM - setenv MODIPSLDIR $LIBOGCM - cd $localdir - if ( `hostname` == rhodes ) then - set NCDFINC=`grep sxnec ../../util/AA_make.gdef| grep NCDF_INC|sed -e "s/^.* =//"` - set NCDFLIB=`grep sxnec ../../util/AA_make.gdef| grep NCDF_LIB|sed -e 's/^.* =//'` - else - if ( `hostname` == nymphea0 ) then - set NCDFINC=`grep fjvpp ../../util/AA_make.gdef| grep NCDF_INC|sed -e "s/^.* =//"` - set NCDFLIB=`grep fjvpp ../../util/AA_make.gdef| grep NCDF_LIB|sed -e 's/^.* =//'` - else if ( `hostname` == mercure ) then - set NCDFINC=`grep sx6nec ../../util/AA_make.gdef| grep NCDF_INC|sed -e "s/^.* =//"` - set NCDFLIB=`grep sx6nec ../../util/AA_make.gdef| grep NCDF_LIB|sed -e 's/^.* =//'` - else if ( `hostname` == brodie ) then - set NCDFINC=`grep sx8brodie ../../util/AA_make.gdef| grep NCDF_INC|sed -e "s/^.* =//"` - set NCDFLIB=`grep sx8brodie ../../util/AA_make.gdef| grep NCDF_LIB|sed -e 's/^.* =//'` - else - echo 'Probleme de definition des variables NCDFINC et NCDFLIB' - endif - endif -else - if ( ! $?LMDGCM ) then - echo You must initialize the variable LMDGCM in your environnement - echo for instance: "setenv LMDGCM /usr/myself/supergcm" in .cshrc - exit - endif - if ( ! $?LIBOGCM ) then - set LIBOGCM=$LMDGCM/libo - endif - if ( ! $?IOIPSLDIR ) then - echo You must initialize the variable IOIPSLDIR in your environnement - echo for instance: "setenv IOIPSLDIR /usr/myself/ioipsl" in .cshrc - exit - else - setenv MODIPSLDIR $IOIPSLDIR - endif - if ( ! $?NCDFLIB ) then - echo You must initialize the variable NCDFLIB in your environnement - echo for instance: "setenv NCDFLIB /usr/myself/netcdf" in .cshrc - exit - endif - if ( ! $?NCDFINC ) then - echo You must initialize the variable NCDFINC in your environnement - echo for instance: "setenv NCDFINC /usr/myself/netcdf" in .cshrc - exit - endif -endif -set model=$LMDGCM -set libo=$LIBOGCM - -######################################################################## -# Les differentes platformes reconnues -######################################################################## - -set HP=0 -set IBM=0 -set SUN=0 -set VPP=0 -set CRAY=0 -set DEC=0 -set LINUX=0 -set NEC=0 -set XNEC=0 -set X6NEC=0 -set X8BRODIE=0 -if ( `uname` == HP-UX ) then - set machine=HP - set HP=1 -else if (`uname` == UNIX_System_V ) then - set machine=VPP - set VPP=1 -else if (`uname` == SunOS ) then - set machine=SUN - set SUN=1 -else if ( `uname` == AIX ) then - set machine=IBM - set IBM=1 -else if ( `uname` == OSF1 ) then - set machine=ALPHA - set DEC=1 -else if ( `uname` == Linux && `hostname` != mercure && `hostname` != brodie ) then - set machine=LINUX - set LINUX=1 -else if ( `hostname` == atlas || `hostname` == axis || `hostname` == etoile ) then - set machine=CRAY - set CRAY=1 -else if ( `uname` == SUPER-UX ) then - set machine=NEC - set NEC=1 -else if ( `hostname` == rhodes) then - set machine=XNEC - set XNEC=1 -else if ( `hostname` == mercure) then - set machine=X6NEC - set X6NEC=1 -else if ( `hostname` == brodie) then - set machine=X8BRODIE - set X8BRODIE=1 -else - echo Vous travaillez sur une machine non prevue par le reglement - exit -endif - -if ( ! -d $libo ) then - mkdir $libo -endif - - -if $VPP then -set netcdf=netcdf_v -else -set netcdf=netcdf -endif -######################################################################## -# Quelques initialisations de variables du shell. -######################################################################## - -set dyn= -set opt_link="" -set adjnt="" -set lcosp="" -set opt_dep="" -set libchimie="" - -set optim="" -set optimbis="" -set optim90="" -set oplink="" - -######################################################################## -# Optimisations par defaut suivant les machines -######################################################################## - -echo "Optimisations par defaut suivant les machines" -set libf=$model/libf -#setenv localdir "LOCAL_DIR=`pwd`" -#setenv localdir "`pwd`" -cd $model -############# -if $CRAY then -############# -# set optim="-Wf'-ei' -dp -Wf'-a static'" - set optimbis=" -DCRAY " - set optim90="-Wp'-P' -DCRAY -p$IOIPSLDIR "'-p$(LIBO) -eiv ' - set optim="$optim90" - if ( $io == "ioipsl" ) then - set oplink="-Wl'-DSTACK=128 -f indef' -L$IOIPSLDIR -lioipsl -L$NCDFLIB -lnetcdf " - else - set oplink="-Wl'-DSTACK=128 -f indef' -L$IOIPSLDIR -L$NCDFLIB -lnetcdf " - endif - set mod_loc_dir=" " - set mod_suffix=" " -################# -else if $SUN then -################# - set optim=" -fast " - set optimbis=" " - set optim90=" -fast -fixed " - set optimtru90=" -fast -free " - if ( $io == "ioipsl" ) then - set opt_link="-lf77compat -L$MODIPSLDIR -lsechiba -lparameters -lstomate -lioipsl -L$NCDFLIB -lnetcdf " - else - set opt_link="-lf77compat -L$MODIPSLDIR -lsechiba -lparameters -lstomate -L$NCDFLIB -lnetcdf " - endif - set mod_loc_dir=$localdir - set mod_suffix=mod -################# -else if $HP then -################# - set optim=" +U77 -O +E1 " - set optimbis=" " -################# -else if $IBM then -################# - set optim=" -O3 -qtune=pwr2 -qarch=pwr2" - set optimbis=" " -################# -else if $VPP then -################# -# set optim="-Dasuxm -On, -g -Ad -Potilax -Eciplume -Si" -# set optimbis=" -Wv,-m3 -Wp,-DVPP -Z $LMDGCM/listage" - set optimbis=" -Wp,-DNC_DOUBLE -Ad -Z $LMDGCM/listage -X9" - set optim90="$optim $optimbis -X9 -w" - set mod_loc_dir=$MODIPSLDIR - set mod_suffix=mod -################# -else if $DEC then -################# - set optim=" " - set optimbis=" " -################# -else if $LINUX then -################# - if ( $FC_LINUX == pgf90 || $FC_LINUX == g95 || $FC_LINUX == gfortran ) then - set optim=" $OPT_LINUX " - set optim90=" $OPT_LINUX " - set optimtru90=" $OPT_LINUX " - else - echo 'compilateur linux non reconnu' - exit - endif - set mod_loc_dir=$MODIPSLDIR - set mod_suffix=mod -################# -else if $NEC then -################# - set optim90=' -clear -C hopt -float0 -ew -P static -Wf,"-pvctl fullmsg noassume "' - set optimtru90=' -clear -f4 -C hopt -float0 -ew -P static -Wf,"-pvctl fullmsg noassume "' - set optim="$optim90" - set optimbis=" " - if ( $io == "ioipsl" ) then - set opt_link=" -C hopt -float0 -ew -P static -L$MODIPSLDIR -lioipsl $NCDFLIB -lnetcdf_i8r8_v " - else - set opt_link=" -C hopt -float0 -ew -P static -L$MODIPSLDIR $NCDFLIB -lnetcdf_i8r8_v " - endif - set mod_loc_dir="." - set mod_suffix="mod" -################# -else if $XNEC then -################# - set optdbl='-dw -Wf\"-A dbl4\"' - set optim90=' -clear -float0 -f3 -Ep -DNC_DOUBLE -dw -Wf\"-A dbl4\" -R5 -Wf,"-pvctl loopcnt=40000 fullmsg noassume "' - set optimtru90=' -clear -f4 -float0 -Ep -DNC_DOUBLE -dw -Wf\"-A dbl4\" -R2 -R3 -R4 -R5 -Wf,"-pvctl loopcnt=40000 fullmsg noassume"' - set optim="$optim90" - set optimbis=" " - set mod_suffix="mod" - set mod_loc_dir="./" -################# -else if $X6NEC then -################# - set optdbl='-dw -Wf\"-A dbl4\"' - set optim90=' -clear -float0 -size_t64 -P stack -Wf "-ptr byte -init stack=nan -init heap=nan" -Ep -DNC_DOUBLE -dw -Wf\"-A dbl4\" -R5 -Wf,"-pvctl loopcnt=40000 fullmsg noassume "' - set optimtru90=' -clear -f4 -float0 -size_t64 -P stack -Wf "-ptr byte -init stack=nan -init heap=nan" -Ep -DNC_DOUBLE -dw -Wf\"-A dbl4\" -R2 -R3 -R4 -R5 -Wf,"-pvctl loopcnt=40000 fullmsg noassume"' - set optim="$optim90" - set optimbis=" " - set mod_suffix="mod" - set mod_loc_dir="./" -################# -else if $X8BRODIE then -################## - set optdbl='-dw -Wf\"-A dbl4\"' -# set optim90='-P stack -Wf,-pvctl res=whole,-A dbl4,-init stack=nan,-init heap=nan,-ptr byte -EP -R5 -float0 -dw -Wf,"-pvctl loopcnt=999999 fullmsg noassume" -I/SX/usr/include' - set optim90='-C vopt -Wf,-pvctl res=whole,-A dbl4,-init stack=nan,-init heap=nan,-ptr byte -EP -DNC_DOUBLE -R5 -float0 -dw -Wf,"-pvctl loopcnt=999999 noassume" -I/SX/usr/include' -# set optim90='-C vsafe -P stack -Wf,-pvctl res=whole,-A dbl4,-ptr byte -EP -R5 -float0 -dw -Wf,"-pvctl loopcnt=999999 fullmsg noassume" -I/SX/usr/include' - set optimtru90="$optim90" - set optim90="$optim90" - set optim="$optim90" - set optimbis=" " - set mod_suffix="mod" - set mod_loc_dir="./" -else - set optim="" - set optimbis=" " -endif - -set nomlib=$machine - -######################################################################## -# lecture des options de mymake -######################################################################## - -top: -if ($#argv > 0) then - switch ($1:q) - - case -h: - -######################################################################## -# Manuel en ligne -######################################################################## -more <> tmp -end -set suf= -foreach i ( `sort tmp | uniq ` ) - set suf=$suf$i -end -if ( ! $IBM ) then - set nomlib="$nomlib$suf" -endif -if ( $DEC ) then - set nomlib=DEC -endif -if ( $IBM ) then - set dim=`echo $dim | sed -en 's/[^0-9]/ /g'` - set dim_=`echo $dim | sed -en 's/[^0-9]/_/g'` -else if ( $SUN ) then - set dim=`echo $dim | sed -e 's/[^0-9]/ /g'` - set dim_=`echo $dim | sed -e 's/[^0-9]/_/g'` -else - set dim_=`echo $dim | sed -e 's/[^0-9]/_/g'` - set dim=`echo $dim | sed -e 's/[^0-9]/ /g'` -endif -set nomlib=${nomlib}${physique}_${dim_}_$grille -## M-A-F nomlib trop long sur CRAY pour ar -if ( $CRAY ) then - set nomlib=F90_${dim_} -endif -if ( $NEC || $XNEC || $X6NEC || $X8BRODIE ) then - set nomlib=F90_${dim_}_'phy'${physique}${FLAG_PARA} -endif -echo calcul de la dimension -set dimc=`echo $dim | wc -w` - -if ( "$dimc" == "2" ) then -set include="$include "'-I$(LIBF)/dyn2d ' -set dimh=$dim -else -set include="$include "'-I$(LIBF)/dyn3d${FLAG_PARA} ' -set dimh=`echo $dim | awk ' { print $1 "." $2 } '` -endif -echo $dimc - -######################################################################## -# path pour les #include -######################################################################## - -if ( $XNEC ) then - set include="$include -I$NCDFINC -I$IOIPSLDIR" -else - set include="$include -I$NCDFINC -I$IOIPSLDIR" -endif -echo $include - -######################################################################## -# Gestion des dimensions du modele. -# on cree ou remplace le fichier des dimensions -######################################################################## - -cd $libf/grid -if ( -f dimensions.h ) then -echo 'ATTENTION: vous etes sans doute en train de compiler le modele par ailleurs' -echo "Attendez que la premiere compilation soit terminee pour relancer la suivante." -echo "Si vous etes sur que vous ne compilez pas le modele par ailleurs," -echo vous pouvez continuer en repondant oui. -echo "Voulez-vous vraiment continuer?" -if ( $< == "oui" ) then -\rm -f $libf/grid/dimensions.h -else -exit -endif -endif - -cd dimension -./makdim $dim -cat $libf/grid/dimensions.h - -cd $LMDGCM -set libo=$libo/$nomlib -if ( ! -d $libo ) then - mkdir $libo - cd $model -endif - -######################################################################## -# Differentes dynamiques (3d, 2d, 1d) -######################################################################## - -set dimension=`echo $dim | wc -w` -echo dimension $dimension -if ( $dimension == 1 ) then -echo pas de dynamique -set dyn="L_DYN= DYN= L_FILTRE= DIRMAIN=phy$physique " -endif -endif -cd $model -if ( $dimension == 3 ) then -cd libf/grid -\rm fxyprim.h -cp -p fxy_${grille}.h fxyprim.h -endif - -###################################################################### -# Traitement special pour le nouveau rayonnement de Laurent Li. -###################################################################### - -#if ( -f $libf/phy$physique/raddim.h ) then -# if ( -f $libf/phy$physique/raddim.$dimh.h ) then -# \rm -f $libf/phy$physique/raddim.h -# cp -p $libf/phy$physique/raddim.$dimh.h $libf/phy$physique/raddim.h -# echo $libf/phy$physique/raddim.$dimh.h -# cat $libf/phy$physique/raddim.$dimh.h -# cat $libf/phy$physique/raddim.h -# else -# echo On peut diminuer la taille de l executable en creant -# echo le fichier $libf/phy$physique/raddim.$dimh.h -# \cp -p $libf/phy$physique/raddim.defaut.h $libf/phy$physique/raddim.h -# endif -#endif - -###################################################################### -# Gestion du filtre qui n'existe qu'en 3d. -###################################################################### - -if ( `expr $dimc \> 2` == 1 ) then - set filtre="FILTRE=$filtre" -else - set filtre="FILTRE= L_FILTRE= " -endif -echo MACRO FILTRE $filtre - -echo $dimc - -######################################################################## -# Avant de lancer le make, on recree le makefile si necessaire -######################################################################## -######################################################################## -# c'est a dire dans 3 cas: -# 1. si la liste des fichiers .F et .h a ete modifiee depuis la -# derniere creation du makefile -# 2. si le fichier contenant cette liste "liste_des_sources" -# n'existe pas. -# 3. Si le makefile n'existe pas. -######################################################################## -########################################## -# On adapte d'abord certains include à F90 -########################################## -########################################## -cd $model -find libf -name '*.[Fh]' -print >! tmp77 -find libf -name '*.[Fh]' -exec egrep -i " *use *ioipsl" {} \; -print >! tmp90 -find libf -name '*.[Fh]90' -print >> tmp90 - -if ( `diff tmp77 liste_des_sources_f77 | wc -w` \ - || `diff tmp90 liste_des_sources_f90 | wc -w` \ - || ! -f makefile \ - || ! -f liste_des_sources_f90 \ - || ! -f liste_des_sources_f77 ) then - echo les fichiers suivants ont ete crees ou detruits - echo ou les fichiers suivants sont passes ou ne sont plus en Fortran 90 - diff liste_des_sources_f77 tmp77 - diff liste_des_sources_f90 tmp90 - \cp tmp77 liste_des_sources_f77 - \cp tmp90 liste_des_sources_f90 - echo On recree le makefile - ./create_make_gcm >! tmp - \mv tmp makefile - echo Nouveau makefile cree. -endif - -######################################################################## -# Execution de la comande make -######################################################################## - -echo PHYSIQUE $phys -echo dynamique $dyn $dimension -echo OPTIM="$optim" $filtre LIBO=$libo $dyn PHYS=$phys DIM=$dimc PROG=$code -echo PATH pour les fichiers INCLUDE $include -echo OPLINK="$oplink" - -################# -if $HP then -################# - set f77='fort77 +OP' - set f90='jensaisrien' - set opt_link="$opt_link -lm" -################# -else if $VPP then -################# - set f77=frt - set f90=$f77 - if ($couple == true) then - set opt_link="-Wg,-c $MODIPSLDIR/liboasis2.4_mpi2.a /usr/lang/mpi2/lib64/libmpi.a /usr/lang/mpi2/lib64/libmp.a -L$MODIPSLDIR -lioipsl /usr/local/lib/lib64/libnetcdf_cc_31.a" - set oplink="-Wl,-t,-P,-dy " - else - set opt_link="-Wg,-c -L$MODIPSLDIR -lioipsl /usr/local/lib/lib64/libnetcdf_cc_31.a" - set oplink="-Wl,-t,-dy " - endif - if ($veget == true) then - set opt_link="$opt_link $link_veget -lioipsl /usr/local/lib/lib64/libnetcdf_cc_31.a" - endif -################# -else if $CRAY then -################# - set f77=f90 - set f90=f90 -################# -else if $LINUX then -################# -# set f77=pgf90 -# set f90=pgf90 - set f77=$FC_LINUX - set f90=$FC_LINUX - if ( $FC_LINUX == 'pgf90' ) then - if ( $io == "ioipsl" ) then - set opt_link=" -L$MODIPSLDIR $link_veget -L$NCDFLIB -lioipsl -lnetcdf " - else - set opt_link=" -L$MODIPSLDIR $link_veget -L$NCDFLIB -lnetcdf " - endif - else if ($FC_LINUX == 'g95' || $FC_LINUX == 'gfortran' ) then - if ( $io == "ioipsl" ) then - set opt_link="-L$MODIPSLDIR $link_veget -lioipsl -L$NCDFLIB -lnetcdf -lioipsl -lnetcdf " - else - set opt_link="-L$MODIPSLDIR $link_veget -lioipsl -L$NCDFLIB -lnetcdf -lnetcdf " - endif - else - set opt_link=" " - endif -################# -else if $SUN then -################# - set f77=f90 - set f90=f90 - if ( $io == "ioipsl" ) then - set opt_link="-lf77compat -L$MODIPSLDIR $link_veget -lioipsl -L$NCDFLIB -lnetcdf " - else - set opt_link="-lf77compat -L$MODIPSLDIR $link_veget -L$NCDFLIB -lnetcdf " - endif -################# -else if $NEC then -################# - set f77=f90 -ftrace - set f90=f90 -ftrace - set opt_link="-L$MODIPSLDIR" - if ($veget == true) then - set opt_link="$opt_link $link_veget" - endif - if ($couple == true) then - set opt_link="$opt_link -lioipsl -loasis2.4_mpi2 -float0 -ew -P static $NCDFLIB " - else - set opt_link="$opt_link -L$MODIPSLDIR -lioipsl -float0 -ew -P static $NCDFLIB " - endif - set mod_loc_dir="./" -################# -else if $XNEC then -################# - set f77="sxmpif90 -ftrace" - set f90="sxmpif90 -ftrace" - if $MODIPSL then - set opt_link="-L$MODIPSLDIR" - if ($veget == true) then - set opt_link="$opt_link $link_veget" - endif - if ($couple == true) then - if ($psmile == true) then - set opt_link="$opt_link -lsxioipsl -float0 $optdbl -P static $NCDFLIB " - else - set opt_link="$opt_link -lsxioipsl -loasis2.4_mpi2 -float0 $optdbl -P static $NCDFLIB " - endif - else - set opt_link="$opt_link -lsxioipsl -float0 $optdbl -P static $NCDFLIB " - endif - else - if ($couple == true) then - set opt_link="-L$MODIPSLDIR" - set opt_link="$opt_link $link_veget -lsxioipsl -loasis2.4_mpi2 -float0 $optdbl -P static $NCDFLIB " - else - set opt_link=" -C hopt -float0 $optdbl -P static -L$MODIPSLDIR $link_veget -lsxioipsl $NCDFLIB " - endif - endif - set mod_loc_dir="./" -################## -else if $X6NEC then -################## - set f77=sxmpif90 - set f90=sxmpif90 - if $MODIPSL then - set opt_link="$opt_link -L$MODIPSLDIR" - if ($veget == true) then - set opt_link="$opt_link $link_veget" - endif - if ($couple == true) then - if ($psmile == true) then - set opt_link="$opt_link -lsxioipsl -float0 -size_t64 $optdbl -P static $NCDFLIB " - else - set opt_link="$opt_link -lsxioipsl -loasis2.4_mpi2 -float0 -size_t64 $optdbl -P static $NCDFLIB " - endif - else - set opt_link="$opt_link -lsxioipsl -float0 -size_t64 $optdbl -P static $NCDFLIB " - endif - else -# set opt_link=" -float0 -size_t64 $optdbl -P static -L$MODIPSLDIR -lsxsechiba -lsxparameters -lsxstomate -lsxioipsl $NCDFLIB " - set opt_link=" $opt_link -float0 -size_t64 $optdbl -P static -L$MODIPSLDIR -lsxioipsl $NCDFLIB " - - endif - set mod_loc_dir="./" -################## -else if $X8BRODIE then -################## - set f77=sxmpif90 - set f90=sxmpif90 - if $MODIPSL then - set opt_link="$opt_link -float0 -Wf,-A dbl4 -L$MODIPSLDIR -lblas" - if ($veget == true) then - set opt_link="$opt_link $link_veget" - endif - if ($couple == true) then - set opt_link="$opt_link -lsxioipsl -float0 $optdbl -P static $NCDFLIB " - else - set opt_link="$opt_link -lsxioipsl -float0 $optdbl -P static $NCDFLIB " - endif - else -# set opt_link=" -float0 $optdbl -P static -L$MODIPSLDIR -lsxsechiba -lsxparameters -lsxstomate -lsxioipsl $NCDFLIB " - set opt_link=" -float0 $optdbl -P static -L$MODIPSLDIR -lsxioipsl $NCDFLIB -lblas" - - endif - set mod_loc_dir="./" -################# -else -################# - set f77=f77 - set f90=f90 -endif - -cd $model - -if $VPP then -set make="gmake RANLIB=ls" -else if $CRAY then -set make="make RANLIB=ls" -else if $NEC then -set make="make RANLIB=ls" -else if $LINUX then -set make="make -k RANLIB=ranlib" -else if $XNEC then -set make="gmake RANLIB=ls" -else if $X6NEC then -set make="gmake RANLIB=ls" -else if $X8BRODIE then -set make="gmake RANLIB=ls" -else -set make="make RANLIB=ranlib" -endif - - - - -# -# etat0_netcdf a besoin d'info de la physique -# A revoir -set include="$include"' -I$(LIBF)/phy'"$physique" -# -# le programme principal create_limit a besoin de l'info du module -# startvar: on met donc libo dans les include pour Nec -set include="$include"' -I$(LIBO)' - - -################################################################# -# Execution de la comande make... ENFIN! -################################################################# - -if $VPP then - set optim90=" $optim90 -Am -M$libo" - set optimtru90="$optim90" - \cp $IOIPSLDIR/*.mod $libo -else if $SUN then - set optim90=" $optim90 -M$libo -M$MODIPSLDIR " - set optimtru90=" $optimtru90 -M$libo -M$MODIPSLDIR " - set optim="$optim90" - \cp $IOIPSLDIR/*.mod $libo -else if $NEC then - set optim90=" $optim90 -I$libo " -else if $XNEC then - set optim90=" $optim90 -I$libo " - set optimtru90=" $optimtru90 -I$libo " -else if $X6NEC then - set optim90=" $optim90 -I$libo " - set optimtru90=" $optimtru90 -I$libo " -else if $X8BRODIE then - set optim90=" $optim90 -I$libo " - set optimtru90=" $optimtru90 -I$libo " -else if $LINUX then - if ( $FC_LINUX == "pgf90" ) then - set optimtru90=" $optimtru90 -module $libo " - set optim90=" $optim90 -module $libo " - else if ( $FC_LINUX == 'g95' ) then - set optimtru90=" $optimtru90 -fmod=$libo " - set optim90=" $optim90 -fmod=$libo " - else if ( $FC_LINUX == 'gfortran' ) then - set optimtru90=" $optimtru90 -M $libo " - set optim90=" $optim90 -M $libo " - endif - set optim="$optim90" - set mod_loc_dir=$libo -# \cp /d3/fairhead/sechiba/ioipsl/*.mod $libo -# \cp $IOIPSLDIR/*.mod $libo -endif - -set link="$f90 $optim90" - -set ar=ar - -if $XNEC then - set link="sxld $opt_link" - set link="$f90 " -# set ar=sxar -else if $X6NEC then - set link="sxld $opt_link" - set link="$f90 -Wl,-hlib_cyclic " -else if $X8BRODIE then - set link="sxld $opt_link" - set link="$f90 -Wl,-hlib_cyclic " -endif - - -cd $localdir - -set source_code=${code}.F -if ( -f $LMDGCM/libf/dyn${dimc}d${FLAG_PARA}/${code}.F90 ) then - set source_code=${code}.F90 -endif - -echo $make -f $LMDGCM/makefile \ -OPTION_DEP="$opt_dep" OPTION_LINK="$opt_link" \ -OPTIM90="$optim90" \ -OPTIMTRU90="$optimtru90" \ -OPTIM="$optim$optimbis" \ -INCLUDE="$include" \ -$filtre \ -LIBO=$libo \ -$dyn \ -$phys \ -DIM=$dimc \ -FLAG_PARA="$FLAG_PARA"\ -L_ADJNT="$adjnt" \ -L_COSP="$lcosp" \ -L_CHIMIE="$libchimie" \ -LOCAL_DIR="$localdir" \ -F77="$f77" \ -F90="$f90" \ -OPLINK="$oplink" \ -LINK="$link" \ -GCM="$LMDGCM" \ -MOD_LOC_DIR=$mod_loc_dir \ -MOD_SUFFIX=$mod_suffix \ -AR=$ar \ -SOURCE=$source_code \ -PROG=$code - -$make -f $LMDGCM/makefile \ -OPTION_DEP="$opt_dep" OPTION_LINK="$opt_link" \ -OPTIM90="$optim90" \ -OPTIMTRU90="$optimtru90" \ -OPTIM="$optim$optimbis" \ -INCLUDE="$include" \ -$filtre \ -LIBO=$libo \ -$dyn \ -$phys \ -DIM=$dimc \ -FLAG_PARA="$FLAG_PARA"\ -L_ADJNT="$adjnt" \ -L_COSP="$lcosp" \ -L_CHIMIE="$libchimie" \ -LOCAL_DIR="$localdir" \ -F77="$f77" \ -F90="$f90" \ -OPLINK="$oplink" \ -LINK="$link" \ -GCM="$LMDGCM" \ -MOD_LOC_DIR=$mod_loc_dir \ -MOD_SUFFIX=$mod_suffix \ -AR=$ar \ -SOURCE=$source_code \ -PROG=$code - -\rm -f $libf/grid/dimensions.h Index: DZ4/trunk/makelmdz_fcm =================================================================== --- /LMDZ4/trunk/makelmdz_fcm (revision 6163) +++ (revision ) @@ -1,411 +1,0 @@ -#!/bin/bash -# $Id$ -# This is a script in Bash. - -# FH : on ne crée plus le fichier arch.mk qui est supposé exister par -# FH : ailleurs. -# FH : ulterieurement, ce fichier sera pré-existant pour une série -# FH : de configurations en versions optimisées et debug qui seront -# FH : liés (ln -s) avec arch.mk en fonction de l'architecture. -# FH : Pour le moment, cette version est en test et on peut créer les -# FH : arch.mk en lançant une première fois makegcm. -# -##set -x -######################################################################## -# options par defaut pour la commande make -######################################################################## - -dim="96x72x19" -physique=lmd -filtre=filtrez -grille=reg -couple=false -veget=false -chimie=false -parallel=none -compil_mod=prod -io=ioipsl -LIBPREFIX="" -fcm_path=none -cosp=false - -LMDGCM=`/bin/pwd` -LIBOGCM=$LMDGCM/libo -LIBFGCM=$LMDGCM/libf -COSP_PATH=$LMDGCM/.void_dir - -######################################################################## -# Quelques initialisations de variables du shell. -######################################################################## - -CPP_KEY="" -INCLUDE="" -LIB="" -adjnt="" -COMPIL_FFLAGS="%PROD_FFLAGS" -PARA_FFLAGS="" -PARA_LD="" -EXT_SRC="" - -######################################################################## -# lecture des options de mymake -######################################################################## - -while (($# > 0)) - do - case $1 in - "-h") cat < .void_file -rm -rf .void_dir -mkdir .void_dir -rm -f arch.path -ln -s arch/arch-${arch}.path ./arch.path -source arch.path - -######################################################################## -# Definition des clefs CPP, des chemins des includes et modules -# et des libraries -######################################################################## - -if [[ "$compil_mod" == "prod" ]] -then - COMPIL_FFLAGS="%PROD_FFLAGS" -elif [[ "$compil_mod" == "dev" ]] -then - COMPIL_FFLAGS="%DEV_FFLAGS" -elif [[ "$compil_mod" == "debug" ]] -then - COMPIL_FFLAGS="%DEBUG_FFLAGS" -fi - -if [[ "$physique" != "nophys" ]] -then - #Default planet type is Earth - CPP_KEY="$CPP_KEY CPP_EARTH" -fi - -if [[ "$chimie" == "INCA" ]] -then - CPP_KEY="$CPP_KEY INCA" - INCLUDE="$INCLUDE -I${INCA_INCDIR}" - LIB="$LIB -L${INCA_LIBDIR} -lchimie" -fi - -if [[ "$couple" != "false" ]] -then - CPP_KEY="$CPP_KEY CPP_COUPLE" - INCLUDE="$INCLUDE -I${OASIS_INCDIR}" - LIB="$LIB -L${OASIS_LIBDIR} -lpsmile.${couple} -lmpp_io" -fi - -if [[ "$parallel" == "mpi" ]] -then - CPP_KEY="$CPP_KEY CPP_PARA CPP_MPI" - PARA_FFLAGS="%MPI_FFLAGS" - PARA_LD="%MPI_LD" -elif [[ "$parallel" == "omp" ]] -then - CPP_KEY="$CPP_KEY CPP_PARA CPP_OMP" - PARA_FFLAGS="%OMP_FFLAGS" - PARA_LD="%OMP_LD" -elif [[ "$parallel" == "mpi_omp" ]] -then - CPP_KEY="$CPP_KEY CPP_PARA CPP_MPI CPP_OMP" - PARA_FFLAGS="%MPI_FFLAGS %OMP_FFLAGS" - PARA_LD="%MPI_LD %OMP_LD" -fi - -if [[ ( "$parallel" == "omp" || "$parallel" == "mpi_omp" ) \ - && "$compil_mod" == "debug" ]] -then - echo "Usually, parallelization with OpenMP requires some optimization." - echo "We suggest switching to \"-dev\"." -fi - -if [[ "$veget" == "true" ]] -then - CPP_KEY="$CPP_KEY CPP_VEGET" - INCLUDE="${INCLUDE} -I${ORCH_INCDIR}" - LIB="${LIB} -L${ORCH_LIBDIR} -l${LIBPREFIX}sechiba -l${LIBPREFIX}parameters -l${LIBPREFIX}stomate -l${LIBPREFIX}parallel -l${LIBPREFIX}orglob" -fi - -if [[ $io == ioipsl ]] -then - CPP_KEY="$CPP_KEY CPP_IOIPSL" - INCLUDE="$INCLUDE -I${IOIPSL_INCDIR}" - LIB="$LIB -L${IOIPSL_LIBDIR} -l${LIBPREFIX}ioipsl" -fi -if [[ "$cosp" == "true" ]] -then - CPP_KEY="$CPP_KEY CPP_COSP" - COSP_PATH="$LIBFGCM/cosp" -# LIB="${LIB} -l${LIBPREFIX}cosp" -fi -INCLUDE="$INCLUDE -I${NETCDF_INCDIR}" -LIB="$LIB -L${NETCDF_LIBDIR} -lnetcdf" - -######################################################################## -# calcul du nombre de dimensions -######################################################################## - - -dim_full=$dim -dim=`echo $dim | sed -e 's/[^0-9]/ /g'` -set $dim -dimc=$# -echo calcul de la dimension -echo dim $dim -echo dimc $dimc - - -######################################################################## -# Gestion des dimensions du modele. -# on cree ou remplace le fichier des dimensions -######################################################################## - -cd $LIBFGCM/grid/dimension -./makdim $dim -cat $LIBFGCM/grid/dimensions.h -cd $LMDGCM - - -######################################################################## -# Differentes dynamiques (3d, 2d, 1d) -######################################################################## - -dimension=`echo $dim | wc -w` -echo dimension $dimension - -if (( $dimension == 3 )) -then - cd $LIBFGCM/grid - \rm fxyprim.h - cp -p fxy_${grille}.h fxyprim.h -else - echo "Probleme dans les dimensions de la dynamique !!" - echo "Non reactive pour l'instant !!!" -fi - -###################################################################### -# Traitement special pour le nouveau rayonnement de Laurent Li. -# ---> YM desactive pour le traitemement en parallele -###################################################################### - -#if [[ -f $libf/phy$physique/raddim.h ]] -#then -# if [[ -f $libf/phy$physique/raddim.$dimh.h ]] -#then -# \rm -f $libf/phy$physique/raddim.h -# cp -p $libf/phy$physique/raddim.$dimh.h $libf/phy$physique/raddim.h -# echo $libf/phy$physique/raddim.$dimh.h -# cat $libf/phy$physique/raddim.h -# else -# echo On peut diminuer la taille de l executable en creant -# echo le fichier $libf/phy$physique/raddim.$dimh.h -# \cp -p $libf/phy$physique/raddim.defaut.h $libf/phy$physique/raddim.h -# fi -#fi - -###################################################################### -# Gestion du filtre qui n'existe qu'en 3d. -###################################################################### - -if (( `expr $dimc \> 2` == 1 )) -then - filtre="FILTRE=$filtre" -else - filtre="FILTRE= L_FILTRE= " -fi -echo MACRO FILTRE $filtre - -echo $dimc - - - -###################################################################### -# Creation du suffixe de la configuration -###################################################################### - - -SUFF_NAME=_${dim_full} -SUFF_NAME=${SUFF_NAME}_phy${physique} - -if [[ "$parallel" != "none" ]] -then - SUFF_NAME=${SUFF_NAME}_para - DYN=dyn${dimc}dpar -else - SUFF_NAME=${SUFF_NAME}_seq - DYN=dyn${dimc}d -fi - -if [[ $veget == "true" ]] -then - SUFF_NAME=${SUFF_NAME}_orch -fi - -if [[ $couple != "false" ]] -then - SUFF_NAME=${SUFF_NAME}_couple -fi - -if [[ $chimie == "INCA" ]] -then - SUFF_NAME=${SUFF_NAME}_inca -fi - -cd $LMDGCM -config_fcm="config.fcm" -rm -f $config_fcm -touch $config_fcm -rm -f bin/${code}${SUFF_NAME}.e -rm -f arch.fcm -rm -f arch.opt - -echo "%ARCH $arch" >> $config_fcm -echo "%INCDIR $INCLUDE" >> $config_fcm -echo "%LIB $LIB" >> $config_fcm -echo "%ROOT_PATH $PWD" >> $config_fcm -echo "%LIBF $LIBFGCM" >> $config_fcm -echo "%LIBO $LIBOGCM" >> $config_fcm -echo "%DYN $DYN" >> $config_fcm -echo "%PHYS phy${physique}" >> $config_fcm -echo "%COSP $COSP_PATH" >> $config_fcm -echo "%CPP_KEY $CPP_KEY" >> $config_fcm -echo "%EXEC $code" >> $config_fcm -echo "%SUFF_NAME $SUFF_NAME" >> $config_fcm -echo "%COMPIL_FFLAGS $COMPIL_FFLAGS" >> $config_fcm -echo "%PARA_FFLAGS $PARA_FFLAGS" >> $config_fcm -echo "%PARA_LD $PARA_LD" >> $config_fcm -echo "%EXT_SRC $EXT_SRC" >> $config_fcm - - - -ln -s arch/arch-${arch}.fcm arch.fcm -if test -f arch/arch-${arch}.opt && [ $compil_mod = "prod" ] - then - ln -s arch/arch-${arch}.opt arch.opt -else - ln -s .void_file arch.opt -fi - - -rm -f $LIBOGCM/${arch}${SUFF_NAME}/.config/fcm.bld.lock -./build_gcm - -rm -rf tmp_src -rm -rf config -ln -s $LIBOGCM/${arch}${SUFF_NAME}/.config config -ln -s $LIBOGCM/${arch}${SUFF_NAME}/.config/tmp tmp_src Index: DZ4/trunk/offline.def =================================================================== --- /LMDZ4/trunk/offline.def (revision 6163) +++ (revision ) @@ -1,12 +1,0 @@ -# -# $Header$ -# -T -4 -T --2. -48.1 -1 -T -6 -2 Index: DZ4/trunk/orchidee.def =================================================================== --- /LMDZ4/trunk/orchidee.def (revision 6163) +++ (revision ) @@ -1,63 +1,0 @@ -# -# -# Parameter file for LMDZ4OR_v2 configuration -# See comments : http://forge.ipsl.jussieu.fr/orchidee/ -# -STOMATE_OK_CO2=TRUE -# STOMATE_OK_STOMATE is not set -# STOMATE_OK_DGVM is not set -# STOMATE_WATCHOUT is not set -SECHIBA_restart_in=default -SECHIBA_rest_out=sechiba_rest.nc -SECHIBA_reset_time=y -# -OUTPUT_FILE=sechiba_out.nc -WRITE_STEP=2592000 -SECHIBA_HISTLEVEL=5 -# -SECHIBA_HISTFILE2 = FALSE -SECHIBA_OUTPUT_FILE2 = sechiba_out_2.nc -WRITE_STEP2 = 86400.0 -SECHIBA_HISTLEVEL2 = 1 -# -STOMATE_OUTPUT_FILE=stomate_history.nc -STOMATE_HIST_DT=10. -STOMATE_HISTLEVEL=0 -SECHIBA_DAY=0.0 -SECHIBA_ZCANOP=0.5 -DT_SLOW=86400. -# IMPOSE_VEG is not set -VEGETATION_FILE=carteveg5km.nc -# VEGETATION_FILE=pft_new.nc -DIFFUCO_LEAFCI=233. -CONDVEG_SNOWA=default -# IMPOSE_AZE is not set -SOILALB_FILE=soils_param.nc -SOILTYPE_FILE=soils_param.nc -ENERBIL_TSURF=280. -HYDROL_SNOW=0.0 -HYDROL_SNOWAGE=0.0 -HYDROL_SNOWICE=0.0 -HYDROL_SNOWICEAGE=0.0 -HYDROL_HDRY=1.0 -HYDROL_HUMR=1.0 -HYDROL_BQSB=default -HYDROL_GQSB=0.0 -HYDROL_DSG=0.0 -HYDROL_DSP=default -HYDROL_QSV=0.0 -HYDROL_OK_HDIFF=n -HYDROL_TAU_HDIFF=1800. -THERMOSOIL_TPRO=280. -RIVER_ROUTING=y -ROUTING_FILE=routing.nc -LAI_MAP=y -LAI_FILE=lai2D.nc -SECHIBA_QSINT=0.02 -ALB_BARE_MODEL = FALSE -PERCENT_THROUGHFALL_PFT = 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30. -RVEG_PFT = .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5 -CDRAG_FROM_GCM = .TRUE. -#LAND_USE=y -#VEGET_YEAR=0 -#VEGET_UPDATE=1Y Index: DZ4/trunk/output.def =================================================================== --- /LMDZ4/trunk/output.def (revision 6163) +++ (revision ) @@ -1,794 +1,0 @@ -######## Abderrahmane le 20 11 08 ######################## -# Niveaux de sorties et nom pour chaque variable dans les# -# fichiers histmth histday histhf histins histLES # -########################################################## -# Surface geop.height -flag_phis = 1, 1, 10, 1, 1 -name_phis = phis -# Grid area -flag_aire = 1, 1, 10, 1, 1 -name_aire = aire -# Surfac ter+lic -flag_contfracATM = 10, 1, 1, 10, 10 -name_contfracATM = contfracATM -# Surface terre OR -flag_contfracOR = 10, 1, 1, 10, 10 -name_contfracOR = contfracOR -# Grid area CONT -flag_aireTER = 10, 10, 1, 10, 10 -name_aireTER = aireTER -# Latent heat flux -flag_flat = 10, 1, 10, 10, 1 -name_flat = flat -# Sea Level Pressure -flag_slp = 1, 1, 1, 10, 1 -name_slp = slp -# Surface Temperature -flag_tsol = 1, 1, 1, 1, 1 -name_tsol = tsol -# Temperature 2m -flag_t2m = 1, 1, 1, 1, 1 -name_t2m = t2m -# Temperature min 2m -flag_t2m_min = 1, 1, 10, 10, 10 -name_t2m_min = t2m_min -# Temperature max 2m -flag_t2m_max = 1, 1, 10, 10, 10 -name_t2m_max = t2m_max -# Temperature ter lic oce sic at 2m -flag_t2m_ter = 10, 4, 10, 10, 10 -name_t2m_ter = t2m_ter -flag_t2m_lic = 10, 4, 10, 10, 10 -name_t2m_lic = t2m_lic -flag_t2m_oce = 10, 4, 10, 10, 10 -name_t2m_oce = t2m_oce -flag_t2m_sic = 10, 4, 10, 10, 10 -name_t2m_sic = t2m_sic -# 10-m wind speed -flag_wind10m = 1, 1, 1, 10, 10 -name_wind10m = wind10m -# 10-m wind speed max -flag_wind10max = 10, 1, 10, 10, 10 -name_wind10max = wind10max -# Sea-ice fraction -flag_sicf = 1, 1, 10, 10, 10 -name_sicf = sicf -# Specific humidity 2m -flag_q2m = 1, 1, 1, 1, 1 -name_q2m = q2m -# 10m zonal wind -flag_u10m = 1, 1, 1, 1, 1 -name_u10m = u10m -# 10m meridional wind -flag_v10m = 1, 1, 1, 1, 1 -name_v10m = v10m -# Surface Pressure -flag_psol = 1, 1, 1, 1, 1 -name_psol = psol -# Surface Air humidity -flag_qsurf = 1, 10, 10, 10, 10 -name_qsurf = qsurf -# 10m zonal wind (ter lic oce sic) -flag_u10m_ter = 10, 4, 10, 10, 10 -name_u10m_ter = u10m_ter -flag_u10m_lic = 10, 4, 10, 10, 10 -name_u10m_lic = u10m_lic -flag_u10m_oce = 10, 4, 10, 10, 10 -name_u10m_oce = u10m_oce -flag_u10m_sic = 10, 4, 10, 10, 10 -name_u10m_sic = u10m_sic -# 10m meridien wind (ter oce ice lic) -flag_v10m_ter = 10, 4, 10, 10, 10 -name_v10m_ter = v10m_ter -flag_v10m_lic = 10, 4, 10, 10, 10 -name_v10m_lic = v10m_lic -flag_v10m_oce = 10, 4, 10, 10, 10 -name_v10m_oce = v10m_oce -flag_v10m_sic = 10, 4, 10, 10, 10 -name_v10m_sic = v10m_sic -# Soil watter content -flag_qsol = 1, 10, 10, 1, 1 -name_qsol = qsol -# Number of dayrain(liq+sol) -flag_ndayrain = 1, 10, 10, 10, 10 -name_ndayrain = ndayrain -# Precip Totale liq+sol -flag_precip = 1, 1, 1, 1, 1 -name_precip = precip -# Large-scale Precip -flag_plul = 1, 1, 1, 1, 10 -name_plul = plul -# Convective Precip -flag_pluc = 1, 1, 1, 1, 10 -name_pluc = pluc -# Snow fall -flag_snow = 1, 1, 10, 1, 10 -name_snow = snow -# Evaporation -flag_evap = 1, 1, 10, 1, 10 -name_evap = evap -# Solar rad. at TOA -flag_tops = 1, 1, 10, 10, 10 -name_tops = tops -# CS Solar rad. at TOA -flag_tops0 = 1, 5, 10, 10, 10 -name_tops0 = tops0 -# IR rad. at TOA -flag_topl = 1, 1, 10, 1, 10 -name_topl = topl -# CR IR rad. at TOA -flag_topl0 = 1, 5, 10, 10, 10 -name_topl0 = topl0 -# SWup at TOA -flag_SWupTOA = 1, 4, 10, 10, 10 -name_SWupTOA = SWupTOA -# CR SWup at TOA -flag_SWupTOAclr = 1, 4, 10, 10, 10 -name_SWupTOAclr = SWupTOAclr -# SWdn at TOA -flag_SWdnTOA = 1, 4, 10, 10, 10 -name_SWdnTOA = SWdnTOA -# CR SWdn at TOA -flag_SWdnTOAclr = 1, 4, 10, 10, 10 -name_SWdnTOAclr = SWdnTOAclr -# SWup at 200hPa -flag_SWup200 = 1, 10, 10, 10, 10 -name_SWup200 = SWup200 -# CR SWup at 200hPa -flag_SWup200clr = 10, 1, 10, 10, 10 -name_SWup200clr = SWup200clr -# SWdn at 200hPa -flag_SWdn200 = 1, 10, 10, 10, 10 -name_SWdn200 = SWdn200 -# CR SWdn at 200hPa -flag_SWdn200clr = 10, 1, 10, 10, 10 -name_SWdn200clr = SWdn200clr -# LWup at 200mb -flag_LWup200 = 1, 10, 10, 10, 10 -name_LWup200 = LWup200 -# CR LWup at 200mb -flag_LWup200clr = 1, 10, 10, 10, 10 -name_LWup200clr = LWup200clr -# LWdn at 200mb -flag_LWdn200 = 1, 10, 10, 10, 10 -name_LWdn200 = LWdn200 -# CR LWdn at 200mb -flag_LWdn200clr = 1, 10, 10, 10, 10 -name_LWdn200clr = LWdn200clr -# Solar rad. at surf -flag_sols = 1, 1, 10, 1, 10 -name_sols = sols -# CR Solar rad. at surf -flag_sols0 = 1, 5, 10, 10, 10 -name_sols0 = sols0 -# IR rad. at surface -flag_soll = 1, 1, 10, 1, 10 -name_soll = soll -# CR IR rad. at surface -flag_soll0 = 1, 5, 10, 10, 10 -name_soll0 = soll0 -# Rayonnement au sol -flag_radsol = 1, 1, 10, 10, 10 -name_radsol = radsol -# SWup at surface -flag_SWupSFC = 1, 4, 10, 10, 10 -name_SWupSFC = SWupSFC -# CR SWup at surface -flag_SWupSFCclr = 1, 4, 10, 10, 10 -name_SWupSFCclr = SWupSFCclr -# SWdn at surface -flag_SWdnSFC = 1, 1, 10, 10, 10 -name_SWdnSFC = SWdnSFC -# CR at surface -flag_SWdnSFCclr = 1, 4, 10, 10, 10 -name_SWdnSFCclr = SWdnSFCclr -# LWup at surface -flag_LWupSFC = 1, 4, 10, 10, 10 -name_LWupSFC = LWupSFC -# CR LWup at surface -flag_LWupSFCclr = 1, 4, 10, 10, 10 -name_LWupSFCclr = LWupSFCclr -# LWdn at surface -flag_LWdnSFC = 1, 4, 10, 10, 10 -name_LWdnSFC = LWdnSFC -# CR LWdn at surface -flag_LWdnSFCclr = 1, 4, 10, 10, 10 -name_LWdnSFCclr = LWdnSFCclr -# Surf. total heat flux -flag_bils = 1, 2, 10, 1, 10 -name_bils = bils -# Sensible heat flux -flag_sens = 1, 1, 10, 1, 1 -name_sens = sens -# Heat flux derivation -flag_fder = 1, 2, 10, 1, 10 -name_fder = fder -# Thermal flux for snow melting -flag_ffonte = 1, 10, 10, 10, 10 -name_ffonte = ffonte -# Ice Calving -flag_fqcalving = 1, 10, 10, 10, 10 -name_fqcalving = fqcalving -# Land ice melt -flag_fqfonte = 1, 10, 10, 10, 10 -name_fqfonte = fqfonte -# Zonal wind stress (ter ice liq oce) -flag_taux_ter = 1, 4, 10, 1, 10 -name_taux_ter = taux_ter -flag_taux_lic = 1, 4, 10, 1, 10 -name_taux_lic = taux_lic -flag_taux_oce = 1, 4, 10, 1, 10 -name_taux_oce = taux_oce -flag_taux_sic = 1, 4, 10, 1, 10 -name_taux_sic = taux_sic -# Meridien wind stress (ter ice liq oce) -flag_tauy_ter = 1, 4, 10, 1, 10 -name_tauy_ter = tauy_ter -flag_tauy_lic = 1, 4, 10, 1, 10 -name_tauy_lic = tauy_lic -flag_tauy_oce = 1, 4, 10, 1, 10 -name_tauy_oce = tauy_oce -flag_tauy_sic = 1, 4, 10, 1, 10 -name_tauy_sic = tauy_sic -# % surface (ter ice liq oce) -flag_pourc_ter = 1, 4, 10, 1, 10 -name_pourc_ter = pourc_ter -flag_pourc_lic = 1, 4, 10, 1, 10 -name_pourc_lic = pourc_lic -flag_pourc_oce = 1, 4, 10, 1, 10 -name_pourc_oce = pourc_oce -flag_pourc_sic = 1, 4, 10, 1, 10 -name_pourc_sic = pourc_sic -# Fraction (ter ice liq oce) -flag_fract_ter = 1, 4, 10, 1, 10 -name_fract_ter = fract_ter -flag_fract_lic = 1, 4, 10, 1, 10 -name_fract_lic = fract_lic -flag_fract_oce = 1, 4, 10, 1, 10 -name_fract_oce = fract_oce -flag_fract_sic = 1, 4, 10, 1, 10 -name_fract_sic = fract_sic -# Surface temperature (ter ice liq oce) -flag_tsol_ter = 1, 4, 10, 1, 10 -name_tsol_ter = tsol_ter -flag_tsol_lic = 1, 4, 10, 1, 10 -name_tsol_lic = tsol_lic -flag_tsol_oce = 1, 4, 10, 1, 10 -name_tsol_oce = tsol_oce -flag_tsol_sic = 1, 4, 10, 1, 10 -name_tsol_sic = tsol_sic -# Sensible heat flux (ter ice liq oce) -flag_sens_ter = 1, 4, 10, 1, 10 -name_sens_ter = sens_ter -flag_sens_lic = 1, 4, 10, 1, 10 -name_sens_lic = sens_lic -flag_sens_oce = 1, 4, 10, 1, 10 -name_sens_oce = sens_oce -flag_sens_sic = 1, 4, 10, 1, 10 -name_sens_sic = sens_sic -# Latent heat flux (ter ice liq oce) -flag_lat_ter = 1, 4, 10, 1, 10 -name_lat_ter = lat_ter -flag_lat_lic = 1, 4, 10, 1, 10 -name_lat_lic = lat_lic -flag_lat_oce = 1, 4, 10, 1, 10 -name_lat_oce = lat_oce -flag_lat_sic = 1, 4, 10, 1, 10 -name_lat_sic = lat_sic -# LW (ter ice liq oce) -flag_flw_ter = 1, 10, 10, 10, 10 -name_flw_ter = flw_ter -flag_flw_lic = 1, 10, 10, 10, 10 -name_flw_lic = flw_lic -flag_flw_oce = 1, 10, 10, 10, 10 -name_flw_oce = flw_oce -flag_flw_sic = 1, 10, 10, 10, 10 -name_flw_sic = flw_sic -# SW (ter ice liq oce) -flag_fsw_ter = 1, 10, 10, 10, 10 -name_fsw_ter = fsw_ter -flag_fsw_lic = 1, 10, 10, 10, 10 -name_fsw_lic = fsw_lic -flag_fsw_oce = 1, 10, 10, 10, 10 -name_fsw_oce = fsw_oce -flag_fsw_sic = 1, 10, 10, 10, 10 -name_fsw_sic = fsw_sic -# Bilan sol (ter ice liq oce) -flag_wbils_ter = 1, 10, 10, 10, 10 -name_wbils_ter = wbils_ter -flag_wbils_lic = 1, 10, 10, 10, 10 -name_wbils_lic = wbils_lic -flag_wbils_oce = 1, 10, 10, 10, 10 -name_wbils_oce = wbils_oce -flag_wbils_sic = 1, 10, 10, 10, 10 -name_wbils_sic = wbils_sic -# Bilan eau (ter ice liq oce) -flag_wbilo_ter = 1, 10, 10, 10, 10 -name_wbilo_ter = wbilo_ter -flag_wbilo_lic = 1, 10, 10, 10, 10 -name_wbilo_lic = wbilo_lic -flag_wbilo_oce = 1, 10, 10, 10, 10 -name_wbilo_oce = wbilo_oce -flag_wbilo_sic = 1, 10, 10, 10, 10 -name_wbilo_sic = wbilo_sic -# Momentum drag coef -flag_cdrm = 1, 10, 10, 1, 10 -name_cdrm = cdrm -# Heat drag coef -flag_cdrh = 1, 10, 10, 1, 10 -name_cdrh = cdrh -# Low-level cloudiness -flag_cldl = 1, 1, 10, 10, 10 -name_cldl = cldl -# Mid-level cloudiness -flag_cldm = 1, 1, 10, 10, 10 -name_cldm = cldm -# High-level cloudiness -flag_cldh = 1, 1, 10, 10, 10 -name_cldh = cldh -# Total cloudiness -flag_cldt = 1, 1, 2, 10, 10 -name_cldt = cldt -# Cloud liquid water path -flag_cldq = 1, 1, 10, 10, 10 -name_cldq = cldq -# Cloud water path -flag_lwp = 1, 5, 10, 10, 10 -name_lwp = lwp -# Cloud ice water path -flag_iwp = 1, 5, 10, 10, 10 -name_iwp = iwp -# Zonal energy transport -flag_ue = 1, 10, 10, 10, 10 -name_ue = ue -# Merid energy transport -flag_ve = 1, 10, 10, 10, 10 -name_ve = ve -# Zonal humidity transport -flag_uq = 1, 10, 10, 10, 10 -name_uq = uq -# Merid humidity transport -flag_vq = 1, 10, 10, 10, 10 -name_vq = vq -# Conv avlbl pot ener -flag_cape = 1, 10, 10, 10, 10 -name_cape = cape -# Cld base pressure -flag_pbase = 1, 10, 10, 10, 10 -name_pbase = pbase -# Cld top pressure -flag_ptop = 1, 4, 10, 10, 10 -name_ptop = ptop -# Cld base mass flux -flag_fbase = 1, 10, 10, 10, 10 -name_fbase = fbase -# Precipitable water -flag_prw = 1, 1, 10, 10, 10 -name_prw = prw -# Boundary Layer Height -flag_s_pblh = 1, 10, 10, 1, 1 -name_s_pblh = pblh -# t at Boundary Layer Height -flag_s_pblt = 1, 10, 10, 1, 1 -name_s_pblt = pblt -# Condensation level -flag_s_lcl = 1, 10, 10, 1, 10 -name_s_lcl = lcl -# Conv avlbl pot enerfor ABL -flag_s_capCL = 1, 10, 10, 1, 10 -name_s_capCL = capCL -# Liq Water in BL -flag_s_oliqCL = 1, 10, 10, 1, 10 -name_s_oliqCL = oliqCL -# Instability criteria(ABL) -flag_s_cteiCL = 1, 10, 10, 1, 1 -name_s_cteiCL = cteiCL -# Exces du thermique -flag_s_therm = 1, 10, 10, 1, 1 -name_s_therm = therm -# deep_cape(HBTM2) -flag_s_trmb1 = 1, 10, 10, 1, 10 -name_s_trmb1 = trmb1 -# inhibition (HBTM2) -flag_s_trmb2 = 1, 10, 10, 1, 10 -name_s_trmb2 = trmb2 -# Point Omega (HBTM2) -flag_s_trmb3 = 1, 10, 10, 1, 10 -name_s_trmb3 = trmb3 -# Bilan au sol sur ocean slab -flag_slab_bils = 1, 1, 10, 10, 10 -name_slab_bils = slab_bils -# ALE BL -flag_ale_bl = 1, 1, 1, 1, 10 -name_ale_bl = ale_bl -# alp_bl -flag_alp_bl = 1, 1, 1, 1, 10 -name_alp_bl = alp_bl -# ale_wk -flag_ale_wk = 1, 1, 1, 1, 10 -name_ale_wk = ale_wk -# alp_wk -flag_alp_wk = 1, 1, 1, 1, 10 -name_alp_wk = alp_wk -# ale -flag_ale = 1, 1, 1, 1, 10 -name_ale = ale -# alp -flag_alp = 1, 1, 1, 1, 10 -name_alp = alp -# Convective INhibition -flag_cin = 1, 1, 1, 1, 10 -name_cin = cin -# WAPE -flag_wape = 1, 1, 1, 1, 10 -name_wape = wape -# u, v w t q et phi aux niveaux 200, 500, 700 et 850 hPa -flag_u850 = 1, 1, 3, 10, 10 -name_u850 = u850 -flag_u700 = 1, 1, 3, 10, 10 -name_u700 = u700 -flag_u500 = 1, 1, 3, 10, 10 -name_u500 = u500 -flag_u200 = 1, 1, 3, 10, 10 -name_u200 = u200 -flag_v850 = 1, 1, 3, 10, 10 -name_v850 = v850 -flag_v700 = 1, 1, 3, 10, 10 -name_v700 = v700 -flag_v500 = 1, 1, 3, 10, 10 -name_v500 = v500 -flag_v200 = 1, 1, 3, 10, 10 -name_v200 = v200 -flag_w850 = 1, 1, 3, 10, 10 -name_w850 = w850 -flag_w700 = 1, 1, 3, 10, 10 -name_w700 = w700 -flag_w500 = 1, 1, 3, 10, 10 -name_w500 = w500 -flag_w200 = 1, 1, 3, 10, 10 -name_w200 = w200 -flag_t850 = 1, 1, 3, 10, 10 -name_t850 = t850 -flag_t700 = 1, 1, 3, 10, 10 -name_t700 = t700 -flag_t500 = 1, 1, 3, 10, 10 -name_t500 = t500 -flag_t200 = 1, 1, 3, 10, 10 -name_t200 = t200 -flag_q850 = 1, 1, 3, 10, 10 -name_q850 = q850 -flag_q700 = 1, 1, 3, 10, 10 -name_q700 = q700 -flag_q500 = 1, 1, 3, 10, 10 -name_q500 = q500 -flag_q200 = 1, 1, 3, 10, 10 -name_q200 = q200 -flag_phi850 = 1, 1, 3, 10, 10 -name_phi850 = phi850 -flag_phi700 = 1, 1, 3, 10, 10 -name_phi700 = phi700 -flag_phi500 = 1, 1, 3, 10, 10 -name_phi500 = phi500 -flag_phi200 = 1, 1, 3, 10, 10 -name_phi200 = phi200 -# Temp mixte oce-sic -flag_t_oce_sic = 1, 10, 10, 10, 10 -name_t_oce_sic = t_oce_sic -# Weak inversion -flag_weakinv = 10, 1, 10, 10, 10 -name_weakinv = weakinv -# dTheta mini -flag_dthmin = 10, 1, 10, 10, 10 -name_dthmin = dthmin -# 10m zonal and meriden wind (ter sic lic oce) -flag_u10_ter = 10, 4, 10, 10, 10 -name_u10_ter = u10_ter -flag_u10_lic = 10, 4, 10, 10, 10 -name_u10_lic = u10_lic -flag_u10_oce = 10, 4, 10, 10, 10 -name_u10_oce = u10_oce -flag_u10_sic = 10, 4, 10, 10, 10 -name_u10_sic = u10_sic -flag_v10_ter = 10, 4, 10, 10, 10 -name_v10_ter = v10_ter -flag_v10_lic = 10, 4, 10, 10, 10 -name_v10_lic = v10_lic -flag_v10_oce = 10, 4, 10, 10, 10 -name_v10_oce = v10_oce -flag_v10_sic = 10, 4, 10, 10, 10 -name_v10_sic = v10_sic -# Cloud optical thickness -flag_cldtau = 10, 5, 10, 10, 10 -name_cldtau = cldtau -# Cloud optical emissivity -flag_cldemi = 10, 5, 10, 10, 10 -name_cldemi = cldemi -# Relative humidity at 2m -flag_rh2m = 10, 5, 10, 10, 10 -name_rh2m = rh2m -# Saturant humidity at 2m -flag_qsat2m = 10, 5, 10, 10, 10 -name_qsat2m = qsat2m -# Surface air potential temperature -flag_tpot = 10, 5, 10, 10, 10 -name_tpot = tpot -# Surface air equivalent potential temperature -flag_tpote = 10, 5, 10, 10, 10 -name_tpote = tpote -# TKE, tke max and tke (ter sic lic oce) -flag_tke = 4, 10, 10, 10, 10 -name_tke = tke -flag_tke_max = 4, 10, 10, 10, 10 -name_tke_max = tke_max -flag_tke_ter = 10, 4, 10, 10, 10 -name_tke_ter = tke_ter -flag_tke_lic = 10, 4, 10, 10, 10 -name_tke_lic = tke_lic -flag_tke_oce = 10, 4, 10, 10, 10 -name_tke_oce = tke_oce -flag_tke_sic = 10, 4, 10, 10, 10 -name_tke_sic = tke_sic -flag_tke_max_ter = 10, 4, 10, 10, 10 -name_tke_max_ter = tke_max_ter -flag_tke_max_lic = 10, 4, 10, 10, 10 -name_tke_max_lic = tke_max_lic -flag_tke_max_oce = 10, 4, 10, 10, 10 -name_tke_max_oce = tke_max_oce -flag_tke_max_sic = 10, 4, 10, 10, 10 -name_tke_max_sic = tke_max_sic -# Kz melange -flag_kz = 4, 10, 10, 10, 10 -name_kz = kz -# Kz max melange -flag_kz_max = 4, 10, 10, 10, 10 -name_kz_max = kz_max -# Sfce net SW radiation OR -flag_SWnetOR = 10, 10, 2, 10, 10 -name_SWnetOR = SWnetOR -# Sfce incident SW radiation OR -flag_SWdownOR = 10, 10, 2, 10, 10 -name_SWdownOR = SWdownOR -# Sfce incident LW radiation OR -flag_LWdownOR = 10, 10, 2, 10, 10 -name_LWdownOR = LWdownOR -# Solid Large-scale Precip -flag_snowl = 10, 1, 10, 10, 10 -name_snowl = snowl -# cape max -flag_cape_max = 10, 1, 10, 10, 10 -name_cape_max = cape_max -# Down. IR rad. at surface -flag_solldown = 10, 1, 10, 1, 10 -name_solldown = solldown -# Boundary-layer dTs(o) -flag_dtsvdfo = 10, 10, 10, 1, 10 -name_dtsvdfo = dtsvdfo -# Boundary-layer dTs(t) -flag_dtsvdft = 10, 10, 10, 1, 10 -name_dtsvdft = dtsvdft -# Boundary-layer dTs(g) -flag_dtsvdfg = 10, 10, 10, 1, 10 -name_dtsvdfg = dtsvdfg -# Boundary-layer dTs(g) -flag_dtsvdfi = 10, 10, 10, 1, 10 -name_dtsvdfi = dtsvdfi -# rugosity -flag_rugs = 10, 10, 10, 1, 1 -name_rugs = rugs -# Cloud liquid water content -flag_lwcon = 2, 5, 10, 10, 1 -name_lwcon = lwcon -# Cloud ice water content -flag_iwcon = 2, 5, 10, 10, 10 -name_iwcon = iwcon -# Air temperature -flag_temp = 2, 3, 4, 1, 1 -name_temp = temp -# Potential air temperature -flag_theta = 2, 3, 4, 1, 1 -name_theta = theta -# Specific humidity -flag_ovap = 2, 3, 4, 1, 1 -name_ovap = ovap -# ? -flag_ovapinit = 2, 3, 4, 1, 1 -name_ovapinit = ovapinit -# ? -flag_wvapp = 2, 10, 10, 10, 10 -name_wvapp = wvapp -# Geopotential height -flag_geop = 2, 3, 10, 1, 1 -name_geop = geop -# Zonal wind -flag_vitu = 2, 3, 4, 1, 1 -name_vitu = vitu -# Meridional wind -flag_vitv = 2, 3, 4, 1, 1 -name_vitv = vitv -# Vertical wind -flag_vitw = 2, 3, 10, 10, 1 -name_vitw = vitw -# Air pressure -flag_pres = 2, 3, 10, 1, 1 -name_pres = pres -# Cloud Fraction -flag_rneb = 2, 5, 10, 10, 1 -name_rneb = rneb -# Convective Cloud Fraction -flag_rnebcon = 2, 5, 10, 10, 1 -name_rnebcon = rnebcon -# Relative humidity -flag_rhum = 2, 10, 10, 10, 10 -name_rhum = rhum -# Ozone concentration -flag_ozone = 2, 10, 10, 10, 10 -name_ozone = ozone -# saturated updraft -flag_upwd = 2, 10, 10, 10, 10 -name_upwd = upwd -# Physics dT -flag_dtphy = 2, 10, 10, 10, 1 -name_dtphy = dtphy -# Physics dq -flag_dqphy = 2, 10, 10, 10, 1 -name_dqphy = dqphy -# Convective precipitation lic and ice -flag_pr_con_l = 2, 10, 10, 10, 10 -name_pr_con_l = pr_con_l -flag_pr_con_i = 2, 10, 10, 10, 10 -name_pr_con_i = pr_con_i -# Large scale precipitation lic and ice -flag_pr_lsc_l = 2, 10, 10, 10, 10 -name_pr_lsc_l = pr_lsc_l -flag_pr_lsc_i = 2, 10, 10, 10, 10 -name_pr_lsc_i = pr_lsc_i -# Albedo surf, Snow age and rugosity (ter sic lic oce) -flag_albe_ter = 3, 4, 10, 1, 10 -name_albe_ter = albe_ter -flag_albe_lic = 3, 4, 10, 1, 10 -name_albe_lic = albe_lic -flag_albe_oce = 3, 4, 10, 1, 10 -name_albe_oce = albe_oce -flag_albe_sic = 3, 4, 10, 1, 10 -name_albe_sic = albe_sic -flag_ages_ter = 3, 10, 10, 10, 10 -name_ages_ter = ages_ter -flag_ages_lic = 3, 10, 10, 10, 10 -name_ages_lic = ages_lic -flag_ages_oce = 3, 10, 10, 10, 10 -name_ages_oce = ages_oce -flag_ages_sic = 3, 10, 10, 10, 10 -name_ages_sic = ages_sic -flag_rugs_ter = 3, 4, 10, 1, 10 -name_rugs_ter = rugs_ter -flag_rugs_lic = 3, 4, 10, 1, 10 -name_rugs_lic = rugs_lic -flag_rugs_oce = 3, 4, 10, 1, 10 -name_rugs_oce = rugs_oce -flag_rugs_sic = 3, 4, 10, 1, 10 -name_rugs_sic = rugs_sic -# Surface albedo -flag_albs = 3, 10, 10, 1, 10 -name_albs = albs -# Surface albedo LW -flag_albslw = 3, 10, 10, 1, 10 -name_albslw = albslw -# Convective Cloud Liquid water content -flag_clwcon = 4, 10, 10, 10, 10 -name_clwcon = clwcon -# undilute adiab updraft -flag_Ma = 4, 10, 10, 10, 10 -name_Ma = Ma -# saturated downdraft -flag_dnwd = 4, 10, 10, 10, 10 -name_dnwd = dnwd -# unsat. downdraft -flag_dnwd0 = 4, 10, 10, 10, 10 -name_dnwd0 = dnwd0 -# Dynamics dT dQ dU dV, ..... -flag_dtdyn = 4, 10, 10, 10, 1 -name_dtdyn = dtdyn -flag_dqdyn = 4, 10, 10, 10, 1 -name_dqdyn = dqdyn -flag_dudyn = 4, 10, 10, 10, 1 -name_dudyn = dudyn -flag_dvdyn = 4, 10, 10, 10, 1 -name_dvdyn = dvdyn -flag_dtcon = 4, 5, 10, 10, 10 -name_dtcon = dtcon -flag_ducon = 4, 10, 10, 10, 10 -name_ducon = ducon -flag_dqcon = 4, 5, 10, 10, 10 -name_dqcon = dqcon -flag_dtwak = 4, 5, 10, 10, 10 -name_dtwak = dtwak -flag_dqwak = 4, 5, 10, 10, 10 -name_dqwak = dqwak -flag_wake_h = 4, 5, 10, 10, 10 -name_wake_h = wake_h -flag_wake_s = 4, 5, 10, 10, 10 -name_wake_s = wake_s -flag_wake_deltat = 4, 5, 10, 10, 10 -name_wake_deltat = wake_deltat -flag_wake_deltaq = 4, 5, 10, 10, 10 -name_wake_deltaq = wake_deltaq -flag_wake_omg = 4, 5, 10, 10, 10 -name_wake_omg = wake_omg -flag_Vprecip = 10, 10, 10, 10, 10 -name_Vprecip = Vprecip -flag_ftd = 4, 5, 10, 10, 10 -name_ftd = ftd -flag_fqd = 4, 5, 10, 10, 10 -name_fqd = fqd -flag_dtlsc = 4, 10, 10, 10, 10 -name_dtlsc = dtlsc -flag_dtlschr = 4, 10, 10, 10, 10 -name_dtlschr = dtlschr -flag_dqlsc = 4, 10, 10, 10, 10 -name_dqlsc = dqlsc -flag_dtvdf = 4, 10, 10, 1, 10 -name_dtvdf = dtvdf -flag_dqvdf = 4, 10, 10, 1, 10 -name_dqvdf = dqvdf -flag_dteva = 4, 10, 10, 10, 10 -name_dteva = dteva -flag_dqeva = 4, 10, 10, 10, 10 -name_dqeva = dqeva -flag_ptconv = 4, 10, 10, 10, 10 -name_ptconv = ptconv -flag_ratqs = 4, 10, 10, 10, 10 -name_ratqs = ratqs -flag_dtthe = 4, 10, 10, 10, 10 -name_dtthe = dtthe -flag_f_th = 4, 10, 10, 10, 10 -name_f_th = f_th -flag_e_th = 4, 10, 10, 10, 10 -name_e_th = e_th -flag_w_th = 4, 10, 10, 10, 10 -name_w_th = w_th -flag_lambda_th = 4, 10, 10, 10, 10 -name_lambda_th = lambda_th -flag_q_th = 4, 10, 10, 10, 10 -name_q_th = q_th -flag_a_th = 4, 10, 10, 10, 10 -name_a_th = a_th -flag_d_th = 4, 10, 10, 10, 10 -name_d_th = d_th -flag_f0_th = 4, 10, 10, 10, 10 -name_f0_th = f0_th -flag_zmax_th = 4, 10, 10, 10, 10 -name_zmax_th = zmax_th -flag_dqthe = 4, 10, 10, 10, 1 -name_dqthe = dqthe -flag_dtajs = 4, 10, 10, 10, 10 -name_dtajs = dtajs -flag_dqajs = 4, 10, 10, 10, 10 -name_dqajs = dqajs -flag_dtswr = 4, 10, 10, 10, 1 -name_dtswr = dtswr -flag_dtsw0 = 4, 10, 10, 10, 10 -name_dtsw0 = dtsw0 -flag_dtlwr = 4, 10, 10, 10, 1 -name_dtlwr = dtlwr -flag_dtlw0 = 4, 10, 10, 10, 10 -name_dtlw0 = dtlw0 -flag_dtec = 4, 10, 10, 10, 10 -name_dtec = dtec -flag_duvdf = 4, 10, 10, 10, 10 -name_duvdf = duvdf -flag_dvvdf = 4, 10, 10, 10, 10 -name_dvvdf = dvvdf -flag_duoro = 4, 10, 10, 10, 10 -name_duoro = duoro -flag_dvoro = 4, 10, 10, 10, 10 -name_dvoro = dvoro -flag_dulif = 4, 10, 10, 10, 10 -name_dulif = dulif -flag_dvlif = 4, 10, 10, 10, 10 -name_dvlif = dvlif -###! Attention a refaire correctement -flag_trac01 = 4, 10, 10, 10, 10 -name_trac01 = trac01 -flag_trac02 = 4, 10, 10, 10, 10 -name_trac02 = trac02 - Index: DZ4/trunk/physiq.def =================================================================== --- /LMDZ4/trunk/physiq.def (revision 6163) +++ (revision ) @@ -1,219 +1,0 @@ -## $Id$ -# -# PARAMETRES ANCIENNEMENT DANS gcm.def -## Cycle diurne ou non -cycle_diurne=y -## Soil Model ou non -soil_model=y -## Choix ou non de New oliq -new_oliq=y -## Activation ou non de la parametrisation de Hines pour la strato -ok_hines=n -## Orodr ou non pour l orographie -ok_orodr=y -## Orolf ou non pour l orographie -ok_orolf=y -## Si = .T. , lecture du fichier limit avec la bonne annee -ok_limitvrai=n -## Nombre d'appels des routines de rayonnements ( par jour) -nbapp_rad=24 -## Flag pour la convection : 1 pour LMD, 2 pour Tiedtke, 3 KE(nvlle version JYG), 30 KE(version IPCC AR4), 4 KE vect -iflag_con=30 -## Facteur multiplication des precip convectives dans KE -cvl_corr=1.0 -## Facteur additif pour l'albedo -pmagic=0.008 -# -# -# -# Parametres fichiers de sortie -# -### type_run = type run par rapport aux fichiers et variables de sortie -# - type_run = CLIM/ENSP (=1) -# - type_run = AMIP/CFMI (=2) -type_run=AMIP -### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc -OK_journe=n -### OK_mensuel= y sortir fichier mensuel histmth.nc, =n pas de fichier histmth.nc -OK_mensuel=y -### OK_instan=y, ecrire sorties "instantannees" (chaque pas de temps de la physique) -OK_instan=n -### OK_hf=y, ecrire sorties hautes frequence histhf.nc, =n pas de fichier histhf.nc -ok_hf=n -# -# Parametres niveau de sorties differents fichiers -# -### lev_histhf=0-4, niveau de sortie fichier "histhf.nc" -# - lev_histhf=0 => pas de sorties histhf.nc -# - lev_histhf=2 => defaut -# - lev_histhf=3 => variables sur niveaux standards -# - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele => fichier. histhf3d.nc -lev_histhf=2 -### lev_histday=0-5, niveau de sortie fichier "histday.nc" -# - lev_histday=0 => pas de sorties lev_histday.nc -# - lev_histday=2 => defaut -# - lev_histday=3 => + champs 3D => F. Lott -# - lev_histday=4 => + champs sous-surfaces -# - lev_histday=5 => + champs F. Aires -lev_histday=2 -### lev_histmth=0-4, niveau de sortie fichier "histmth.nc" -# - lev_histmth=0 => pas de sorties lev_histmth.nc -# - lev_histmth=2 => defaut -# - lev_histmth=3 => albedo, rugosite sous-surfaces -# - lev_histmth=4 => champs tendances 3d -lev_histmth=2 -### ecrit_hf = frequence ecriture fichier histhf.nc en jours -ecrit_hf=0.250 -### ecrit_day = frequence ecriture fichier histday.nc en jours -ecrit_day=1. -### ecrit_mth = frequence ecriture fichier histmth.nc en jours -ecrit_mth=30 -### freqin_isccp = frequence input en secondes du simulateur ISCCP -freq_ISCCP=10800. -### freqout_isccp = frequence output en jours du simulateur ISCCP -ecrit_ISCCP=30 -### niveau du diagnostique de conservation d energie -if_ebil=0 -# -# parametres KE -# -### epmax = Efficacite precipitation maximale -epmax = .999 -### ok_adj_ema = ?? pas utilise -ok_adj_ema = n -### iflag_clw Flag calcul eau liquide -# - iflag_clw=0 : qcond_incld(i,l) = em_qcondc(l) -# - iflag_clw=1 : qcond_incld(i,l) = em_qcond(l) -# - iflag_clw=2 : eau liquide diagnostique en fonction de la Precip -iflag_clw = 1 -# -# parametres nuages -# -### cld_lc_lsc contenu en eau liquide des nuages large-scale (fisrtilp) -cld_lc_lsc = 4.16e-4 -### cld_lc_con contenu en eau liquide des nuages convectifs (fisrtilp) -cld_lc_con = 4.16e-4 -### cld_tau_lsc cte de temps utilisee pour eliminer l eau large-scale (fisrtilp) -cld_tau_lsc = 1800. -### cld_tau_con cte de temps utilisee pour eliminer l eau convective (fisrtilp) -cld_tau_con = 1800. -### ffallv_lsc cte utilisee dans calcul vitesse de chute cristaux de glace large-scale (fisrtilp) -ffallv_lsc = 0.5 -### ffallv_lsc cte utilisee dans calcul vitesse de chute cristaux de glace convectifs (fisrtilp) -ffallv_con = 0.5 -### coef_eva coef evaporation precips eau/glace (fisrtilp/fisrtilp_tr?/conlmd?) -coef_eva = 2.e-5 -### reevap_ice reevaporation de toute la precip dans la couche du dessous pour la glace (fisrtilp) -reevap_ice = y -### iflag_cldcon flag pour calculer ratqsc=F(ratqsbas,fact_cldcon,q_seri) (physiq) -# - iflag_cldcon<=-1 diag. rain_Tiedtke -# - iflag_cldcon=1, ratqsc=ratqsbas+fact_cldcon*(q_seri(1)-q_seri(k))/q_seri(k) -# - iflag_cldcon=1/2, ratqs=max(ratqs,ratqsc) -# - iflag_cldcon=3, ratqs=ratqss -iflag_cldcon = 3 -### iflag_pdf : flag calcul distribution sous-maille de l eau et des nuages -# - iflag_pdf=0, version ratqs, -# - iflag_pdf=1, calcul eau condensee, fraction nuageuse, eau nuageuse a partir -# - des PDFS de Sandrine Bony -iflag_pdf = 1 -### fact_cldcon constante calcul ratqsc (voir iflag_cldcon) et proprietes nuages convectifs, clwcon0 (physiq.F) -fact_cldcon = 1. -### facttemps= facteur de relaxation de ratqs (iflag_cldcon=1/2) et rnebcon (iflag_cldcon=3) -facttemps = 0. -## ok_newmicro =y appel newmicro , =n appel nuage (calcul epaisseur optique et emmissivite des nuages) -ok_newmicro = y -### iflag_ratqs=0 correspond a la version IPCC AR4 -iflag_ratqs=0 -### ratqsbas ratqs en bas si iflag_cldcon=1 -ratqsbas = 0.005 -### ratqshaut ratqs en haut pour ratqss "stables" -ratqshaut = 0.33 -### rad_froid rayon cristaux des nuages de glace "froids" -rad_froid = 35 -### rad_chau1 rayon goutelettes d eau chauds", en haut: k=4-klev -rad_chau1=12 -### rad_chau2 rayon goutelettes d eau chauds", en bas: k=1-3 -rad_chau2=11 -# -# Coefficient et parametres sur les drags -# -f_cdrag_ter=1. -f_cdrag_stable=1. -f_cdrag_oce=0.8 -f_rugoro=0. -### ksta_ter coef.diffusion minimale sur terre/sic/lic -ksta_ter=1.e-7 -### cdmmax = cdrag maximum pour le moment -cdmmax = 2.5E-3 -### cdhmax = cdrag maximum pour l energie -cdhmax = 2.0E-3 -# -# Parametres "orbitaux/ ere geologique" -# -### R_ecc = Excentricite -R_ecc = 0.016715 -### R_peri = Equinoxe -R_peri = 102.7 -### R_incl = Inclinaison -R_incl = 23.441 -### solaire = Constante solaire -solaire = 1361. -# -# Taux gaz a effet de serre -# -### co2_ppm = taux CO2 en ppm -co2_ppm = 367. -### CH4_ppb = taux CH4 en ppb -CH4_ppb = 1760. -### N2O_ppb = taux N2O en ppb -N2O_ppb = 316. -### CFC11_ppt = taux CFC11 en ppt -CFC11_ppt = 741.2 -### CFC12_ppt = taux CFC12 en ppt -CFC12_ppt = 535. -# -# Parametres effets directs/indirects des "aerosols" -# -### ok_ade=y/n flag Aerosol direct effect -ok_ade=n -### ok_aie=y/n flag Aerosol indirect effect -ok_aie=n -### aer_type = Aerosol variation type : actuel / preind / scenario / annuel -aer_type=actuel -### type of coupled aerosol =1 (default) =2 => bc only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol -flag_aerosol=6 -### bl95_b0 = Parameter in CDNC-maer link (Boucher&Lohmann 1995) -bl95_b0=1.7 -### bl95_b1 = Parameter in CDNC-maer link (Boucher&Lohmann 1995) -bl95_b1=0.2 -# -# Parametre de lecture de l'ozone -# -# Allowed values are 0, 1 and 2 -# 0: do not read an ozone climatology -# 1: read a single ozone climatology that will be used day and night -# 2: read two ozone climatologies, the average day and night climatology and the daylight climatology -read_climoz=0 -# -# Parametres simulateur COSP (CFMIP Observational Simulator Package) -# -### ok_cosp=y/n flag simulateur COSP -ok_cosp=n -## freq_COSP = frequence d'appel de COSP en secondes -freq_COSP=10800. -# -# Parametres simulateur ISCCP -# -### ok_isccp=y/n flag simulateur ISCCP -ok_isccp=n -### top_height = flag choix calcul nuages par le simulateur en utilisant -# - les donnees IR et/ou VIS et l algorithme ISCCP-D1 -# - top_height = 1 -> algo IR-VIS -# - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)" -# - top_height = 3 -> algo IR -top_height = 1 -### overlap = Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP -# - overlap=1 Max overlap -# - overlap=2 Random overlap -# - overlap=3 Max/Random overlap -overlap = 3 Index: DZ4/trunk/run.def =================================================================== --- /LMDZ4/trunk/run.def (revision 6163) +++ (revision ) @@ -1,47 +1,0 @@ -# $Id$ -# -## Fichier de configuration general -## -INCLUDEDEF=physiq.def -INCLUDEDEF=gcm.def -INCLUDEDEF=orchidee.def -INCLUDEDEF=output.def -## Type de calendrier utilise -## valeur possible: earth_360d (defaut), earth_365d, earth_366d -calend=earth_360d -## Jour de l'etat initial ( = 350 si 20 Decembre ,par expl. ,comme ici ) -dayref=1 -## Annee de l'etat initial ( avec 4 chiffres ) -anneeref=1980 -## Nombre de jours d'integration -nday=5 -## Remise a zero de la date initiale -raz_date=0 -## periode de sortie des variables de controle (en pas) -iconser=240 -## sorties instantanees dans la dynamique (fichiers dyn_hist.nc and co.) -ok_dyn_ins=n -## periode d'ecriture des sorties instantanees dans la dynamique -## (en pas dynamiques) -iecri=960 -## sorties de valeurs moyennes dans la dynamique (fichiers dyn_hist_ave.nc and co.) -ok_dyn_ave=n -## periode de stockage des moyennes dans la dynamique et dans dynzon -periodav=30 -## flag de sortie dynzon -ok_dynzon=n -## activation du calcul d equilibrage de charge -adjust=n -## activation du filtre fft -use_filtre_fft=n -## niveau d'impression de controle -prt_level=1 -## -## Informations sur la configuration utilisee -## -### type_ocean = force / slab /couple -type_ocean=force -### version_ocean = nemo / opa8 -version_ocean=nemo -### VEGET= y si ORCHIDEE, =n si bucket -VEGET=n Index: DZ4/trunk/traceur.def =================================================================== --- /LMDZ4/trunk/traceur.def (revision 6163) +++ (revision ) @@ -1,5 +1,0 @@ -4 -14 14 H2Ov -10 10 H2Ol -10 10 RN -10 10 PB