[Dart-dev] DART/branches Revision: 11702

[dart at ucar.edu]

thoar at ucar.edu
2017-06-07 21:01:46 -0600 (Wed, 07 Jun 2017)
331
These files reproduce what the CICE results Yongfei generated with 
with an earlier version of the rma_cice branch (which we pruned a while ago).
The results were generated on yellowstone. 
Have yet to test Cheyenne.
Only tested the 'simple squeeze' method in dart_to_cice.
The test case is in /glade/p/image/DART_test_cases/cice




Modified: DART/branches/rma_trunk/models/cice/dart_to_cice.f90
===================================================================
--- DART/branches/rma_trunk/models/cice/dart_to_cice.f90	2017-06-07 21:45:41 UTC (rev 11701)
+++ DART/branches/rma_trunk/models/cice/dart_to_cice.f90	2017-06-08 03:01:46 UTC (rev 11702)
@@ -33,8 +33,13 @@
 !------------------------------------------------------------------
 
 character(len=256) :: dart_to_cice_input_file = 'dart_restart.nc'
+character(len=256) :: original_cice_input_file = 'cice_restart.nc'
+character(len=256) :: previous_cice_input_file = 'pre_restart.nc'
 character(len=128) :: balance_method = 'simple_squeeze'
-namelist /dart_to_cice_nml/ dart_to_cice_input_file, balance_method
+namelist /dart_to_cice_nml/ dart_to_cice_input_file, & 
+                            original_cice_input_file, &
+                            previous_cice_input_file, &
+                            balance_method
 
 character(len=512) :: string1, string2, msgstring
 character(len=15)  :: varname
@@ -42,19 +47,87 @@
 
 integer :: Nx, Ny
 integer :: Ncat   ! number of categories in ice-thickness dist
+integer, parameter :: Nilyr = 8   ! number of layers in ice, hardwired
+integer, parameter :: Nslyr = 8   ! number of layers in snow, hardwired
 
+real(r8), allocatable :: aicen_original(:,:,:)
+real(r8), allocatable :: vicen_original(:,:,:)
+real(r8), allocatable :: vsnon_original(:,:,:)
+real(r8), allocatable :: aice_original(:,:)
+!real(r8), allocatable :: vice_original(:,:)
+!real(r8), allocatable :: vsno_original(:,:)
+real(r8), allocatable :: hicen_original(:,:,:)
+real(r8), allocatable :: hsnon_original(:,:,:)
+
+real(r8), allocatable :: aicen_pre(:,:,:)
+real(r8), allocatable :: vicen_pre(:,:,:)
+real(r8), allocatable :: vsnon_pre(:,:,:)
+real(r8), allocatable :: aice_pre(:,:)
+!real(r8), allocatable :: vice_pre(:,:)
+!real(r8), allocatable :: vsno_pre(:,:)
+
 real(r8), allocatable :: aicen(:,:,:)
 real(r8), allocatable :: vicen(:,:,:)
 real(r8), allocatable :: vsnon(:,:,:)
-real(r8), allocatable ::  aice(:,:)
+real(r8), allocatable :: Tsfcn(:,:,:)
+real(r8), allocatable :: sice001(:,:,:)
+real(r8), allocatable :: sice002(:,:,:)
+real(r8), allocatable :: sice003(:,:,:)
+real(r8), allocatable :: sice004(:,:,:)
+real(r8), allocatable :: sice005(:,:,:)
+real(r8), allocatable :: sice006(:,:,:)
+real(r8), allocatable :: sice007(:,:,:)
+real(r8), allocatable :: sice008(:,:,:)
+real(r8), allocatable :: qice001(:,:,:)
+real(r8), allocatable :: qice002(:,:,:)
+real(r8), allocatable :: qice003(:,:,:)
+real(r8), allocatable :: qice004(:,:,:)
+real(r8), allocatable :: qice005(:,:,:)
+real(r8), allocatable :: qice006(:,:,:)
+real(r8), allocatable :: qice007(:,:,:)
+real(r8), allocatable :: qice008(:,:,:)
+real(r8), allocatable :: qsno001(:,:,:)
+real(r8), allocatable :: qsno002(:,:,:)
+real(r8), allocatable :: qsno003(:,:,:)
+real(r8), allocatable :: aice(:,:)
+real(r8), allocatable :: vice(:,:)
+real(r8), allocatable :: vsno(:,:)
 
-integer  :: i, j, k
-integer  :: VarID, ncid, iunit, io, ndims
+real(r8), allocatable :: aice_temp(:,:)
+real(r8), allocatable :: increment_aice(:,:)
+!real(r8), allocatable :: increment_vice(:,:)
+!real(r8), allocatable :: increment_vsno(:,:)
+
+real(r8), allocatable :: tendency_aice(:,:)
+!real(r8), allocatable :: tendency_vice(:,:)
+!real(r8), allocatable :: tendency_vsno(:,:)
+real(r8), allocatable :: tendency_aicen(:,:,:)
+!real(r8), allocatable :: tendency_vicen(:,:,:)
+!real(r8), allocatable :: tendency_vsnon(:,:,:)
+
+real(r8) :: R, weight_aicen     !, weight_vicen, weight_vsnon
+
+integer  :: i, j, k, n
+integer  :: VarID, ncid,ncid2, iunit, io, ndims
 real(r8) :: squeeze
 
-! If nothing needs to be rebalanced, no need to write to file.
-logical  :: update_restart = .false.
 
+real(r8), parameter :: &      ! from ice_therm_vertical.F90
+     phi_init = 0.75_r8, &    ! initial liquid fraction of frazil ice 
+     dSin0_frazil = 3.0_r8    ! bulk salinity reduction of newly formed frazil 
+
+real(r8), parameter :: &      ! from pop constants.F90
+     sss = 34.7_r8  !  use ocn_ref_salinity since we dont know truth
+
+real(r8), parameter :: &