Skip to content

ACCESS-CM3 windstress coupling #48

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
kieranricardo merged 9 commits into from
Aug 12, 2026
Merged

ACCESS-CM3 windstress coupling #48

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
add4ba1
change field names and remove unused fields
kieranricardo Mar 16, 2026
ac759db
add in sea ice wind stress, delete unused atmosphere state fields
kieranricardo Mar 30, 2026
6613241
fix import
kieranricardo Mar 30, 2026
fcbe445
scale wind stress by ice fraction
kieranricardo Apr 14, 2026
3412483
update wind stress name
kieranricardo May 4, 2026
7b1121a
change icepack to head of access-cm3-Icepack-1.5.1
kieranricardo Jul 13, 2026
d21b801
add missing halo update
kieranricardo Jul 14, 2026
a19b0e5
angle sign fix
kieranricardo Jul 14, 2026
708d590
add HALO updates
kieranricardo Aug 12, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 4 additions & 0 deletions cicecore/cicedyn/dynamics/ice_dyn_eap.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -264,6 +264,10 @@ subroutine eap (dt)
file=__FILE__, line=__LINE__)

if (.not. calc_strair) then
call ice_HaloUpdate (strax, halo_info, &
field_loc_center, field_type_vector)
call ice_HaloUpdate (stray, halo_info, &
field_loc_center, field_type_vector)
call grid_average_X2Y('F', strax, grid_atm_dynu, strairxU, 'U')
call grid_average_X2Y('F', stray, grid_atm_dynv, strairyU, 'U')
else
Expand Down
4 changes: 4 additions & 0 deletions cicecore/cicedyn/dynamics/ice_dyn_evp.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -461,6 +461,10 @@ subroutine evp (dt)
file=__FILE__, line=__LINE__)

if (.not. calc_strair) then
call ice_HaloUpdate (strax, halo_info, &
field_loc_center, field_type_vector)
call ice_HaloUpdate (stray, halo_info, &
field_loc_center, field_type_vector)
call grid_average_X2Y('F', strax, grid_atm_dynu, strairxU, 'U')
call grid_average_X2Y('F', stray, grid_atm_dynv, strairyU, 'U')
else
Expand Down
4 changes: 4 additions & 0 deletions cicecore/cicedyn/dynamics/ice_dyn_vp.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -318,6 +318,10 @@ subroutine implicit_solver (dt)
file=__FILE__, line=__LINE__)

