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