Changeset 1146 for LMDZ4/trunk/makegcm

Timestamp:

Apr 9, 2009, 12:11:35 PM (15 years ago)

Author:

Laurent Fairhead

Message:

Réintegration dans le tronc des modifications issues de la branche LMDZ-dev
comprises entre la révision 1074 et 1145
Validation: une simulation de 1 jour en séquentiel sur PC donne les mêmes
résultats entre la trunk et la dev
LF

Location:

LMDZ4/trunk

Files:

• . (modified) (1 prop)
• makegcm (modified) (16 diffs)

LMDZ4/trunk/makegcm

@@ -11 +11 @@
 set phys="PHYS=$physique"
 set include='-I$(LIBF)/grid -I$(LIBF)/bibio -I$(LIBF)/filtrez -I. '
-set ntrac = 4
 set filtre=filtrez
 set grille=reg
@@ -38 +37 @@
 #setenv LIBOGCM $LMDGCM/libo
 set INCALIB=../INCA3/config/lib
-#set LMDGCM="`pwd`"
-#setenv LIBOGCM $LMDGCM/libo
+set LMDGCM="`pwd`"
+setenv LIBOGCM $LMDGCM/libo
 #set LMDGCM=/d4/fairhead/V4/
 #setenv LIBOGCM $LMDGCM/libo
@@ -48 +47 @@
 #setenv NCDFINC /distrib/local/netcdf/pgi_32bits/include
 #setenv NCDFLIB /distrib/local/netcdf/pgi_32bits/lib/
-#setenv IOIPSLDIR /d4/fairhead/Install_LMDZ4_FH/LMDZ4preV4/modipsl/lib/
-#setenv MODIPSLDIR /d4/fairhead/Install_LMDZ4_FH/LMDZ4preV4/modipsl/lib/
-#setenv NCDFINC /d4/fairhead/Install_LMDZ4_FH/LMDZ4preV4/netcdf-3.6.1/include
-#setenv NCDFLIB /d4/fairhead/Install_LMDZ4_FH/LMDZ4preV4/netcdf-3.6.1/lib
+setenv IOIPSLDIR /d4/fairhead/Install_LMDZ4_FH/LMDZ4preV4/modipsl/lib/
+setenv MODIPSLDIR /d4/fairhead/Install_LMDZ4_FH/LMDZ4preV4/modipsl/lib/
+setenv NCDFINC /d4/fairhead/Install_LMDZ4_FH/LMDZ4preV4/netcdf-3.6.1/include
+setenv NCDFLIB /d4/fairhead/Install_LMDZ4_FH/LMDZ4preV4/netcdf-3.6.1/lib
 
 
@@ -207 +206 @@
    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 "
+   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=" "
@@ -217 +220 @@
    set optim90=" -fast -fixed "
    set optimtru90=" -fast -free "
-   set opt_link="-lf77compat -L$MODIPSLDIR -lsechiba -lparameters -lstomate -lioipsl -L$NCDFLIB -lnetcdf "
+   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
@@ -271 +278 @@
    set optim="$optim90"
    set optimbis=" "
-   set opt_link=" -C hopt -float0 -ew -P static -L$MODIPSLDIR -lioipsl  $NCDFLIB -lnetcdf_i8r8_v "
+   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"
@@ -365 +376 @@
 -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 
+             L'effet des options -d est d'ecraser le fichier 
              $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
-             $LMDGCM/libf/grid/dimension/dimensions.im.jm.lm.tntrac
+             qui contient sous forme de 3 PARAMETER FORTRAN les 3 dimensions
+             de la grille horizontale im, jm et verticale lm, par un nouveu fichier
+             $LMDGCM/libf/grid/dimension/dimensions.im.jm.lm
              Si ce fichier n'existe pas encore, il est cree par le script
              $LMDGCM/libf/grid/dimension/makdim
@@ -472 +477 @@
         set parallel="$2" ; shift ; shift ; goto top
   
-     case -t
-        set ntrac=$2 ; shift ; shift ; goto top
-
      case -include
         set include="$include -I$2" ; shift ; shift ; goto top
@@ -543 +545 @@
    set phys="L_PHY= LIBPHY="
 else
-   set cppflags="$cppflags -DCPP_PHYS"
+   #Default planet type is Earth
+   set cppflags="$cppflags -DCPP_EARTH"
 endif
 
@@ -581 +584 @@
 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
 
@@ -631 +615 @@
    set dim=`echo $dim | sed -e 's/[^0-9]/ /g'`
 endif
-set nomlib=${nomlib}${physique}_${dim_}_t${ntrac}_$grille
+set nomlib=${nomlib}${physique}_${dim_}_$grille
 ## M-A-F nomlib trop long sur CRAY pour ar
 if ( $CRAY ) then
-    set nomlib=F90_${dim_}_t${ntrac}
+    set nomlib=F90_${dim_}
 endif
 if ( $NEC || $XNEC || $X6NEC || $X8BRODIE ) then
-    set nomlib=F90_${dim_}_t${ntrac}_'phy'${physique}${FLAG_PARA}
+    set nomlib=F90_${dim_}_'phy'${physique}${FLAG_PARA}
 endif
 echo calcul de la dimension
@@ -664 +648 @@
 ########################################################################
 # Gestion des dimensions du modele.
-# on cree ou remplace le fichier des dimensions/nombre de traceur
+# on cree ou remplace le fichier des dimensions
 ########################################################################
 
@@ -682 +666 @@
 
 cd dimension
-./makdim $ntrac $dim
+./makdim $dim
 cat $libf/grid/dimensions.h
 
@@ -714 +698 @@
 ######################################################################
 
-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
+#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
 
 ######################################################################
@@ -822 +806 @@
    set f90=$FC_LINUX
    if ( $FC_LINUX == 'pgf90' ) then
-   set opt_link=" -L$MODIPSLDIR $link_veget -L$NCDFLIB -lioipsl -lnetcdf "
+     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') then
-     set opt_link="-L$MODIPSLDIR $link_veget -lioipsl -L$NCDFLIB -lnetcdf -lioipsl -lnetcdf "
+     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=" "
@@ -833 +825 @@
    set f77=f90
    set f90=f90
-   set opt_link="-lf77compat -L$MODIPSLDIR $link_veget -lioipsl -L$NCDFLIB -lnetcdf "
+   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