Context Navigation

runmeso @ 235

Last change on this file since 235 was 235, checked in by aslmd, 14 years ago
MESOSCALE: corrected runmeso and gcm initialization so that it could be used on GNOME ifort cluster and virtually everywhere. changed choice numbering and added choices to initialize but not execute (good for job submission). Also corrected what was broken in the two previous commits on the version of the mesoscale model with old physics
Property svn:executable set to ``*
File size: 21.8 KB

Line
1	#! /bin/bash
2
3	#########################################################
4	# runmeso : #
5	# utility to combine compilation + preprocessing + run #
6	#########################################################
7	###############################
8	# Author: A. Spiga #
9	# New version : December 2008 #
10	# Modified : 2009-2010 #
11	# New physics : 2011 #
12	###############################
13
14	# TODO: - GCM compile
15	# - NETCDF
16	# - (modif registry)
17	# - quick restart
18
19
20	###########################################
21	###########################################
22	###########################################
23
24	changeregis=0 ## 1 if changed registry
25	#changeregis=1 ## 1 if changed registry
26	i_want_to_compile=1
27
28	scenario=""
29	#scenario="storm"
30
31	if [[ "${MESO}" = "" ]]
32	then
33	echo "Please initialize the variable MESO in your environnement (or modify scripts)"
34	exit
35	fi
36	####################################
37	# A FEW DEFINITIONS YOU CAN MODIFY #
38	####################################
39	main_folder=$MESO
40	meso_folder=${main_folder}'/LMD_MM_MARS'
41	calendar=${main_folder}'/LMD_MM_MARS/SIMU/calendar'
42	simu_folder=$PWD
43
44	################
45	# WELCOME USER #
46	################
47	echo '**********************************************'
48	echo '** LMD Martian Mesoscale Model launcher **'
49	echo '**********************************************'
50	#echo '** Hello '$(whoami)' on '$(hostname)' !'
51	#echo '** ['$(uname -a)']'
52	echo '** User >>> '$(whoami)
53	echo '** Hostname >>> '$(hostname)
54	echo '** Hardware >>> '$(uname -m)
55	echo '** Processor >>> '$(uname -p)
56	echo '**********************************************'
57
58	######## NEWPHYS ???
59	if [[ ! ( -f "run.def" ) ]]
60	then
61	newphys=0
62	else
63	echo '** New LMD physics >>> YES [run.def is here]'
64	echo '**********************************************'
65	newphys=1
66	if [[ ! ("${scenario}" == "") ]]
67	then
68	echo '** Scenario >>>' ${scenario}
69	echo '**********************************************'
70	fi
71	fi
72
73	##########################
74	# GET INFO FROM NAMELIST #
75	##########################
76	#echo '***********************'
77	#echo '** Your settings **'
78	#echo '***********************'
79	echo '** Simulation directory >>> '${simu_folder}
80	if [[ ! ( -f "namelist.input" ) ]]
81	then
82	echo 'no namelist.input file !'
83	exit
84	fi
85	#
86	# READ INFO NEEDED FOR INIT
87	#
88	year=$(more ${simu_folder}/namelist.input \| grep start_year \| awk '{print $3}' \| sed s+','+''+g)
89	month=$(more ${simu_folder}/namelist.input \| grep start_month \| awk '{print $3}' \| sed s+','+''+g)
90	day=$(more ${simu_folder}/namelist.input \| grep start_day \| awk '{print $3}' \| sed s+','+''+g)
91	hour=$(more ${simu_folder}/namelist.input \| grep start_hour \| awk '{print $3}' \| sed s+','+''+g)
92	wrf_date=${year}'-'${month}'-'${day}'_00:00:00'
93	daygcm=$(more ${calendar} \| grep ${wrf_date} \| awk '{print $2}')
94	lsubs=$(more ${calendar} \| grep ${wrf_date} \| awk '{print $3}')
95	eyear=$(more ${simu_folder}/namelist.input \| grep end_year \| awk '{print $3}' \| sed s+','+''+g)
96	emonth=$(more ${simu_folder}/namelist.input \| grep end_month \| awk '{print $3}' \| sed s+','+''+g)
97	eday=$(more ${simu_folder}/namelist.input \| grep end_day \| awk '{print $3}' \| sed s+','+''+g)
98	ehour=$(more ${simu_folder}/namelist.input \| grep end_hour \| awk '{print $3}' \| sed s+','+''+g)
99	ewrf_date=${eyear}'-'${emonth}'-'${eday}'_00:00:00'
100	edaygcm=$(more ${calendar} \| grep ${ewrf_date} \| awk '{print $2}')
101	elsubs=$(more ${calendar} \| grep ${ewrf_date} \| awk '{print $3}')
102	echo '** Start date >>> '${year}'-'${month}'-'${day}" [sol="${daygcm}", Ls="${lsubs}"]"
103	echo '** End date >>> '${eyear}'-'${emonth}'-'${eday}" [sol="${edaygcm}", Ls="${elsubs}"]"
104	#
105	# READ INFO NEEDED FOR DOMAIN
106	#
107	interval_seconds="$(more ${simu_folder}/namelist.input \| grep -e '[^!]interval_seconds =')"
108	if [[ "${interval_seconds}" = "" ]]
109	then
110	interval_seconds="interval_seconds = 3700 !! (p2) Frequency of large-scale fields update (s)"
111	fi
112	#
113	max_dom="$(more ${simu_folder}/namelist.input \| grep -e '[^!]max_dom =')"
114	if [[ "${max_dom}" = "" ]]
115	then
116	max_dom="max_dom = 1 !! (r)(n)(p2) Total number of domains"
117	fi
118	#
119	parent_id="$(more ${simu_folder}/namelist.input \| grep -e '[^!]parent_id =')"
120	if [[ "${parent_id}" = "" ]]
121	then
122	parent_id="parent_id = 0, !! (n)(p2) Associated parent domain"
123	fi
124	#
125	i_parent_start="$(more ${simu_folder}/namelist.input \| grep -e '[^!]i_parent_start =')"
126	if [[ "${i_parent_start}" = "" ]]
127	then
128	i_parent_start="i_parent_start = 0, !! (n)(p2) x-position of the bottom-left nest corner"
129	fi
130	#
131	j_parent_start="$(more ${simu_folder}/namelist.input \| grep -e '[^!]j_parent_start =')"
132	if [[ "${j_parent_start}" = "" ]]
133	then
134	j_parent_start="j_parent_start = 0, !! (n)(p2) y-position of the bottom-left nest corner"
135	fi
136	#
137	parent_grid_ratio="$(more ${simu_folder}/namelist.input \| grep -e '[^!]parent_grid_ratio =')"
138	if [[ "${parent_grid_ratio}" = "" ]]
139	then
140	parent_grid_ratio="parent_grid_ratio = 1, !! (n)(p2) Ratio of horizontal resolution parent/nest"
141	fi
142	#
143	e_we="$(more ${simu_folder}/namelist.input \| grep -e '[^!]e_we =')"
144	e_sn="$(more ${simu_folder}/namelist.input \| grep -e '[^!]e_sn =')"
145	#
146	dx="$(more ${simu_folder}/namelist.input \| grep -e '[^!]dx =' \| awk '{print $3}' \| sed s+','+''+g)"
147	dy="$(more ${simu_folder}/namelist.input \| grep -e '[^!]dy =' \| awk '{print $3}' \| sed s+','+''+g)"
148
149
150	#
151	# READ INFO NEEDED FOR COMPILATION
152	#
153	nx=$(more ${simu_folder}/namelist.input \| grep -e '[^!]e_we =' \| awk '{print $3}' \| sed s+','+''+g)
154	ny=$(more ${simu_folder}/namelist.input \| grep -e '[^!]e_sn =' \| awk '{print $3}' \| sed s+','+''+g)
155	nz=$(more ${simu_folder}/namelist.input \| grep -e '[^!]e_vert =' \| awk '{print $3}' \| sed s+','+''+g)
156	nest=$(more ${simu_folder}/namelist.input \| grep -e '[^!]max_dom =' \| awk '{print $3}' \| sed s+','+''+g)
157	if [[ "${nest}" = "" ]]
158	then
159	nest='1'
160	fi
161	numcallphys=$(ls callphys*def \| wc -l)
162	if [[ ${numcallphys} -lt ${nest} ]]
163	then
164	echo '**********'
165	echo You should have ${nest} callphys.def' file(s)' !!
166	exit
167	fi
168	mars_scenario=$(more ${simu_folder}/namelist.input \| grep mars \| awk '{print $3}' \| sed s+','+''+g)
169	case ${mars_scenario} in
170	0) tracers='1' ;;
171	1) tracers='2' ;;
172	2) tracers='1' ;;
173	3) tracers='2' ;;
174	11) tracers='4' ;;
175	*) tracers='1' ;;
176	esac
177	echo '** Domain(s) >>> '${nest}
178	echo '** Tracer(s) >>> '${tracers}
179
180	############ SCATTERERS SCATTERERS ######
181	if [[ ${newphys} -eq 1 ]]
182	then
183	#### CHANGE THIS IF YOU ADD NEW SCENARIOS
184	case ${mars_scenario} in
185	0) scat=1 ;;
186	1) echo "!! How many scatterers are you using : 1,2,..." ; read scat ;;
187	2) scat=1 ;;
188	3) scat=1 ;;
189	11) echo "!! How many scatterers are you using : 1,2,..." ; read scat ;;
190	20) scat=1 ;;
191	*) scat=1 ;;
192	esac
193	echo '** Scatterers >>> '${scat}
194	else
195	#### USELESS here BUT COMPATIBILITY
196	scat=1
197	fi
198	#########################################
199
200	echo '** Grid points >>> '${nx}" x "${ny}" x "${nz}
201	echo '** Horizontal resolution >>> '${dx}
202
203	#echo "OK ? ('no' to exit, any key if OK)" ; read wellok
204	# if [[ "${wellok}" = "no" ]]
205	# then
206	# exit
207	# fi
208
209	###########
210	# ACTIONS #
211	###########
212	echo '**********************************************'
213	echo "** The runmeso script does it for you !"
214	echo '**********************************************'
215	echo "** Step 0: compile the model if not done yet"
216	echo "** Step 1: run GCM and prepare inputs"
217	echo "** Step 2: interpolate on domain horizontal grid"
218	echo "** Step 3: interpolate on chosen vertical levels"
219	echo "** Step 4: run the model"
220	echo "************************* Please make a choice"
221	echo " 1 --> Perform steps 0 1 2 3 4 [ -1: step 0 skipped]"
222	echo " 2 --> Perform steps 0 2 3 4 [ -2: step 0 skipped]"
223	echo " 3 --> Perform steps 0 3 4 [ -3: step 0 skipped]"
224	echo " 4 --> Perform steps 0 4 [ -4: step 0 skipped]"
225	echo "11 --> Perform steps 0 1 2 3 [-11: step 0 skipped]"
226	echo "12 --> Perform steps 0 2 3 [-12: step 0 skipped]"
227	echo "13 --> Perform steps 0 3 [-13: step 0 skipped]"
228	echo " 0 --> Perform step 0 only"
229	echo "99 --> Oops ! Exit now please."
230	echo '**********************************************'
231	echo "What is your choice? " ; read step
232	if [[ ${step} -lt 0 ]]
233	then
234	step=$(expr 0 - ${step})
235	i_want_to_compile=0
236	echo ".... ok, I skip the compiling & linking step. you'll just be asked the compiler you previously defined."
237	fi
238	case ${step} in
239	99) exit ;;
240	0) gcmstep=0 ; wpsstep=0 ; realstep=0 ; wrfstep=0 ;;
241	1) gcmstep=1 ; wpsstep=1 ; realstep=1 ; wrfstep=1 ;;
242	2) gcmstep=0 ; wpsstep=1 ; realstep=1 ; wrfstep=1 ;;
243	3) gcmstep=0 ; wpsstep=0 ; realstep=1 ; wrfstep=1 ;;
244	4) gcmstep=0 ; wpsstep=0 ; realstep=0 ; wrfstep=1 ;;
245	11) gcmstep=1 ; wpsstep=1 ; realstep=1 ; wrfstep=0 ;;
246	12) gcmstep=0 ; wpsstep=1 ; realstep=1 ; wrfstep=0 ;;
247	13) gcmstep=0 ; wpsstep=0 ; realstep=1 ; wrfstep=0 ;;
248	) echo "* I guess you do not know what to do. Let's start from the beginning."
249	gcmstep=1 ; wpsstep=1 ; realstep=1 ; wrfstep=1 ;;
250	esac
251
252
253	###################################
254	# LINK THE EXECUTABLES OR COMPILE #
255	###################################
256	echo '**********************************************'
257	echo "** Compiling option "
258	echo "** <1> pgf90"
259	echo "** <2> g95"
260	echo "** <3> mpi+pgf90"
261	echo "** <4> ifort"
262	echo "** <5> mpi+ifort"
263	echo "** .... more in makemeso ...."
264	echo '**********************************************'
265	echo "Your choice ? (1-5)" ; read compilo
266	if [[ "${compilo}" == "" \|\| ${compilo} -lt 1 \|\| ${compilo} -gt 5 ]]
267	then
268	echo "** well, let's try g95"
269	compilo=2
270	fi
271	if [[ ${compilo} -eq 3 \|\| ${compilo} -eq 5 ]]
272	then
273	echo How many processors ? e.g. 1, 2, 4, 6, 8, 12, 16, 20, 24, 32, 64, 128 ; read numproc
274	else
275	numproc=1
276	fi
277	#if [[ ${i_want_to_compile} -ge 1 ]]
278	#then
279	cd ${meso_folder}
280	echo ${compilo} > tmp
281	if [[ ${compilo} -eq 3 ]]
282	then
283	echo ${numproc} >> tmp
284	fi
285	echo ${nest} >> tmp
286
287	if [[ ${newphys} -eq 0 ]]
288	then
289	./makemeso -d < tmp > /dev/null
290	else
291	if [[ "${scenario}" == "" ]]
292	then
293	./makemeso -d -p < tmp > /dev/null
294	else
295	./makemeso -d -p -s ${scenario} < tmp > /dev/null
296	fi
297	fi
298	\rm tmp 2> /dev/null
299	work_folder=${meso_folder}'/'$(more what_folder)
300	echo 'work_folder : '${work_folder}
301	\rm what_folder 2> /dev/null
302
303	##################### NB: define the domain if was not done yet
304	if [[ -f "${simu_folder}/namelist.wps" ]]
305	then
306	happy='y'
307	else
308	happy=''
309	fi
310	if [[ ${step} -le 2 ]]
311	then
312	while [[ "${happy}" != "y" ]]
313	do
314	#vi ${simu_folder}/namelist.input
315	if [[ ! ( -d "${work_folder}/WPS" ) ]]
316	then
317	echo 'no WPS. you need to install it for your architecture'
318	exit
319	fi
320	\rm ${simu_folder}/namelist.wps 2> /dev/null
321	echo '**********************************************'
322	echo " !!! NO namelist.wps FILE FOUND !!!"
323	echo '**********************************************'
324	echo 'Please define the domain (not for nests for the moment) : '
325	echo " - 1/4 Projection: polar lambert mercator ? previous is "${proj} ; read proj
326	echo " - 2/4 Data source (multiple if nests): 'gcm' '32ppd' '64ppd' '64ppd_noHRti' '32ppd_HRalb' ? previous is "${data_source} ; read data_source
327	echo " - 3/4 North Latitude of domain centre ? previous is "${lat} ; read lat
328	echo " - 4/4 East Longitude of domain centre ? previous is "${lon} ; read lon
329	#####
330	sed s+"__INTERVAL_SECONDS__"+"${interval_seconds}"+g ${meso_folder}/SIMU/namelist.wps \| sed s+"__MAX_DOM__"+"${max_dom}"+g \| sed s+"__PARENT_ID__"+"${parent_id}"+g \| sed s+"__PARENT_GRID_RATIO__"+"${parent_grid_ratio}"+g \| sed s+"__I_PARENT_START__"+"${i_parent_start}"+g \| sed s+"__J_PARENT_START__"+"${j_parent_start}"+g \| sed s+"__E_WE__"+"${e_we}"+g \| sed s+"__E_SN__"+"${e_sn}"+g \| sed s+"__DX__"+"${dx}"+g \| sed s+"__DY__"+"${dy}"+g \| sed s+"__PROJ__"+"${proj}"+g \| sed s+"__DATA_SOURCE__"+"${data_source}"+g \| sed s+"__LAT__"+"${lat}"+g \| sed s+"__LON__"+"${lon}"+g \| sed s+'0000-00-00_00:00:00'+${year}'-'${month}'-'${day}'_'${hour}':00:00'+g \| sed s+'1111-11-11_11:11:11'+${eyear}'-'${emonth}'-'${eday}'_'${ehour}':00:00'+g > ${simu_folder}/namelist.wps
331	####
332	ln -sf ${simu_folder}/namelist.wps ${work_folder}/WPS/namelist.wps
333	#vi ${work_folder}/WPS/namelist.wps
334	echo '**********************************************'
335	echo OK, thank you. Please wait... ncview will pop up soon !
336	cd ${work_folder}/WPS/
337	\rm geo_em.d0*.nc 2> /dev/null
338	./geogrid.exe > /dev/null 2> /dev/null
339	ncview geo_em.d01.nc > /dev/null 2> /dev/null &
340	if [[ -f "geo_em.d02.nc" ]]
341	then
342	ncview geo_em.d02.nc > /dev/null 2> /dev/null &
343	fi
344	if [[ -f "geo_em.d03.nc" ]]
345	then
346	ncview geo_em.d03.nc > /dev/null 2> /dev/null &
347	fi
348	cd ${simu_folder}/
349	echo '**********************************************'
350	echo 'Are you happy with the results ? type y if yes, any key if no' ; read happy
351	done
352	fi
353	#####################
354	if [[ ${i_want_to_compile} -ge 1 ]]
355	then
356	cd ${simu_folder}
357	## corresponding executables
358	wrfexe=${work_folder}/wrf_x${nx}_y${ny}_z${nz}_d${nest}_t${tracers}_p${numproc}_s${scat}.exe
359	realexe=${work_folder}/real_x${nx}_y${ny}_z${nz}_d${nest}_t${tracers}_p${numproc}_s${scat}.exe
360	if [[ -f ${wrfexe} ]]
361	then
362	echo '**********************************************'
363	echo " !!! wrf.exe exists in ${work_folder}"
364	echo '**********************************************'
365	echo "Recompile ? ('yes' to accept, any key to refuse)" ; read reply
366	else
367	reply='yes'
368	fi
369	echo '********************************************** OK thanks. Now please wait...'
370	echo ' '
371	if [[ "${reply}" == "yes" ]]
372	then
373	#####################
374	echo "** "Begin compilation process...
375	\rm ${wrfexe} 2> /dev/null
376	\rm ${realexe} 2> /dev/null
377	#
378	# CREATE makemeso_input
379	#
380	\rm ${simu_folder}/makemeso_input 2> /dev/null
381	touch ${simu_folder}/makemeso_input
382	echo ${compilo} >> ${simu_folder}/makemeso_input
383	if [[ "${compilo}" = "3" \|\| "${compilo}" = "5" ]]
384	then
385	echo ${numproc} >> ${simu_folder}/makemeso_input
386	fi
387	echo ${nx} >> ${simu_folder}/makemeso_input
388	echo ${ny} >> ${simu_folder}/makemeso_input
389	echo ${nz} >> ${simu_folder}/makemeso_input
390	echo ${nest} >> ${simu_folder}/makemeso_input
391	echo ${tracers} >> ${simu_folder}/makemeso_input
392	if [[ ${newphys} -eq 1 ]]
393	then
394	echo ${scat} >> ${simu_folder}/makemeso_input
395	#echo 'n' >> ${simu_folder}/makemeso_input
396	fi
397	#
398	# COMPILE THE MODEL
399	#
400	cd ${meso_folder}
401	if [[ "${changeregis}" = "1" ]]
402	then
403	cd ${work_folder}/WRFV2/Registry
404	\rm Registry 2> /dev/null
405	./Registry.bash
406	cd ../../..
407	fi
408
409	##******
410	if [[ ${newphys} -eq 0 ]]
411	then
412	./makemeso -r < ${simu_folder}/makemeso_input
413	else
414	if [[ "${scenario}" == "" ]]
415	then
416	./makemeso -r -p < ${simu_folder}/makemeso_input
417	else
418	./makemeso -r -p -s ${scenario} < ${simu_folder}/makemeso_input
419	fi
420	fi
421	##******
422
423
424	if [[ -f ${wrfexe} ]]
425	then
426	echo "** "seems OK ... linking executables
427	ln -sf ${wrfexe} ${simu_folder}/wrf.exe
428	ln -sf ${realexe} ${simu_folder}/real.exe
429	chmod a+x ${realexe} ${wrfexe} ## sait-on jamais
430	cd $simu_folder
431	else
432	echo no ${wrfexe}
433	echo problem with the compilation
434	exit
435	fi
436	else
437	#echo simply linking existing executables...
438	ln -sf ${wrfexe} ${simu_folder}/wrf.exe
439	ln -sf ${realexe} ${simu_folder}/real.exe
440	chmod a+x ${realexe} ${wrfexe} ## sait-on jamais
441	fi
442	fi
443
444	######################################
445	# 1. RUN GCM FOR INIT/BDY CONDITIONS #
446	######################################
447	if [ ${gcmstep} -eq 1 ]
448	then
449	echo '************************************'
450	echo '** Step 1: Prepare GCM inputs **'
451	echo '************************************'
452	#
453	# run GCM
454	#
455	if [[ ${newphys} -eq 0 ]]
456	then
457	cd ${main_folder}'/LMDZ.MARS/'
458	else
459	cd ${main_folder}'/LMDZ.MARS.new/'
460	fi
461	if [[ ! ( -f "gcm.e" ) ]]
462	then
463	echo "Oops. There is no gcm.e, you shall compile the GCM first."
464	exit
465	fi
466	cd myGCM
467	echo ${daygcm} \| ./launch_gcm
468	#
469	# prepare files for WPS
470	#
471	cd ${work_folder}
472	if [[ -f ${work_folder}/PREP_MARS/create_readmeteo.exe && -f ${work_folder}/PREP_MARS/readmeteo.exe ]]
473	then
474	### This is now generic so PREP_MARS could be moved out of DIRCOMPs files
475	cd ${work_folder}/PREP_MARS/
476	if [[ ${compilo} -eq 1 \|\| ${compilo} -eq 3 ]]
477	then
478	ln -sf compile_pgf compile
479	fi
480	if [[ ${compilo} -eq 2 ]]
481	then
482	ln -sf compile_g95 compile
483	fi
484	if [[ ${compilo} -eq 4 \|\| ${compilo} -eq 5 ]]
485	then
486	ln -sf compile_ifort compile
487	fi
488	if [[ ${compilo} -gt 5 ]]
489	then
490	echo "please adapt compile in PREP_MARS"
491	fi
492	./compile_and_exec
493	#echo 1 \| ./create_readmeteo.exe
494	#./readmeteo.exe < readmeteo.def
495	else
496	echo please compile the PREP_MARS utilities in ${work_folder}
497	exit
498	fi
499	fi
500
501	###################################
502	# 2. RUN WRF PREPROCESSING SYSTEM #
503	###################################
504	if [ ${wpsstep} -eq 1 ]
505	then
506	echo '*****************************************'
507	echo '**** Step 2: Perform WPS interpolations *'
508	echo '*****************************************'
509	if [[ -f ${work_folder}/WPS/geogrid.exe && -f ${work_folder}/WPS/metgrid.exe ]]
510	then
511	\rm -rf ${simu_folder}/met_em 2> /dev/null
512	mkdir ${simu_folder}/met_em
513	cd ${work_folder}/WPS/
514	\rm -rf ${work_folder}/WPS/WRFFEED/current 2> /dev/null
515	mkdir ${work_folder}/WPS/WRFFEED/current
516	ln -sf ${simu_folder}/namelist.wps ${work_folder}/WPS/namelist.wps
517	#mv namelist.wps namelist.wps.tmp
518	#sed s+'0000-00-00_00:00:00'+${year}'-'${month}'-'${day}'_'${hour}':00:00'+g namelist.wps.tmp \| sed s+'1111-11-11_11:11:11'+${eyear}'-'${emonth}'-'${eday}'_'${ehour}':00:00'+g > namelist.wps
519	./geogrid.exe
520	./metgrid.exe
521	#mv namelist.wps ${simu_folder}
522	#mv namelist.wps.tmp namelist.wps
523	namedata=$(date +simu_%d%h%Y_%H%M)
524	mv ${work_folder}/WPS/WRFFEED/current ${work_folder}'/WPS/WRFFEED/'${namedata}
525	echo met_em files in $namedata
526	cd ${simu_folder}/met_em
527	ln -sf ${work_folder}/WPS/WRFFEED/${namedata}/* .
528	cd ${simu_folder}
529	else
530	echo please compile the WPS utilities in ${work_folder}
531	exit
532	fi
533	fi
534
535	###################
536	# 3. RUN REAL.EXE #
537	###################
538	if [ ${realstep} -eq 1 ]
539	then
540	echo '*******************************'
541	echo '**** Step 3: Execute real.exe *'
542	echo '*******************************'
543	cd ${simu_folder}
544	if [[ ! -d met_em ]]
545	then
546	ls
547	echo "PB: no met_em files:"
548	echo "mkdir met_em folder then link met_em** files !"
549	exit
550	fi
551	cd ${simu_folder}/met_em
552	ln -sf ${simu_folder}/namelist.input .
553	ln -sf ${simu_folder}/real.exe .
554	ln -sf ${simu_folder}/levels .
555	rm -rf ${simu_folder}/met_em/wrfbdy* 2> /dev/null
556	rm -rf ${simu_folder}/met_em/wrfinput* 2> /dev/null
557	if [[ ${compilo} -eq 3 \|\| ${compilo} -eq 5 ]]
558	then
559	rm -rf ${simu_folder}/met_em/rsl.* 2> /dev/null
560	echo "check in "${simu_folder}"/met_em/rsl.out.0000"
561	./real.exe
562	tail -n 4 ${simu_folder}"/met_em/rsl.out.0000"
563	else
564	rm -rf ${simu_folder}/met_em/log* 2> /dev/null
565	echo "check in "${simu_folder}"/met_em/log_real"
566	real.exe > log_real
567	tail -n 4 ${simu_folder}"/met_em/log_real"
568	fi
569	if [[ -f ${simu_folder}/met_em/wrfinput_d01 && -f ${simu_folder}/met_em/wrfbdy_d01 ]]
570	then
571	echo ok... link ini and bdy conditions
572	cd ${simu_folder}
573	ln -sf ${simu_folder}/met_em/wrf_d0 .
574	else
575	echo mmm... there was a problem, see log files in met_em
576	exit
577	fi
578	fi
579
580	##################
581	# 4. RUN WRF.EXE #
582	##################
583	if [ ${wrfstep} -eq 1 ]
584	then
585	xeyes &
586	echo '******************************'
587	echo '**** Step 4: Execute wrf.exe *'
588	echo '******************************'
589	cd ${simu_folder}
590	rm -rf ${simu_folder}/wrfout*
591	if [[ ${compilo} -eq 3 ]]
592	then
593	rm -rf ${simu_folder}/rsl.* 2> /dev/null
594	# ${WHERE_MPI}/mpirun -host smith -np ${numproc} wrf.exe < /dev/null &
595	${WHERE_MPI}/mpirun -np ${numproc} wrf.exe < /dev/null &
596	else
597	rm -rf ${simu_folder}/log* 2> /dev/null
598	nohup wrf.exe > log_wrf &
599	fi
600	echo OK ! Simulation should be running now !
601	echo Check with the command: "ps -e -o '%u %p %C %x %c' \| grep wrf.exe"
602	echo '******************************'
603	echo '**** End.'
604
605	fi
606	exit
607
608
609

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

Download in other formats:

Original Format