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

Jump to bottom

Update icepack_write_parameters #436

Merged
merged 1 commit into from
Apr 18, 2023
Merged

Update icepack_write_parameters #436

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
336 changes: 169 additions & 167 deletions columnphysics/icepack_parameters.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1785,187 +1785,189 @@ 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
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Is this removing tabs?

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

No tabs, just shifting spaces.

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Oh, this is just shifting whitespace?

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,*) " min_bgc = ",min_bgc
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,*) " 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,*) " 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,*) " 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,*) " 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,*) " 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,*) " 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
4 changes: 3 additions & 1 deletion configuration/driver/icedrv_InitMod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -12,7 +12,7 @@ module icedrv_InitMod
use icepack_intfc, only: icepack_warnings_flush, icepack_warnings_aborted
use icepack_intfc, only: icepack_query_parameters, icepack_query_tracer_flags
use icepack_intfc, only: icepack_write_tracer_flags, icepack_write_tracer_indices
use icepack_intfc, only: icepack_write_tracer_sizes
use icepack_intfc, only: icepack_write_tracer_sizes, icepack_write_parameters
use icedrv_system, only: icedrv_system_abort

implicit none
Expand Down Expand Up @@ -120,6 +120,8 @@ subroutine icedrv_initialize
if (restart) &
call init_shortwave ! initialize radiative transfer

call icepack_write_parameters(nu_diag)

istep = istep + 1 ! update time step counters
istep1 = istep1 + 1
time = time + dt ! determine the time and date
Expand Down