Project

General

Profile

DML_0025_INT2LM.sh

Matthias Demuzere, 11/02/2015 08:28 AM

 
1
#!/bin/bash
2
#
3

    
4
EXTPAR=DML_0025_EXTPAR.nc
5
EXTPARDIR=/media/Niels_1/TUTORIAL/archive/
6
OUTPUTDIR=/media/Niels_1/TUTORIAL/data/ERAINT_DML_0025_INTLM/
7
WORKDIR=$OUTPUTDIR/log        # set working directory
8
INPUTDIR='/media/Niels_1/TUTORIAL/data/ANT_022_CCLM/out01/'
9
 
10
#INT2LM_BIN=/home/rcs/software/gcc/int2lm_131101_2.00_clm3/int2lm_src/tstint2lm
11

    
12
INT2LM_BIN=/home/rcs/software/intel/cclm-sp_1.5/src/int2lm/int2lm.rcs
13
#ulimit -s unlimited
14

    
15
NPROCX=1
16
NPROCY=1
17
NPROCXY=1
18

    
19
#INT2LM_BIN=/home/rcs/software/intel/cclm-sp_1.5/src/int2lm/int2lm.rcs
20
#ulimit -s unlimited
21

    
22
if [ -d $WORKDIR ] 
23
then
24
rm -R $WORKDIR
25
fi
26

    
27
mkdir -p $DATAINT2LM
28
mkdir -p $WORKDIR
29

    
30
cd $WORKDIR
31

    
32
### FOR DML: AR event only
33
start_date=2011020100
34
end_date=2011021700
35
ydate_ini="2011021700"
36
ydate_first=$ydate_ini
37

    
38
#################################################
39
# Namelist
40
#################################################
41
cat > INPUT <<**
42
 &CONTRL
43
  ydate_ini = "${start_date}", 
44
  hstart = 0.0, hstop = 200000.0, hincbound = 3,
45
  linitial = .TRUE.,  lboundaries = .TRUE., 
46
  lfilter_oro =.TRUE.,
47
  nprocx = ${NPROCX}, nprocy = ${NPROCY},
48
  l_cressman = .FALSE.,l_bicub_spl = .FALSE., lbdclim = .TRUE.,
49
  ltime_mean = .TRUE., luvcor = .TRUE. , lreorder = .TRUE.,
50
  lmulti_layer_in = .TRUE., lmulti_layer_lm = .TRUE.,
51
  yinput_model = "COSMO",
52
  itype_w_so_rel = 1,
53
!  itype_t_cl = 1, 
54
  itype_rootdp = 4,
55
  itype_aerosol = 2,
56
  lprog_qi = .TRUE.,
57
  lforest = .TRUE.,
58
  llake = .FALSE.,
59
  llbc_smooth = .FALSE., nlbc_smooth = 10,
60
  lsso = .FALSE.,
61
  nincwait = 5, nmaxwait = 20,
62
  luse_t_skin = .TRUE.,
63
  idbg_level = 11,
64
  itype_ndvi = 0,
65
  itype_calendar = 0,
66
!  itype_albedo = 2,
67
 /END
68
 &GRID_IN
69
  pcontrol_fi = 30000.,
70
  ie_in_tot = 296, je_in_tot = 296, ke_in_tot = 40,
71
  startlat_in_tot = -32.56, startlon_in_tot = -32.56,
72
  endlat_in_tot = 32.56, endlon_in_tot = 32.56,
73
  pollat_in = 6.08, pollon_in = 13.08,
74
  dlon_in=0.22, dlat_in=0.22,
75
!  lushift_in = .FALSE., .FALSE., lvshift_in = .FALSE., .FALSE.,
76
  ke_soil_in = 9,   
77
  czml_soil_in = 0.005, 0.025, 0.07, 0.16, 0.34, 0.70, 1.42, 2.86, 5.74, 11.5,
78
  polgam_in = 180.,
79
!  east_add_in=1,  ! Neceassary having pole in the domain
80
!  south_add_in=1, ! Neceassary having pole in the domain
81
!  west_add_in=1,  ! Neceassary having pole in the domain
82
 /END
83
 &LMGRID
84
  ielm_tot = 300,
85
  jelm_tot = 300,
86
  kelm_tot = 40,
87
  pollat = 18.05, pollon = 23.35, polgam = 180., 
88
  dlon = 0.025, dlat = 0.025, 
89
  startlat_tot = -0.95,
90
  startlon_tot = -3.75,
91
  ke_soil_lm = 9,
92
  czml_soil_lm = 0.005, 0.025, 0.07, 0.16, 0.34, 0.70, 1.42, 2.86, 5.74, 11.5,
93
  czvw_so_lm = 0.75, 0.75, 0.75, 0.75, 0.75, 0.75, 0.75, 0.75, 0.75, 0.75,
94
  irefatm = 2,
95
  ivctype = 2,
96
 /END
97
 &DATABASE
98
 /END
99
 &DATA
100
  ie_ext = 304, je_ext = 304,
101
  ylmext_lfn = "${EXTPAR}",
102
  ylmext_cat = "${EXTPARDIR}",
103
  yinext_lfn = "lffd${start_date}c.nc",
104
  yinext_cat = "${INPUTDIR}",
105
  yin_cat = "${INPUTDIR}",
106
!  yinput_type = 'analysis',
107
  ytunit_in = 'd',
108
  ytunit_out = 'd',
109
  ylm_cat = "${OUTPUTDIR}",
110
  yinext_form_read = "ncdf",
111
  ylmext_form_read = 'ncdf',
112
  yin_form_read = "ncdf",
113
  ylm_form_write = 'ncdf',
114
  nprocess_ini = 131, nprocess_bd = 132,
115
 /END
116
 &PRICTR
117
  igp_tot = 36,
118
  jgp_tot = 30,
119
  lchkin = .TRUE., lchkout = .TRUE.,
120
 /END
121
**
122

    
123
echo ----- start INT2CLM
124

    
125
mpirun -np ${NPROCXY} $INT2LM_BIN
126

    
127
echo ----- INT2CLM finished
128