source: BOL/Multi_atlas/METRICS/make_mask.sh @ 4515

Last change on this file since 4515 was 3684, checked in by idelkadi, 4 years ago

Repository under svn of a first version of Multiatlas diagnostics for LMDZ. This version is adapted to be able to run a LMDZ multiatlas on an individual account on the ciclad machine of the IPSL. In this version, the parts to be modified are identified so as to subsequently adapt it to other machines.
This version is still under development.

  • Property svn:executable set to *
File size: 1.3 KB
Line 
1if [ 0 = 0 ] ; then
2for mask in pourc_ter pourc_oce pourc_sic pourc_lic ; do
3cdo timavg -selvar,$mask /prodigfs/ipslfs/dods/fabric/lmdz/SE/ORIG/NPv5.4_SE_1982_1989_1M_histmth.nc $mask.nc
4done
5cdo add pourc_oce.nc pourc_sic.nc pourc_ocesic.nc
6cdo remapcon,tas.nc pourc_ocesic.nc tmp.nc
7fi
8
9cat <<eod>| tmp.jnl
10use tas.nc
11use tmp.nc
12let lat=y*(pourc_oce*0+1)
13let lon=x*(pourc_oce*0+1)
14reg/l=1
15! shade un
16!let maskcircA=if ( abs(lat+55) lt 15 and pourc_oce gt 80 ) then 1
17let maskcircA=if ( abs(lat+50) le 10 and pourc_oce gt 80 ) then 1
18let maskoce=if ( pourc_oce gt 80 ) then 1
19let maskter=if ( pourc_oce lt 20 ) then 1
20let maskNH=if ( lat gt 0 ) then 1
21let maskSH=if ( lat lt 0 ) then 1
22let maskNAtl=if ( lat ge 45 and lat le 60 and lon ge 310 and lon le 350 ) then 1
23let maskArct=if ( lat ge 70 and pourc_oce ge 95 ) then 1
24let maskSibe=if ( lat ge 40 and lon ge 50 and lon le 180 and pourc_oce le 5 ) then 1
25let maskUSA=if ( lat ge 30 and lat le 55 and lon ge 250 and lon le 310 and pourc_oce le 5 ) then 1
26save/clobber/file=maskcircA.nc maskcircA
27save/clobber/file=maskoce.nc maskoce
28save/clobber/file=maskter.nc maskter
29save/clobber/file=maskNH.nc maskNH
30save/clobber/file=maskSH.nc maskSH
31save/clobber/file=maskNAtl.nc maskNAtl
32save/clobber/file=maskArct.nc maskArct
33save/clobber/file=maskSibe.nc maskSibe
34save/clobber/file=maskUSA.nc maskUSA
35quit
36eod
37
38ferret -script tmp.jnl
Note: See TracBrowser for help on using the repository browser.