source: trunk/MESOSCALE/LMD_MM_MARS/makemeso @ 1937

Last change on this file since 1937 was 1937, checked in by aslmd, 7 years ago

MESOSCALE: makemeso: made compil for LES compatible with ciclad env and Occigen env i.e. netcdff is put forward by configure in Occigen but not in ciclad

  • Property svn:executable set to *
File size: 38.9 KB
Line 
1#! /bin/bash
2
3
4##################################
5# makemeso allows you to compile #
6# WRF with modified LMD physics  #
7##################################
8
9##############################
10# Author: A. Spiga           #
11# New version : October 2008 #
12# Last update : January 2009 #
13#               November  09 #
14#               January   11 #
15##############################
16
17###############################
18#### Type makemeso -h for help
19###############################
20
21
22
23##############################################################################################
24##############################################################################################
25##############################################################################################
26##############################################################################################
27
28donotcompile=0
29donotallow=0
30config='real'
31donotcompilephys=0
32justphys=0
33debug=0
34fresh_start=0
35phys="mars_lmd"
36scenario=""
37from_scratch=0
38
39while getopts "drc:njhgp:fs:xe" options; do
40  case $options in
41   d ) donotcompile=1;;       ## just to check the compile folder
42   r ) donotallow=1;;         ## allow only known config (useful with 'makemeso < last')
43   c ) config="${OPTARG}";;   ## idealized cases
44   n ) donotcompilephys=1;;   ## do not recompile physics
45   j ) justphys=1;;           ## just compile LMD physics
46   g ) debug=1;;              ## debug mode
47   p ) phys="${OPTARG}";;     ## choose physics
48   f ) fresh_start=1;;        ## a fresh start
49   e ) from_scratch=1;;       ## a fresh start with a completely new folder
50   s ) scenario="${OPTARG}";; ## a specific scenario (with precompiling flags)
51   x ) donotcompile=1;phys="void_lmd_new";donotcompilephys=1;; ## a case with no LMD physics included
52   h ) echo "
53# Use:   
54#
55# makemeso                   ## basic use (real-case configuration)
56#
57# makemeso -d                ## no compilation, just check the name of the compile folder
58#
59# makemeso -c ideal          ## idealized mode (convective cell, mountain wave, etc...)
60# makemeso -c les            ## large-eddy simulations mode based on WRFV3
61#
62# makemeso -n                ## do not recompile LMD physics (must have been compiled before)
63#
64# makemeso < last            ## basic use + skip questions [! script must have been executed at least once]
65# makemeso -r < last         ## basic use + skip questions + only known config
66# makemeso -nr < last        ## basic use + skip questions + only known config + no LMD phys recompile
67#
68# makemeso -j                ## just compile the LMD physics
69#
70# makemeso -g                ## debug mode
71#
72# makemeso -h                ## display options
73#
74# makemeso -p mars_lmd_new   ## choose physics
75#             venus_lmd_new
76#
77# makemeso -f                ## fresh start [clean -a]
78#
79# makemeso -e                ## a completely new recompile with erasing the WHOLE compiling folder
80#
81# makemeso -s DUSTSTORM      ## a specific scenario (with precompiling flags)
82#
83# makemeso -x                ## a case with no LMD physics included
84" ; exit ;;
85  esac
86done   
87
88#--------------
89# talk w/ user
90#--------------
91  echo '****************************************'
92  echo ' LMD Mesoscale Model Compiler. Welcome.'
93  echo '****************************************'
94  if [[ "${phys}" == *"new"* ]]
95  then
96    echo '*********** PHYSICS: ' ${phys}
97    echo '****************************************'
98  fi
99  # computer
100  uname -a | grep x86_64 > /dev/null
101  if [[ "$?" == 0 ]]
102  then
103    machine='64'
104  else
105    machine='32'
106  fi
107  # compiler
108  echo "Supported compiler options are "
109  echo "       <1> pgf90"
110  echo "       <2> g95"
111  echo "       <3> pgf90 + mpi"
112  echo "       <4> ifort"
113  echo "       <5> ifort + mpi"
114  echo "       <7> xlf + mpi (IBM AIX)"
115  echo "       <8> gfortran"
116  echo "Your choice ?" ; read reply
117  case ${reply} in
118    1) compilo='pgf' ; numproc=1 ;;
119    2) compilo='g95' ; numproc=1 ;;
120    3) compilo='mpi'
121       if [[ "${WHERE_MPI}" = "" ]]
122       then
123         echo Please initialize the variable WHERE_MPI in your environnement
124         exit
125       fi ;;
126    4) compilo='ifort' ; numproc=1 ;;
127    5) compilo='mpifort' ;;
128    7) compilo='mpixlf' ;;
129    8) compilo='gnu' ; numproc=1 ;;
130    *) echo not supported by this script ; exit ;;
131  esac
132  # number of processors
133  if [[ "${phys}" == *"new"* ]]
134  then
135   numproc=999
136  else
137   case ${reply} in
138    3) echo How many processors ? e.g. 1, 2, 4, 6, 8, 12, 16, 20, 24, 32, 64, 128 ; read numproc ;;
139    5) echo How many processors ? e.g. 1, 2, 4, 6, 8, 12, 16, 20, 24, 32, 64, 128 ; read numproc ;;
140    7) echo How many processors ? e.g. 1, 2, 4, 6, 8, 12, 16, 20, 24, 32, 64, 128 ; read numproc ;;
141   esac
142  fi
143
144  # dimensions
145  if [ ${donotcompile} -eq 0 ]
146  then   
147  # ngrid: now dynamically set in newphys 03/2014
148   if [[ "${phys}" == *"new"* ]]
149   then
150      lon=999
151      lat=999
152      level=999
153   else
154      echo Grid points in longitude ? ; read lon
155      echo Grid points in latitude ? ; read lat
156      echo Number of vertical levels ? ; read level
157   fi
158  fi
159  echo Number of domains ? ; read dom
160
161###PB lecture dom si < last
162
163    case ${dom} in
164       1) single='_single' ;;
165       *) single='_nest' ;;
166    esac
167      if [ "${compilo}" = "mpi" -o "${compilo}" = "mpifort" -o "${compilo}" = "mpixlf" ]
168      then
169           single=''
170      fi
171
172testflag=''
173####
174#testflag='_test'
175####
176
177    conf_wrf="${phys}_${config}${scenario}_${compilo}_${machine}${single}${testflag}"
178
179    if [ ${debug} -eq 1 ]
180    then
181     conf_wrf="debug_${conf_wrf}"
182    fi
183
184    \rm what_folder 2> /dev/null
185    echo ${conf_wrf} > what_folder
186         if [ ${donotcompile} -eq 0 ]
187         then
188            ##\rm what_folder 2> /dev/null
189            ##echo '**********************'
190            ##echo '*** Your folder is ...'
191            ##echo '**********************'
192            ##echo ${conf_wrf} | tee what_folder
193            ##cat what_folder
194            ##echo ${conf_wrf} > what_folder
195            ##echo ${reply} > what_compilo
196            ##echo ${numproc} > what_numproc
197            #if [[ "${phys}" != "nophys_" ]]
198            #then
199            #  exit
200            #fi
201         #else
202            # tracers: now dynamically set in newphys 09/2013
203            if [[ "${phys}" == *"new"* ]]
204            then
205              tra=999
206            else
207              echo Number of tracers ? ; read tra
208              if [ ${tra} -eq 0 ]
209              then
210                tra=1   
211              fi
212            fi
213         fi
214  # 'from scratch' case
215  if [ ${from_scratch} -eq 1 ]
216  then
217      echo "***** I ERASE THE FOLDER "${conf_wrf}
218      \rm -rf ${conf_wrf}
219  fi
220  # folder
221
222  mkdir ${conf_wrf} 2> /dev/null
223  if [[ "$?" == 0 ]]
224  then
225    echo new folder ... link sources
226    ###############################################
227     echo ${phys}
228     ./SRC/SCRIPTS/copy_model -p ${phys} -c ${config}
229     mv zeWRFV2 ${conf_wrf}/WRFV2
230    ###############################################
231     if [[ "${phys}" == *"new"* ]]
232     then
233      cd ${conf_wrf}/WRFV2
234        if [[ "${config}" != "les" ]]
235        then
236           cp ../../SRC/WRFV2/Registry/Registry.EM.newphys Registry/Registry.EM
237        else
238           cp -L ../../SRC/LES/WRFV2/Registry/Registry.EM.newphys Registry/Registry.EM
239        fi
240      cd phys
241      ln -sf module_lmd_driver.F.new module_lmd_driver.F
242      cd ..
243
244     else
245      cd ${conf_wrf}/WRFV2
246      cd Registry ; ./Registry.bash ; cd ..
247     fi
248          #### sparadrap consequent a l'utilisation de copy_model pour les liens
249          #### -- car alors il manque fftpack
250          if [[ "${config}" == "les" ]]
251          then
252            cp ../../SRC/LES/correcfft ./
253            ./correcfft
254            \rm correcfft
255          fi
256  else
257    cd ${conf_wrf}/WRFV2
258    if [ ${fresh_start} -eq 1 ]
259    then
260      echo '*** FRESH START, I clean everything'
261      ### clean the whole place 
262      clean -a > /dev/null 2> /dev/null
263      ### get a possibly modified registry !
264      if [[ "${phys}" == *"new"* ]]
265      then
266          cp ../../SRC/WRFV2/Registry/Registry.EM.newphys Registry/Registry.EM
267      else
268          cp ../../SRC/WRFV2/Registry/Registry.EM Registry/
269      fi
270      if [[ "${config}" == "les" ]]
271      then
272        cp ../../SRC/LES/WRFV2/Registry/Registry.EM Registry/
273      fi
274      ### ensure the model will be recompiled from scratch with correct registry
275      cd Registry ; mv Registry Registry.bak ; Registry.bash ; cd ..
276#     answer='y' ## a voir... pas si sur
277#    else
278#            if [[ "${phys}" != "nophys_" ]]
279#            then
280#      echo Did you modify anything in the Registry or clean ? y for yes, any key for no ; read answer
281#      case ${answer} in
282#        y) cd Registry ; mv Registry Registry.bak ; Registry.bash ; cd .. ;;
283#        *) answer='no' ;;
284#      esac
285#            fi
286    fi
287  fi
288
289  # summary
290  echo '**********************'
291  echo '*** Your folder is ...'
292  echo '**********************'
293  echo ${conf_wrf}
294  echo '**********************'
295  echo '****************************************'
296  echo so ...
297  echo your computer is ${machine} bits
298  echo ${compilo} is your compiler
299       if [[ "${compilo}" = "mpi" ]]
300       then
301         echo MPICH is in ${WHERE_MPI}
302       fi
303  echo you have ${dom} domains
304  if [[ "${phys}" == *"new"* ]]
305  then
306   echo number of processors can be changed without recompiling
307   echo x points can be changed without recompiling
308   echo y points can be changed without recompiling
309   echo z points can be changed without recompiling
310   echo number of tracers can be changed without recompiling
311  else
312   echo you will use ${numproc} processors
313   echo you have ${lon} x points
314   echo '         '${lat} y points
315   echo '         '${level} z points
316   echo '         '${tra} tracers
317  fi
318  echo '****************************************'
319  # save answer
320  \rm last 2> /dev/null
321  touch last
322  echo ${reply} >> last
323  if [[ "${phys}" != *"new"* ]]
324  then
325    if [ "${compilo}" = "mpi" -o "${compilo}" = "mpifort" -o "${compilo}" = "mpixlf" ]
326    then
327     echo ${numproc} >> last
328    fi
329    echo ${lon} >> last
330    echo ${lat} >> last
331    echo ${level} >> last
332  fi
333  echo ${dom} >> last
334  if [[ "${phys}" != **"new"* ]]
335  then
336   echo ${tra} >> last
337  fi
338
339#------------
340# log files
341#------------
342  \rm log_compile 2> /dev/null
343  \rm log_error 2> /dev/null
344
345#-----------------
346# configure WRF
347#-----------------
348  mv configure.wrf configure.wrf.bak 2> /dev/null
349  ########################
350  conf_wrf_forall="${compilo}_${machine}${single}${testflag}"
351
352         #### in this case options are different because for LES, WRFV3 is used
353         if [[ "${config}" == "les" ]]
354         then
355         conf_wrf_forall="${config}${compilo}_${machine}${single}${testflag}"
356         fi
357
358  case ${conf_wrf_forall} in
359
360     #######LES-specific (WRFV3-based) -- previously in 'makeles'
361     #######
362
363     lesmpi_64)        ## MPI (dm) 64 bits [PS: remplacer 3 par 4 pour openMP]
364                       \rm conf > /dev/null 2> /dev/null ; touch conf ; echo 3 >> conf ; echo 1 >> conf
365                       \rm configure.wrf > /dev/null 2> /dev/null ; ./configure < conf > /dev/null 2> /dev/null
366                       sed s+"-lnetcdf"+"-lnetcdf -L../$phys/libo -llmd"+g   configure.wrf > yeah ; mv -f yeah configure.wrf
367                       if [[ "$(hostname)" == "ciclad1.ipsl.jussieu.fr" ]]
368                       then
369            echo "SPECIFIC CHANGES FOR CICLAD CLUSTER. EDIT makemeso IF YOU ENCOUNTER PROBLEMS."
370            #sed s+"mpif90 -f90=$(SFC)"+"/usr/lib64/openmpi/1.4.2-pgf/bin/mpif90"+g                      configure.wrf > yeah ; mv -f yeah configure.wrf
371            #sed s+"mpicc -cc=$(SCC)"+"/usr/lib64/openmpi/1.4.2-gfortran/bin/mpicc -DMPI2_SUPPORT"+g     configure.wrf > yeah ; mv -f yeah configure.wrf
372            sed s+"mpif90 -f90=\$(SFC)"+"/usr/lib64/openmpi/1.4.3-pgfgcc/bin/mpif90"+g                   configure.wrf > yeah ; mv -f yeah configure.wrf
373            sed s+"mpicc -cc=\$(SCC)"+"/usr/lib64/openmpi/1.4.3-pgfgcc/bin/mpicc -DMPI2_SUPPORT"+g       configure.wrf > yeah ; mv -f yeah configure.wrf
374                       else
375            #sed s+"mpif90"+"$WHERE_MPI/mpif90"+g                 configure.wrf > yeah ; mv -f yeah configure.wrf
376            #sed s+"mpicc"+"$WHERE_MPI/mpicc -DMPI2_SUPPORT"+g    configure.wrf > yeah ; mv -f yeah configure.wrf
377            #### OPENMPI
378            sed s+"mpif90 -f90=\$(SFC)"+"$WHERE_MPI/mpif90"+g                 configure.wrf > yeah ; mv -f yeah configure.wrf
379            sed s+"mpicc -cc=\$(SCC)"+"$WHERE_MPI/mpicc -DMPI2_SUPPORT"+g    configure.wrf > yeah ; mv -f yeah configure.wrf
380                       fi
381                       sed s+"-fastsse"+" "+g                               configure.wrf > yeah ; mv -f yeah configure.wrf
382                          ### pas forcement necessaire ici mais OK
383                          sed s+"-llmd"+"-llmd $NETCDF/lib/libnetcdf.a"+g     configure.wrf > yeah ; mv -f yeah configure.wrf
384                       \rm conf > /dev/null ;;
385
386     lesmpifort_64)    ## MPI (dm) 64 bits IFORT
387                       \rm conf > /dev/null 2> /dev/null ; touch conf ; echo 7 >> conf ; echo 1 >> conf
388                       \rm configure.wrf > /dev/null 2> /dev/null ; ./configure < conf > /dev/null 2> /dev/null
389                       if [[ "${phys}" != *"void"* ]]
390                       then
391                         sed s+"-lnetcdff"+"TOZTOZ"+g configure.wrf > yeah ; mv -f yeah configure.wrf
392                         # ... in case there is netcdff save it in TOZTOZ so that it is not impacted in the following sed
393                         sed s+"-lnetcdf"+"-lnetcdf -L../$phys/libo -llmd $NETCDF/lib/libnetcdf.a"+g   configure.wrf > yeah ; mv -f yeah configure.wrf
394                         sed s+"TOZTOZ"+"-lnetcdff"+g configure.wrf > yeah ; mv -f yeah configure.wrf
395                       else
396                         sed s+"-lnetcdf"+"-lnetcdf $NETCDF/lib/libnetcdf.a"+g   configure.wrf > yeah ; mv -f yeah configure.wrf
397                       fi
398                       sed s+"mpif90"+"$WHERE_MPI/mpif90"+g                 configure.wrf > yeah ; mv -f yeah configure.wrf
399                       sed s+"mpicc"+"$WHERE_MPI/mpicc -DMPI2_SUPPORT"+g    configure.wrf > yeah ; mv -f yeah configure.wrf
400                       sed s+"-fastsse"+" "+g                               configure.wrf > yeah ; mv -f yeah configure.wrf
401                       sed s+"O3"+"O3 -mcmodel=large -shared-intel"+g configure.wrf > yeah ; mv -f yeah configure.wrf
402                       #sed s+"-lnetcdf -lnetcdff"+"-lnetcdf -lnetcdff -L../mars_lmd/libo -llmd"+g     configure.wrf > yeah ; mv -f yeah configure.wrf
403                       if [ ${debug} -ne 0 ]   # not working for xlf!
404                             then
405                             echo 'DEBUG DEBUG DEBUG DEBUG'
406                             sed s+"# -g"+" -g -no-ftz -traceback -ftrapuv -fp-stack-check -check"+g configure.wrf > yeah
407                             mv -f yeah configure.wrf
408                       fi
409                       \rm conf > /dev/null ;;
410
411     #######
412     #######LES-specific (WRFV3-based) -- previously in 'makeles'   
413
414                     # GFORTRAN, 64 bits, no nesting     
415     gnu_64_single)    cd arch ; ln -sf ../configure.defaults.gfortran64 configure.defaults ; cd ..
416                       echo 15 | configure > log_compile 2> log_error
417                       sed -f physics.sed configure.wrf > yeah ; mv -f yeah configure.wrf 
418                       sed s+"-L../$phys/libo -llmd"+"-L../$phys/libo -llmd $NETCDF/lib/libnetcdf.a"+g configure.wrf > yeah ; mv -f yeah configure.wrf ;;
419                     # PGF90, 32 bits, no nesting
420     pgf_32_single)    echo 1 | configure > log_compile 2> log_error
421                       sed -f physics.sed configure.wrf > yeah ; mv -f yeah configure.wrf ;;
422                     # PGF90, 32 bits, nesting
423     pgf_32_nest)      echo 2 | configure > log_compile 2> log_error
424                       sed -f physics.sed configure.wrf > yeah ; mv -f yeah configure.wrf ;;
425                     # PGF90, 64 bits, no nesting
426     pgf_64_single)    echo 1 | configure > log_compile 2> log_error
427                       #### the following line is necessary because readtesassim has a line "use netcdf" in it!
428                       #sed s+"-L../$phys/libo -llmd"+"-L../$phys/libo -llmd $NETCDF/lib/libnetcdf.a"+g configure.wrf > yeah ; mv -f yeah configure.wrf
429                       sed -f physics.sed configure.wrf > yeah ; mv -f yeah configure.wrf ;;
430                     # PGF90, 64 bits, nesting
431     pgf_64_nest)      cd arch ; ln -sf ../configure.defaults.fix64nest configure.defaults ; cd ..
432                       echo 4 | configure > log_compile 2> log_error
433                       sed -f physics.sed configure.wrf > yeah ; mv -f yeah configure.wrf ;;
434                     # G95, 32 bits, no nesting
435     g95_32_single)    echo 13 | configure > log_compile 2> log_error
436                       sed -f physics.sed configure.wrf > yeah ; mv -f yeah configure.wrf ;;
437                     # G95, 64 bits, no nesting  ### we modify configure.defaults to add x86_64 next to g95 in the header comment
438     g95_64_single)    cd arch ; sed s/"PC Linux i486 i586 i686, g95 compiler"/"PC Linux i486 i586 i686 x86_64, g95 compiler"/g configure.defaults > yeahyeahyeah ; rm configure.defaults ; mv yeahyeahyeah configure.defaults ; cd ..
439                       echo 14 | configure > log_compile 2> log_error
440                         ##### marche pas ....
441                         #sed s+"-fno-second-underscore"+""+g configure.wrf > yeah ; mv -f yeah configure.wrf
442                         #sed -f physics.sed configure.wrf > yeah ; mv -f yeah configure.wrf
443                         #sed s+"-L../$phys/libo -llmd"+"-L../$phys/libo -llmd $NETCDF/lib/libnetcdf.a"+g configure.wrf > yeah ; mv -f yeah configure.wrf ;;
444                       sed -f physics.sed configure.wrf > yeah ; mv -f yeah configure.wrf ;;
445                       #sed s+"=       g95"+"=       /home/physastro/aspiga/mysoft/g95/g95-install/bin/x86_64-unknown-linux-gnu-g95"+g configure.wrf > yeah ; mv -f yeah configure.wrf;;
446                     # IFORT, 64 bits, no nesting
447     ifort_64_single)  echo 5 | configure > log_compile 2> log_error
448                       sed s+"-DIFORT_KLUDGE"+" "+g configure.wrf > yeah ; mv -f yeah configure.wrf
449                       if [ ${debug} -ne 0 ]   # not working for xlf!
450                             then
451                             echo 'DEBUG DEBUG DEBUG DEBUG'
452                             sed s+"# -g"+" -g -no-ftz -traceback -ftrapuv -fp-stack-check -check"+g configure.wrf > yeah
453                             mv -f yeah configure.wrf
454                       fi
455                         ## !!! NETCDF must be defined
456                         sed s+"lio_grib_share"+"lio_grib_share -L../$phys/libo -llmd $NETCDF/lib/libnetcdf.a"+g configure.wrf > yeah
457                         mv -f yeah configure.wrf ;;
458     mpifort_64)     # MPI+IFORT, 64 bits, no nesting / nesting
459                       echo 9 | ./configure > log_compile 2> log_error
460                       sed -f physics.sed configure.wrf > yeah ; mv -f yeah configure.wrf
461                       sed s+"-f90=ifort"+" "+g configure.wrf > yeah ; mv -f yeah configure.wrf
462                       sed s+"-cc=icc"+" "+g configure.wrf > yeah ; mv -f yeah configure.wrf
463                       sed s+"-DIFORT_KLUDGE"+" "+g configure.wrf > yeah ; mv -f yeah configure.wrf
464                       sed s+"-O1"+"-O1 -mcmodel=medium -shared-intel"+g configure.wrf > yeah ; mv -f yeah configure.wrf
465                       sed s+"-O2"+"-O2 -mcmodel=medium -shared-intel"+g configure.wrf > yeah ; mv -f yeah configure.wrf
466                       sed s+"-O3"+"-O3 -mcmodel=medium -shared-intel"+g configure.wrf > yeah ; mv -f yeah configure.wrf
467                       if [ ${debug} -ne 0 ]   # not working for xlf!
468                             then
469                             echo 'DEBUG DEBUG DEBUG DEBUG'
470                             sed s+"#-g"+" -g -no-ftz -traceback -ftrapuv -fp-stack-check -check"+g configure.wrf > yeah
471                             mv -f yeah configure.wrf
472                       fi
473                       sed s+"O3"+"O2"+g configure.wrf > yeah ; mv -f yeah configure.wrf
474                                        #sed s+"O3"+"O0"+g configure.wrf > yeah ; mv -f yeah configure.wrf  ### pour compilation rapide
475                                        #-w -ftz -align all -fno-alias -fp-model precise >>> options indiquees dans WRF3
476                                        #http://software.intel.com/en-us/articles/performance-tools-for-software-developers-building-wrf-with-the-intel-compilers/
477                            #### NB: RSL is replaced by RSL_LITE, better memory capacity, less seg fault
478                            sed s+"../external/RSL/RSL"+"../external/RSL_LITE"+g configure.wrf > yeah ; mv -f yeah configure.wrf
479                            sed s+"-DWRF_RSL_IO"+" "+g configure.wrf > yeah ; mv -f yeah configure.wrf
480                            sed s+"-DRSL"+" "+g configure.wrf > yeah ; mv -f yeah configure.wrf
481                            sed s+"librsl.a"+"librsl_lite.a"+g configure.wrf > yeah ; mv -f yeah configure.wrf
482                            sed s+"../external/RSL/gen"+"../external/RSL_LITE/gen"+g configure.wrf > yeah ; mv -f yeah configure.wrf
483                            sed s+"../external/RSL/module"+"../external/RSL_LITE/module"+g configure.wrf > yeah ; mv -f yeah configure.wrf
484                            sed s+"-lrsl"+"-lrsl_lite"+g configure.wrf > yeah ; mv -f yeah configure.wrf
485                            sed s+"linux"+" "+g configure.wrf > yeah ; mv -f yeah configure.wrf
486                         ### necessary even if mpi-selector is used (no need in makegcm though)
487                         sed s+"mpif90"+"${WHERE_MPI}/mpif90"+g configure.wrf | sed s+"mpicc"+"${WHERE_MPI}/mpicc"+g > yeah ; mv -f yeah configure.wrf
488                         sed s+"mpicc"+"mpicc -DMPI2_SUPPORT"+g configure.wrf > yeah ; mv -f yeah configure.wrf
489                         ## !!! NETCDF must be defined
490                         if [[ "${phys}" != *"void"* ]]
491                         then
492                           sed s+"-L../$phys/libo -llmd"+"-L../$phys/libo -llmd $NETCDF/lib/libnetcdf.a"+g configure.wrf > yeah
493                         else
494                           sed s+"-lesmf_time"+"-lesmf_time $NETCDF/lib/libnetcdf.a"+g configure.wrf > yeah
495                         fi
496                         mv -f yeah configure.wrf
497                         #### POUR LE TRAITEMENT PARTICULIERS des NESTS sur iDATAPLEX [cf. module_lmd_driver]
498                         sed s+"ARCHFLAGS       =       "+"ARCHFLAGS       =       -DSPECIAL_NEST_SAVE "+g configure.wrf > yeah
499                         mv -f yeah configure.wrf ;;
500                     # MPICH, 64 bits, no nesting / nesting
501     mpi_64)           cd arch ; ln -sf ../configure.defaults.fix64nest configure.defaults ; cd ..
502                       echo 3 | configure > log_compile 2> log_error
503                       sed -f physics.sed configure.wrf > yeah ; mv -f yeah configure.wrf
504                       ### the following line is necessary because readtesassim has a line "use netcdf" in it!
505                       sed s+"-L../$phys/libo -llmd"+"-L../$phys/libo -llmd $NETCDF/lib/libnetcdf.a"+g configure.wrf > yeah ; mv -f yeah configure.wrf
506                       sed s+"mpif90"+"${WHERE_MPI}/mpif90"+g configure.wrf | sed s+"mpicc"+"${WHERE_MPI}/mpicc"+g | sed s+"mpicc"+"mpicc -DMPI2_SUPPORT"+g > yeah
507                       #sed s+"mpif90"+"${WHERE_MPI}/mpif90"+g configure.wrf | sed s+"mpicc"+"${WHERE_MPI}/mpicc"+g | sed s+"mpicc"+"mpicc -DMPI2_SUPPORT -w -O3"+g > yeah
508                                #| sed s+"fastsse"+"fast"+g > yeah
509                       mv -f yeah configure.wrf ;;
510                       #    #sed s+"-fastsse"+"-O2 -Munroll -Mcache_align"+g configure.wrf > yeah
511                       #    ##sed s+"-fastsse"+"-O2 -fpic"+g configure.wrf > yeah ## marche pas, ILM problem
512                       #    ##sed s+"-fastsse"+"-mcmodel=medium -Mlarge_arrays"+g configure.wrf > yeah ## marche pas, ILM problem
513                       #sed s+"-fastsse"+"-O3 -tp=nehalem-64 -fastsse -Mscalarsse -Mvect=sse -Mipa=fast"+g configure.wrf > yeah
514                       #sed s+"-fastsse"+"-O2 -Munroll -Mcache_align"+g configure.wrf > yeah
515                       #sed s+"-fastsse"+""+g configure.wrf > yeah
516                       #sed s+"-fastsse"+"-fast"+g configure.wrf > yeah
517                       #sed s+"-fastsse"+"-Mvect=noaltcode -Msmartalloc -Mprefetch=distance:8 -Mfprelaxed"+g configure.wrf > yeah
518                       #mv -f yeah configure.wrf ;;
519                     # MPICH, 64 bits, OK with periodic BC but no positive definite
520     mpi_64_test)      cd arch ; ln -sf ../configure.defaults.fix64nest configure.defaults ; cd ..
521                       echo 2 | configure > log_compile 2> log_error
522                       sed -f physics.sed configure.wrf > yeah ; mv -f yeah configure.wrf
523                       sed s+"mpif90"+"${WHERE_MPI}/mpif90"+g configure.wrf | sed s+"mpicc"+"${WHERE_MPI}/mpicc"+g | sed s+"mpicc"+"mpicc -DMPI2_SUPPORT"+g > yeah
524                       mv -f yeah configure.wrf ;;
525                     # XLF+MPICH on AIX machine (64 bits) with nesting, NB: the 32bits notice is not to be considered
526     mpixlf_32)        cd arch ; ln -sf ../configure.defaults.fix64nest configure.defaults ; cd ..
527                       echo 3 | configure > log_compile 2> log_error
528                       sed -f physics.sed configure.wrf > yeah
529                       mv -f yeah configure.wrf ;;
530                     # ANYTHING ELSE
531     *)                echo NO PRESETS ...
532                       if [ ${donotallow} -eq 0 ] 
533                       then
534                         configure
535                       else
536                         # problem when an input file is used
537                         echo 'please cd to '$PWD' and type ./configure' 
538                         exit   
539                       fi ;;
540  esac
541  ########################
542
543
544  ### here a case structure would be great
545
546  if [[ "${phys}" == *"new"* ]]
547  then
548     if [[ "${config}" == "les" ]]  ### LES is different because of WRFV3
549     then
550       sed s+"ARCH_LOCAL      =       "+"ARCH_LOCAL      =       -DNEWPHYS "+g configure.wrf > yeah
551     else
552       sed s+"ARCHFLAGS       =       "+"ARCHFLAGS       =       -DNEWPHYS "+g configure.wrf > yeah
553     fi
554     mv -f yeah configure.wrf
555  fi
556
557        ### SCENARIO SCENARIO SCENARIO
558        if [[ ! ("${scenario}" == "") ]]   ### not supported with LES for the moment?
559        then
560          # precompiling flag in WRF part
561          sed s+"ARCHFLAGS       =       "+"ARCHFLAGS       =       -D${scenario} "+g configure.wrf > yeah
562          mv -f yeah configure.wrf
563          # precompiling flag in LMD part
564          cd $phys
565          cp -Lrf makegcm_${compilo} yeah
566          \rm makegcm_${compilo}
567          sed s+"-DMESOSCALE"+"-D${scenario} -DMESOSCALE"+ yeah > makegcm_${compilo}
568          chmod 755 makegcm_${compilo}
569          \rm yeah
570          cd ..
571        fi
572
573  ################
574  #if [[ "${phys}" == "nophys_" ]]
575  #then
576  #   if [[ "${config}" == "les" ]]  ### LES is different because of WRFV3
577  #   then
578  #     sed s+"ARCH_LOCAL      =       "+"ARCH_LOCAL      =       -DNOPHYS "+g configure.wrf > yeah ; mv -f yeah configure.wrf
579  #     sed s+"-L../$phys/libo -llmd"+""+g configure.wrf > yeah ; mv -f yeah configure.wrf
580  #   else
581  #     ### not tested yet but should be working 
582  #     echo CAUTION CAUTION CAUTION NOT FULLY TESTED
583  #     sed s+"ARCHFLAGS       =       "+"ARCHFLAGS       =       -DNOPHYS "+g configure.wrf > yeah ; mv -f yeah configure.wrf
584  #     sed s+"-L../$phys/libo -llmd"+""+g configure.wrf > yeah ; mv -f yeah configure.wrf
585  #   fi
586  #fi
587  ################
588
589  if [[ "${phys}" == *"venus"* ]]
590  then
591    sed s+"LIB_BUNDLED     = "+"LIB_BUNDLED     =    /usr/lib64/liblapack.so  \ "+g configure.wrf > yeah ; mv -f yeah configure.wrf
592  fi
593
594  if [ ${debug} -ne 0 ]   # not working for xlf!
595  then
596    echo 'DEBUG DEBUG DEBUG DEBUG'
597    sed s+"# -g"+" -g -no-ftz -traceback -ftrapuv -fp-stack-check -check"+g configure.wrf > yeah
598    #sed s+"# -g"+"-g -Ktrap=fp -Mbounds"+g configure.wrf > yeah
599    #sed s+"# -g"+"-g"+g configure.wrf > yeah
600    mv -f yeah configure.wrf
601  fi
602
603  if [[ "${phys}" == *"prescribed"* ]]
604  then
605     donotcompilephys=1
606  fi
607
608  # some compiler and environment mess things up when using -C
609  sed s+"-C"+""+g configure.wrf > yeah ; mv -f yeah configure.wrf   
610
611#exit
612##################################################
613# compile physics
614##################################################
615if [ ${donotcompilephys} -eq 0 ] #|| [[ "${phys}" != *"prescribed"* ]]
616then
617
618 cd $phys/
619
620 if [[ "${phys}" == *"new"* ]]
621 then
622  divx=1
623  divy=1
624 else
625  # required size
626  #----------------
627  case ${numproc} in
628      1) divx=1 ; divy=1  ;;
629      2) divx=1 ; divy=2  ;;
630      4) divx=2 ; divy=2  ;;
631      6) divx=2 ; divy=3  ;;
632      8) divx=2 ; divy=4  ;;
633     12) divx=3 ; divy=4  ;;
634     16) divx=4 ; divy=4  ;;
635     20) divx=4 ; divy=5  ;;
636     24) divx=4 ; divy=6  ;;
637     32) divx=4 ; divy=8  ;;
638     64) divx=8 ; divy=8  ;;
639    128) divx=8 ; divy=16 ;;
640      *) echo not yet supported ; exit ;;
641  esac
642 fi
643  physx=$(expr ${lon} - 1)
644  restex=$(expr ${physx} \% ${divx})   
645  physy=$(expr ${lat} - 1)
646  restey=$(expr ${physy} \% ${divy})
647  if [[ ${restex} != 0 || ${restey} != 0 ]]
648  then
649     echo nx-1 shall be divided by ${divx}
650     echo ny-1 shall be divided by ${divy}
651     exit
652  fi
653  physx=$(expr ${physx} \/ ${divx})
654  physy=$(expr ${physy} \/ ${divy})
655  physz=$(expr ${level} - 1)
656
657  # change this if you change num_soil_layers in WRF
658  # -- default is 10
659  if [[ "${phys}" == *"new"* ]]
660  then
661     soilsize=18 ## nouvelle physique. ne sert a rien ici. voir comsoil.
662  else
663     soilsize=10
664  fi
665
666  # GCM environment variables
667  #--------------------------
668  export LMDGCM=$PWD
669  export LIBOGCM=$PWD/libo
670
671  # generate the appropriate dimphys
672  #---------------------------------
673  if [[ "${phys}" == *"mars"* ]]
674  then
675    cd libf/phymars
676  elif [[ "${phys}" == *"venus"* ]]
677  then
678    cd libf/phyvenus
679  fi
680  \rm dimphys.h 2> /dev/null
681#  if [[ "${phys}" == "newphys_" ]]
682#  then
683#    ln -sf ../../meso_dimphys.h_ref . ## ne sert plus !!
684#  else
685  if [[ "${phys}" == "mars_lmd" ]]
686  then
687    physize=$(expr ${physx} \* ${physy})
688    sed s/--xsize--/${physx}/g meso_dimphys.h_ref | sed s/--ysize--/${physy}/g | sed s/--physize--/${physize}/g | sed s/--zsize--/${physz}/g | sed s/--soilsize--/${soilsize}/g > dimphys.h
689    head -15 dimphys.h
690  fi
691     #### TEST new new phys
692     #if [[ "${phys}" == "newphys_" ]]
693     #then
694     #  touch gr_fi_dyn.F.lien
695     #    ln -sf ../dyn3d/gr_fi_dyn.F .  ## dommage, a corriger
696     ################################
697     ################################
698     ################################
699     #fi
700     ####
701     ####
702  cd ../..
703
704  # prepare for nesting
705  #---------------------
706  cd libf
707  duplicate${dom} 2> /dev/null
708  cd ..
709
710  # compile physics
711  #--------------------------
712  \rm libf/grid/dimensions.h 2> /dev/null
713  \rm -rf libo/* 2> /dev/null
714  echo 1. compiling LMD physics ...
715  echo compilation info in:
716  echo $PWD/libo/log_compile_phys
717
718  if [ ${debug} -ne 0 ] 
719  then
720    if [[ "${phys}" == *"mars"* ]]
721    then
722      echo 'DEBUG DEBUG DEBUG DEBUG'
723      nohup makegcm_${compilo} -debug -full -t ${tra} -p mars -d ${physz} testphys1d | tee libo/log_compile_phys | grep 'warnings' >> libo/log_compile_phys
724    elif [[ "${phys}" == *"venus"* ]]
725    then
726      echo 'DEBUG DEBUG DEBUG DEBUG'
727      nohup ./makelmdz -debug -full -p venus -d 51 rcm1d -arch linux-ifort -cpp MESOSCALE  | tee libo/log_compile_phys | grep 'warnings' >> libo/log_compile_phys # on /u/ and /planeto/
728      nohup ./makelmdz --debug p venus -d 51 rcm1d -arch linux-ifort -cpp MESOSCALE  | tee libo/log_compile_phys | grep 'warnings' >> libo/log_compile_phys
729    fi
730  else
731  if [[ "${phys}" == *"mars"* ]]
732    then
733      nohup makegcm_${compilo} -t ${tra} -p mars -d ${physz} testphys1d | tee libo/log_compile_phys | grep 'warnings' >> libo/log_compile_phys
734    elif [[ "${phys}" == *"venus"* ]]
735    then
736      nohup ./makelmdz -full -p venus -d 51 rcm1d -arch linux-ifort_meso -cpp MESOSCALE  | tee libo/log_compile_phys | grep 'warnings' >> libo/log_compile_phys # on /u/ and /planeto/
737      nohup ./makelmdz -p venus -d 51 rcm1d -arch linux-ifort_meso -cpp MESOSCALE  | tee libo/log_compile_phys | grep 'warnings' >> libo/log_compile_phys
738      #nohup ./makelmdz -full -p venus -d 51 rcm1d -arch CICLADifort -cpp MESOSCALE  | tee libo/log_compile_phys | grep 'warnings' >> libo/log_compile_phys # on CICLAD
739      #nohup ./makelmdz -p venus -d 51 rcm1d -arch CICLADifort -cpp MESOSCALE  | tee libo/log_compile_phys | grep 'warnings' >> libo/log_compile_phys
740      #nohup ./makelmdz -full -p venus -d 51 rcm1d -arch X64_MESU -cpp MESOSCALE  | tee libo/log_compile_phys | grep 'warnings' >> libo/log_compile_phys #on MESU
741      #nohup ./makelmdz -p venus -d 51 rcm1d -arch X64_MESU -cpp MESOSCALE  | tee libo/log_compile_phys | grep 'warnings' >> libo/log_compile_phys
742      #nohup makelmdz -full -p venus -d 51 rcm1d -arch X64_OCCIGEN -cpp MESOSCALE  | tee libo/log_compile_phys | grep 'warnings' >> libo/log_compile_phys #on OCCIGEN
743      #nohup makelmdz -p venus -d 51 rcm1d -arch X64_OCCIGEN -cpp MESOSCALE  | tee libo/log_compile_phys | grep 'warnings' >> libo/log_compile_phys
744    fi
745  fi
746  echo ... done
747        # clean the duplicate routines
748        cd libf
749        duplicate1 2> /dev/null
750        cd ..
751
752  # merge LMD executables in one lib
753  #--------------------------------------
754  cd libo
755  mkdir temp
756  if [[ "${phys}" == *"mars"* ]]
757  then
758    cp -f LINUX*/*.a temp
759  else #VENUS
760    cp -f linux*/*.a temp/ # on /u/ and /planeto/
761    #cp -f /home/mlefevre/modipsl/lib/libioipsl.a temp/ #on CICLAD
762    #cp -f X64*/*.a temp #on MESU and OCCIGEN
763    cp -f ../ioipsl/libioipsl.a temp/
764  fi
765  cd temp
766  if [[ "${phys}" == *"mars"* ]]
767   then   
768   ar x libbibio.a
769   ar x libphymars.a
770   ar x libaeronomars.a
771  else
772   ar x libphyvenus.a
773   ar x libmisc.a
774   ar x libphy_common.a
775   ar x libioipsl.a
776  fi
777  \rm *.a
778  ar r liblmd.a *
779  cp -f liblmd.a ..
780  cd ..
781  \rm -r temp
782  nm liblmd.a > liblmd_content
783  # finish merge
784  cd ..
785  if [[ "${phys}" == *"new"* ]]
786  then
787    if [[ "${phys}" == *"mars"* ]]
788    then
789      # include modules for interfacing
790      mv  *.mod ../inc/
791    elif [[ "${phys}" == *"venus"* ]] 
792    then
793      mv  ./libo/linux*/*.mod ../inc/ # on /u/ and /planeto/
794      #mv  ./libo/CICLAD*/*.mod ../inc/ # on CICLAD
795      #mv  ./libo/X64*/*.mod ../inc/ # on MESU and OCCIGEN
796      cp -f ./ioipsl/*.mod ../inc/
797    fi
798  fi
799
800
801  # save a copy
802  #--------------
803  cp -f libo/liblmd.a libo/liblmd.a_${lon}_${lat}_${level}_${dom}_${tra}
804  echo '****************************************'
805
806  # ok
807  #----
808  cd ..
809
810  # manage nest includes in module_lmd_driver.F
811  #----
812  cp call_meso_inifis$dom.inc call_meso_inifis.inc
813  cp call_meso_physiq$dom.inc call_meso_physiq.inc
814
815
816fi
817##################################################
818# END compile physics
819##################################################
820
821#------------------
822# compile WRF
823#------------------
824if [ ${justphys} -eq 0 ]
825then
826
827  echo 2. compiling WRF dynamical core ...
828
829  #if [[ "${phys}" == "nophys_" ]]
830  #then
831  # echo 'NO LMD PHYSICS included'
832  #else
833  # if [[ ! ( -f "call_meso_physiq.inc" ) ]]
834  # then
835  #  echo 'did you compile the physics ? no call_meso_physiq.inc !'
836  #  exit
837  # fi
838  #fi
839
840  # be sure to compile with the most recent physics
841  touch phys/module_lmd_driver.F
842
843  # talk to user
844  echo '>>> compiling ... this may be long ... <<<'
845  echo check progress in:
846  echo $PWD/log_compile
847  echo check possible errors in:
848  echo $PWD/log_error
849 
850  # compile ...
851  echo '>>> YOUR CONFIG IS : '${config}
852  case ${config} in
853    'real')         ./compile em_real > log_compile 2> log_error
854                    # save executables 
855                    cd main
856                      if [[ -f real.exe ]]
857                      then
858                        echo 'Looks good ! real.exe is here...'
859                      fi
860                      if [[ -f wrf.exe ]]
861                      then
862                        echo 'Looks good ! wrf.exe is here...'
863                      fi
864                    cp -f real.exe ../../real_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.exe
865                    cp -f wrf.exe ../../wrf_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.exe
866                    cd ..  ;;
867    'ideal')        #mkdir 'test/em_quarter_ss' 2> /dev/null
868                    #echo '>>> YOUR CONFIG IS : '${config}
869                    ./compile em_quarter_ss > log_compile 2> log_error
870                    # save executables 
871                    cd main
872                      if [[ -f ideal.exe ]]
873                      then
874                        echo 'Looks good ! ideal.exe is here...'
875                      fi
876                      if [[ -f wrf.exe ]]
877                      then
878                        echo 'Looks good ! wrf.exe is here...'
879                      fi
880                    cp -f ideal.exe ../../ideal_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.exe
881                    cp -f wrf.exe ../../wrf_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.exe
882                    cd ..  ;;
883    'les')          #echo '>>> YOUR CONFIG IS : '${config}
884                    mkdir 'test/em_les' 2> /dev/null                   
885                    ./compile em_les > log_compile 2> log_error
886                    # save executables 
887                    cd main
888                      if [[ -f ideal.exe ]]
889                      then
890                        echo 'Looks good ! ideal.exe is here...'
891                      fi
892                      if [[ -f wrf.exe ]]
893                      then
894                        echo 'Looks good ! wrf.exe is here...'
895                      fi
896                    cp -f ideal.exe ../../ideal_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.exe
897                    cp -f wrf.exe ../../wrf_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.exe
898                    cd ..  ;;
899    *)              echo not supported... please use ; echo ideal les ; exit ;;
900  esac
901  echo '*******last lines from log_error*********'
902  tail -n 20 log_error
903fi
904
905  # the end
906  echo '****************************************'
907  echo 'done.'
908  echo '****************************************'
909  cp last ../makemeso_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}
910  mv last ../../
911  svn info ../../ > ../makemeso_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.svn.info
912
913  # link latest compiled folder as LATEST
914  cd ../../
915  rm -rf LATEST
916  ln -sf ${conf_wrf} LATEST
917
918  # link latest compiled exec
919  cd LATEST
920  case ${config} in
921    'real')         ln -sf real_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.exe real.exe ;;
922    'ideal')        ln -sf ideal_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.exe ideal.exe ;;
923    'les')          ln -sf ideal_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.exe ideal.exe ;;
924  esac
925  ln -sf wrf_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.exe wrf.exe
926  cd ..
927
928  # add here specific messages
929  if [[ "${dom}" != "1" ]]
930  then
931  nest=$(expr ${lon} + 4)
932  echo NB: in namelist.input, please set:
933  echo '
934  max_dom = '$dom'
935  s_we    = 1,1,
936  e_we    = '$lon','$nest',
937  s_sn    = 1,1,
938  e_sn    = '$lat','$nest',
939  s_vert  = 1,1,
940  e_vert  = '$level','$level','
941  fi
942
943
Note: See TracBrowser for help on using the repository browser.