===================================================================
Cube83
same as cube81 but
ALLOW_CLIMSST_RELAXATION to Reynold's daily 1/4-deg AVHRR analysis
(tauThetaClimRelax = 1296000.
relax daily SST within 15 days, too weak constraint?)
===================================================================
Cube82
same as cube81 but
new values for data.exf (now
exf_albedo = 0.16,
previous
exf_albedo = 0.15)
===================================================================
Cube81
same as cube79 but with
new values for data.seaice (compromised value from optimal Arctic and
Antarctic regios);
new values for data.exf (now
exf_albedo = 0.15,
exf_scal_BulkCdn = 1.015,
while previous
exf_albedo = 0.1507,
exf_scal_BulkCdn = 1.0185);
DIFFKR_2_20_1_lat6070_cube81, which has diffKrS/T=3.55e-7 in Arctic;
GMRedi turned off
===================================================================
Cube80
same as cube79 but with
SEAICE_BREAKS_VERIFICATION (SST relaxation depends on SIarea)
ALLOW_CLIMSST_RELAXATION to Reynold's daily 1/4-deg AVHRR analysis
ALLOW_CLIMSSS_RELAXATION to monthly 1-deg WOA05 analysis
===================================================================
Cube79
same as cube78 (blend3) but with
useRelativeWind=.TRUE. and SEAICE_gamma_t=259200
===================================================================
Cube78
blend #3
===================================================================
Cube77
same as cube75 but with SEAICE_SALINITY defined
SEAICE_salinity = 0.3,
and with era40_ecmwf_blend/cube59_GPCP forcing
===================================================================
Cube76
same as cube77 but 1979-2006 integration
===================================================================
Cube75
blend #2
===================================================================
Cube74
same as cube71 but initialized on January 1, 1989 from WGHC-January
climatology, where seasonal cycle is from WOA01.
===================================================================
Cube73
same as cube71 but initialized on January 1, 1989.
===================================================================
Cube72
same as cube70 but with useGMRedi=.TRUE.
GM_background_K=100, GM_Kmin_horiz=98, GM_taper_scheme='ldd97'
===================================================================
Cube71
same as cube66 but with useGMRedi=.FALSE.
hydrogThetaFile='PHC_THETA_JAN_510x6x510x50.bin'
===================================================================
Cube70
same as cube66 but with GM_background_K=150, GM_Kmin_horiz=50
Bug: PHCT is in-situ rather than potential temperature
THETA output was converted back to potential after integration
but expected to impact regions of deep convection.
===================================================================
Cube69
same as cube65 but with blending wind fields (cube49)
===================================================================
Cube68
same as cube65 but with blending spfh2m field (cube49)
===================================================================
Cube67
same as cube65 but with GPCP precipitation (cube59)
===================================================================
cube66
same as cube67 but with
SEAICE_MULTICATEGORY, SEAICE_VARIABLE_FREEZING_POINT,
rhonil=1027.5, [temp/salt]AdvScheme=7, convertFW2Salt=-1, LSR_ERROR=1e-6
useRealFreshWaterFlux=.TRUE., runoff_lat_inc=179*1,
SEAICE_availHeatFrac=.5, HeffFile=heff_file8, HEFF_MAX=10
useGMRedi=.TRUE., GM_background_K=200, GM_taper_scheme='ldd97'
bathyFile=GEBCO_510x6x510_ver06_dig.bin
hydrogSaltFile ='PHCS_JAN_510x6x510x50.bin'
hydrogThetaFile='PHCT_JAN_510x6x510x50.bin'
Bug: PHCT is in-situ rather than potential temperature
THETA output was converted back to potential after integration
but expected to impact regions of deep convection.
===================================================================
Cube65
same as cube64 but with diffKr as function of (y,z) instead of const.
Bryan and Lewis 1979 vertical diffusivity parameters
with latitudinal dependence:
#define ALLOW_BL79_LAT_VARY in CPP_OPTIONS.h
diffKrFile = 'DIFFKR_2_20_1_lat6070', in data
===================================================================
Cube64
same as cube55 but with
new ERA40_full_resolution (320x160) forcing
from /skylla/data1/data/era40/era40_fullres
===================================================================
Cube63
(crashed on 92/11/2)
same as cube61 but with
a new blending of wind:
inside -90S---50N: cube61 winds, outside 60N: NCEP_RGAU
between: linear interp
see
/host/skylla/data1/data/blend_forcing/cube63_wind
/host/skylla/data1/data/blend_forcing/cube61_forcing
===================================================================
Cube62
same as cube61 but with
a new blending of wind:
inside -50S---50N: cube61 winds, outside 60N/S: NCEP_RGAU
between: linear interp
See
/host/skylla/data1/data/blend_focing/cube62_wind
/host/skylla/data1/data/blend_focing/cube61_forcing
===================================================================
Cube61
same as cube55 but with
a new blend forcing, which has the same weight as cube49 forcing field
but has corrections in 1): corrected shift of ERA40 time step
2): corrected CORE precipiation and radiation fields
see
/host/skylla/data1/data/blend_focing/cube61_forcing
/host/skylla/data1/data/ncep/ncep_core/*corrected*
===================================================================
Cube60
same as cube55 but with
wind field derived from SSM/I 4xdaily (360x181)
see /host/skylla/data1/data/SSMI_3.0: download.sh and generating code
===================================================================
Cube59
same as cube55 but with
precipitation derived from GPCP 5-day (144x73),
see /host/skylla/data1/data/GPCP: download.sh and generating code
===================================================================
Cube58
(crashed on 93/01/18)
same as cube55 but with
& ht = 10.000 d0 ,
& hq = 10.000 d0
in exf_constants.h
===================================================================
Cube57
same as cube55 but with
SEAICE_dryIceAlb = 0.89311,
SEAICE_drySnowAlb = 0.98491,
===================================================================
Cube56
same as cube55 but with
SEAICE_wetIceAlb = 0.80131
SEAICE_wetSnowAlb = 0.84211
===================================================================
Cube55
same as cube50 but with SEAICEuseFlooding=.TRUE. which
converts snow to sea ice if flooded.
===================================================================
Cube54
same as cube50 but with local, North Atlantic Mode Water Region
perturbation to humidity (aqhfile='cube54_spfh2m/RG_spfh2m')
===================================================================
Cube53
same as cube51 but with exf_scal_BulkCdn = 0.9
===================================================================
Cube52
same as cube51 but with exf_albedo = 0.16
===================================================================
Cube51
same as cube50 but with bottomDragQuadratic = 0.0012
(but data.kpp same as cube47)
===================================================================
Cube50
same as cube48 but with SEAICEadvSnow=.true.
(but data.kpp same as cube49)
Crashed June 2002 due to 1000+ m thick, oscillating hsnow
near Antarctic coastline.
Restarted from 263016 (Jan. 1, 2002) with 20 m hsnow cap in
pickup_seaice.0000263016.data for ix=321:323; iy=364; of face 6.
Had to restart a second time from 273888 after caping snow at 12 m.
Note that same integration but with data.kpp same as cube48
crashed after one year, but nothing abnormal was noticed in
simulation leading up to crash, so probably a fluke.
===================================================================
Cube49
first blend integration, includes monthly Arctic Ocean runoff
and latest growth.F bug fixes but not snow advection
===================================================================
Cube48
same as cube49 but with more realistic initial sea-ice conditions
from a PSC estimate (skylla:/data4/cube/matlab/mk_pickup48.m)
(but data.kpp same as cube47)
===================================================================
Cube47
same as cube43
but including Peter Winsor's new monthly Arctic Ocean runoff
climatology runoff-360x180x12.bin
===================================================================
Cube46
same as cube43 but with "exf_inscal_aqh=1.1" in data.exf_hr
new data.diagnostics
crashes right after 05-Mar-1997 06:00:00
===================================================================
Cube45
same as cube43 but for forcing dataset.
The dataset is mixed NCEP adjusted by GODAE v3.i85 and NCEP.
Wind forcing fields are adjusted from GODAE(see cube44)
while other forcing fields are NCEP_RGAU (same as cube43)
Also same EXE file as cube43, without compilation
===================================================================
Cube44
same as cube43 but for forcing dataset.
The dataset is NCEP adjusted by ECCO-GODAE version 3, iter85
atemp/aqh/precip/uwind/vwind/dsw, dsw is adjusted indirectly
through swflux while dlw is unadjusted.
Also same EXE file as cube43, without compilation
===================================================================
Cube43
same as cube40 but with bug fixes in growth.F,
i.e., treatment of EVAP, HSNOW, and AREA(i,j,2)
===================================================================
Cube42
same as cube40 but for forcing dataset.
Wind forcing fields are taken from CORE(see cube41)
while other forcing fields are NCEP_RGAU (same as cube40)
Also same EXE file as cube40, without compilation
(as sea-ice package has possible modification then)
running on /nobackup2b/menemenl/cube42/MITgcm/verification/global_ocean.cs32x15/run
===================================================================
Cube41
same as cube40 but for forcing dataset.
The dataset is CORE, which is taken from
http://data1.gfdl.noaa.gov/nomads/forms/mom4/CORE/CIAF_1p0.html
Corrected Inter-Annual Forcing (CIAF_1p0)
and necessary change is made:
4Xdaily {Air temp, specific humidity, velocity}: changed to include leap yr
daily {downward shortwave,longwave radiation}: changed to 4xdaily and leap-yr
monthly {Precipitation}: changed to 4Xdaily and leap-yr
The temperature/humidity are specified at 10m height as opposed to 2m.
Also without specifying net shortwave radiation (swflux)
/nobackup2a/menemenl/cube41/MITgcm/verification/global_ocean.cs32x15/run
There was a problem with reference height. See Chris' e-mail dated 11/20/2006
6:58 PM: NCEP winds and CORE winds, air temperature, and humidity are
referenced to 10 m while NCEP air temperature and humidity are referenced to 2
m above the surface. Cube58 is an attempt to integrate with ht=hq=10 but
cube58 crashed on 93/01/18.
===================================================================
Cube40
same as cube38 but without specifying net shortwave radiation
swflux for open water, instead it is computed from downward shortwave.
this may cause considerably different summer sea ice extent.
/nobackup2b/menemenl/cube40/MITgcm/verification/global_ocean.cs32x15/run
===================================================================
Cube39
same as cube38 but surface air temperature is increased by 0.5 deg K
exf_offset_atemp=273.65 instead of 273.15 in data.exf
/nobackup2b/menemenl/cube39/MITgcm/verification/global_ocean.cs32x15/run
===================================================================
Cube38
same as cube37 but starting from WOCE climatology WGHC_(T/S)_510x6x510.bin
/nobackup2a/menemenl/cube38/MITgcm/verification/global_ocean.cs32x15/run
===================================================================
Cube37
same as cube36 (cube 31 with June 10th06 code) but
in SEAICE_OPTIONS.h
#define SEAICE_CGRID
in data.seaice
# SEAICEadvScheme = 33
(because there was a bug with implementation of SEAICEadvScheme=33 on CGRID)
An noted that following time stepd have bad AREAtave, also SST,
and possibly other fields:
0000054936 %02-Feb-1994
0000055008 %03-Feb-1994
0000055080 %04-Feb-1994
===================================================================
Cube36
cube 31 with newer code (as of June 10th '06)
in SEAICE_OPTIONS.h
#define SEAICE_TEST_ICE_STRESS_1
in data.seaice
SEAICEadvScheme = 33
pickup.0000000216.cube27
Integration is in
columbia2:/nobackup/menemenl/trm/MITgcm/verification/global_ocean.cs32x15/run36
It crashed on niter=9468, added viscA4Grid=0.01 on niter=8712
crashed again on niter=9468 (11 may 1992)
===================================================================
Cube35: cube 31 with
data.kpp: &KPP_PARM01: Riinfty = 0.701
revert to pickup.0000000216.cube27
===================================================================
Cube34: cube 31 with
data: &PARM03: tauSaltClimRelax = 7684760.,
yr: 2002, mon: 10, day: 1
===================================================================
Cube33: cube 31 with
data.kpp: &KPP_PARM01: Ricr = 0.3
yr: 2002, mon: 3, day: 1
===================================================================
Cube32: cube 31 with
data.kpp: &KPP_PARM01: Riinfty = 0.75
!incomplete
yr: 1996, mon: 3, day: 1
===================================================================
Cube31: same as Cube28 but data.seaice has following
parameters default values:
SEAICE_drag = 0.002 _d 0
OCEAN_drag = 0.001 _d 0
SEAICE_waterDrag = 5.5 _d 0
SEAICE_strength = 2.75 _d +04
This integration runs faster, close to 20 hours per year
on 216 CPUs.
Started with pickup.0000000216.cube31
===================================================================
Cube30: same as cube29 but with
diffKrT=1.0e-5,
diffKrS=1.0e-5,
===================================================================
Cube29: same as cube28 but with
SEAICE_strength = 2.75e4
===================================================================
cube28: same as cube26 but with
viscC4Leith=1.5,
viscC4Leithd=1.5,
highOrderVorticity = .TRUE.
useAreaViscLength=.TRUE.
===================================================================
cube27: same as cube26 but with
viscC4Leith=1,
viscC4Leithd=1,
useAreaViscLength=.TRUE.,
SadournyCoriolis=.TRUE.,
useAbsVorticity=.TRUE.,
===================================================================
cube26 : same as cube25 with bug-fix to exf_bulkformulae
====================================================================
cube25 : same as cube24 but
viscar is increased back to 1e-3
completed on c15:
c15:menemenl/cube/MITgcm/verification/global_ocean.cs32x15/cube25
====================================================================
cube24 : same as cube22 but
limit ice surface temperature to -40C
MIN_ATEMP = -40
MIN_TICE = -40
====================================================================
cube23 : same as cube22 but with reduced viscosity
viscA4GridMin=5.e-5,
viscA4GridMax=1.e-3,
viscC4Leith=0.05,
deltaT=200,
after month 1, visca4 is of order 5e-10 and fairly homogeneous
and advcfl_w=6e-2, so switch to
viscA4GridMin=1.e-5,
after month 2, advcfl_w=1e-1, so switch to
viscA4GridMin=1.e-6,
viscC4Leith=0.02,
after month 3, advcfl_w=1e-1, 2e6<visca4<4e10 so switch to:
viscA4GridMin=1.e-7,
viscA4GridMax=1.e-2,
viscC4Leith=0.001,
above parameters crash so switch to
viscC4Leith=0.005,
above parameters crash after six months so switch to
viscA4GridMin=5.e-7,
viscA4GridMax=5.e-3,
viscC4Leith=0.005,
above parameters crash:
(PID.TID 0000.0001) %MON advcfl_wvel_max = 1.6543807544341E-01
(PID.TID 0000.0001) %MON advcfl_wvel_max = 1.0639473729765E-01
(PID.TID 0000.0001) %MON advcfl_wvel_max = 2.0187775095744E-01
(PID.TID 0000.0001) %MON advcfl_wvel_max = 2.6024594944192E-01
(PID.TID 0000.0001) %MON advcfl_wvel_max = 6.4082003867612E-01
so switch to:
# viscA4GridMin=5.e-7,
# viscA4GridMax=5.e-3,
viscC4Leith=0.01,
after month 9 switch back to
viscC4Leith=0.005,
fails again, so switch to
viscC4Leith=0.0075,
started on c9:
c9:menemenl/cube/MITgcm/verification/global_ocean.cs32x15/run
continuing on c20
c20:menemenl/cube/MITgcm/verification/global_ocean.cs32x15/run
continuing on c16
c16:menemenl/cube/MITgcm/verification/global_ocean.cs32x15/run
missing files
THETA 2004 month 7, 8 , 9
note that SALT is continued by SALTanom starting january 2001
====================================================================
cube22 : same as cube19 but starting from end of cube19
and with eddy diagnostics turned on
SEAICE_strength = 3.20e4
SEAICE_drag = 0.002
OCEAN_drag = 0.0002
SEAICE_waterDrag = 1
c20:menemenl/cube/MITgcm/verification/global_ocean.cs32x15/cube22
also on lou:/u/menemenl/cube/cube22
====================================================================
on c16 there is a 2-year integration, repeat of cube22
years 2003-2004 with hbl/mld diags for ME's project
/nobackup16/menemenl/cube/MITgcm/verification/global_ocean.cs32x15/cube22HblMld
being copied to
lou:/ecco/menemenl/cube22/cube22HblMld
====================================================================
cube21 : same as cube20 but starting from end of cube20
SEAICE_strength = 3.20e4
SEAICE_drag = 0.002
OCEAN_drag = 0.0005
SEAICE_waterDrag = 2.75
====================================================================
cube20diags :
repeat of cube20 with eddy diagnostics by Ed Hill on
columbia15:/nobackup15/menemenl/cube_eh3/MITgcm/verification/global_ocean.cs32x15/run
Note that out.2003 and out.2004 contain 6-hourly ETAN and PHIBOT output.
====================================================================
cube20 : same as cube16 but starting from Levitus
SEAICE_strength = 3.20e4
SEAICE_drag = 0.002
OCEAN_drag = 0.0005
SEAICE_waterDrag = 2.75
====================================================================
cube19 : same as cube16 but with further reduced ocean/water drag
SEAICE_strength = 3.20e4
SEAICE_drag = 0.002
OCEAN_drag = 0.0002
SEAICE_waterDrag = 1
====================================================================
cube18 : same as cube15 but with SEAICE_TEST_ICE_STRESS_1 and
with reduced SEAICE_drag and ocean/water drag
SEAICE_strength = 3.20e4
SEAICE_drag = 0.0002
OCEAN_drag = 0.0002
SEAICE_waterDrag = 1
====================================================================
cube17 : same as cube15 but with undef SEAICE_ALLOW_DYNAMICS
Crashes after 54 days with a very strange problem: see
c20:global_ocean.cs32x15/run17/out.1992/icecrash1.jpg
c20:menemenl/cube/MITgcm/verification/global_ocean.cs32x15/run17
====================================================================
cube16 : same as cube15 but with reduced ocean/water drag
SEAICE_strength = 3.20e4
SEAICE_drag = 0.002
OCEAN_drag = 0.0005
SEAICE_waterDrag = 2.75
c18:menemenl/cube/MITgcm/verification/global_ocean.cs32x15/run16
transferred output to lou
fails in 2003, day 120, time step 297864
====================================================================
cube15 : same as cube13 but with increased p*
SEAICE_strength = 3.20e4
SEAICE_drag = 0.002
OCEAN_drag = 0.001
SEAICE_waterDrag = 5.5
====================================================================
cube14 : same as cube13 but with
SEAICE_TEST_ICE_STRESS_1
and default seaice parameter values
1992-1994, integrates OK, but produces even larger Antarctic polynyas
====================================================================
cube13 : same as cube12, but all albedos increased by 0.1
# SEAICE_dryIceAlb = 0.75 _d 0
# SEAICE_wetIceAlb = 0.66 _d 0
# SEAICE_drySnowAlb = 0.84 _d 0
# SEAICE_wetSnowAlb = 0.7 _d 0
# SEAICE_waterAlbedo = 0.15
#
SEAICE_dryIceAlb = 0.85
SEAICE_wetIceAlb = 0.76
SEAICE_drySnowAlb = 0.94
SEAICE_wetSnowAlb = 0.8
SEAICE_waterAlbedo = 0.25
#
cube13 has crashed in 2001
so starting in 2001 switched to:
viscA4GridMin=3.e-4,
viscA4GridMax=4.e-3,
viscC4Leith=0.2,
====================================================================
cube12 : same as c6ncp10 but
started from end of 2002 c6ncp10
and viscC4Leith=0.2,
cube12 has crashed in 1999
so starting in 1999 switched to:
viscA4GridMin=3.e-4,
viscA4GridMax=4.e-3,
viscC4Leith=0.2,
====================================================================
c6ncp10 : SEAICE_freeze = -1.8
====================================================================
c6era01: same as cube 9 but with ERA surface boundary conditions.
====================================================================
Cube9: same as cube 6 but with Leith biharmonic viscosity.
====================================================================
Cube6: same as cube 5 but with yet smaller horizontal biharmonic
viscosity: 5e10 m^4/s.
====================================================================
Cube5: same as cube 3 but shortwave radiation and albedo have been
adjusted for a more realistic sea-ice extent.
====================================================================
Cube3: same as cube 2 but with half the horizontal biharmonic
viscosity: 2.4e11 m^4/s.
====================================================================
Cube2: 1992-2002, with sea-ice dynamics, surface salinity relaxation,
and lower vertical viscosity and diffusivity.
====================================================================
cube1: 1992-2002, no sea-ice dynamics, no surface salinity
relaxation. Horizontal biharmonic viscosity: 4.7e11 m^4/s.
====================================================================
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
failed tests:
=============
failed13a: same as cube13 but with
LSR_ERROR = 0.0001
SEAICE_freeze = -1.96
crashed after 51 days integration
half or 1/4 time step do not solve problem
notice spurious large anomalies in sflux and hflux from sea-ice model
=============
failed13b: same as cube13, with ice-modified stress:
SEAICE_TEST_ICE_STRESS_1
crashes after 183 days
/nobackup13/menemenl/cube/MITgcm/verification/global_ocean.cs32x15/failed13b
=============
failed13c: same as failed13b but with
LSR_ERROR = 0.0001
SEAICE_freeze = -1.96
crashes after 183 days
=============
failed13d: same as failed13c ???
crashes after 219 days
=============
failed16a : same as cube15 but with reduced SEAICE_drag
SEAICE_strength = 3.20e4
SEAICE_drag = 0.001
OCEAN_drag = 0.001
SEAICE_waterDrag = 5.5
fails at day 450
deleted
=============
failed17a : same as cube15 but with SEAICE_TEST_ICE_STRESS_1 and
SEAICE_strength = 3.20e4
SEAICE_drag = 0.001
OCEAN_drag = 0.001
SEAICE_waterDrag = 5.5
viscA4GridMin=3.e-4,
viscA4GridMax=4.e-3,
viscC4Leith=0.2,
crashed in 1997, day 2031, step 146232
--------------------------------------
restarted from beginning of 1997 with
viscA4GridMin=4.e-4,
viscA4GridMax=4.e-3,
viscC4Leith=0.2,
crashed in 1997, day 2061, step 148392
--------------------------------------
restarted from beginning of 1997 with
viscA4GridMin=5.e-4,
viscA4GridMax=5.e-3,
viscC4Leith=0.2,
crashed in 1997, day 2031, step 146232
--------------------------------------
restarted from beginning of 1997 with
viscA4GridMin=4.e-4,
viscA4GridMax=4.e-3,
viscC4Leith=0.25,
crashed in 1997, day 2031, step 146232
=============
failed17b : same as cube15 but with SEAICEuseDYNAMICS = .FALSE.
it crashed because I forgot to undef SEAICE_TEST_ICE_STRESS_1
=============
failed18a : same as cube17 but with reduced ocean and water drag
SEAICE_strength = 3.20e4
SEAICE_drag = 0.001
OCEAN_drag = 0.0005
SEAICE_waterDrag = 2.75
crashes after 219 days
=============
failed18b : same as cube17 but with increased SEAICE_drag and reduced ocean/water drag
SEAICE_strength = 3.20e4
SEAICE_drag = 0.002
OCEAN_drag = 0.0005
SEAICE_waterDrag = 2.75
fails at day 183
=============
failed18c : same as cube17 but with reduced SEAICE_drag and ocean/water drag
SEAICE_strength = 3.20e4
SEAICE_drag = 0.0005
OCEAN_drag = 0.0005
SEAICE_waterDrag = 2.75
crashed on day 468
====================================================================
ssh lou
cd cube
cp cube02/1992/HEFFtave.0000014688.data ../h2
cp cube03/out.1992/HEFFtave.0000014688.data ../h3
cp cube04/out.1992/HEFFtave.0000014688.data ../h4
cp cube05/out.1992/HEFFtave.0000014688.data ../h5
cp cube06/1992/HEFFtave.0000014688.data ../h6
cp cube09/1992/HEFFtave.0000014688.data ../h9
cp c6era01/1992/HEFFtave.0000014688.data ../h10
cp c6ncp10/1992/HEFFtave.0000014688.data ../h11
cp cube12/out.1992/HEFFtave.0000014688.data ../h12
cp cube13/out.1992/HEFFtave.0000014688.data ../h13
tl=6;
runs=[12:18 116 118];
h=zeros(510,510,length(runs)); n=0;
for i=runs, disp(i)
n=n+1;
tmp=readbin(['h' int2str(i)],[510 6 510]);
h(:,:,n)=squeeze(tmp(:,tl,:));
end
clf reset, orient landscape, wysiwyg
for i=1:length(runs)
subplot(3,3,i),mypcolor(h(:,:,i)');caxis([0 .01]),colorbar
title(runs(i))
end
%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%
ssh lou
cd cube
dmget cube12/out.1993/HEFFtave.0000041472.data
dmget cube13/out.1993/HEFFtave.0000041472.data
\cp cube12/out.1993/HEFFtave.0000041472.data ../h12.41472
\cp cube13/out.1993/HEFFtave.0000041472.data ../h13.41472
cd ..
ml
tl=6;
runs=[12 13 15 17];
h=zeros(510,510,length(runs)); n=0;
for i=runs, disp(i)
n=n+1;
tmp=readbin(['h' int2str(i) '.41472'],[510 6 510]);
h(:,:,n)=squeeze(tmp(:,tl,:));
end
clf reset, orient tall, wysiwyg
for i=1:length(runs)
subplot(2,2,i),mypcolor(h(:,:,i)');caxis([0 .01]),colorbar
title(runs(i))
end