Context Navigation

makemeso @ 155

Last change on this file since 155 was 153, checked in by aslmd, 13 years ago
MESOSCALE: modified makemeso and runmeso to comply with changes performed in r148 about scatterers in newphys. a new question is now asked if you run makemeso -p or if you run runmeso with run.def in the folder. note that runmeso only asks the question for scenarios mars=1 or mars=11 where water cycle is on, otherwise by default number of scatterers is 1 (i.e. only dust).
Property svn:executable set to ``*
File size: 30.5 KB

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

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

Download in other formats:

Original Format