Context Navigation

makemeso @ 1242

Last change on this file since 1242 was 1236, checked in by aslmd, 11 years ago
MESOSCALE. A necessary complement commit to r1234 where a upgraded interface making use of modules was proposed. Completed the new formulation for module_lmd_driver for newphys with improved interface with both ini/bdy conditions and physical parameterizations. Changed the Registry accordingly. Finished changes about I/O with the LMD physics (see LMDZ.MARS/README). Made all those changes compatible for old interface, and LES runs (checked with test cases), as well as old input files. Changed makemeso to account for full flexibility on changin nx ny ntracers nproc with newphys. Cleaned the now obsolete bits of code used in LMD physics shared with the GCM. ----- Everything is now ready to properly code both restart runs and nesting for mesoscale runs with new physics.
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_inc/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