Changeset 148 for trunk/LMDZ.MARS/makegcm

Timestamp:

Jun 8, 2011, 3:26:57 PM (14 years ago)

Author:

emillour

Message:

Mars GCM:

Significant update on how the number of scatterers is managed:

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

also removed some superfluous files (tmp, tmp77, ...).

EM

File:

• trunk/LMDZ.MARS/makegcm (modified) (8 diffs)

trunk/LMDZ.MARS/makegcm

@@ -13 +13 @@
 set grille=reg
 set dyntype="dyn"
+set scatterers="1"
 ########################################################################
 # path a changer contenant les sources et les objets du modele
@@ -290 +291 @@
              $LMDGCM/libf/grid/dimension/makdim
 
+-s nscat   Number of radiatively active scatterers
+
 -p PHYS    Selects the set of physical parameterizations
            you want to compile the model with.
@@ -346 +349 @@
      case -p
         set physique="$2" ; set phys="PHYS=$physique" ; shift ; shift ; goto top
+
+     case -s
+        set scatterers="$2" ; shift ; goto top
 
      case -g
@@ -459 +465 @@
 # build final name of libraries directory:
 if ( "$dyntype" == "olddyn" ) then
-  set nomlib=${nomlib}_${dim_}_t${ntrac}_${grille}_olddyn
+  set nomlib=${nomlib}_${dim_}_t${ntrac}_s${scatterers}_${grille}_olddyn
 else
-  set nomlib=${nomlib}_${dim_}_t${ntrac}_${grille}
+  set nomlib=${nomlib}_${dim_}_t${ntrac}_s${scatterers}_${grille}
 endif
 
@@ -507 +513 @@
 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
+  echo "WARNING: you are probably already compiling the model (perhaps"
+  echo " elsewhere). Wait until the first compilation ends before "
+  echo " compiling a different configuration."
+  echo "If you are sure that you are not already compiling, then you"
+  echo " may continue this compilation by answering yes."
+  echo "Do you want to continue?"
+  if ( $< == "yes" ) then
     #remove old dimensions.h file
     \rm -f $libf/grid/dimensions.h
+    #remove old scatterers.h file
+    \rm -f $libf/phymars/scatterers.h
   else
     exit
@@ -525 +534 @@
 # echo contents of dimensions.h to standard output
 cat $libf/grid/dimensions.h
+
+cd $libf/phymars/scatterers
+# Build the appropriate 'scatterers.h' file
+./make_scatterers $scatterers
+# echo contents of scatterers.h to standard output
+cat $libf/phymars/scatterers.h
 
 cd $LMDGCM
@@ -681 +696 @@
 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
-# A revoir
 set include="$include"" -I$libf/phy$physique"
 
@@ -812 +784 @@
 PROG=$code
 
-\rm $libf/grid/dimensions.h
+# cleanup, remove dimensions.h and scatterers.h
+\rm -f $libf/grid/dimensions.h
+\rm -f $libf/phymars/scatterers.h