Subversion Repositories lagranto.ecmwf

Rev

Rev 21 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 21 Rev 27
Line 37... Line 37...
37 
\rm -f P${date}_q
 37 
\rm -f P${date}_q
38 
cdo -f nc -t ecmwf copy -invertlat ${grbdir}/an${date}_q P${date}_q
 38 
cdo -f nc -t ecmwf copy -invertlat ${grbdir}/an${date}_q P${date}_q
39 
 
 39 
 
40 
# Convert an${date}_ps
 40 
# Convert an${date}_ps
41 
\rm -f P${date}_ps
 41 
\rm -f P${date}_ps
42 
cdo -f nc -t ecmwf copy -invertlat ${grbdir}/an${date}_ps P${date}_ps
 42 
cdo -f nc -t ecmwf copy -invertlat ${grbdir}/an${date}_ps P${date}_ps_scratch
43 
ncap2   -O -s 'PS=0.01f*exp(LNSP)' P${date}_ps  P${date}_ps
 43 
#ncap2   -O -s 'PS=0.01f*exp(LNSP)' P${date}_ps  P${date}_ps
- 
 
 44 
cdo -O -f nc -b 64 -L expr,'PS=0.01*exp(LNSP)' P${date}_ps_scratch P${date}_ps
44 
 
 45 
 
45 
# Merge all files
 46 
# Merge all files
46 
\rm -f P${date}
 47 
\rm -f P${date}
47 
cdo -f nc merge P${date}_tuvw P${date}_q P${date}_ps P${date}
 48 
cdo -f nc merge P${date}_tuvw P${date}_q P${date}_ps P${date}
48 
\rm -f P${date}_tuvw
 49 
\rm -f P${date}_tuvw
49 
\rm -f P${date}_q
 50 
\rm -f P${date}_q
50 
\rm -f P${date}_ps
 51 
\rm -f P${date}_ps
- 
 
 52 
\rm -f P${date}_ps_scratch
51 
 
 53 
 
52 
# Proceed to next date
 54 
# Proceed to next date
53 
set date=`newtime ${date} ${timestep}`
55 
set date=`newtime ${date} ${timestep}`
54 
if ( "${date}" != "${finaldate}" ) goto loop
 56 
if ( "${date}" != "${finaldate}" ) goto loop
55 
 
 57