===================================================================
Cube84
same as cube81 but
useRelativeWind=.FALSE.
===================================================================
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
|
