Context Navigation

makemeso @ 1243

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

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format