if (.not. calc_strair) then
call ice_HaloUpdate (strax, halo_info, &
field_loc_center, field_type_vector)
call ice_HaloUpdate (stray, halo_info, &
field_loc_center, field_type_vector)
call grid_average_X2Y('F', strax, grid_atm_dynu, strairxU, 'U')
call grid_average_X2Y('F', stray, grid_atm_dynv, strairyU, 'U')
else
Expand Down
2 changes: 1 addition & 1 deletion cicecore/cicedyn/infrastructure/ice_grid.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2094,7 +2094,7 @@ subroutine mom_grid_rotation_angle(lon_cnr, lat_cnr, lon_cen, angle)
+ lon_adj
enddo ; enddo
lon_scale = cos(p25*(lat_cnr(I,J) + lat_cnr(I+1,J+1) + lat_cnr(I+1,J) + lat_cnr(I,J+1)))
angle(i,j) = atan2(lon_scale*((lonB(1,2) - lonB(2,1) + lonB(2,2) - lonB(1,1))), &
angle(i,j) = -atan2(lon_scale*((lonB(1,2) - lonB(2,1) + lonB(2,2) - lonB(1,1))), &
(lat_cnr(I,J+1) - lat_cnr(I+1,J) + lat_cnr(I+1,J+1) - lat_cnr(I,J)) )
enddo
enddo
Expand Down
205 changes: 21 additions & 184 deletions cicecore/drivers/access/cmeps/ice_import_export.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -20,8 +20,7 @@ module ice_import_export
use ice_flux_bgc , only : faero_atm, faero_ocn
use ice_flux_bgc , only : fiso_atm, fiso_ocn, fiso_evap
use ice_flux_bgc , only : Qa_iso, Qref_iso, HDO_ocn, H2_18O_ocn, H2_16O_ocn
use ice_flux , only : fresh_ai, fsalt_ai, zlvl, uatm, vatm, potT, Tair, Qa
use ice_flux , only : rhoa, swvdr, swvdf, swidr, swidf, flw, frain
use ice_flux , only : fresh_ai, fsalt_ai, flw, frain, strax, stray
use ice_flux , only : fsnow, uocn, vocn, sst, ss_tltx, ss_tlty, frzmlt
use ice_flux , only : send_i2x_per_cat
use ice_flux , only : sss, Tf, wind, fsw
Expand Down Expand Up @@ -187,21 +186,10 @@ subroutine ice_advertise_fields(gcomp, importState, exportState, flds_scalar_nam
end if

! from atmosphere
call fldlist_add(fldsToIce_num, fldsToIce, 'Sa_z' )
call fldlist_add(fldsToIce_num, fldsToIce, 'Sa_u' )
call fldlist_add(fldsToIce_num, fldsToIce, 'Sa_v' )
call fldlist_add(fldsToIce_num, fldsToIce, 'Sa_shum' )
call fldlist_add(fldsToIce_num, fldsToIce, 'Sa_tbot' )
call fldlist_add(fldsToIce_num, fldsToIce, 'Sa_pbot' )
call fldlist_add(fldsToIce_num, fldsToIce, 'Faxa_swvdr' )
call fldlist_add(fldsToIce_num, fldsToIce, 'Faxa_swvdf' )
call fldlist_add(fldsToIce_num, fldsToIce, 'Faxa_swndr' )
call fldlist_add(fldsToIce_num, fldsToIce, 'Faxa_swndf' )
call fldlist_add(fldsToIce_num, fldsToIce, 'Faxa_lwdn' )
call fldlist_add(fldsToIce_num, fldsToIce, 'Faxa_rain' )
call fldlist_add(fldsToIce_num, fldsToIce, 'Faxa_snow' )
call fldlist_add(fldsToIce_num, fldsToIce, 'Sa_ptem' ) !cesm
call fldlist_add(fldsToIce_num, fldsToIce, 'Sa_dens' ) !cesm
call fldlist_add(fldsToIce_num, fldsToIce, 'Faia_taux' )
call fldlist_add(fldsToIce_num, fldsToIce, 'Faia_tauy' )

! the following are advertised but might not be connected if they are not present
! in the cmeps esmFldsExchange_xxx_mod.F90 that is model specific
Expand Down Expand Up @@ -260,15 +248,6 @@ subroutine ice_advertise_fields(gcomp, importState, exportState, flds_scalar_nam
call fldlist_add(fldsFrIce_num, fldsFrIce, 'Si_floediam' )
end if

! ice/atm fluxes computed by ice
call fldlist_add(fldsFrIce_num, fldsFrIce, 'Faii_taux' )
call fldlist_add(fldsFrIce_num, fldsFrIce, 'Faii_tauy' )
call fldlist_add(fldsFrIce_num, fldsFrIce, 'Faii_lat' )
call fldlist_add(fldsFrIce_num, fldsFrIce, 'Faii_sen' )
call fldlist_add(fldsFrIce_num, fldsFrIce, 'Faii_lwup' )
call fldlist_add(fldsFrIce_num, fldsFrIce, 'Faii_evap' )
call fldlist_add(fldsFrIce_num, fldsFrIce, 'Faii_swnet' )

! ice/ocn fluxes computed by ice
call fldlist_add(fldsFrIce_num, fldsFrIce, 'Fioi_melth' )
call fldlist_add(fldsFrIce_num, fldsFrIce, 'Fioi_swpen' )
Expand Down Expand Up @@ -452,7 +431,7 @@ subroutine ice_import( importState, rc )
integer , intent(out) :: rc

! local variables
integer,parameter :: nflds=16
integer,parameter :: nflds=5
integer,parameter :: nfldv=6
integer :: i, j, iblk, n, k
integer :: ilo, ihi, jlo, jhi !beginning and end of physical domain
Expand Down Expand Up @@ -518,63 +497,19 @@ subroutine ice_import( importState, rc )
call state_getimport(importState, 'So_s', output=aflds, index=2, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

! import atm states

call state_getimport(importState, 'Sa_z', output=aflds, index=3, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

if (State_FldChk(importState, 'Sa_ptem') .and. State_fldchk(importState, 'Sa_dens')) then
call state_getimport(importState, 'Sa_ptem', output=aflds, index=4, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return
call state_getimport(importState, 'Sa_dens', output=aflds, index=5, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return
else if (State_FldChk(importState, 'Sa_pbot')) then
call state_getimport(importState, 'Sa_pbot', output=aflds, index=6, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return
else
call abort_ice(trim(subname)//&
": ERROR either Sa_ptem and Sa_dens OR Sa_pbot must be in import state")
end if

call state_getimport(importState, 'Sa_tbot', output=aflds, index=7, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

call state_getimport(importState, 'Sa_shum', output=aflds, index=8, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

! import ocn/ice fluxes

call state_getimport(importState, 'Fioo_q', output=aflds, index=9, &
call state_getimport(importState, 'Fioo_q', output=aflds, index=3, &
areacor=med2mod_areacor, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

! import atm fluxes

call state_getimport(importState, 'Faxa_swvdr', output=aflds, index=10, &
areacor=med2mod_areacor, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

call state_getimport(importState, 'Faxa_swndr', output=aflds, index=11, &
call state_getimport(importState, 'Faxa_rain', output=aflds, index=4, &
areacor=med2mod_areacor, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

call state_getimport(importState, 'Faxa_swvdf', output=aflds, index=12, &
areacor=med2mod_areacor, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

call state_getimport(importState, 'Faxa_swndf', output=aflds, index=13, &
areacor=med2mod_areacor, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

call state_getimport(importState, 'Faxa_lwdn', output=aflds, index=14, &
areacor=med2mod_areacor, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

call state_getimport(importState, 'Faxa_rain', output=aflds, index=15, &
areacor=med2mod_areacor, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

call state_getimport(importState, 'Faxa_snow', output=aflds, index=16, &
call state_getimport(importState, 'Faxa_snow', output=aflds, index=5, &
areacor=med2mod_areacor, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

Expand Down Expand Up @@ -618,37 +553,16 @@ subroutine ice_import( importState, rc )
end do
!$OMP END PARALLEL DO

! flatn_f = - Foxx_evap(i,j,cat,k) * Lsub !! latent heat -
! fcondtopn_f = botmelt
! fsurfn_f (:,:,cat,:) = topmelt(:,:,cat,:) + botmelt(:,:,cat,:)
! if (um_tsfice(i,j,cat,k) > 0.0) then
! trcrn(i,j,nt_Tsfc,cat,k) = 0.0
! else if (um_tsfice(i,j,cat,k) < -60.0) then
! trcrn(i,j,nt_Tsfc,cat,k) = -60.0
! else
! trcrn(i,j,nt_Tsfc,cat,k) = um_tsfice(i,j,cat,k)
! endif

!$OMP PARALLEL DO PRIVATE(iblk,i,j)
do iblk = 1, nblocks
do j = 1,ny_block
do i = 1,nx_block
sst (i,j,iblk) = aflds(i,j, 1,iblk)
sss (i,j,iblk) = aflds(i,j, 2,iblk)
zlvl (i,j,iblk) = aflds(i,j, 3,iblk)
! see below for 4,5,6
Tair (i,j,iblk) = aflds(i,j, 7,iblk)
Qa (i,j,iblk) = aflds(i,j, 8,iblk)
frzmlt (i,j,iblk) = aflds(i,j, 9,iblk)
swvdr(i,j,iblk) = 0.0 ! aflds(i,j,10,iblk)
swidr(i,j,iblk) = 0.0 !aflds(i,j,11,iblk)
swvdf(i,j,iblk) = 0.0 !aflds(i,j,12,iblk)
swidf(i,j,iblk) = 0.0 !aflds(i,j,13,iblk)
flw (i,j,iblk) = 0.0 ! aflds(i,j,14,iblk)
frain(i,j,iblk) = aflds(i,j,15,iblk)
fsnow(i,j,iblk) = aflds(i,j,16,iblk)
! strax !! windstress - already handled, come back to this
! stray !! windstress
frzmlt (i,j,iblk) = aflds(i,j, 3,iblk)
frain(i,j,iblk) = aflds(i,j,4,iblk)
fsnow(i,j,iblk) = aflds(i,j,5,iblk)

end do
end do
Expand Down Expand Up @@ -678,40 +592,6 @@ subroutine ice_import( importState, rc )
end do
end if

if ( State_fldChk(importState, 'Sa_ptem') .and. State_fldchk(importState,'Sa_dens')) then
!$OMP PARALLEL DO PRIVATE(iblk,i,j)
do iblk = 1, nblocks
do j = 1,ny_block
do i = 1,nx_block
potT (i,j,iblk) = aflds(i,j, 4,iblk)
rhoa (i,j,iblk) = aflds(i,j, 5,iblk)
end do
end do
end do
!$OMP END PARALLEL DO
else if (State_fldChk(importState, 'Sa_pbot')) then
!$OMP PARALLEL DO PRIVATE(iblk,i,j)
do iblk = 1, nblocks
do j = 1,ny_block
do i = 1,nx_block
inst_pres_height_lowest = aflds(i,j,6,iblk)
if (inst_pres_height_lowest > 0.0_ESMF_KIND_R8) then
potT (i,j,iblk) = Tair(i,j,iblk) * (100000._ESMF_KIND_R8/inst_pres_height_lowest)**0.286_ESMF_KIND_R8
else
potT (i,j,iblk) = 0.0_ESMF_KIND_R8
end if
if (Tair(i,j,iblk) /= 0._ESMF_KIND_R8) then
rhoa(i,j,iblk) = inst_pres_height_lowest / &
(287.058_ESMF_KIND_R8*(1._ESMF_KIND_R8+0.608_ESMF_KIND_R8*Qa(i,j,iblk))*Tair(i,j,iblk))
else
rhoa(i,j,iblk) = 1.2_ESMF_KIND_R8
endif
end do !i
end do !j
end do !iblk
!$OMP END PARALLEL DO
end if

deallocate(aflds)
allocate(aflds(nx_block,ny_block,nfldv,nblocks))
aflds = c0
Expand All @@ -722,10 +602,9 @@ subroutine ice_import( importState, rc )
if (ChkErr(rc,__LINE__,u_FILE_u)) return
call state_getimport(importState, 'So_v', output=aflds, index=2, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

call state_getimport(importState, 'Sa_u', output=aflds, index=3, rc=rc)
call state_getimport(importState, 'Faia_taux', output=aflds, index=3, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return
call state_getimport(importState, 'Sa_v', output=aflds, index=4, rc=rc)
call state_getimport(importState, 'Faia_tauy', output=aflds, index=4, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

call state_getimport(importState, 'So_dhdx', output=aflds, index=5, rc=rc)
Expand All @@ -746,8 +625,8 @@ subroutine ice_import( importState, rc )
do i = 1,nx_block
uocn (i,j,iblk) = aflds(i,j, 1,iblk)
vocn (i,j,iblk) = aflds(i,j, 2,iblk)
uatm (i,j,iblk) = aflds(i,j, 3,iblk)
vatm (i,j,iblk) = aflds(i,j, 4,iblk)
strax (i,j,iblk) = aflds(i,j, 3,iblk) * aice(i, j, iblk)
stray (i,j,iblk) = aflds(i,j, 4,iblk) * aice(i, j, iblk)
ss_tltx(i,j,iblk) = aflds(i,j, 5,iblk)
ss_tlty(i,j,iblk) = aflds(i,j, 6,iblk)
enddo !i
Expand Down Expand Up @@ -930,16 +809,13 @@ subroutine ice_import( importState, rc )
do i = 1, nx_block

! atmosphere
workx = uatm(i,j,iblk) ! wind velocity, m/s
worky = vatm(i,j,iblk)
uatm (i,j,iblk) = workx*cos(ANGLET(i,j,iblk)) & ! convert to POP grid
+ worky*sin(ANGLET(i,j,iblk)) ! note uatm, vatm, wind
vatm (i,j,iblk) = worky*cos(ANGLET(i,j,iblk)) & ! are on the T-grid here
workx = strax(i,j,iblk) ! wind stress
worky = stray(i,j,iblk)
strax (i,j,iblk) = workx*cos(ANGLET(i,j,iblk)) & ! convert to POP grid
+ worky*sin(ANGLET(i,j,iblk)) ! note strax and stray
stray (i,j,iblk) = worky*cos(ANGLET(i,j,iblk)) & ! are on the T-grid here
- workx*sin(ANGLET(i,j,iblk))

wind (i,j,iblk) = sqrt(uatm(i,j,iblk)**2 + vatm(i,j,iblk)**2)
fsw (i,j,iblk) = swvdr(i,j,iblk) + swvdf(i,j,iblk) &
+ swidr(i,j,iblk) + swidf(i,j,iblk)
enddo
enddo
enddo
Expand Down Expand Up @@ -1212,45 +1088,6 @@ subroutine ice_export(importState, exportState, cpl_dt, rc )
if (ChkErr(rc,__LINE__,u_FILE_u)) return
end if

! ------
! ice/atm fluxes computed by ice
! ------

! Zonal air/ice stress
call state_setexport(exportState, 'Faii_taux' , input=tauxa, lmask=tmask, ifrac=ailohi, &
areacor=mod2med_areacor, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

! Meridional air/ice stress
call state_setexport(exportState, 'Faii_tauy' , input=tauya, lmask=tmask, ifrac=ailohi, &
areacor=mod2med_areacor, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

! Latent heat flux (atm into ice)
call state_setexport(exportState, 'Faii_lat' , input=flat, lmask=tmask, ifrac=ailohi, &
areacor=mod2med_areacor, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

! Sensible heat flux (atm into ice)
call state_setexport(exportState, 'Faii_sen' , input=fsens, lmask=tmask, ifrac=ailohi, &
areacor=mod2med_areacor, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

! longwave outgoing (upward), average over ice fraction only
call state_setexport(exportState, 'Faii_lwup' , input=flwout, lmask=tmask, ifrac=ailohi, &
areacor=mod2med_areacor, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

! Evaporative water flux (kg/m^2/s)
call state_setexport(exportState, 'Faii_evap' , input=evap, lmask=tmask, ifrac=ailohi, &
areacor=mod2med_areacor, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

! Shortwave flux absorbed in ice and ocean (W/m^2)
call state_setexport(exportState, 'Faii_swnet' , input=fswabs, lmask=tmask, ifrac=ailohi, &
areacor=mod2med_areacor, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

! ------
! ice/ocn fluxes computed by ice
! ------
Expand Down Expand Up @@ -2135,7 +1972,7 @@ subroutine ice_advertise_fields_access_export(gcomp, importState, exportState, f
call fldlist_add(fldsFrIce_num , fldsFrIce, 'Si_topk', ungridded_lbound=1, ungridded_ubound=ncat) ! from ice flux: keffn_top
call fldlist_add(fldsFrIce_num , fldsFrIce, 'Si_pndf_n', ungridded_lbound=1, ungridded_ubound=ncat) ! from icepack_shorwave: apeffn
call fldlist_add(fldsFrIce_num , fldsFrIce, 'Si_pndt_n', ungridded_lbound=1, ungridded_ubound=ncat) ! from ice state field: trcrn
call fldlist_add(fldsFrIce_num , fldsFrIce, 'sstfrz')
call fldlist_add(fldsFrIce_num , fldsFrIce, 'Si_Tf')

write (tmpString, *) ncat
call ESMF_LogWrite("CICE number of ice categories: " // trim(tmpString))
Expand Down Expand Up @@ -2193,7 +2030,7 @@ subroutine ice_export_access(importState, exportState, ailohi, cpl_dt, rc)
call state_setexport(exportState, 'Si_vice_n', input=vicen , lmask=tmask, ifrac=ailohi, rc=rc, index=n, ungridded_index=n)
end do

call state_setexport(exportState, 'sstfrz', input=Tf , lmask=tmask, ifrac=ailohi, rc=rc)
call state_setexport(exportState, 'Si_Tf', input=Tf , lmask=tmask, ifrac=ailohi, rc=rc)

! To conserve pond areas, scale by ice fractions before mapping. Unscale in the atmosphere after mapping
pndfn_scaled(:,:,:,:) = apeffn(:,:,:,:) * aicen(:,:,:,:)
Expand Down
2 changes: 1 addition & 1 deletion icepack
Open in desktop
Submodule icepack updated 33 files
+7 −13 .zenodo.json
+ − COPYRIGHT.pdf
+2 −1 columnphysics/icepack_algae.F90
+11 −12 columnphysics/icepack_flux.F90
+15 −12 columnphysics/icepack_fsd.F90
+1 −1 columnphysics/icepack_intfc.F90
+6 −22 columnphysics/icepack_itd.F90
+40 −71 columnphysics/icepack_parameters.F90
+63 −63 columnphysics/icepack_shortwave.F90
+0 −3 columnphysics/icepack_snow.F90
+57 −64 columnphysics/icepack_therm_bl99.F90
+10 −3 columnphysics/icepack_therm_itd.F90
+34 −38 columnphysics/icepack_therm_vertical.F90
+13 −19 columnphysics/icepack_wavefracspec.F90
+4 −2 columnphysics/icepack_zbgc_shared.F90
+1 −1 columnphysics/version.txt
+1 −1 configuration/driver/icedrv_MAIN.F90
+0 −1 configuration/driver/icedrv_flux.F90
+8 −8 configuration/driver/icedrv_forcing.F90
+12 −34 configuration/driver/icedrv_init.F90
+13 −12 configuration/driver/icedrv_step.F90
+0 −3 configuration/scripts/icepack_in
+0 −1 configuration/scripts/options/set_nml.fsd1
+0 −1 configuration/scripts/options/set_nml.fsd12
+4 −3 configuration/scripts/tests/poll_queue.csh
+3 −3 doc/source/conf.py
+4 −6 doc/source/icepack_index.rst
+1 −1 doc/source/intro/copyright.rst
+16 −30 doc/source/science_guide/sg_thermo.rst
+45 −48 doc/source/user_guide/interfaces.include
+3 −4 doc/source/user_guide/lg_overview.rst
+0 −5 doc/source/user_guide/ug_case_settings.rst
+1 −5 doc/source/user_guide/ug_running.rst
Loading