Changeset 671 for trunk/LMDZ.GENERIC/makegcm_g95

Timestamp:

May 25, 2012, 11:11:42 AM (13 years ago)

Author:

emillour

Message:

Generic GCM:

• Significant update on how the number of scatterers is managed: Instead of having to manualy change 'nearkind' in radinc_h.F90, the number of scatterers must now be set when compiling, using makegcm "makegcm -s 1" for one scatterer or "makegcm -s 2" for two (e.g. dust and water ice), default behaviour (ie not specifying -s #) is -s 1 Modified phystd/radinc_h.F90 , added directory phystd/scatterers with script make_scatterers , and adapted makegcm* scripts.

EM

File:

• trunk/LMDZ.GENERIC/makegcm_g95 (modified) (15 diffs)

trunk/LMDZ.GENERIC/makegcm_g95

@@ -6 +6 @@
 ########################################################################
 set dim="64x48x32"
-set physique=mars
+set physique=std
 set phys="PHYS=$physique"
-set include='-I$(LIBF)/grid -I$(LIBF)/bibio -I. -I$(LIBF)/aeronomars'
+set include='-I$(LIBF)/grid -I$(LIBF)/bibio -I.'
 set ntrac = 1
 set filtre=filtrez
 set grille=reg
 set dyntype="dyn"
+set bands="32x36"
+set scatterers="1"
 ########################################################################
 # path a changer contenant les sources et les objets du modele
@@ -25 +27 @@
 setenv LIBOGCM $LMDGCM/libo
 # NetCDF
-# 64 bit machines
-#  setenv NCDFLIB /donnees/emlmd/netcdf64_g95/netcdf-3.6.1/lib
-#  setenv NCDFINC /donnees/emlmd/netcdf64_g95/netcdf-3.6.1/include
+# 64 bit machines what is below is OK on LMD machines
+  setenv NCDFLIB /donnees/emlmd/netcdf64-4.0.1_g95/lib
+  setenv NCDFINC /donnees/emlmd/netcdf64-4.0.1_g95/include
 #else
 #  setenv NCDFLIB /distrib/local/netcdf/g95_32bits/lib
@@ -180 +182 @@
 else if $LINUX then
 # Ehouarn 'g95' compiler
-   set optim90=" -O2 -fstatic -fzero -Wall "
-   set optimtru90=" -O2 -fstatic -fzero -Wall "
+   set optim90="-O3 -fstatic -funroll-loops "
+   set optimtru90="-O3 -fstatic -funroll-loops "
 #   set opt_link=" -Mfree -L/usr/local/pgi/linux86/lib -lpgf90 -lpgftnrtl -lpghpf -lpghpf2 -L$NCDFLIB -lnetcdf -Wl,-Bstatic "
 # Ehouarn: there is no /usr/local/pgi/linux86/lib, but it doesn't matter
@@ -241 +243 @@
 
 
-
-
-Par default, la commande makegcm:
----------------------------------
-
-1. compile une serie de sous programmes se trouvant dans des sous-repertoires
-de $LMDGCM/libf.
-Les sous programmes sont ensuite stokes sur dans des librairies FORTRAN
-sur $LIBOGCM.
-
-2. Ensuite, makegcm compile le programme prog.f se trouvant par default sur
-$LMDGCM/libf/dyn3d et effectue le lien avec l ensemble des librairies.
-
-La variable '$LMDGCM' doit etre initialisee dans votre .cshrc ou en dur
-dans la comande makegcm.
-
-La commande makegcm est faite pour permettre de gerer en parallele des
-versions differentes du modele, compilees avec des options de compilation
-et des dimensions differentes sans avoir a chaque fois a recompiler tout
-le modele.
-
-Les librairies FORTRAN sont stoquees sur le directory $LIBOGCM.
+The makegcm script:
+-------------------
+
+1. compiles a series of subroutines located in the $LMDGCM/libf
+ sub-directories.
+ The objects are then stored in the libraries in $LIBOGCM.
+
+2. then, makegcm compiles program prog.f located by default in
+$LMDGCM/libf/dyn3d and makes the link with the libraries.
+
+Environment Variables '$LMDGCM' and '$LIBOGCM'
+ must be set as environment variables or directly
+ in the makegcm file.
+
+The makegcm command is used to control the different versions of the model
+ in parallel, compiled using the compilation options 
+ and the various dimensions, without having to recompile the whole model.
+
+The FORTRAN libraries are stored in directory $LIBOGCM.
 
 
@@ -268 +267 @@
 --------
 
-Les options suivantes peuvent etre definies soit par defaut en editant le
-"script" makegcm, soit en interactif:
-
--d imxjmxlm  ou im, jm, et lm sont resp. le nombre de longitudes, latitudes
-             et couches verticales. 
-
--t ntrac   selectionne le nombre de traceur advectes par la dynamique.
-           Dans les versions courantes du modele terrestre on a par exemple
-           ntrac=2 pour l'eau vapeur et liquide
-
-             L'effet des options -d et -t est d'ecraser le fichier 
+The following options can either be defined by default by editing the
+makegcm "script", or in interactive mode:
+
+-d imxjmxlm  where im, jm, and lm are the number of longitudes,
+             latitudes and vertical layers respectively.
+
+-t ntrac   Selects the number of tracers present in the model
+
+             Options -d and -t overwrite file 
              $LMDGCM/libf/grid/dimensions.h
-             qui contient sous forme de 4 PARAMETER FORTRAN les 3 dimensions
-             de la grille horizontale im, jm, lm plus le nombre de traceurs
-             advectes passivement par la dynamique ntrac, par un nouveu fichier
+             which contains the 3 dimensions of the
+             horizontal grid 
+             im, jm, lm plus the number of tracers passively advected
+             by the dynamics ntrac,
+             in 4 PARAMETER FORTRAN format 
+             with a new file:
              $LMDGCM/libf/grid/dimension/dimensions.im.jm.lm.tntrac
-             Si ce fichier n'existe pas encore, il est cree par le script
+             If the file does not exist already
+             it is created by the script
              $LMDGCM/libf/grid/dimension/makdim
 
--p PHYS    pour selectionner le jeu de parametrisations physiques avec
-           lequel on veut compiler le modele.
-           Le modele sera alors compile en prenant les sources des
-           parametrisations physiques dans le repertoire:
+-s nscat   Number of radiatively active scatterers
+
+-p PHYS    Selects the set of physical parameterizations
+           you want to compile the model with.
+           The model is then compiled using the physical
+           parameterization sources in directory:
             $LMDGCM/libf/phyPHYS
 
--g grille  selectionne le type de grille qu'on veut utiliser.
-           L'effet de cette option est d'ecraser le fichier
+-g grille  Selects the grid type.
+           This option overwrites file
            $LMDGCM/libf/grid/fxyprim.h
-           avec le fichier
+           with file
            $LMDGCM/libf/grid/fxy_grille.h
-           grille peut prendre les valeurs:
-           1. reg pour la grille reguliere
-           2. sin pour avoir des points equidistants en sinus de la latitude
-           3. new pour pouvoir zoomer sur une partie du globe
-
--O "optimisation fortran" ou les optimisations fortran sont les options de la
-            commande f77
+           the grid can take the following values:
+           1. reg - the regular grid
+           2. sin - to obtain equidistant points in terms of sin(latitude)
+           3. new - to zoom into a part of the globe
+
+-O "compilation options" set of fortran compilation options to use
 
 -include path
-           Dans le cas ou on a dans des sous programmes des fichiers 
-           #include (cpp) qui se trouve sur des repertoires non references
-           par defaut
-
--adjnt     Pour compiler la l'adjoint du code dynamique
+           Used if the subroutines contain #include files (ccp) that 
+           are located in directories that are not referenced by default.
+
+-adjnt     Compiles the adjoint model to the dynamical code.
 
 -olddyn    To compile GCM with "old dynamics"
 
--filtre  filtre
-           Pour choisir le filtre en longitude dans les regions polaires.
-           "filtre" correspond au nom d'un repertoire se trouvant sur
-           $LMDGCM/libf. Le filtre standard du modele est "filtrez" qui peut
-           etre utilise aussi bien pour une grille reguliere que pour une 
-           grille zoomee en longitude.
+-filtre  filter
+           To select the longitudinal filter in the polar regions.
+           "filter" corresponds to the name of a directory located in
+           $LMDGCM/libf. The standard filter for the model is "filtrez"
+           which can be used for a regular grid and for a  
+           grid with longitudinal zoom.
 
 -link "-Ldir1 -lfile1 -Ldir2 -lfile2 ..."
-          Pour rajouter un lien avec les librairies FORTRAN
-          libfile1.a, libfile2.a ... se trouvant respectivement sur les
-          repertoires dir1, dir2 ...
-          Si dirn est un repertoire dont le chemin est automatique (comme
-          par exemple /usr/lib ...) il n'est pas besoin de specifier -Ldirn.
-
-Auteur: Frederic Hourdin  (hourdin@lmd.jussieu.fr)
+           Adds a link to FORTRAN libraries
+           libfile1.a, libfile2.a ... 
+           located in directories dir1, dir2 ...respectively
+           If dirn is a directory with an automatic path 
+           (/usr/lib ... for example) 
+           there is no need to specify  -Ldirn.
+
 eod
 exit
@@ -346 +347 @@
         set physique="$2" ; set phys="PHYS=$physique" ; shift ; shift ; goto top
 
+     case -s
+        set scatterers="$2" ; shift ; goto top
+
      case -g
         set grille="$2" ; shift ; shift ; goto top
 
+    case -b
+        set bands=$2 ; shift ; shift ; goto top
+           
      case -t
         set ntrac=$2 ; shift ; shift ; goto top
@@ -377 +384 @@
            set optim90="$optim90"" -G1 "
         else if $LINUX then
-           set optim90="-g -Mbounds "
+           ## for g95
+           set optim="-g -Wall -ftrace=full -fbounds-check"
+           set optim90="-g -Wall -ftrace=full -fbounds-check"
+           set optimtru90="-g -Wall -ftrace=full -fbounds-check"
         else 
            echo "pas d option debug predefinie pour cette machine"
@@ -397 +407 @@
 if ( "$physique" == 'nophys' ) then
    set phys="L_PHY= LIBPHY="
-endif
-
-########################################################################
-# choix du nombre de traceur par defaut si il n'a pas ete choisi,
-# suivant la physique
-########################################################################
-
-if ( $ntrac == 0  ) then
-    if ( "$physique" == 'nophys' ) then
-        set ntrac=1
-    else if ( "$physique" == 'lmd' ) then
-        set ntrac=2
-    else if ( "$physique" == 'lmd_test_li' ) then
-        set ntrac=2
-    else if ( "$physique" == 'ec' ) then
-        set ntrac=1
-    else
-        set ntrac = 1
-    endif
 endif
 
@@ -439 +430 @@
 
 # dimension
-
 echo "dimension avant sed $dim"
 if ( $IBM ) then
@@ -448 +438 @@
    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'`
      set dim=`echo $dim | sed -e 's/[^0-9]/ /g'`
      set dim_=`echo $dim | sed -e 's/[^0-9]/_/g'`
 endif
 
+# bands
+echo "bands avant sed $bands"
+if ( $IBM ) then
+   set bands=`echo $bands | sed -en 's/[^0-9]/ /g'`
+   set bands_=`echo $bands | sed -en 's/[^0-9]/_/g'`
+else if ( $SUN || $XNEC ) then
+   set bands=`echo $bands | sed -e 's/[^0-9]/ /g'` 
+   set bands_=`echo $bands | sed -e 's/[^0-9]/_/g'`
+else
+     set bands=`echo $bands | sed -e 's/[^0-9]/ /g'`
+     set bands_=`echo $bands | sed -e 's/[^0-9]/_/g'`
+endif
+
+
 # build final name of libraries directory:
 if ( "$dyntype" == "olddyn" ) then
-  set nomlib=${nomlib}_${dim_}_t${ntrac}_${grille}_olddyn
-else
-  set nomlib=${nomlib}_${dim_}_t${ntrac}_${grille}
-endif
+  set nomlib=${nomlib}_${dim_}_t${ntrac}_s${scatterers}_${grille}_olddyn
+else
+  set nomlib=${nomlib}_${dim_}_t${ntrac}_s${scatterers}_${grille}
+endif
+
+# Append number of bands to nomlib (new for universal model)
+set nomlib=${nomlib}_${bands_}
 
 # Append 'physique' type, if it is not mars, to nomlib
@@ -504 +509 @@
 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
-    #remove old dimensions.h file
-    \rm $libf/grid/dimensions.h
+  echo "WARNING: you are already compiling the model somewhere else"
+  echo "Wait until the first compilation is finished before starting."
+  echo "If you are sure that you are not compiling elsewhere, you can"
+  echo "type [yes] to continue."
+  echo "Do you want to continue?"
+  if ( $< == "yes" ) then
+    #remove old dimensions.h , bands.h and scatterers.h files
+    \rm -f $libf/grid/dimensions.h
+    \rm -f $libf/grid/bands.h
+    \rm -f $libf/phy${physique}/scatterers.h
   else
     exit
@@ -522 +529 @@
 # echo contents of dimensions.h to standard output
 cat $libf/grid/dimensions.h
+
+# Build the appropriate 'bands.h' file
+makbands $bands
+# echo contents of bands.h to standard output
+cat $libf/grid/bands.h
+
+# Build the appropriate 'scatterers.h' file
+cd $libf/phy$physique/scatterers
+./make_scatterers $scatterers
+# echo contents of scatterers.h to standard output
+cat $libf/phy$physique/scatterers.h
 
 cd $LMDGCM
@@ -549 +567 @@
   \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 $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
 
@@ -680 +680 @@
 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 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
-
-#
-# Encore un test temporaire: probleme de compilation sur VPP
-# l'optimisation de startvar se passe mal
-#
-if ($code == 'create_etat0' && $VPP) 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 optim90="$optim90"" -X9 -w"
-   set opt_link="-L$NCDFLIB -lnetcdf"
-   touch $LMDGCM/libf/dyn3d/startvar.F
-   touch $LMDGCM/libf/dyn3d/etat0_netcdf.F
-endif
-
 #
 # etat0_netcdf a besoin d'info de la physique
@@ -815 +775 @@
 PROG=$code
 
-\rm $libf/grid/dimensions.h
+\rm -f $libf/grid/dimensions.h
+\rm -f $libf/grid/bands.h
+\rm -f $libf/phy$physique/scatterers.h