Parent Directory
| 
Revision Log
| 
Revision Graph
| Links to HEAD: | (view) (download) (annotate) |
| Links to MAIN: | (view) (download) (annotate) |
| Sticky Tag: |
comment out filling of diagnostics SIuwind,SIvwind,SIfu,SIfv
add a store for hsnow, necessary after adding flag to add snow mass to seaiceMassC/U/V
make new code taf-able. There are still recomputations left that I do not understand, so bottom drag code is not included by default
- ECCO_OPTIONS.h is needed when including ecco_cost.h, ecco.h - AUTODIFF_OPTIONS.h is needed when including tamc.h, tamc_keys.h - CTRL_OPTIONS.h is needed when including ctrl.h, etc - pkg/seaice/seaice_cost*.F : clean up CPP brackets - SEAICE_SIZE.h : replace ALLOW_AUTODIFF_TAMC with ALLOW_AUTODIFF to avoid needing AUTODIFF_OPTIONS.h anytime SEAICE_SIZE.h is included (it seems that THSICE_SIZE.h, PTRACERS_SIZE.h have the same issue...)
- move most of "PART 1" of S/R seaice_growth into a new and separate
S/R seaice_reg_ridge that will both regularize the state variable after
advection call the actual ridging routine;
- requires to make 5 fields global fields: d_HEFFbyNEG,d_HSNWbyNEG,
d_AREAbyRLX,d_HEFFbyRLX,saltFluxAdjust;
- remove ridging code from seaice_advdiff
- clean up remaining part 1 of seaice_growth in an attempt to make it
more readable (this is probably subjective)
- finally: add new ridging routine for ITD code
- introduce tile exchange for HEFFITD, AREAITD, and HSNOWITD; - move calls to DIAGNOSTICS_FILL for SIheffN and SIareaN from seaice_model.F to seaice_diagnostics_state.F
avoid un-used variables
moving calls to SEAICE_ITD_REDIST and SEAICE_ITD_SUM from seaice_model to seaice_advdiff and seaice_growth because HEFF needs to be updated within these two routines
add uice/vice as uc/vc to the parameter list of s/r seaice_advdiff, similar to thsice_advdiff. Note that for the B-grid case the routine parameters uc/vc are overwritten.
allows to compile without pkg/exf
avoid unused variables
- removing all "ToM" comments - removing all SEAICE_DEBUG lines, which were introduced with ITD development
- implement new sequence of calls for thsice+seaice:
previously: ice-Dyn,ice-Advect,ice-Thermo(thsice)
new sequence: ice-Thermo(thsice),ice-Dyn,ice-Advect
- allows (with temporary CPP option "#define OLD_THSICE_CALL_SEQUENCE"
in CPP_OPTIONS.h) to recover old sequence
- remove legacy branch code. - retire SEAICE_GROWTH_LEGACY and SEAICE_CAP_HEFF accordingly.
remove variable "divergence" and associated code; contained artificial divergence rate for test of ITD code in 1-D test case; never meant to be permanently introduced
remove tabs ; remove un-used variables.
Step 1 of merging ice-thickness distribution (ITD) code from MITgcm_contrib/torge/itd/code/ to main repository (author: Torge Martin)
Prepare adjoint of SITRACER
- SItracer : added run time params, and allow coupling to ocn salinity.
new run time params in SEAICE_TRACER.h
--------------------------------------
SItrNumInUse - Number of tracers that are in use (must be less that SItrMaxNum)
SItrFromOcean0 - for new ice from ocean freeze - constant tracer value
SItrFromOceanFrac - for new ice from ocean freeze - fraction of ocean tracer
SItrFromFlood0 - for new ice from snow flood - constant tracer value
SItrFromFloodFrac - for new ice from snow flood - fraction of ocean tracer
SItrExpand0 - for ice cover thermo. exans. - constant tracer value
new run time params in SEAICE_PARAMS.h
--------------------------------------
SEAICE_salinityTracer :: use SItracer to exchange and trace ocean salt in ice
SEAICE_age Tracer :: use SItracer to trace the age of ice
simple propagation of SItrNumInUse in
-------------------------------------
seaice_diagnostics_init.F seaice_diagnostics_state.F
seaice_model.F seaice_advdiff.F
seaice_monitor.F seaice_output.F
isolate saltPlumeFlux array exch into a specific S/R called from do_oceanic_phys.F
- removal of the old way of seaice age tracer, which is now replaced by particular cases of SITRACER. - retired params : SEAICEadvAge, SEAICEadvSchAge, SEAICEdiffKhAge, IceAgeTrFile. - added to SITRACER : IceAgeTrFile, check pickups, monitor, output.
forgot his case in previous changes: needs to fill areaNm1 & hEffNm1 for S/R SEAICE_GROWTH (case SEAICE_GROWTH_LEGACY).
- remove YNEG, RIVER, TMIX from common block; - remove also areaNm1 & hEffNm1 when SEAICE_GROWTH_LEGACY is undef
avoid unused variables
- remove this patholigical case treatment in evolution branch, where it is un-necessary and breaks conservation.
Messing around with LSR adjoint
rename former OBCS_APPLY_UVICE to OBCS_ADJUST_UVICE
Add flag DISABLE_SEAICE_GROWTH - just for An!
- introducing ALLOW_SITRACER and seaice_tracer_phys.F to handle generic seaice tracer. For now it covers, and was tested for, salinity and age (work in progress). - introducing siEps (1e-5, parameter, defined in SEAICE_PARAMS.h).
comment out calls to S/R OBCS_COPY_TRACER (no longer needed with the the use of interior mask in advect & diffusion S/R)
Renaming of compile time flag SEAICE_SALINITY to SEAICE_VARIABLE_SALINITY and SEAICE_salinity to SIsalFRAC. Added logical checks for the use of the retired variable names in seaice_readparms and seaice_check. ---------------------------------------------------------------------- Modified Files: pkg/autodiff/autodiff_restore.F pkg/autodiff/autodiff_store.F pkg/obcs/obcs_apply_seaice.F pkg/seaice/SEAICE.h pkg/seaice/SEAICE_OPTIONS.h pkg/seaice/SEAICE_PARAMS.h pkg/seaice/seaice_ad_check_lev1_dir.h pkg/seaice/seaice_advdiff.F pkg/seaice/seaice_check.F pkg/seaice/seaice_diagnostics_state.F pkg/seaice/seaice_growth.F pkg/seaice/seaice_init_varia.F pkg/seaice/seaice_model.F pkg/seaice/seaice_monitor.F pkg/seaice/seaice_output.F pkg/seaice/seaice_read_pickup.F pkg/seaice/seaice_readparms.F pkg/seaice/seaice_write_pickup.F ----------------------------------------------------------------------
avoid un-used variables (placed within proper #ifdef)
Start a seaice tracer-like code but for now confinded to ice age
remove SEAICE_ALLOW_TD_IF pieces of code
switch some test on debugLevel value to debugMode test
call OBCS_COPY_TRACER for advected seaice variable (AREA, HEFF, HSNOW, HSALT and ICEAGE) (should fix Pb when using high-order scheme with OB near tile edge)
- move areaMin, hiceMin, areaMax to common blocks. - sort out diagnostics in seaice_growth.F. - practical approximations in seaice_growth.F adjoint. - reset TICE to celsius2K when no ice is present.
fix seaice model adjoint -- to be followed by updates in ajoint model test results by jmc
Update store directives for EVP
uses the right type of S/R for diagnostic filling
remove un-used variable
Changed initialisation and re-initialisations of seaiceMaskU/V
add "#ifdef SEAICE_CGRID / #endif" arround seaiceMaskU,V ; (this fixes offline_exf_seaice AD)
o initializations of *nm1 to break dependencies o adjust some store directives
third and step of replacing 3D versions of UICE,VICE,HEFF,AREA by 2D versions. Keep your fingers crossed.
o adding JM's doc/diags_changes.txt: a very useful summary that I refer to all the time o adding pkg/seaice diagnostics for atmospheric and oceanic heat flux and for atmospheric and oceanic ice growth
move diagnostic of Atmospheric Fresh-Water flux with other forcing diag (since it is not a state variable)
Adding stressDivergenceX/Y to store.
comment out un-used variables
fix wrong type of EXCH for "saltPlumeFlux" + change macros (EXCH & GLOBAL_SUM/MAX) sufix _R4/_R8 to _RS/_RL when applied to _RS/_RL variable
Update tamc.h for single-prec comlev option
added sea ice age tracer; #define SEAICE_AGE in SEAICE_PARAMS.h
Add Ian Fenty sea-ice thermodynamics modifications. Currently this is in separate routines to enable alternative testing. (relevant routines are ..._if.F) CPP switch is via #define SEAICE_ALLOW_TD_IF I.F. code doesnt contain latest SEAICE_SALINITY code
o Update stores after recent IF (DIFFERENT_MULTIPLE ... o re-initialise to break some dependencies (not sure if all correct) o turn some scalares into tw-dim.
o pkg/obcs: replace obcs_apply_area/heff/hsalt/hsnow with obcs_apply_seaice
o pkg/seaice: moved DIFFERENT_MULTIPLE from seaice_model to seaice_dynsolver http://forge.csail.mit.edu/pipermail/mitgcm-devel/2008-January/003173.html
Moved seaice tracer obcs after call to seaice_growth.
pkg/seaice: removed some unused halos
o pkg/seaice: removed SEAICE_FFIELDS.h and seaice_get_forcing.F seaice forcing fields can now be read only through pkg/exf
added check, readparms, and more diagnostice to pkg/salt_plume also changed package flag from useSaltPlume to useSALT_PLUME
added pkg/salt_plume by gathering code, which were previously spread around various files in model/inc and model/src results remain numerically identical to before, as a first step towards adding more options, etc., to this package
one more bug fix: wrong routine arguments
bug fixes - thx Manfredi
Added open boundary conditions capability for seaice UICE, and VICE UICE and VICE are reset at the edges after calling the solver and before advection/diffusion. Needs testing.
Committing An Nguyen's modifications to SALT_PLUME code. This includes addition of a saltPlumeFlux array to FFIELDS and of routine plumefrac.F
adding ice salinity HSALT as a prognostic variable
o pkg/seaice: added computation of saltFlux in seaice_growth; For time being, sea ice salinity is assumed constant contribution of snow flooding to freshening of sea ice is neglected. SEAICE_salinity is set in SEAICE_PARM01 of data.seaice and it defaults to 0.0 in order to maintain backward compatibility. Typical values for SEAICE_salinity are 4 to 10 g/kg.
Another flag was needed.
Bracket diagnostics calls according to package.
move diagnosing "non state variable" from seaice_diagnostics_fill to seaice_model
Seaice_model suddenly needs much more storing. Don't have time to check what has changed (again).
remove ALLOW_SEAICE from exf pkg files and add #define ALLOW_RUNOFF in SEAICE_OPTIONS.h
- fix CPP comments (was giving warnings). - fix EXCH of uwind,vwind for CS-grid : it does change the results of my cs32-grid test, but might have even larger impact on Lat-Lon-Cap grid.
fix previous modifs so that it compiles.
some parameter clarifications and a little bit of spring cleaning
add all needed EXF #ifdef and define them in SEAICE_OPTIONS.h if not using EXF
remove code which deals with advection of pkg/thsice variables
-remove seaice_exch and seaice_exch_uv -remove superfluous exchange of theta
add debug calls.
Updating seaice adjoint, step 1 (everything, except SEAICE_EVP).
overhaul of thermodynamics - growth.F and groatb.F are replaced by new routine seaice_growth.F - budget.F is replaced by two new routines seaice_budget_ocean/ice.F - move a few global fields out of SEAICE.h into seaice_growth.F and make them 2D (FICE/QNETO/ ...) - remove FO (it is the same as QNETO) - introduce a few local fields to avoid modifying external fields such as atemp, etc. lab_sea does not change, but hopefully it will be easier for Patrick do the adjoint
o add flooding algorithm and advection of HSNOW; both are turned off by default and do not fix any of the observed problems.
implement sea-ice loading (when using "RealFreshWaterFlux")
o modifications and additions necessary for coupling thsice to seaice
- major change in seaice_advdiff to make advection and diffusion
more transparent
- two new routines for mapping and computing advection and
diffusion, makes seaice_diffusion obsolete
One more STORE change.
Avoid recomputation call to seaice_advdiff to circumvene TAF bug (1.8.42)
First stab at SEAICE_CGRID.
- add c-grid version of the dynamic solver:
+ three new routines that are the c-grid substitute lsr, dynsolver
and ostres
+ put a few fields that were local to dynsolver into global common
blocks, so that I can move the computation of stresses etc into
seaice_lsr (saves coding but may break the adjoint; Patrick, I am
sorry!).
- replace more hardwired parameters by runtime parameters
- add ice masks that mask the rhs of the implicit solvers where there
is no ice (commented out in seaice_dynsolver, because i am not sure
if this works properly), eventually this should replace the clipping
of ice velocities in seaice_dynsolver to +/-40cm/s.
- correct units of turning angles (now they are really in degrees) - move computation of GWATX/Y out of main routine into dynsolver - initialize GWATX/Y in seaice_init
add a few new advection schemes to seaice: ENUM_UPWIND_1RST, ENUM_DST2, ENUM_FLUX_LIMIT, ENUM_DST3, ENUM_DST3_FLUX_LIMIT Default is still the old one
remove "baseTime" (no used) from arg. list of DIFF_BASE_MULTIPLE and rename it: DIFFERENT_MULTIPLE
use baseTime as time origin ; DIFF_BASE_MULTIPLE replaces DIFFERENT_MULTIPLE
o added seaice_summary.F and removed obsolete ALLOW_SEAICE's from pkg/seaice
o require forcing_In_AB=.FALSE. and remove surfaceTendencyTice from pkg/seaice - modified verification/lab_sea and natl_box_adjoint accordingly
Fooling around with ApiReference for pkg/kpp, pkg/sbo, and pkg/seaice
removed a debug print statement
- added SEAICE_deltaTtherm and SEAICE_deltaTdyn in order to reduce cost
of pkg/seaice. For time being, SEAICE_deltaTtherm = deltaTtracer
and SEAICE_deltaTdyn must be integer multiple of SEAICE_deltaTtherm
- changed proxy for geostrophic velocity and updated
verification/lab_sea/results/* accordingly
- changed proxy for geostrophic velocity and updated
verification/lab_sea/results/* accordingly
50 -> 50.0 for compiler compatibility
- changed proxy for geostrophic velocity and updated
verification/lab_sea/results/output.txt accordingly
o moved SEAICE_DO_DIAGS to S/R do_the_model_io o bug fix for pkg/seaice pickups
move "call seaice_write_pickup" to packages_write_pickup with other pkgs
added "exch_uv" exchanges to pkg/seaice
o modifications to make FREEZE flux visible to pkg/kpp
- moved surfaceTendencyTice from pkg/seaice to main code
- FREEZE moved to FORWARD_STEP
- subroutine FREEZE now limits only surface temperature
this means new output.txt for global_ocean.90x40x15,
global_ocean.cs32x15, and global_with_exf, but note
that results for these three experiments remain
bit-identical to before if allowFreezing=.FALSE.
o added surface flux output variables to TIMEAVE_STATVARS
o time-averaged output for pkg/ptracers
o merging from ecco-branch o updated/extended some store directives for seaice
o undid all of the cp51 checkin pending some ongoing code cleanups and discussion
o modifications to make FREEZE flux visible to pkg/kpp
- moved surfaceTendencyTice from pkg/seaice to main code
- FREEZE & EXTERNAL_FORCING_SURF moved to FORWARD_STEP
- subroutine FREEZE now limits only surface temperature
(this means new output.txt for global_ocean.90x40x15,
global_ocean.cs32x15, and global_with_exf)
o added surface flux output variables to TIMEAVE_STATVARS
o first check-in for the "branch-genmake2" merge
o verification suite as run on shelley (gcc 3.2.2):
Wed Oct 8 23:42:29 EDT 2003
T S U V
G D M c m s m s m s m s
E p a R g m m e . m m e . m m e . m m e .
N n k u 2 i a a d i a a d i a a d i a a d
2 d e n d n x n . n x n . n x n . n x n .
OPTFILE=NONE
Y Y Y Y 13 16 16 16 0 16 16 16 16 16 16 16 16 13 12 0 0 pass adjustment.128x64x1
Y Y Y Y 16 16 16 16 0 16 16 16 16 16 16 0 0 16 16 0 0 pass adjustment.cs-32x32x1
Y Y Y Y 16 16 16 16 0 16 16 16 16 16 16 22 0 16 16 22 0 pass adjust_nlfs.cs-32x32x1
Y Y Y Y -- 13 13 16 16 13 13 13 13 16 16 16 16 16 16 16 16 N/O advect_cs
Y Y Y Y -- 22 16 16 16 16 16 16 13 16 16 16 16 16 16 16 16 N/O advect_xy
Y Y Y Y -- 13 16 13 16 16 16 16 16 16 16 22 16 16 16 16 16 N/O advect_xz
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass aim.5l_cs
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 16 16 16 16 13 16 pass aim.5l_Equatorial_Channel
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 13 16 16 13 13 16 pass aim.5l_LatLon
Y Y Y Y 13 16 16 16 16 16 16 16 16 16 13 12 13 13 16 13 16 pass exp0
Y Y Y Y 14 16 16 16 16 16 16 16 22 16 16 16 13 16 16 22 16 pass exp1
Y Y Y Y 13 13 16 13 16 16 16 16 16 13 13 16 16 13 13 13 13 pass exp2
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass exp4
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 22 16 16 16 22 16 pass exp5
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass front_relax
Y Y Y Y 14 16 16 13 13 16 16 13 13 16 13 13 16 12 13 13 16 pass global_ocean.90x40x15
Y Y Y Y 10 16 16 13 13 16 13 16 16 13 13 13 13 16 16 13 16 FAIL global_ocean.cs32x15
Y Y Y Y 6 11 12 13 13 12 13 16 13 9 9 9 9 10 9 9 11 FAIL global_ocean_pressure
Y Y Y Y 14 16 16 13 16 16 16 13 13 13 13 13 16 12 16 13 16 pass global_with_exf
Y Y Y Y 14 16 16 16 16 16 16 16 16 11 13 22 13 16 16 9 16 pass hs94.128x64x5
Y Y Y Y 13 16 16 16 16 16 16 16 16 11 16 16 16 13 16 22 13 pass hs94.1x64x5
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 13 13 16 16 22 13 pass hs94.cs-32x32x5
Y Y Y Y 10 10 16 13 13 16 16 16 22 16 13 13 13 13 13 22 13 FAIL ideal_2D_oce
Y Y Y Y 8 16 16 16 16 16 16 16 16 13 13 8 16 16 16 16 16 FAIL internal_wave
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 13 22 13 13 13 22 16 pass inverted_barometer
Y Y Y Y 12 16 16 16 16 16 16 16 16 16 13 12 13 13 13 13 13 FAIL lab_sea
Y Y Y Y 11 16 16 16 16 16 16 16 13 13 13 12 13 16 13 12 13 FAIL natl_box
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass plume_on_slope
Y Y Y Y 13 16 16 16 16 13 16 16 16 16 16 16 16 13 16 16 16 pass solid-body.cs-32x32x1
Merging for c51 vs. e34
checkpoint50g_post
o merged with release1_p17 (pkg/seaice and verification/lab_sea)
- added SEAICE_MULTILEVEL for 8-category sea-ice thermodynamics
- LSR sea-ice dynamic solver moved to SouthWest B-grid location and
made the default because of faster convergence than ADI
checkpoint50c_pre Merging from release1_p13: o bug fix for pkg/seaice dynamic solver o Added SEAICE_initialHEFF to pkg/seaice
Merging from release1_p12: o Modifications for using pkg/exf with pkg/seaice - improved description of the various forcing configurations - added basic radiation bulk formulae to pkg/exf - units/sign fix for evap computation in exf_getffields.F - updated verification/global_with_exf/results/output.txt o Added pkg/sbo for computing IERS Special Bureau for the Oceans (SBO) core products, including oceanic mass, center-of-mass, angular, and bottom pressure (see pkg/sbo/README.sbo). o Lower bound for viscosity/diffusivity in pkg/kpp/kpp_routines.F to avoid negative values in shallow regions. - updated verification/natl_box/results/output.txt - updated verification/lab_sea/results/output.txt o MPI gather, scatter: eesupp/src/gather_2d.F and scatter_2d.F o Added useSingleCpuIO option (see PARAMS.h). o Updated useSingleCpuIO option in mdsio_writefield.F to work with multi-field files, e.g., for single-file pickup. o pkg/seaice: - bug fix in growth.F: QNET for no shortwave case - added HeffFile for specifying initial sea-ice thickness - changed SEAICE_EXTERNAL_FLUXES wind stress implementation o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf, kpp_transport_t.F, forward_step.F, and the_main_loop.F o pkg/seaice: - adjoint-friendly modifications - added a SEAICE_WRITE_PICKUP at end of the_model_main.F
checkpoint47f_post
Merging from release1_p10:
o modifications for using pkg/exf with pkg/seaice
- pkg/seaice CPP options SEAICE_EXTERNAL_FORCING
and SEAICE_EXTERNAL_FLUXES
- pkg/exf CPP options EXF_READ_EVAP and
EXF_NO_BULK_COMPUTATIONS
- usage examples are Experiments 8 and 9 in
verification/lab_sea/README
- verification/lab_sea default experiment now uses
pkg/gmredi, pkg/kpp, pkg/seaice, and pkg/exf
checkpoint47b_post
Merging from release1_p9:
o pkg/seaice
- removed GOTO's and added taf directives
- double precision constants to reduce the g77 (Linux)
to F77 (SGI) differences reported in release1_p8
o tools/genmake
- added SGI options
o verification/testscript
- updated to that of checkpoint47a_post
o verification/global_ocean.90x40x15/input/eedata
- modified for SGI f77 compatibility
o verification/lab_sea
- added description of sea-ice model
- added missing matlab routines
- added test of thermodynamics parallelization
Modified Files:
doc/tag-index pkg/seaice/SEAICE_FFIELDS.h
pkg/seaice/SEAICE_PARAMS.h pkg/seaice/adi.F
pkg/seaice/advect.F pkg/seaice/budget.F pkg/seaice/diffus.F
pkg/seaice/dynsolver.F pkg/seaice/groatb.F pkg/seaice/growth.F
pkg/seaice/lsr.F pkg/seaice/ostres.F
pkg/seaice/seaice_do_diags.F pkg/seaice/seaice_get_forcing.F
pkg/seaice/seaice_init.F pkg/seaice/seaice_model.F
pkg/seaice/seaice_readparms.F tools/genmake
verification/global_ocean.90x40x15/input/eedata
verification/lab_sea/README
verification/lab_sea/matlab/lookat_exp1.m
verification/lab_sea/matlab/lookat_exp2.m
verification/lab_sea/matlab/lookat_exp3.m
verification/lab_sea/matlab/lookat_exp4.m
verification/lab_sea/matlab/lookat_exp5.m
verification/lab_sea/matlab/lookat_exp6.m
verification/lab_sea/results/AREAtave.0000000010.data
verification/lab_sea/results/HEFFtave.0000000010.data
verification/lab_sea/results/UICEtave.0000000010.data
verification/lab_sea/results/VICEtave.0000000010.data
verification/lab_sea/results/output.txt
Added Files:
verification/lab_sea/seaice.ps
verification/lab_sea/matlab/lookat_exp7.m
verification/lab_sea/matlab/mmax.m
verification/lab_sea/matlab/mypcolor.m
verification/lab_sea/matlab/myquiver.m
verification/lab_sea/matlab/readbin.m
verification/lab_sea/matlab/wysiwyg.m
Removed Files:
verification/lab_sea/code/KPP_OPTIONS.h
Merging from release1_p8 branch: o New package: pkg/seaice Sea ice model by D. Menemenlis (JPL) and Jinlun Zhang (Seattle). The sea-ice code is based on Hibler (1979-1980). Two sea-ice dynamic solvers, ADI and LSR, are included. In addition to computing prognostic sea-ice variables and diagnosing the forcing/external data fields that drive the ocean model, SEAICE_MODEL also sets theta to the freezing point under sea-ice. The implied surface heat flux is then stored in variable surfaceTendencyTice, which is needed by KPP package (kpp_calc.F and kpp_transport_t.F) to diagnose surface buoyancy fluxes and for the non-local transport term. Because this call precedes model thermodynamics, temperature under sea-ice may not be "exactly" at the freezing point by the time theta is dumped or time-averaged.
file seaice_model.F was initially added on branch release1.
This form allows you to request diffs between any two revisions of this file. For each of the two "sides" of the diff, select a symbolic revision name using the selection box, or choose 'Use Text Field' and enter a numeric revision.
| ViewVC Help | |
| Powered by ViewVC 1.1.22 |