Skip to content
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

Merge Consortium Icepack main #b0283b287f April 18, 2023 #23

Merged
merged 3 commits into from
Apr 21, 2023
Merged

Merge Consortium Icepack main #b0283b287f April 18, 2023 #23

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
45bb78e
Fix bug in recent fsd mods to support optional wlat argument. (#435)
apcraig Apr 18, 2023
b0283b2
Update icepack_write_parameters and add call in icepack driver. (#436)
apcraig Apr 18, 2023
e5a95e4
Merge branch 'main' of https://github.com/cice-consortium/icepack int…
apcraig Apr 18, 2023
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
338 changes: 170 additions & 168 deletions columnphysics/icepack_parameters.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1809,190 +1809,192 @@ subroutine icepack_write_parameters(iounit)

character(len=*),parameter :: subname='(icepack_write_parameters)'

write(iounit,*) ""
write(iounit,*) subname
write(iounit,*) " rhos = ",rhos
write(iounit,*) " rhoi = ",rhoi
write(iounit,*) " rhow = ",rhow
write(iounit,*) " cp_air = ",cp_air
write(iounit,*) " rhos = ",rhos
write(iounit,*) " rhoi = ",rhoi
write(iounit,*) " rhow = ",rhow
write(iounit,*) " cp_air = ",cp_air
write(iounit,*) " emissivity = ",emissivity
write(iounit,*) " floediam = ",floediam
write(iounit,*) " hfrazilmin = ",hfrazilmin
write(iounit,*) " cp_ice = ",cp_ice
write(iounit,*) " cp_ocn = ",cp_ocn
write(iounit,*) " depressT = ",depressT
write(iounit,*) " dragio = ",dragio
write(iounit,*) " calc_dragio = ",calc_dragio
write(iounit,*) " cp_ice = ",cp_ice
write(iounit,*) " cp_ocn = ",cp_ocn
write(iounit,*) " depressT = ",depressT
write(iounit,*) " dragio = ",dragio
write(iounit,*) " calc_dragio= ",calc_dragio
write(iounit,*) " iceruf_ocn = ",iceruf_ocn
write(iounit,*) " thickness_ocn_layer1 = ",thickness_ocn_layer1
write(iounit,*) " albocn = ",albocn
write(iounit,*) " gravit = ",gravit
write(iounit,*) " albocn = ",albocn
write(iounit,*) " gravit = ",gravit
write(iounit,*) " viscosity_dyn = ",viscosity_dyn
write(iounit,*) " Tocnfrz = ",Tocnfrz
write(iounit,*) " rhofresh = ",rhofresh
write(iounit,*) " zvir = ",zvir
write(iounit,*) " vonkar = ",vonkar
write(iounit,*) " cp_wv = ",cp_wv
write(iounit,*) " Tocnfrz = ",Tocnfrz
write(iounit,*) " rhofresh = ",rhofresh
write(iounit,*) " zvir = ",zvir
write(iounit,*) " vonkar = ",vonkar
write(iounit,*) " cp_wv = ",cp_wv
write(iounit,*) " stefan_boltzmann = ",stefan_boltzmann
write(iounit,*) " Tffresh = ",Tffresh
write(iounit,*) " Lsub = ",Lsub
write(iounit,*) " Lvap = ",Lvap
write(iounit,*) " Timelt = ",Timelt
write(iounit,*) " Tsmelt = ",Tsmelt
write(iounit,*) " Tffresh = ",Tffresh
write(iounit,*) " Lsub = ",Lsub
write(iounit,*) " Lvap = ",Lvap
write(iounit,*) " Timelt = ",Timelt
write(iounit,*) " Tsmelt = ",Tsmelt
write(iounit,*) " ice_ref_salinity = ",ice_ref_salinity
write(iounit,*) " iceruf = ",iceruf
write(iounit,*) " Cf = ",Cf
write(iounit,*) " Pstar = ",Pstar
write(iounit,*) " Cstar = ",Cstar
write(iounit,*) " kappav = ",kappav
write(iounit,*) " kice = ",kice
write(iounit,*) " ksno = ",ksno
write(iounit,*) " zref = ",zref
write(iounit,*) " hs_min = ",hs_min
write(iounit,*) " snowpatch = ",snowpatch
write(iounit,*) " rhosi = ",rhosi
write(iounit,*) " sk_l = ",sk_l
write(iounit,*) " saltmax = ",saltmax
write(iounit,*) " phi_init = ",phi_init
write(iounit,*) " min_salin = ",min_salin
write(iounit,*) " salt_loss = ",salt_loss
write(iounit,*) " iceruf = ",iceruf
write(iounit,*) " Cf = ",Cf
write(iounit,*) " Pstar = ",Pstar
write(iounit,*) " Cstar = ",Cstar
write(iounit,*) " kappav = ",kappav
write(iounit,*) " kice = ",kice
write(iounit,*) " ksno = ",ksno
write(iounit,*) " zref = ",zref
write(iounit,*) " hs_min = ",hs_min
write(iounit,*) " snowpatch = ",snowpatch
write(iounit,*) " rhosi = ",rhosi
write(iounit,*) " sk_l = ",sk_l
write(iounit,*) " saltmax = ",saltmax
write(iounit,*) " phi_init = ",phi_init
write(iounit,*) " min_salin = ",min_salin
write(iounit,*) " salt_loss = ",salt_loss
write(iounit,*) " Tliquidus_max = ",Tliquidus_max
write(iounit,*) " min_bgc = ",min_bgc
write(iounit,*) " min_bgc = ",min_bgc
write(iounit,*) " dSin0_frazil = ",dSin0_frazil
write(iounit,*) " hi_ssl = ",hi_ssl
write(iounit,*) " hs_ssl = ",hs_ssl
write(iounit,*) " awtvdr = ",awtvdr
write(iounit,*) " awtidr = ",awtidr
write(iounit,*) " awtvdf = ",awtvdf
write(iounit,*) " awtidf = ",awtidf
write(iounit,*) " qqqice = ",qqqice
write(iounit,*) " TTTice = ",TTTice
write(iounit,*) " qqqocn = ",qqqocn
write(iounit,*) " TTTocn = ",TTTocn
write(iounit,*) " argcheck = ",argcheck
write(iounit,*) " puny = ",puny
write(iounit,*) " bignum = ",bignum
write(iounit,*) " secday = ",secday
write(iounit,*) " pi = ",pi
write(iounit,*) " pih = ",pih
write(iounit,*) " piq = ",piq
write(iounit,*) " pi2 = ",pi2
write(iounit,*) " hi_ssl = ",hi_ssl
write(iounit,*) " hs_ssl = ",hs_ssl
write(iounit,*) " awtvdr = ",awtvdr
write(iounit,*) " awtidr = ",awtidr
write(iounit,*) " awtvdf = ",awtvdf
write(iounit,*) " awtidf = ",awtidf
write(iounit,*) " qqqice = ",qqqice
write(iounit,*) " TTTice = ",TTTice
write(iounit,*) " qqqocn = ",qqqocn
write(iounit,*) " TTTocn = ",TTTocn
write(iounit,*) " argcheck = ",trim(argcheck)
write(iounit,*) " puny = ",puny
write(iounit,*) " bignum = ",bignum
write(iounit,*) " secday = ",secday
write(iounit,*) " pi = ",pi
write(iounit,*) " pih = ",pih
write(iounit,*) " piq = ",piq
write(iounit,*) " pi2 = ",pi2
write(iounit,*) " rad_to_deg = ",rad_to_deg
write(iounit,*) " Lfresh = ",Lfresh
write(iounit,*) " cprho = ",cprho
write(iounit,*) " Cp = ",Cp
write(iounit,*) " ktherm = ", ktherm
write(iounit,*) " conduct = ", conduct
write(iounit,*) " fbot_xfer_type = ", fbot_xfer_type
write(iounit,*) " calc_Tsfc = ", calc_Tsfc
write(iounit,*) " update_ocn_f = ", update_ocn_f
write(iounit,*) " dts_b = ", dts_b
write(iounit,*) " ustar_min = ", ustar_min
write(iounit,*) " hi_min = ", hi_min
write(iounit,*) " a_rapid_mode = ", a_rapid_mode
write(iounit,*) " Rac_rapid_mode = ", Rac_rapid_mode
write(iounit,*) " Lfresh = ",Lfresh
write(iounit,*) " cprho = ",cprho
write(iounit,*) " Cp = ",Cp
write(iounit,*) " ktherm = ", ktherm
write(iounit,*) " conduct = ", trim(conduct)
write(iounit,*) " fbot_xfer_type = ", trim(fbot_xfer_type)
write(iounit,*) " calc_Tsfc = ", calc_Tsfc
write(iounit,*) " update_ocn_f = ", update_ocn_f
write(iounit,*) " dts_b = ", dts_b
write(iounit,*) " ustar_min = ", ustar_min
write(iounit,*) " hi_min = ", hi_min
write(iounit,*) " a_rapid_mode = ", a_rapid_mode
write(iounit,*) " Rac_rapid_mode = ", Rac_rapid_mode
write(iounit,*) " aspect_rapid_mode = ", aspect_rapid_mode
write(iounit,*) " dSdt_slow_mode = ", dSdt_slow_mode
write(iounit,*) " phi_c_slow_mode = ", phi_c_slow_mode
write(iounit,*) " phi_i_mushy = ", phi_i_mushy
write(iounit,*) " shortwave = ", shortwave
write(iounit,*) " albedo_type = ", albedo_type
write(iounit,*) " albicev = ", albicev
write(iounit,*) " albicei = ", albicei
write(iounit,*) " albsnowv = ", albsnowv
write(iounit,*) " albsnowi = ", albsnowi
write(iounit,*) " ahmax = ", ahmax
write(iounit,*) " R_ice = ", R_ice
write(iounit,*) " R_pnd = ", R_pnd
write(iounit,*) " R_snw = ", R_snw
write(iounit,*) " dT_mlt = ", dT_mlt
write(iounit,*) " rsnw_mlt = ", rsnw_mlt
write(iounit,*) " kalg = ", kalg
write(iounit,*) " kstrength = ", kstrength
write(iounit,*) " krdg_partic = ", krdg_partic
write(iounit,*) " krdg_redist = ", krdg_redist
write(iounit,*) " mu_rdg = ", mu_rdg
write(iounit,*) " atmbndy = ", atmbndy
write(iounit,*) " calc_strair = ", calc_strair
write(iounit,*) " formdrag = ", formdrag
write(iounit,*) " highfreq = ", highfreq
write(iounit,*) " natmiter = ", natmiter
write(iounit,*) " atmiter_conv = ", atmiter_conv
write(iounit,*) " tfrz_option = ", tfrz_option
write(iounit,*) " saltflux_option = ", saltflux_option
write(iounit,*) " kitd = ", kitd
write(iounit,*) " kcatbound = ", kcatbound
write(iounit,*) " floeshape = ", floeshape
write(iounit,*) " wave_spec = ", wave_spec
write(iounit,*) " wave_spec_type= ", wave_spec_type
write(iounit,*) " nfreq = ", nfreq
write(iounit,*) " hs0 = ", hs0
write(iounit,*) " frzpnd = ", frzpnd
write(iounit,*) " dpscale = ", dpscale
write(iounit,*) " rfracmin = ", rfracmin
write(iounit,*) " rfracmax = ", rfracmax
write(iounit,*) " pndaspect = ", pndaspect
write(iounit,*) " hs1 = ", hs1
write(iounit,*) " hp1 = ", hp1
write(iounit,*) " snwredist = ", snwredist
write(iounit,*) " snw_aging_table = ", snw_aging_table
write(iounit,*) " snwgrain = ", snwgrain
write(iounit,*) " dSdt_slow_mode = ", dSdt_slow_mode
write(iounit,*) " phi_c_slow_mode = ", phi_c_slow_mode
write(iounit,*) " phi_i_mushy= ", phi_i_mushy
write(iounit,*) " shortwave = ", trim(shortwave)
write(iounit,*) " albedo_type= ", trim(albedo_type)
write(iounit,*) " albicev = ", albicev
write(iounit,*) " albicei = ", albicei
write(iounit,*) " albsnowv = ", albsnowv
write(iounit,*) " albsnowi = ", albsnowi
write(iounit,*) " ahmax = ", ahmax
write(iounit,*) " R_ice = ", R_ice
write(iounit,*) " R_pnd = ", R_pnd
write(iounit,*) " R_snw = ", R_snw
write(iounit,*) " dT_mlt = ", dT_mlt
write(iounit,*) " rsnw_mlt = ", rsnw_mlt
write(iounit,*) " kalg = ", kalg
write(iounit,*) " kstrength = ", kstrength
write(iounit,*) " krdg_partic= ", krdg_partic
write(iounit,*) " krdg_redist= ", krdg_redist
write(iounit,*) " mu_rdg = ", mu_rdg
write(iounit,*) " atmbndy = ", trim(atmbndy)
write(iounit,*) " calc_strair= ", calc_strair
write(iounit,*) " formdrag = ", formdrag
write(iounit,*) " highfreq = ", highfreq
write(iounit,*) " natmiter = ", natmiter
write(iounit,*) " atmiter_conv = ", atmiter_conv
write(iounit,*) " tfrz_option= ", trim(tfrz_option)
write(iounit,*) " saltflux_option = ", trim(saltflux_option)
write(iounit,*) " kitd = ", kitd
write(iounit,*) " kcatbound = ", kcatbound
write(iounit,*) " floeshape = ", floeshape
write(iounit,*) " wave_spec = ", wave_spec
write(iounit,*) " wave_spec_type = ", trim(wave_spec_type)
write(iounit,*) " nfreq = ", nfreq
write(iounit,*) " hs0 = ", hs0
write(iounit,*) " frzpnd = ", trim(frzpnd)
write(iounit,*) " dpscale = ", dpscale
write(iounit,*) " rfracmin = ", rfracmin
write(iounit,*) " rfracmax = ", rfracmax
write(iounit,*) " pndaspect = ", pndaspect
write(iounit,*) " hs1 = ", hs1
write(iounit,*) " hp1 = ", hp1
write(iounit,*) " snwredist = ", trim(snwredist)
write(iounit,*) " snw_aging_table = ", trim(snw_aging_table)
write(iounit,*) " snwgrain = ", snwgrain
write(iounit,*) " use_smliq_pnd = ", use_smliq_pnd
write(iounit,*) " rsnw_fall = ", rsnw_fall
write(iounit,*) " rsnw_tmax = ", rsnw_tmax
write(iounit,*) " rhosnew = ", rhosnew
write(iounit,*) " rhosmin = ", rhosmin
write(iounit,*) " rhosmax = ", rhosmax
write(iounit,*) " windmin = ", windmin
write(iounit,*) " drhosdwind = ", drhosdwind
write(iounit,*) " snwlvlfac = ", snwlvlfac
write(iounit,*) " isnw_T = ", isnw_T
write(iounit,*) " isnw_Tgrd = ", isnw_Tgrd
write(iounit,*) " isnw_rhos = ", isnw_rhos
write(iounit,*) " snowage_rhos = ", snowage_rhos(1)
write(iounit,*) " snowage_Tgrd = ", snowage_Tgrd(1)
write(iounit,*) " snowage_T = ", snowage_T(1)
write(iounit,*) " snowage_tau = ", snowage_tau(1,1,1)
write(iounit,*) " snowage_kappa = ", snowage_kappa(1,1,1)
write(iounit,*) " snowage_drdt0 = ", snowage_drdt0(1,1,1)
write(iounit,*) " rsnw_fall = ", rsnw_fall
write(iounit,*) " rsnw_tmax = ", rsnw_tmax
write(iounit,*) " rhosnew = ", rhosnew
write(iounit,*) " rhosmin = ", rhosmin
write(iounit,*) " rhosmax = ", rhosmax
write(iounit,*) " windmin = ", windmin
write(iounit,*) " drhosdwind = ", drhosdwind
write(iounit,*) " snwlvlfac = ", snwlvlfac
write(iounit,*) " isnw_T = ", isnw_T
write(iounit,*) " isnw_Tgrd = ", isnw_Tgrd
write(iounit,*) " isnw_rhos = ", isnw_rhos
! write(iounit,*) " snowage_rhos = ", snowage_rhos(1)
! write(iounit,*) " snowage_Tgrd = ", snowage_Tgrd(1)
! write(iounit,*) " snowage_T = ", snowage_T(1)
! write(iounit,*) " snowage_tau = ", snowage_tau(1,1,1)
! write(iounit,*) " snowage_kappa = ", snowage_kappa(1,1,1)
! write(iounit,*) " snowage_drdt0 = ", snowage_drdt0(1,1,1)
write(iounit,*) " snw_ssp_table = ", trim(snw_ssp_table)
write(iounit,*) " bgc_flux_type = ", bgc_flux_type
write(iounit,*) " z_tracers = ", z_tracers
write(iounit,*) " scale_bgc = ", scale_bgc
write(iounit,*) " solve_zbgc = ", solve_zbgc
write(iounit,*) " dEdd_algae = ", dEdd_algae
write(iounit,*) " modal_aero = ", modal_aero
write(iounit,*) " bgc_flux_type = ", trim(bgc_flux_type)
write(iounit,*) " z_tracers = ", z_tracers
write(iounit,*) " scale_bgc = ", scale_bgc
write(iounit,*) " solve_zbgc = ", solve_zbgc
write(iounit,*) " dEdd_algae = ", dEdd_algae
write(iounit,*) " modal_aero = ", modal_aero
write(iounit,*) " conserv_check = ", conserv_check
write(iounit,*) " skl_bgc = ", skl_bgc
write(iounit,*) " solve_zsal = ", solve_zsal
write(iounit,*) " grid_o = ", grid_o
write(iounit,*) " l_sk = ", l_sk
write(iounit,*) " initbio_frac = ", initbio_frac
write(iounit,*) " grid_oS = ", grid_oS
write(iounit,*) " l_skS = ", l_skS
write(iounit,*) " phi_snow = ", phi_snow
write(iounit,*) " fr_resp = ", fr_resp
write(iounit,*) " algal_vel = ", algal_vel
write(iounit,*) " R_dFe2dust = ", R_dFe2dust
write(iounit,*) " dustFe_sol = ", dustFe_sol
write(iounit,*) " T_max = ", T_max
write(iounit,*) " fsal = ", fsal
write(iounit,*) " op_dep_min = ", op_dep_min
write(iounit,*) " fr_graze_s = ", fr_graze_s
write(iounit,*) " fr_graze_e = ", fr_graze_e
write(iounit,*) " fr_mort2min = ", fr_mort2min
write(iounit,*) " fr_dFe = ", fr_dFe
write(iounit,*) " k_nitrif = ", k_nitrif
write(iounit,*) " t_iron_conv = ", t_iron_conv
write(iounit,*) " max_loss = ", max_loss
write(iounit,*) " max_dfe_doc1 = ", max_dfe_doc1
write(iounit,*) " fr_resp_s = ", fr_resp_s
write(iounit,*) " y_sk_DMS = ", y_sk_DMS
write(iounit,*) " t_sk_conv = ", t_sk_conv
write(iounit,*) " t_sk_ox = ", t_sk_ox
write(iounit,*) " frazil_scav = ", frazil_scav
write(iounit,*) " sw_redist = ", sw_redist
write(iounit,*) " sw_frac = ", sw_frac
write(iounit,*) " sw_dtemp = ", sw_dtemp
write(iounit,*) " skl_bgc = ", skl_bgc
write(iounit,*) " solve_zsal = ", solve_zsal
write(iounit,*) " grid_o = ", grid_o
write(iounit,*) " l_sk = ", l_sk
write(iounit,*) " initbio_frac = ", initbio_frac
write(iounit,*) " grid_oS = ", grid_oS
write(iounit,*) " l_skS = ", l_skS
write(iounit,*) " phi_snow = ", phi_snow
write(iounit,*) " fr_resp = ", fr_resp
write(iounit,*) " algal_vel = ", algal_vel
write(iounit,*) " R_dFe2dust = ", R_dFe2dust
write(iounit,*) " dustFe_sol = ", dustFe_sol
write(iounit,*) " T_max = ", T_max
write(iounit,*) " fsal = ", fsal
write(iounit,*) " op_dep_min = ", op_dep_min
write(iounit,*) " fr_graze_s = ", fr_graze_s
write(iounit,*) " fr_graze_e = ", fr_graze_e
write(iounit,*) " fr_mort2min= ", fr_mort2min
write(iounit,*) " fr_dFe = ", fr_dFe
write(iounit,*) " k_nitrif = ", k_nitrif
write(iounit,*) " t_iron_conv= ", t_iron_conv
write(iounit,*) " max_loss = ", max_loss
write(iounit,*) " max_dfe_doc1 = ", max_dfe_doc1
write(iounit,*) " fr_resp_s = ", fr_resp_s
write(iounit,*) " y_sk_DMS = ", y_sk_DMS
write(iounit,*) " t_sk_conv = ", t_sk_conv
write(iounit,*) " t_sk_ox = ", t_sk_ox
write(iounit,*) " frazil_scav= ", frazil_scav
write(iounit,*) " sw_redist = ", sw_redist
write(iounit,*) " sw_frac = ", sw_frac
write(iounit,*) " sw_dtemp = ", sw_dtemp
write(iounit,*) ""

end subroutine icepack_write_parameters

Expand Down
Loading