diff --git a/genie-biogem/src/fortran/biogem_data_netCDF.f90 b/genie-biogem/src/fortran/biogem_data_netCDF.f90 index acacb8260..ac90ff5d9 100644 --- a/genie-biogem/src/fortran/biogem_data_netCDF.f90 +++ b/genie-biogem/src/fortran/biogem_data_netCDF.f90 @@ -647,7 +647,7 @@ SUBROUTINE sub_save_netcdf_3d_USER() & trim(loc_unitsname),const_real_zero,const_real_zero) call sub_putvar3d_g('misc_rCdtoCa',loc_iou,n_i,n_j,n_k,loc_ntrec,loc_ijk(:,:,:),loc_mask) end IF - IF (ocn_select(io_Cd) .AND. ocn_select(io_PO4) .AND. ctrl_data_save_slice_diag) THEN + IF (ocn_select(io_Cd) .AND. ocn_select(io_PO4)) THEN loc_unitsname = 'nmol kg-1 (umol kg-1)-1' loc_ijk(:,:,:) = const_real_null DO i=1,n_i @@ -663,7 +663,7 @@ SUBROUTINE sub_save_netcdf_3d_USER() & trim(loc_unitsname),const_real_zero,const_real_zero) call sub_putvar3d_g('misc_rCdtoPO4',loc_iou,n_i,n_j,n_k,loc_ntrec,loc_ijk(:,:,:),loc_mask) end IF - IF (ocn_select(io_Cd) .AND. ocn_select(io_PO4) .AND. ocn_select(io_Ca) .AND. ctrl_data_save_slice_diag) THEN + IF (ocn_select(io_Cd) .AND. ocn_select(io_PO4) .AND. ocn_select(io_Ca)) THEN loc_unitsname = 'umol kg-1 (mmol kg-1)-1' loc_ijk(:,:,:) = const_real_null DO i=1,n_i @@ -1364,6 +1364,26 @@ SUBROUTINE sub_save_netcdf_2d_USER() & n_i,n_j,loc_ntrec,loc_isij(is,:,:),loc_mask_surf) end do end if + !----------------------------------------------------------------------- + ! nutrient availablity diagnostics + !----------------------------------------------------------------------- + if (ctrl_data_save_slice_bio .AND. ctrl_data_save_slice_diag_bio) then + if ( ocn_select(io_PO4) .AND. ocn_select(io_SiO2) ) then + loc_unitsname = 'n/a' + loc_ij(:,:) = const_real_null + DO i=1,n_i + DO j=1,n_j + IF (n_k >= goldstein_k1(i,j)) THEN + loc_ij(i,j) = int_ocn_timeslice(io_SiO2,i,j,n_k)/int_t_timeslice - & + & par_bio_red_POP_PON*int_ocn_timeslice(io_PO4,i,j,n_k)/int_t_timeslice + end IF + END DO + END DO + call sub_adddef_netcdf(loc_iou,3,'misc_sur_SiSTAR','Si Star', & + & trim(loc_unitsname),const_real_zero,const_real_zero) + call sub_putvar2d('misc_sur_SiSTAR',loc_iou,n_i,n_j,loc_ntrec,loc_ij(:,:),loc_mask_surf) + end if + end if ! ### INSERT CODE TO SAVE ADDITIONAL 2-D DATA FIELDS ######################################################################### ! ! ! ############################################################################################################################ ! @@ -1578,7 +1598,7 @@ SUBROUTINE sub_save_netcdf_2d() ! ! save diagnostics data !---------------------------------------------------------------- - If (ctrl_data_save_slice_diag .AND. flag_rokgem) then + If (ctrl_data_save_slice_diag_geochem .AND. flag_rokgem) then loc_unitsname = 'mol kg-1 yr-1' DO ib=1,n_diag_bio loc_ij(:,:) = int_diag_bio_timeslice(ib,:,:)/int_t_timeslice diff --git a/genie-ecogem/data/input/20170811.p0055c.8P8Z.pal.1x_SST.dat b/genie-ecogem/data/input/20170811.p0055c.8P8Z.pal.1x_SST.dat new file mode 100644 index 000000000..d310e0173 --- /dev/null +++ b/genie-ecogem/data/input/20170811.p0055c.8P8Z.pal.1x_SST.dat @@ -0,0 +1,36 @@ +1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +282.71 282.19 281.65 281.54 281.36 281.46 281.52 281.68 282.09 282.62 1e-05 280.32 280.82 281.44 281.67 281.55 281.61 281.79 282.18 282.45 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 278.73 278.42 280.65 281.46 281.91 282.71 +287.18 286.62 286.08 285.84 285.23 284.76 284.44 284.47 284.88 284.67 284.66 282.68 283.09 286.29 286.5 285.41 284.69 284.24 284.34 284.57 284.41 284.24 284.02 283.86 283.86 284.02 284.45 282.86 281.92 282.59 283 1e-05 284.88 285.36 285.18 286.23 +290.13 289.45 289.03 288.04 287.21 286.47 285.94 285.7 285.77 286.04 286.71 1e-05 1e-05 288.75 289.4 288.29 287.29 286.52 286.46 286.81 286.65 286.35 286.1 285.97 285.93 286.11 286.92 1e-05 1e-05 1e-05 1e-05 1e-05 288.81 289.25 289.16 290.3 +292.14 291.03 290.44 289.79 289.12 288.45 287.94 287.67 287.61 287.79 288.41 1e-05 288.45 291.18 291.96 290.6 289.29 288.72 288.67 289.28 289.63 289.45 289.15 288.78 288.48 288.48 289.2 1e-05 1e-05 1e-05 1e-05 1e-05 291.87 292.3 292.63 293.12 +294.1 292.95 292.3 291.7 291.08 290.44 289.95 289.65 289.5 289.62 290.27 1e-05 290.31 293.82 294.52 293.06 291.54 290.94 290.63 291.32 292.87 292.58 291.95 291.37 290.91 290.8 291.36 1e-05 1e-05 1e-05 1e-05 1e-05 294.36 295.01 295.41 295.13 +295.83 294.9 294.22 293.63 293.02 292.43 291.95 291.63 291.43 291.56 292.3 1e-05 1e-05 294.42 295.92 295.63 293.92 293.28 1e-05 293.02 295.63 295.66 294.52 293.79 293.27 293.2 293.54 294.01 1e-05 1e-05 1e-05 1e-05 296.3 296.94 296.94 296.6 +297.58 296.93 296.38 295.9 295.36 294.78 294.29 293.95 293.77 294.01 294.49 1e-05 1e-05 294.49 296.05 296.8 296.57 295.81 1e-05 1e-05 298.08 298.41 297.55 296.68 296.29 296.23 296.36 296.29 296.34 1e-05 1e-05 296.74 298.01 298.69 298.42 298.11 +299.91 299.57 299.24 299.01 298.58 298 297.42 297.08 297.03 297.42 295.65 1e-05 1e-05 1e-05 297.08 297.3 298.16 296.48 1e-05 1e-05 298.4 299.01 299.08 298.46 298.52 298.67 298.52 298.4 298.28 298.11 1e-05 298.91 300.36 300.95 300.56 300.2 +301.53 300.84 300.33 299.91 299.55 299.18 298.82 298.46 297.81 297.09 296.05 1e-05 1e-05 1e-05 298.04 297.82 298.34 296.71 1e-05 1e-05 298.42 298.87 298.82 298.51 298.4 298.47 298.6 298.86 299.28 299.77 1e-05 301.21 302.85 305.06 304.91 302.79 +301.64 301.16 300.6 300.03 299.5 299.02 298.62 298.17 297.6 297.07 296.57 1e-05 1e-05 1e-05 298.51 298.31 298.77 297.26 1e-05 1e-05 298.45 299.3 299.32 298.98 298.93 298.96 299.13 299.38 299.74 300.22 300.96 301.95 302.5 302.5 302.48 302.09 +302.21 301.61 300.88 300.2 299.59 299.08 298.62 298.16 297.7 297.34 297.15 1e-05 1e-05 1e-05 1e-05 298.88 299.15 298 1e-05 1e-05 1e-05 299.79 299.86 299.46 299.54 299.61 299.79 300.01 300.33 300.77 301.4 302.17 302.63 302.59 302.64 302.57 +302.6 301.81 301.02 300.32 299.71 299.19 298.73 298.35 298.01 297.74 297.74 1e-05 1e-05 1e-05 1e-05 299.4 299.57 298.84 1e-05 1e-05 1e-05 300.24 300.44 299.99 300.22 300.29 300.45 300.68 301.02 301.55 302.25 303.04 303.4 303.43 303.53 303.25 +302.39 301.6 300.88 300.25 299.7 299.24 298.87 298.61 298.42 298.27 1e-05 1e-05 1e-05 1e-05 1e-05 299.75 299.92 299.76 1e-05 1e-05 300.98 300.55 301.25 301.16 300.74 300.83 301.07 301.38 301.82 302.56 303.53 304.4 304.6 304.5 304.03 303.22 +301.98 301.25 300.63 300.1 299.65 299.28 299.01 298.9 298.96 299.12 1e-05 1e-05 1e-05 1e-05 1e-05 300.2 300.26 300.76 1e-05 1e-05 301.42 300.87 301.88 302.66 300.97 300.92 301.67 302.2 302.75 303.78 305.22 306.23 305.58 304.79 303.79 302.81 +301.77 301.11 300.54 300.04 299.62 299.29 299.07 299.06 299.47 300.23 1e-05 1e-05 1e-05 1e-05 1e-05 300.76 300.46 1e-05 1e-05 1e-05 1e-05 301.27 302.04 1e-05 301.75 301.05 301.9 303.05 303.46 304.5 306.12 306.24 305.27 304.37 303.4 302.53 +301.49 300.97 300.53 300.18 299.91 299.72 299.65 299.78 300.3 301.5 1e-05 1e-05 1e-05 1e-05 299.76 300.85 300.91 1e-05 1e-05 1e-05 1e-05 301.57 301.86 1e-05 304.05 302.22 302.57 304.01 303.53 304.17 305.49 305.59 304.61 303.71 302.84 302.1 +301.59 301.24 300.99 300.84 300.77 300.77 300.83 301.03 301.61 302.78 1e-05 1e-05 1e-05 1e-05 299.54 301.58 301.17 1e-05 1e-05 1e-05 1e-05 301.95 301.92 1e-05 1e-05 303.47 303.11 304.25 303.38 304.27 304.73 304.95 304.1 303.35 302.65 302.06 +301.57 301.24 301.04 300.97 301 301.16 301.44 301.76 302.29 303.18 1e-05 1e-05 1e-05 299.38 299.73 302.57 1e-05 1e-05 1e-05 1e-05 1e-05 302.3 301.89 1e-05 302.76 303.8 303.46 303.5 1e-05 305.04 304.15 304.55 303.84 303.2 302.58 302.02 +301.12 300.74 300.5 300.39 300.44 300.69 301.26 301.94 302.1 301.89 301.75 1e-05 299 298.46 300.04 303.22 1e-05 1e-05 1e-05 1e-05 1e-05 302.2 301.9 1e-05 303.12 303.87 303.61 303.55 1e-05 303.93 303.86 303.88 303.29 302.73 302.16 301.61 +300.75 300.32 299.99 299.79 299.72 299.87 300.3 301.03 301.43 300.85 299.83 298.65 298.64 299.28 299.67 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 301.68 301.64 303.72 303.7 303.87 303.31 1e-05 1e-05 303.19 303.53 303.35 302.78 302.26 301.75 301.24 +300.47 300.08 299.76 299.53 299.4 299.41 299.67 300.22 300.98 300.36 299.52 299.08 299.09 299.35 299.58 300.3 1e-05 1e-05 1e-05 1e-05 1e-05 300.9 300.93 303.66 303.81 302.92 302.41 1e-05 304.17 303.01 303.14 302.81 302.33 301.85 301.37 300.91 +300.18 299.79 299.43 299.12 298.89 298.8 298.88 299.42 300.71 299.73 299.2 299.12 298.95 298.86 298.77 298.91 1e-05 1e-05 1e-05 1e-05 1e-05 300 300.19 302.64 302.58 301.7 1e-05 1e-05 303.47 302.78 302.7 302.27 301.81 301.37 300.95 300.56 +299.95 299.58 299.15 298.75 298.42 298.22 298.2 298.58 1e-05 299.38 299.1 298.96 298.62 298.34 298 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 299.23 299.57 301.53 1e-05 1e-05 1e-05 1e-05 302.4 302.35 302.36 301.84 301.38 300.96 300.59 300.27 +299.75 299.41 298.91 298.4 297.96 297.67 297.57 297.76 1e-05 298.36 298.82 298.83 298.45 297.95 297.41 1e-05 1e-05 1e-05 297.84 297.47 1e-05 298.75 298.81 1e-05 1e-05 1e-05 1e-05 1e-05 301.57 301.67 302.17 301.58 301.06 300.62 300.26 299.99 +299.6 299.25 298.69 298.07 297.52 297.12 296.87 296.8 1e-05 297.16 298.11 298.81 298.6 297.8 296.92 296.31 1e-05 1e-05 297.03 297.32 297.81 298.05 297.67 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 300.89 301.92 301.74 301.16 300.66 300.2 299.86 +299.77 299.26 298.49 297.78 297.21 296.72 296.11 295.73 1e-05 296.18 297.39 298.8 298.81 298.06 296.82 295.88 295.54 1e-05 296.84 297.45 297.54 297.23 296.45 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 300.01 301.15 301.67 301.7 301.31 300.68 300.13 +300.21 300.13 298.7 297.85 297.5 296.73 295.27 1e-05 1e-05 295.48 296.22 296.5 296.56 297.34 297.85 295.88 295.21 295.88 297.73 297.69 297.31 296.84 295.52 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 299.25 299.85 300.05 299.69 299.4 299.24 299.4 +296.79 296.95 296.72 296.36 296.64 295.8 294 1e-05 1e-05 1e-05 1e-05 294.25 294.63 294.8 294.91 295.76 1e-05 295.1 296.04 295.05 295.07 1e-05 294.75 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 298.6 298.53 298.08 297.55 297.05 296.79 +294.71 294.56 294.2 293.63 293.13 292.81 292.2 1e-05 1e-05 1e-05 1e-05 1e-05 293.4 293.17 292.38 292.64 292.62 292.81 1e-05 293.47 293.5 292.65 292.3 291.89 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 297.35 297.21 297.08 296.24 295.5 294.98 +293.14 292.71 292.03 291.23 290.6 290.33 290.36 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 289.54 289.99 1e-05 1e-05 1e-05 291.27 290.98 290.44 289.6 289.18 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 295.05 294.43 295.11 294.97 294.16 293.56 +291.98 291.12 290.11 289.1 288.43 288.17 288.42 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 286.09 287.09 287.81 1e-05 1e-05 1e-05 1e-05 288.04 287.35 286.29 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 291.4 291.83 291.98 292.54 292.58 +289.15 289.24 288.1 286.91 286.31 286.29 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 283.84 284.33 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 282.88 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 288 288.44 288.13 288.37 288.74 +284.77 284.78 284.43 284 284.12 284.93 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 278.46 278.74 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 284.71 284.42 284.38 284.56 +282.11 282.47 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 274.88 275.07 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 281.55 281.86 +1e-05 1e-05 1e-05 1e-05 1e-05 273.83 273.79 273.74 273.7 273.69 273.66 273.64 273.6 273.56 273.55 273.55 273.56 273.56 273.59 273.6 273.63 273.67 273.69 273.72 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 diff --git a/genie-ecogem/data/input/20170811.p0055c.8P8Z.pal.3x_SST.dat b/genie-ecogem/data/input/20170811.p0055c.8P8Z.pal.3x_SST.dat new file mode 100644 index 000000000..9a610f44f --- /dev/null +++ b/genie-ecogem/data/input/20170811.p0055c.8P8Z.pal.3x_SST.dat @@ -0,0 +1,36 @@ +1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +287.06 286.55 286.05 285.94 285.77 285.86 285.91 286.04 286.42 286.92 1e-05 284.92 285.4 286.05 286.18 285.96 285.98 286.13 286.51 286.79 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 283.46 283.52 285.31 285.98 286.37 287.16 +291.45 290.91 290.38 290.15 289.54 289.09 288.77 288.76 289.13 288.97 289.04 287.23 287.51 290.67 290.88 289.73 288.97 288.5 288.61 288.81 288.63 288.47 288.31 288.16 288.15 288.32 288.8 287.32 286.5 287.06 287.45 1e-05 289.85 290.28 289.9 290.68 +294.4 293.71 293.29 292.33 291.47 290.77 290.24 289.99 290.03 290.29 290.93 1e-05 1e-05 293.02 293.66 292.49 291.5 290.75 290.7 291.04 290.9 290.71 290.47 290.29 290.22 290.38 291.15 1e-05 1e-05 1e-05 1e-05 1e-05 293.48 293.89 293.63 294.6 +296.42 295.3 294.7 294.04 293.33 292.67 292.16 291.89 291.82 291.98 292.57 1e-05 292.77 295.35 296.09 294.71 293.41 292.86 292.86 293.48 293.82 293.66 293.38 293.01 292.74 292.71 293.34 1e-05 1e-05 1e-05 1e-05 1e-05 296.35 296.72 297.02 297.44 +298.31 297.17 296.51 295.9 295.24 294.6 294.1 293.79 293.63 293.73 294.34 1e-05 294.56 297.9 298.57 297.1 295.6 295.01 294.76 295.47 296.97 296.69 296.08 295.49 295.05 294.92 295.41 1e-05 1e-05 1e-05 1e-05 1e-05 298.69 299.33 299.68 299.36 +300.02 299.08 298.4 297.77 297.13 296.52 296.01 295.68 295.46 295.56 296.26 1e-05 1e-05 298.55 299.95 299.6 297.88 297.25 1e-05 297.14 299.63 299.66 298.54 297.82 297.3 297.19 297.51 297.97 1e-05 1e-05 1e-05 1e-05 300.54 301.19 301.18 300.81 +301.72 301.04 300.46 299.93 299.36 298.75 298.23 297.88 297.67 297.87 298.33 1e-05 1e-05 298.68 300.14 300.78 300.45 299.66 1e-05 1e-05 302.06 302.33 301.45 300.59 300.17 300.09 300.2 300.15 300.25 1e-05 1e-05 300.97 302.2 302.9 302.63 302.28 +303.96 303.57 303.18 302.88 302.41 301.79 301.18 300.83 300.75 301.11 299.45 1e-05 1e-05 1e-05 301.16 301.31 301.98 300.32 1e-05 1e-05 302.44 302.97 302.97 302.32 302.33 302.44 302.29 302.21 302.16 302.07 1e-05 303.05 304.47 305.06 304.65 304.29 +305.52 304.79 304.24 303.76 303.34 302.93 302.53 302.15 301.5 300.81 299.86 1e-05 1e-05 1e-05 302.12 301.83 302.17 300.59 1e-05 1e-05 302.52 302.89 302.78 302.44 302.29 302.32 302.45 302.73 303.18 303.72 1e-05 305.26 306.89 309.02 308.83 306.74 +305.67 305.14 304.53 303.91 303.33 302.82 302.37 301.91 301.35 300.84 300.43 1e-05 1e-05 1e-05 302.61 302.32 302.62 301.16 1e-05 1e-05 302.6 303.37 303.32 302.96 302.88 302.86 303.04 303.32 303.71 304.24 305.02 306.02 306.61 306.61 306.57 306.15 +306.25 305.59 304.83 304.1 303.46 302.91 302.42 301.94 301.49 301.17 301.05 1e-05 1e-05 1e-05 1e-05 302.88 303.02 301.93 1e-05 1e-05 1e-05 303.89 303.89 303.47 303.54 303.56 303.75 304.01 304.37 304.86 305.51 306.3 306.79 306.76 306.77 306.65 +306.65 305.82 304.99 304.25 303.6 303.05 302.56 302.17 301.85 301.63 301.7 1e-05 1e-05 1e-05 1e-05 303.41 303.47 302.79 1e-05 1e-05 1e-05 304.38 304.49 304.03 304.24 304.29 304.45 304.73 305.12 305.71 306.41 307.2 307.59 307.62 307.67 307.35 +306.48 305.63 304.87 304.21 303.62 303.13 302.73 302.47 302.3 302.2 1e-05 1e-05 1e-05 1e-05 1e-05 303.78 303.85 303.73 1e-05 1e-05 305.21 304.74 305.32 305.2 304.81 304.89 305.11 305.47 305.97 306.74 307.7 308.55 308.77 308.66 308.19 307.35 +306.1 305.33 304.66 304.09 303.6 303.2 302.92 302.8 302.86 303.07 1e-05 1e-05 1e-05 1e-05 1e-05 304.21 304.21 304.73 1e-05 1e-05 305.68 305.09 305.98 306.65 305.08 305.05 305.75 306.32 306.93 307.98 309.36 310.34 309.73 308.98 307.99 306.98 +305.92 305.22 304.6 304.07 303.62 303.26 303.02 303.01 303.39 304.14 1e-05 1e-05 1e-05 1e-05 1e-05 304.75 304.42 1e-05 1e-05 1e-05 1e-05 305.56 306.2 1e-05 305.85 305.23 306.01 307.17 307.65 308.73 310.28 310.39 309.48 308.6 307.63 306.72 +305.64 305.08 304.6 304.21 303.91 303.7 303.6 303.72 304.21 305.36 1e-05 1e-05 1e-05 1e-05 303.92 304.84 304.9 1e-05 1e-05 1e-05 1e-05 305.91 306.07 1e-05 307.98 306.37 306.69 308.09 307.75 308.45 309.73 309.8 308.87 307.98 307.09 306.3 +305.73 305.33 305.03 304.84 304.73 304.7 304.74 304.9 305.46 306.58 1e-05 1e-05 1e-05 1e-05 303.69 305.54 305.18 1e-05 1e-05 1e-05 1e-05 306.3 306.15 1e-05 1e-05 307.54 307.22 308.37 307.73 308.63 309.05 309.22 308.4 307.63 306.9 306.26 +305.7 305.33 305.08 304.96 304.95 305.07 305.32 305.61 306.12 306.99 1e-05 1e-05 1e-05 303.66 303.87 306.47 1e-05 1e-05 1e-05 1e-05 1e-05 306.6 306.15 1e-05 306.81 307.83 307.55 307.69 1e-05 309.39 308.52 308.85 308.14 307.48 306.81 306.2 +305.25 304.83 304.54 304.41 304.42 304.64 305.15 305.8 305.94 305.77 305.69 1e-05 303.27 302.7 304.07 307.07 1e-05 1e-05 1e-05 1e-05 1e-05 306.51 306.14 1e-05 307.11 307.89 307.69 307.78 1e-05 308.37 308.25 308.21 307.61 307.01 306.4 305.79 +304.88 304.41 304.05 303.82 303.73 303.85 304.24 304.93 305.31 304.81 303.9 302.79 302.79 303.37 303.62 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 305.99 305.84 307.65 307.62 307.87 307.41 1e-05 1e-05 307.67 307.91 307.69 307.1 306.53 305.97 305.41 +304.59 304.15 303.8 303.55 303.4 303.38 303.62 304.12 304.85 304.31 303.57 303.16 303.16 303.35 303.54 304.26 1e-05 1e-05 1e-05 1e-05 1e-05 305.19 305.12 307.55 307.7 306.95 306.55 1e-05 308.54 307.46 307.53 307.16 306.64 306.12 305.58 305.07 +304.25 303.82 303.43 303.1 302.85 302.74 302.81 303.31 304.55 303.66 303.2 303.11 302.93 302.8 302.72 302.91 1e-05 1e-05 1e-05 1e-05 1e-05 304.26 304.36 306.53 306.51 305.79 1e-05 1e-05 307.82 307.19 307.07 306.6 306.1 305.62 305.13 304.68 +303.98 303.56 303.1 302.68 302.34 302.12 302.09 302.45 1e-05 303.29 303.06 302.89 302.53 302.23 301.92 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 303.41 303.68 305.43 1e-05 1e-05 1e-05 1e-05 306.76 306.72 306.69 306.14 305.63 305.16 304.73 304.35 +303.72 303.34 302.82 302.28 301.83 301.52 301.4 301.61 1e-05 302.28 302.74 302.71 302.3 301.79 301.29 1e-05 1e-05 1e-05 302.09 301.72 1e-05 302.85 302.86 1e-05 1e-05 1e-05 1e-05 1e-05 305.94 306.03 306.46 305.83 305.27 304.77 304.35 304.02 +303.52 303.14 302.56 301.91 301.34 300.92 300.66 300.62 1e-05 301.12 302.02 302.64 302.39 301.58 300.76 300.23 1e-05 1e-05 301.2 301.45 301.88 302.07 301.67 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 305.26 306.19 305.93 305.29 304.73 304.22 303.82 +303.63 303.09 302.31 301.58 300.98 300.46 299.85 299.5 1e-05 300.16 301.29 302.58 302.54 301.78 300.59 299.74 299.49 1e-05 300.89 301.42 301.5 301.18 300.42 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 304.38 305.39 305.83 305.77 305.29 304.61 304.01 +304 303.87 302.45 301.59 301.2 300.41 299 1e-05 1e-05 299.46 300.14 300.33 300.35 301.06 301.52 299.65 299.11 299.76 301.59 301.54 301.13 300.69 299.43 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 303.56 304.09 304.26 303.82 303.46 303.23 303.25 +300.73 300.82 300.54 300.14 300.36 299.49 297.76 1e-05 1e-05 1e-05 1e-05 298.14 298.47 298.6 298.68 299.49 1e-05 298.99 299.92 298.94 298.91 1e-05 298.59 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 302.84 302.73 302.25 301.68 301.13 300.81 +298.74 298.53 298.11 297.51 296.98 296.63 296.02 1e-05 1e-05 1e-05 1e-05 1e-05 297.26 297.01 296.23 296.47 296.48 296.74 1e-05 297.39 297.36 296.57 296.17 295.79 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 301.56 301.41 301.25 300.39 299.62 299.05 +297.21 296.73 296.02 295.18 294.55 294.25 294.27 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 293.49 293.89 1e-05 1e-05 1e-05 295.3 294.91 294.29 293.54 293.18 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 299.29 298.73 299.33 299.14 298.31 297.67 +296.06 295.17 294.15 293.11 292.46 292.16 292.38 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 290.15 291.09 291.72 1e-05 1e-05 1e-05 1e-05 292.1 291.35 290.34 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 295.74 296.14 296.24 296.71 296.7 +293.33 293.34 292.18 291 290.37 290.33 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 287.89 288.4 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 287.01 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 292.37 292.8 292.5 292.67 293.04 +289.07 288.98 288.57 288.15 288.24 289.03 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 282.68 282.96 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 289.08 288.83 288.79 288.92 +286.41 286.73 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 278.11 278.74 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 285.93 286.22 +1e-05 1e-05 1e-05 1e-05 1e-05 275.11 275.09 275 274.95 274.95 274.93 274.89 274.8 274.71 274.69 274.69 274.72 274.71 274.75 274.75 274.77 274.8 274.73 274.85 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 diff --git a/genie-ecogem/data/input/20170811.worjh2.PO4.8P8Z.pal.1x_SST.dat b/genie-ecogem/data/input/20170811.worjh2.PO4.8P8Z.pal.1x_SST.dat new file mode 100644 index 000000000..6ed5f7069 --- /dev/null +++ b/genie-ecogem/data/input/20170811.worjh2.PO4.8P8Z.pal.1x_SST.dat @@ -0,0 +1,36 @@ +1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 271.51 271.54 271.67 271.99 272.5 272.59 272.56 272.03 271.7 271.57 271.52 271.51 271.54 271.62 271.7 271.89 272.27 272.41 271.84 271.47 271.35 271.32 271.31 271.32 271.31 271.31 271.31 271.32 271.31 271.31 271.33 271.34 271.34 271.36 1e-05 +274.15 274.35 274.98 275.72 276.37 276.88 276.74 276.3 275.52 274.8 274.55 274.65 274.64 274.75 274.91 275.04 275.43 276.09 276.41 275.55 274.9 274.09 273.45 273.05 272.82 272.57 272.58 272.67 272.87 272.82 272.55 272.84 273.09 273.04 273.37 273.95 +279.74 280.33 281.1 281.87 282.7 283.34 283.32 282.63 281.5 281.16 281.32 281.73 281.86 281.88 281.8 281.71 281.86 282.13 282.96 282.47 281.12 279.33 278.69 278.5 278.58 278.44 278.58 278.67 278.77 278.8 278.42 278.54 278.77 278.83 278.99 279.26 +282.43 283.03 283.78 284.46 285.38 286.21 286.55 286.04 285.06 284.71 284.85 285.26 285.25 285.07 284.9 284.68 284.46 284.39 285.4 1e-05 282.28 281.46 281.24 281.26 281.47 281.54 281.7 281.95 282.26 282.46 282.22 281.88 282.13 281.89 281.8 282 +284.85 285.47 286.28 287.04 287.95 289.29 290.02 289.95 289.01 288.77 288.56 288.55 288.25 287.86 287.48 287.05 286.69 286.55 287.69 1e-05 284.53 284.15 284.19 284.28 284.42 284.3 284.39 284.7 285.6 286.4 286.01 285.04 285.06 284.64 284.4 284.47 +287.22 287.96 289 289.88 289.91 291.83 292.07 291.47 291.75 292 291.47 291.13 290.66 290.15 289.63 289.1 288.73 288.65 290.08 1e-05 287.86 287.82 287.64 287.3 287.23 286.91 286.74 287.06 288.57 290.5 289.79 288.46 287.99 287.39 286.88 286.84 +289.59 290.76 291.83 292.76 1e-05 294.34 294.45 293.73 294.05 294.48 293.97 293.76 293.41 292.74 292.07 291.5 291.04 290.58 293.2 1e-05 291.01 291.94 291.38 290.35 289.71 289.3 288.95 289.97 290.86 293.44 293.29 291.82 290.92 290.11 289.44 289.33 +291.44 294.86 295.62 1e-05 1e-05 295.99 296.32 295.76 295.79 295.72 295.39 296.4 297.43 297.38 295.95 294.63 293.56 291.97 296.94 1e-05 1e-05 295.27 295.28 294.06 292.84 291.86 290.68 295.69 1e-05 294.69 294.3 294.36 294.3 292.91 292.18 291.88 +292.73 296.35 1e-05 1e-05 1e-05 296.91 296.73 296.43 296.35 296.2 295.96 296.12 296.24 295.93 295.43 295.17 294.19 292.85 295.43 1e-05 1e-05 295.27 295.55 296.48 295.98 293.04 291.43 293.64 1e-05 294.87 294.58 294.33 294.08 293.79 293.61 293.38 +293.86 297.04 1e-05 1e-05 1e-05 297.52 297.22 296.96 296.88 296.63 296.53 296.54 296.29 295.89 295.43 294.91 294.31 293.68 295.37 1e-05 1e-05 295.41 295.87 295.52 294.7 293.17 292.2 294.17 1e-05 295.54 295.57 295.21 294.75 294.4 294.21 294 +295.14 297.96 1e-05 1e-05 1e-05 297.89 297.89 297.68 297.58 297.43 297.31 297.06 296.7 296.21 295.69 295.19 294.73 294.55 296.1 1e-05 1e-05 1e-05 296.82 295.73 294.62 293.65 293 1e-05 1e-05 296.42 296.83 296 295.5 295.19 294.99 294.94 +296.4 298.43 1e-05 1e-05 298.74 298.47 298.59 298.63 298.58 298.41 298.12 297.69 297.22 296.63 296.04 295.6 295.33 295.31 296.92 1e-05 1e-05 1e-05 297.83 296.07 295.1 294.47 294.1 1e-05 1e-05 1e-05 297.95 296.77 296.27 296.07 295.96 296.01 +297.61 299.4 1e-05 1e-05 299.21 299.28 299.45 299.46 299.41 299.24 298.82 298.25 297.57 296.97 296.47 296.1 295.94 296 1e-05 1e-05 1e-05 1e-05 298.2 296.41 295.61 295.32 295.56 1e-05 1e-05 1e-05 298.22 297.59 297.14 297.01 297 297.13 +298.87 300.23 1e-05 1e-05 300.43 300.57 300.65 300.54 300.37 299.94 299.3 298.58 297.88 297.33 296.88 296.58 296.56 297.01 1e-05 1e-05 1e-05 1e-05 298.52 296.86 296.24 296.28 297.26 1e-05 1e-05 1e-05 298.09 298.11 298.26 298.07 298.17 298.37 +300.43 1e-05 301.86 302.52 302.07 302.18 302.28 301.74 301.17 300.22 299.37 298.61 298.02 297.49 297.07 296.89 297.19 298.61 1e-05 1e-05 1e-05 1e-05 298.82 297.21 296.68 297.09 298.6 1e-05 1e-05 1e-05 298.76 299.11 299.54 299.42 299.67 299.76 +302.1 1e-05 302.6 303.92 303.52 303.86 303.45 302.27 301.14 300.01 299.19 298.51 298.09 297.75 297.55 297.54 297.91 299.83 1e-05 1e-05 1e-05 1e-05 299.04 297.61 297.36 297.79 299.64 1e-05 1e-05 1e-05 299.84 300.41 301.1 300.91 300.89 301.41 +1e-05 1e-05 303.56 304.37 304.56 304.5 303.48 301.99 300.84 300 299.38 298.82 298.51 298.29 298.4 298.53 298.82 300.56 1e-05 1e-05 1e-05 299.91 298.83 298.46 298.36 298.56 301.04 1e-05 1e-05 1e-05 300.63 301.8 302.13 301.29 301.02 301.77 +1e-05 1e-05 303.82 303.64 304.59 303.62 302.52 301.35 300.43 299.76 299.27 298.73 298.54 298.37 298.53 298.78 299.05 300.54 1e-05 1e-05 1e-05 299.74 299.11 299.13 298.95 299.19 301.42 1e-05 1e-05 1e-05 300.6 301.96 301.7 300.83 300.52 302.13 +1e-05 304.09 304.43 303.24 302.84 301.74 301 300.23 299.53 298.91 298.49 298.21 298.1 298.18 298.5 299.01 299.53 300.87 1e-05 1e-05 299.95 299.26 299.03 299.64 300.12 300.12 301.65 1e-05 1e-05 1e-05 299.95 300.94 300.66 300.22 300.14 302.22 +1e-05 302.65 303.11 302.02 301.21 300.42 299.9 299.31 298.72 298.26 297.95 297.92 297.96 298.31 299 299.66 299.97 300.23 1e-05 1e-05 298.81 298.3 298.49 299.56 301.17 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 299.8 299.41 299.99 300.09 302.37 +1e-05 301.75 301.38 300.58 299.99 299.47 299.03 298.55 297.97 297.57 297.25 297.13 297.22 297.47 298.57 299.55 299.34 298.33 1e-05 298.26 297.77 297.56 297.53 298.2 301.41 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 298.61 298.38 299.55 299.81 302.77 +1e-05 300.62 300.15 299.42 299.02 298.55 298.13 297.77 297.24 296.76 296.38 296.15 296.13 296.37 297.53 300.48 301.63 1e-05 296.72 297.04 297.13 296.86 296.57 296.69 300.09 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 297.56 297.63 299.01 299.16 301.2 +1e-05 298.37 298.78 298.54 298.2 297.8 297.37 296.96 296.53 296.18 295.71 295.29 295.13 295.28 296.72 302.45 1e-05 297.17 296.76 296.68 296.53 296.2 295.74 295.46 298.35 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 296.91 297.17 1e-05 298.49 300.32 +1e-05 297.13 297.49 297.88 297.59 297.15 296.69 296.2 295.75 295.44 295.01 294.47 294.18 294.22 295.77 1e-05 296.77 296.23 296.11 296.16 296.02 295.55 295 294.47 296.64 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 296.76 297.53 1e-05 1e-05 1e-05 +1e-05 1e-05 296.79 297.43 297.36 296.9 296.44 295.66 295.06 294.69 294.18 293.68 293.31 293.1 294.73 1e-05 296.16 295.45 295.69 295.82 295.64 295.1 294.36 293.6 295.46 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 298.3 299.05 1e-05 1e-05 1e-05 +1e-05 1e-05 295.59 297.88 298.02 297.71 297.79 295.95 295.06 294.88 294.1 293.15 292.25 291.83 294.05 1e-05 295.75 295.41 296.3 296.27 296.07 295.44 294.2 292.79 293.82 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 1e-05 294.77 295.58 295.17 294.36 293.64 293.07 292.77 292.76 293.4 293.5 291.14 290.52 292.26 1e-05 1e-05 1e-05 295.1 294.68 294.57 295.23 294.37 292.07 291.66 1e-05 1e-05 293.7 292.88 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 1e-05 1e-05 294.09 293.33 292.22 291.31 290.57 290.04 289.83 289.86 290.2 289.27 289.73 1e-05 1e-05 1e-05 1e-05 292.4 291.86 291.72 291.63 291.62 290.86 290.11 293.27 292.1 291.96 291.56 291.77 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 1e-05 1e-05 291.71 290.87 289.94 289.05 288.43 287.93 287.59 287.42 287.48 287.33 288.85 1e-05 1e-05 1e-05 1e-05 290.06 289.56 289.35 289.22 289.05 288.8 288.53 1e-05 290.16 290.92 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 1e-05 1e-05 288.58 287.92 287.12 286.5 286.07 285.7 285.4 285.07 284.93 285.16 286.52 1e-05 1e-05 1e-05 1e-05 1e-05 287.58 287.48 287.02 286.73 286.41 286.6 287.55 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 1e-05 1e-05 1e-05 284.92 284.15 283.61 283.66 283.48 283.12 282.62 282.41 282.79 284.56 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 285.05 284.27 284.35 284.25 284.35 285.36 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 1e-05 1e-05 1e-05 281.59 281.66 281.56 281.67 281.24 280.55 279.8 279.79 281.14 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 282.29 281.53 281.7 281.84 282.23 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 277.59 277.52 278.12 279.32 277.52 277.16 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 278.93 278.84 279.07 279.12 278.87 276.28 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 273.46 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 274.07 274.32 1e-05 1e-05 273.9 273.9 273.91 273.92 273.86 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +272.07 272.04 272 271.96 271.92 271.87 271.83 271.8 271.78 271.78 271.77 271.77 271.78 271.8 271.81 271.82 271.84 271.86 271.9 271.9 1e-05 1e-05 1e-05 1e-05 272 272.05 272.1 272.12 272.1 272.11 272.13 272.14 272.13 272.13 272.11 272.09 diff --git a/genie-ecogem/data/input/20170811.worjh2.PO4.8P8Z.pal.3x_SST.dat b/genie-ecogem/data/input/20170811.worjh2.PO4.8P8Z.pal.3x_SST.dat new file mode 100644 index 000000000..f3a953f2f --- /dev/null +++ b/genie-ecogem/data/input/20170811.worjh2.PO4.8P8Z.pal.3x_SST.dat @@ -0,0 +1,36 @@ +1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 275.05 275.58 276.15 276.61 276.93 276.89 276.65 276.17 275.59 275.12 274.86 274.76 274.8 274.92 275.12 275.55 276.03 276.16 275.74 275.11 274.27 273.47 273.03 272.61 272.24 271.89 271.82 271.89 272.1 273.75 274.46 274.84 275.2 275.5 1e-05 +279.52 279.77 280.33 280.98 281.59 282.07 281.88 281.37 280.47 279.64 279.36 279.36 279.33 279.42 279.62 279.85 280.31 280.91 281.12 280.4 279.73 278.91 278.44 278.39 278.32 277.97 277.65 277.45 277.41 277.4 277.44 278.02 278.42 278.66 279.01 279.4 +284.67 285.2 285.78 286.42 287.15 287.74 287.72 287.03 285.85 285.51 285.58 285.86 285.98 286.01 286.08 286.01 286.14 286.47 287.2 286.75 285.46 283.78 283.26 283.23 283.34 283.23 283.29 283.31 283.38 283.53 283.47 283.65 283.94 284.01 284.07 284.27 +287.03 287.53 288.17 288.75 289.55 290.35 290.66 290.15 289.19 288.85 288.89 289.18 289.22 289.06 288.91 288.74 288.59 288.51 289.38 1e-05 286.56 285.81 285.69 285.77 285.95 285.93 286.01 286.21 286.51 286.82 286.85 286.66 286.87 286.7 286.56 286.69 +289.23 289.77 290.5 291.17 291.99 293.24 293.88 293.89 292.99 292.68 292.4 292.34 292.07 291.75 291.42 291.02 290.68 290.52 291.53 1e-05 288.73 288.45 288.53 288.61 288.7 288.55 288.55 288.81 289.67 290.57 290.35 289.56 289.56 289.18 288.91 288.93 +291.39 292.07 293.06 293.9 293.9 295.62 295.85 295.37 295.62 295.76 295.26 294.89 294.45 293.99 293.51 292.99 292.62 292.49 293.8 1e-05 291.97 292.05 291.85 291.48 291.34 290.97 290.79 291.03 292.52 294.43 293.94 292.73 292.27 291.67 291.17 291.09 +293.58 294.7 295.75 296.67 1e-05 298.1 298.21 297.58 297.85 298.21 297.64 297.43 297.11 296.49 295.83 295.27 294.78 294.3 296.79 1e-05 294.98 295.79 295.36 294.31 293.65 293.21 292.87 293.8 294.69 297.28 297.18 295.85 294.97 294.17 293.51 293.37 +295.33 298.64 299.4 1e-05 1e-05 299.78 300.1 299.58 299.57 299.45 299.08 300.03 301.02 300.97 299.54 298.22 297.15 295.62 300.46 1e-05 1e-05 299.08 298.99 297.78 296.59 295.62 294.47 299.36 1e-05 298.51 298.18 298.27 298.16 296.79 296.07 295.76 +296.6 300.11 1e-05 1e-05 1e-05 300.74 300.58 300.27 300.15 299.97 299.69 299.82 299.89 299.56 299.04 298.74 297.77 296.5 298.99 1e-05 1e-05 299.13 299.38 300.16 299.59 296.74 295.21 297.34 1e-05 298.73 298.49 298.25 297.99 297.67 297.46 297.23 +297.74 300.84 1e-05 1e-05 1e-05 301.4 301.11 300.84 300.73 300.46 300.31 300.29 300.01 299.57 299.07 298.53 297.93 297.32 298.99 1e-05 1e-05 299.35 299.73 299.27 298.41 296.9 295.98 297.89 1e-05 299.4 299.47 299.15 298.7 298.33 298.11 297.89 +299.04 301.81 1e-05 1e-05 1e-05 301.83 301.82 301.59 301.47 301.29 301.12 300.84 300.44 299.92 299.36 298.84 298.37 298.2 299.75 1e-05 1e-05 1e-05 300.66 299.51 298.36 297.39 296.79 1e-05 1e-05 300.28 300.7 299.95 299.48 299.16 298.94 298.87 +300.32 302.31 1e-05 1e-05 302.68 302.45 302.56 302.55 302.48 302.27 301.96 301.48 300.98 300.36 299.74 299.27 298.98 298.99 300.57 1e-05 1e-05 1e-05 301.66 299.88 298.85 298.2 297.87 1e-05 1e-05 1e-05 301.79 300.72 300.26 300.06 299.94 299.98 +301.55 303.29 1e-05 1e-05 303.16 303.28 303.42 303.42 303.34 303.13 302.67 302.06 301.35 300.71 300.18 299.78 299.61 299.72 1e-05 1e-05 1e-05 1e-05 302.07 300.26 299.38 299.05 299.29 1e-05 1e-05 1e-05 302.09 301.55 301.15 301.02 301 301.12 +302.81 304.12 1e-05 1e-05 304.35 304.52 304.59 304.47 304.28 303.83 303.17 302.42 301.68 301.1 300.62 300.29 300.25 300.71 1e-05 1e-05 1e-05 1e-05 302.43 300.75 300.03 300 300.94 1e-05 1e-05 1e-05 302.01 302.07 302.22 302.05 302.15 302.34 +304.34 1e-05 305.66 306.32 305.93 306.05 306.15 305.63 305.08 304.13 303.26 302.47 301.85 301.29 300.84 300.63 300.9 302.26 1e-05 1e-05 1e-05 1e-05 302.76 301.14 300.53 300.84 302.25 1e-05 1e-05 1e-05 302.68 303.02 303.45 303.34 303.59 303.68 +305.94 1e-05 306.38 307.66 307.31 307.66 307.28 306.14 305.05 303.95 303.09 302.38 301.93 301.54 301.31 301.27 301.61 303.46 1e-05 1e-05 1e-05 1e-05 302.97 301.54 301.2 301.55 303.26 1e-05 1e-05 1e-05 303.71 304.25 304.92 304.75 304.75 305.26 +1e-05 1e-05 307.32 308.12 308.32 308.28 307.31 305.87 304.76 303.91 303.25 302.66 302.31 302.05 302.11 302.21 302.49 304.18 1e-05 1e-05 1e-05 303.85 302.79 302.34 302.14 302.3 304.64 1e-05 1e-05 1e-05 304.46 305.57 305.89 305.11 304.88 305.6 +1e-05 1e-05 307.58 307.42 308.37 307.45 306.38 305.26 304.35 303.65 303.12 302.55 302.32 302.12 302.22 302.45 302.72 304.14 1e-05 1e-05 1e-05 303.65 303.02 302.95 302.7 302.91 305.02 1e-05 1e-05 1e-05 304.42 305.71 305.46 304.66 304.42 305.95 +1e-05 307.9 308.21 307.07 306.7 305.66 304.95 304.18 303.47 302.82 302.35 302.03 301.87 301.92 302.18 302.65 303.17 304.47 1e-05 1e-05 303.78 303.17 302.91 303.4 303.8 303.79 305.22 1e-05 1e-05 1e-05 303.82 304.73 304.46 304.08 304.05 306.03 +1e-05 306.5 306.94 305.9 305.15 304.41 303.9 303.29 302.67 302.16 301.8 301.71 301.7 302 302.62 303.25 303.58 303.87 1e-05 1e-05 302.68 302.23 302.35 303.29 304.86 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 303.62 303.25 303.83 303.98 306.15 +1e-05 305.64 305.3 304.54 303.98 303.48 303.02 302.51 301.9 301.44 301.08 300.9 300.94 301.14 302.16 303.11 302.99 302.12 1e-05 302.1 301.66 301.46 301.37 301.94 305.05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 302.49 302.24 303.38 303.69 306.5 +1e-05 304.53 304.11 303.41 303.02 302.55 302.11 301.69 301.13 300.59 300.17 299.91 299.84 300.04 301.12 303.94 305.18 1e-05 300.65 300.91 300.99 300.72 300.38 300.45 303.7 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 301.47 301.49 302.84 303.04 304.98 +1e-05 302.39 302.77 302.52 302.18 301.76 301.3 300.84 300.36 299.98 299.46 299.02 298.82 298.92 300.29 305.88 1e-05 300.95 300.62 300.51 300.34 300 299.51 299.22 301.96 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 300.8 301.01 1e-05 302.33 304.11 +1e-05 301.14 301.48 301.84 301.53 301.06 300.56 300.03 299.54 299.19 298.72 298.16 297.84 297.83 299.31 1e-05 300.44 300.01 299.93 299.95 299.77 299.29 298.71 298.19 300.29 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 300.58 301.3 1e-05 1e-05 1e-05 +1e-05 1e-05 300.74 301.34 301.24 300.74 300.24 299.44 298.81 298.4 297.87 297.33 296.93 296.7 298.25 1e-05 299.81 299.21 299.46 299.56 299.33 298.76 298 297.29 299.1 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 301.97 302.7 1e-05 1e-05 1e-05 +1e-05 1e-05 299.5 301.68 301.74 301.46 301.47 299.63 298.73 298.5 297.71 296.75 295.85 295.42 297.5 1e-05 299.36 299.1 299.97 299.93 299.68 298.96 297.74 296.43 297.49 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 1e-05 298.66 299.42 299.03 298.19 297.47 296.86 296.53 296.46 297 297.01 294.72 294.12 295.73 1e-05 1e-05 1e-05 298.76 298.33 298.15 298.64 297.73 295.71 295.37 1e-05 1e-05 297.27 296.5 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 1e-05 1e-05 297.86 297.18 296.09 295.2 294.45 293.92 293.64 293.6 293.77 292.91 293.34 1e-05 1e-05 1e-05 1e-05 296.12 295.58 295.38 295.19 295.03 294.47 293.82 296.77 295.65 295.56 295.18 295.4 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 1e-05 1e-05 295.46 294.66 293.81 292.96 292.34 291.82 291.45 291.23 291.24 291.04 292.45 1e-05 1e-05 1e-05 1e-05 293.82 293.34 293.09 292.87 292.52 292.37 292.21 1e-05 293.76 294.51 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 1e-05 1e-05 292.34 291.75 291.03 290.44 290.02 289.63 289.31 288.95 288.76 288.93 290.15 1e-05 1e-05 1e-05 1e-05 1e-05 291.37 291.27 290.71 290.17 289.82 290.21 291.11 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 1e-05 1e-05 1e-05 288.79 288.08 287.59 287.62 287.44 287.07 286.56 286.32 286.64 288.27 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 288.99 288.09 287.81 287.48 287.8 289.03 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 1e-05 1e-05 1e-05 285.63 285.6 285.5 285.64 285.17 284.5 283.83 283.79 285.03 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 286.31 285.51 285.37 284.7 285.08 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 281.65 281.61 282.16 283.37 281.62 281.3 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 282.98 282.89 282.78 279.75 279.91 279.72 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 277.62 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 277.65 278.11 1e-05 1e-05 277.71 277.11 277.15 277.17 277.2 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 1e-05 +272.79 272.76 272.73 272.7 272.66 272.63 272.61 272.58 272.56 272.54 272.54 272.53 272.54 272.55 272.56 272.58 272.59 272.61 272.62 272.53 1e-05 1e-05 1e-05 1e-05 272.69 272.76 272.84 272.91 272.86 272.83 272.83 272.83 272.84 272.86 272.85 272.82 diff --git a/genie-main/configs/muffin.CBE.p0055c.BASES.config b/genie-main/configs/muffin.CBE.p0055c.BASES.config new file mode 100644 index 000000000..dcd2307ea --- /dev/null +++ b/genie-main/configs/muffin.CBE.p0055c.BASES.config @@ -0,0 +1,150 @@ +# ******************************************************************* +# CONFIGURATION TEMPLATE (EMBM + GOLDSTEIN + GOLDSTEIN SEA-ICE + ATCHEM + BIOGEM) +# ******************************************************************* + +# ******************************************************************* +# GENIE COMPONENT SELECTION +# ******************************************************************* +ma_flag_ebatmos=.TRUE. +ma_flag_goldsteinocean=.TRUE. +ma_flag_goldsteinseaice=.TRUE. +ma_flag_biogem=.TRUE. +ma_flag_atchem=.TRUE. +ma_flag_ecogem=.TRUE. +# ******************************************************************* + +# ******************************************************************* +# GRID CONFIGURATION +# ******************************************************************* +# Topography +ma_fname_topo='p0055c' +# Grid resolution of climate components +GENIENXOPTS='$(DEFINE)GENIENX=36' +GENIENYOPTS='$(DEFINE)GENIENY=36' +GENIENLOPTS='$(DEFINE)GENIENL=1' +GOLDSTEINNLONSOPTS='$(DEFINE)GOLDSTEINNLONS=36' +GOLDSTEINNLATSOPTS='$(DEFINE)GOLDSTEINNLATS=36' +GOLDSTEINNLEVSOPTS='$(DEFINE)GOLDSTEINNLEVS=16' +# Topography for climate components +ea_topo='p0055c' +go_topo='p0055c' +gs_topo='p0055c' +# Wind forcing for climate components +ea_taux_u="p0055c_taux_u.interp" +ea_tauy_u="p0055c_tauy_u.interp" +ea_taux_v="p0055c_taux_v.interp" +ea_tauy_v="p0055c_tauy_v.interp" +ea_adv_u="p0055c_uncep.silo" +ea_adv_v="p0055c_vncep.silo" +#assumed longitudinal offset of the grid +gm_par_grid_lon_offset=-180.0 +# ******************************************************************* + +# ******************************************************************* +# TRACER CONFIGURATION +# ******************************************************************* +# Set number of tracers +GOLDSTEINNTRACSOPTS='$(DEFINE)GOLDSTEINNTRACS=14' +# Set selected tracers +gm_atm_select_3=.true. +gm_atm_select_4=.true. +gm_atm_select_6=.true. +gm_atm_select_16=.true. +gm_ocn_select_3=.true. +gm_ocn_select_4=.true. +gm_ocn_select_8=.true. +gm_ocn_select_10=.true. +gm_ocn_select_12=.true. +gm_ocn_select_15=.true. +gm_ocn_select_16=.true. +gm_ocn_select_20=.true. +gm_ocn_select_35=.true. +gm_ocn_select_38=.true. +gm_ocn_select_40=.true. +gm_ocn_select_50=.true. +gm_sed_select_3=.true. +gm_sed_select_4=.true. +gm_sed_select_8=.true. +gm_sed_select_14=.true. +gm_sed_select_15=.true. +gm_sed_select_22=.true. +gm_sed_select_32=.true. +gm_sed_select_33=.true. +gm_sed_select_34=.true. +gm_sed_select_36=.true. +gm_sed_select_87=.true. +# Set tracer initial values +ac_atm_init_3=278.0E-06 +ac_atm_init_4=-6.5 +ac_atm_init_6=0.2095 +bg_ocn_init_3=2.244E-03 +bg_ocn_init_4=0.4 +bg_ocn_init_8=2.159E-06 +bg_ocn_init_10=1.696E-04 +bg_ocn_init_12=1.700E-03 +bg_ocn_init_35=1.822E-02 +bg_ocn_init_38=1.500E-02 +bg_ocn_init_40=0.0 +bg_ocn_init_50=2.989E-02 +# ******************************************************************* + +# ******************************************************************* +# PHYSICAL CLIMATE CONFIGURATION #1: main goin parameters +# ******************************************************************* +# kraus-turner mixed layer scheme on (1) or off (0) +go_imld = 1 +# temp0 +go_10=10.0 +# temp1 +go_11=10.0 +# rel +go_12=0.9000000 +# scl_tau / SclTau +go_13=1.531013488769531300 +# ocean diffusivites iso (or horiz) / OcnHorizDiff +go_14=1494.438354492187500000 +# ocean diffusivites dia (or vert) / OcnVertDiff +go_15=0.000025363247914356 +# inverse minimum drag in days / InvDrag +go_16=2.710164785385131800 +# scl_tau (should be same as GOLDSTEIN's value) / SclTau +ea_11=1.531013488769531300 +# atm. diff. amp. for T / AtmDiffAmpT +ea_12=5204945.000000000000000000 +# atm. diff. amp. for q / AtmDiffAmpQ +ea_13=1173269.250000000000000000 +# dist'n width / AtmDiffWidthT +ea_14=1.410347938537597700 +# dist'n slope / AtmDiffSlopeT +ea_15=0.090003050863742828 +# atm. advection factor for T_z / AtmAdvectCoeffTm +ea_16=0.001037851092405617 +# atm. advection factor for T_m / AtmAdvectCoeffQm +ea_17=0.0000000E+00 +# atm. advection factor for q_z / AtmAdvectCoeffQz +ea_18=0.164652019739151000 +# atm. advection factor for q_m / AtmAdvectCoeffQz +ea_19=0.164652019739151000 +# SclFWF +ea_28=0.0 +# sea-ice eddy diffusivity / SeaiceDiff +gs_11=3573.718017578125000000 +# solar constant reduced by 0.46% for end Paleocene +ma_genie_solar_constant=1361.7 +# adjusted planetary albedo +ea_albedop_offs=0.200 +ea_albedop_amp=0.260 +ea_albedop_skew=0.0 +ea_albedop_skewp=0 +ea_albedop_mod2=-0.000 +ea_albedop_mod4=0.000 +ea_albedop_mod6=0.250 +# minus 1 PSU +go_saln0=33.9 +# set CO2-climate feedback +ea_36=y +# ******************************************************************* + +# ******************************************************************* +# USER-APPENDED OPTIONS FOLLOW ... +# ******************************************************************* diff --git a/genie-main/configs/muffin.CBE.worjh2.BASES.config b/genie-main/configs/muffin.CBE.worjh2.BASES.config new file mode 100644 index 000000000..1ce25f58d --- /dev/null +++ b/genie-main/configs/muffin.CBE.worjh2.BASES.config @@ -0,0 +1,137 @@ +# ******************************************************************* +# CONFIGURATION TEMPLATE (EMBM + GOLDSTEIN + GOLDSTEIN SEA-ICE + ATCHEM + BIOGEM + ECOGEM) +# ******************************************************************* + +# ******************************************************************* +# GENIE COMPONENT SELECTION +# ******************************************************************* +ma_flag_ebatmos=.TRUE. +ma_flag_goldsteinocean=.TRUE. +ma_flag_goldsteinseaice=.TRUE. +ma_flag_atchem=.TRUE. +ma_flag_biogem=.TRUE. +ma_flag_ecogem=.TRUE. +# ******************************************************************* + +# ******************************************************************* +# GRID CONFIGURATION +# ******************************************************************* +# Topography +ma_fname_topo='worjh2' +# Grid resolution of climate components +GENIENXOPTS='-DGENIENX=36' +GENIENYOPTS='-DGENIENY=36' +GENIENLOPTS='-DGENIENL=1' +GENIENXOPTS='$(DEFINE)GENIENX=36' +GENIENYOPTS='$(DEFINE)GENIENY=36' +GOLDSTEINNLONSOPTS='$(DEFINE)GOLDSTEINNLONS=36' +GOLDSTEINNLATSOPTS='$(DEFINE)GOLDSTEINNLATS=36' +GOLDSTEINNLEVSOPTS='$(DEFINE)GOLDSTEINNLEVS=16' +# Topography for climate components +ea_topo=worjh2 +go_topo=worjh2 +gs_topo=worjh2 +# ******************************************************************* + +# ******************************************************************* +# TRACER CONFIGURATION +# ******************************************************************* +# Set number of tracers +GOLDSTEINNTRACSOPTS='$(DEFINE)GOLDSTEINNTRACS=14' +# Set selected tracers +gm_atm_select_3=.true. +gm_atm_select_4=.true. +gm_atm_select_6=.true. +gm_atm_select_16=.true. +gm_ocn_select_3=.true. +gm_ocn_select_4=.true. +gm_ocn_select_8=.true. +gm_ocn_select_10=.true. +gm_ocn_select_12=.true. +gm_ocn_select_15=.true. +gm_ocn_select_16=.true. +gm_ocn_select_20=.true. +gm_ocn_select_35=.true. +gm_ocn_select_38=.true. +gm_ocn_select_40=.true. +gm_ocn_select_50=.true. +gm_sed_select_3=.true. +gm_sed_select_4=.true. +gm_sed_select_8=.true. +gm_sed_select_14=.true. +gm_sed_select_15=.true. +gm_sed_select_22=.true. +gm_sed_select_32=.true. +gm_sed_select_33=.true. +gm_sed_select_34=.true. +gm_sed_select_36=.true. +# Set tracer initial values +ac_atm_init_3=278.0E-06 +ac_atm_init_4=-6.5 +ac_atm_init_6=0.2095 +bg_ocn_init_3=2.244E-03 +bg_ocn_init_4=0.4 +bg_ocn_init_8=2.159E-06 +bg_ocn_init_10=1.696E-04 +bg_ocn_init_12=2.363E-03 +bg_ocn_init_35=1.025E-02 +bg_ocn_init_38=1.500E-02 +bg_ocn_init_40=0.0 +bg_ocn_init_50=5.282E-02 +# ******************************************************************* + +# ******************************************************************* +# PHYSICAL CLIMATE CONFIGURATION #1: main goin parameters +# ******************************************************************* +# temp0 +go_10=0.0 +# temp1 +go_11=0.0 +# rel +go_12=0.9000000 +# scl_tau / SclTau +go_13=1.2 +# ocean diffusivites iso (or horiz) / OcnHorizDiff +go_14=1494.438354492187500000 +# ocean diffusivites dia (or vert) / OcnVertDiff +go_15=0.000025363247914356 +# inverse minimum drag in days / InvDrag +go_16=2.710164785385131800 +# scl_tau (should be same as GOLDSTEIN's value) / SclTau +ea_11=1.2 +# atm. diff. amp. for T / AtmDiffAmpT +ea_12=5204945.000000000000000000 +# atm. diff. amp. for q / AtmDiffAmpQ +ea_13=1173269.250000000000000000 +# dist'n width / AtmDiffWidthT +ea_14=1.410347938537597700 +# dist'n slope / AtmDiffSlopeT +ea_15=0.090003050863742828 +# atm. advection factor for T_z / AtmAdvectCoeffTm +ea_16=0.001037851092405617 +# atm. advection factor for T_m / AtmAdvectCoeffQm +ea_17=0.0000000E+00 +# atm. advection factor for q_z / AtmAdvectCoeffQz +ea_18=0.164652019739151000 +# atm. advection factor for q_m / AtmAdvectCoeffQz +ea_19=0.164652019739151000 +## extra1a +#ea_25=-2.1228021E-02 +## extra1b +#ea_26=0.1202921 +## extra1c +#ea_27=0.1273681 +# SclFWF +ea_28=0.8 +# sea-ice eddy diffusivity / SeaiceDiff +gs_11=3573.718017578125000000 +#diffusivity scaling factor +ea_diffa_scl=0.25 +#grid point distance over which scalar is applied (j direction) +ea_diffa_len=3 +# ******************************************************************* + +# ******************************************************************* +# USER-APPENDED OPTIONS FOLLOW ... +# ******************************************************************* + diff --git a/genie-main/configs/muffin.CBE.worjh2.BASESFeTDTL.config b/genie-main/configs/muffin.CBE.worjh2.BASESFeTDTL.config new file mode 100644 index 000000000..e675159e0 --- /dev/null +++ b/genie-main/configs/muffin.CBE.worjh2.BASESFeTDTL.config @@ -0,0 +1,150 @@ +# ******************************************************************* +# CONFIGURATION TEMPLATE (EMBM + GOLDSTEIN + GOLDSTEIN SEA-ICE + ATCHEM + BIOGEM + ECOGEM) +# ******************************************************************* + +# ******************************************************************* +# GENIE COMPONENT SELECTION +# ******************************************************************* +ma_flag_ebatmos=.TRUE. +ma_flag_goldsteinocean=.TRUE. +ma_flag_goldsteinseaice=.TRUE. +ma_flag_atchem=.TRUE. +ma_flag_biogem=.TRUE. +ma_flag_ecogem=.TRUE. +# ******************************************************************* + +# ******************************************************************* +# GRID CONFIGURATION +# ******************************************************************* +# Topography +ma_fname_topo='worjh2' +# Grid resolution of climate components +GENIENXOPTS='-DGENIENX=36' +GENIENYOPTS='-DGENIENY=36' +GENIENLOPTS='-DGENIENL=1' +GENIENXOPTS='$(DEFINE)GENIENX=36' +GENIENYOPTS='$(DEFINE)GENIENY=36' +GOLDSTEINNLONSOPTS='$(DEFINE)GOLDSTEINNLONS=36' +GOLDSTEINNLATSOPTS='$(DEFINE)GOLDSTEINNLATS=36' +GOLDSTEINNLEVSOPTS='$(DEFINE)GOLDSTEINNLEVS=16' +# Topography for climate components +ea_topo=worjh2 +go_topo=worjh2 +gs_topo=worjh2 +# ******************************************************************* + +# ******************************************************************* +# TRACER CONFIGURATION +# ******************************************************************* +# Set number of tracers +GOLDSTEINNTRACSOPTS='$(DEFINE)GOLDSTEINNTRACS=17' +# Select tracers +gm_atm_select_3=.true. +gm_atm_select_4=.true. +gm_atm_select_6=.true. +gm_atm_select_16=.true. +gm_ocn_select_3=.true. +gm_ocn_select_4=.true. +gm_ocn_select_8=.true. +gm_ocn_select_10=.true. +gm_ocn_select_12=.true. +gm_ocn_select_15=.true. +gm_ocn_select_16=.true. +gm_ocn_select_20=.true. +gm_ocn_select_22=.true. +gm_ocn_select_90=.true. +gm_ocn_select_42=.true. +gm_ocn_select_35=.true. +gm_ocn_select_38=.true. +gm_ocn_select_40=.true. +gm_ocn_select_50=.true. +gm_sed_select_3=.true. +gm_sed_select_4=.true. +gm_sed_select_8=.true. +gm_sed_select_10=.true. +gm_sed_select_13=.true. +gm_sed_select_14=.true. +gm_sed_select_15=.true. +gm_sed_select_21=.true. +gm_sed_select_22=.true. +gm_sed_select_25=.true. +gm_sed_select_32=.true. +gm_sed_select_33=.true. +gm_sed_select_34=.true. +gm_sed_select_36=.true. +# Set tracer initial values +ac_atm_init_3=278.0E-06 +ac_atm_init_4=-6.5 +ac_atm_init_6=0.2095 +bg_ocn_init_3=2.244E-03 +bg_ocn_init_4=0.4 +bg_ocn_init_8=2.159E-06 +bg_ocn_init_10=1.696E-04 +bg_ocn_init_12=2.363E-03 +bg_ocn_init_15=0.0 +bg_ocn_init_16=0.0 +bg_ocn_init_20=0.0 +bg_ocn_init_22=0.0 +bg_ocn_init_90=0.650E-09 +bg_ocn_init_42=1.000E-09 +bg_ocn_init_35=1.025E-02 +bg_ocn_init_38=2.916E-02 +bg_ocn_init_40=0.0 +bg_ocn_init_50=5.282E-02 +# ******************************************************************* + +# ******************************************************************* +# PHYSICAL CLIMATE CONFIGURATION #1: main goin parameters +# ******************************************************************* +# temp0 +go_10=0.0 +# temp1 +go_11=0.0 +# rel +go_12=0.9000000 +# scl_tau / SclTau +go_13=1.2 +# ocean diffusivites iso (or horiz) / OcnHorizDiff +go_14=1494.438354492187500000 +# ocean diffusivites dia (or vert) / OcnVertDiff +go_15=0.000025363247914356 +# inverse minimum drag in days / InvDrag +go_16=2.710164785385131800 +# scl_tau (should be same as GOLDSTEIN's value) / SclTau +ea_11=1.2 +# atm. diff. amp. for T / AtmDiffAmpT +ea_12=5204945.000000000000000000 +# atm. diff. amp. for q / AtmDiffAmpQ +ea_13=1173269.250000000000000000 +# dist'n width / AtmDiffWidthT +ea_14=1.410347938537597700 +# dist'n slope / AtmDiffSlopeT +ea_15=0.090003050863742828 +# atm. advection factor for T_z / AtmAdvectCoeffTm +ea_16=0.001037851092405617 +# atm. advection factor for T_m / AtmAdvectCoeffQm +ea_17=0.0000000E+00 +# atm. advection factor for q_z / AtmAdvectCoeffQz +ea_18=0.164652019739151000 +# atm. advection factor for q_m / AtmAdvectCoeffQz +ea_19=0.164652019739151000 +## extra1a +#ea_25=-2.1228021E-02 +## extra1b +#ea_26=0.1202921 +## extra1c +#ea_27=0.1273681 +# SclFWF +ea_28=0.8 +# sea-ice eddy diffusivity / SeaiceDiff +gs_11=3573.718017578125000000 +#diffusivity scaling factor +ea_diffa_scl=0.25 +#grid point distance over which scalar is applied (j direction) +ea_diffa_len=3 +# ******************************************************************* + +# ******************************************************************* +# USER-APPENDED OPTIONS FOLLOW ... +# ******************************************************************* + diff --git a/genie-main/src/xml-config/xml/definition.xml b/genie-main/src/xml-config/xml/definition.xml index ee8e9d35b..8fc3e3421 100644 --- a/genie-main/src/xml-config/xml/definition.xml +++ b/genie-main/src/xml-config/xml/definition.xml @@ -6975,6 +6975,10 @@ .TRUE. short circuit atmosphere by not taking CO2 directly, instead have less DIC put into ocean. + + .FALSE. + scale (global) weathering with runoff + Global_avg weathering scheme ID string ('Global_avg','GKWM',or 'GEM-CO2') diff --git a/genie-rokgem/src/fortran/rokgem_box.f90 b/genie-rokgem/src/fortran/rokgem_box.f90 index 673a83ead..b5ed38432 100644 --- a/genie-rokgem/src/fortran/rokgem_box.f90 +++ b/genie-rokgem/src/fortran/rokgem_box.f90 @@ -461,6 +461,7 @@ SUBROUTINE sub_glob_avg_weath(dum_dts,dum_sfcatm1,dum_runoff,dum_photo,dum_respv REAL :: loc_R0 REAL :: loc_maxR REAL :: loc_minR + REAL :: loc_totR REAL :: loc_avP REAL :: loc_P(n_i,n_j) REAL :: loc_P0 @@ -597,13 +598,13 @@ SUBROUTINE sub_glob_avg_weath(dum_dts,dum_sfcatm1,dum_runoff,dum_photo,dum_respv ! calculate mean runoff (mm s-1) ! for equal area grid: - loc_R = 0.0 + loc_totR = 0.0 loc_maxR = 0.0 loc_minR = 0.0 DO i=1,n_i DO j=1,n_j m = landmask(i,j) * dum_runoff(i,j) - loc_R = loc_R + m + loc_totR = loc_totR + m IF ((m.GT.loc_maxR).AND.(landmask(i,j).EQ.1)) THEN loc_maxR = m ENDIF @@ -612,7 +613,7 @@ SUBROUTINE sub_glob_avg_weath(dum_dts,dum_sfcatm1,dum_runoff,dum_photo,dum_respv ENDIF END DO END DO - loc_R = loc_R/nlandcells + loc_R = loc_totR/nlandcells ! convert atm pCO2 to ppm DO i=1,n_i @@ -1086,31 +1087,29 @@ SUBROUTINE sub_glob_avg_weath(dum_dts,dum_sfcatm1,dum_runoff,dum_photo,dum_respv ! ######################################################################################################################### ! ! Spread out atmosphere variables' fluxes onto land - DO k=1,n_atm - IF(k.gt.2) THEN - loc_force_flux_weather_a_percell(k) = loc_force_flux_weather_a(k)/nlandcells - loc_force_flux_weather_a_land(k,:,:) = landmask(:,:) * loc_force_flux_weather_a_percell(k) - end IF + DO k=3,n_atm + loc_force_flux_weather_a_percell(k) = loc_force_flux_weather_a(k)/nlandcells + loc_force_flux_weather_a_land(k,:,:) = landmask(:,:) * loc_force_flux_weather_a_percell(k) END DO ! no need to route to the atmosphere - just take it straight from the cells above the land (assuming same grid) ! convert from Mol/yr to Mol/sec/m^2 and put it into passing array (only take variable altered here - pCO2) - dum_sfxatm1(ia_PCO2,:,:) = loc_force_flux_weather_a_land(ia_PCO2,:,:)/(phys_rok(ipr_A,:,:)*conv_yr_s) + dum_sfxatm1(ia_PCO2,:,:) = loc_force_flux_weather_a_land(ia_PCO2,:,:)/(phys_rok(ipr_A,:,:)*conv_yr_s) dum_sfxatm1(ia_PCO2_13C,:,:) = loc_force_flux_weather_a_land(ia_PCO2_13C,:,:)/(phys_rok(ipr_A,:,:)*conv_yr_s) ! Spread out ocean variables' fluxes onto land - DO k=1,n_ocn - IF(k.gt.2) THEN - loc_force_flux_weather_o_percell(k) = loc_force_flux_weather_o(k)/nlandcells + DO k=3,n_ocn + IF (opt_weather_runoff) THEN + loc_force_flux_weather_o_percell(k) = loc_force_flux_weather_o(k)/loc_totR + loc_force_flux_weather_o_land(k,:,:) = landmask(:,:) * dum_runoff(:,:) * loc_force_flux_weather_o_percell(k) + else + loc_force_flux_weather_o_percell(k) = loc_force_flux_weather_o(k)/nlandcells loc_force_flux_weather_o_land(k,:,:) = landmask(:,:) * loc_force_flux_weather_o_percell(k) - end IF + end if END DO ! route it into the coastal ocean cells (to pass to biogem in coupled model) and save the output to file - DO k=1,n_ocn - IF(k.gt.2) THEN -!!$ IF((k.EQ.io_ALK).OR.(k.EQ.io_DIC).OR.(k.EQ.io_Ca).OR.(k.EQ.io_DIC_13C).OR.(k.EQ.io_DIC_14C)) THEN - CALL sub_coastal_output( loc_force_flux_weather_o_land(k,:,:), & - & runoff_drainto(:,:,:),runoff_detail(:,:), & - & loc_force_flux_weather_o_ocean(k,:,:) ) - ENDIF + DO k=3,n_ocn + CALL sub_coastal_output( loc_force_flux_weather_o_land(k,:,:), & + & runoff_drainto(:,:,:),runoff_detail(:,:), & + & loc_force_flux_weather_o_ocean(k,:,:) ) END DO ! convert from Mol/yr to Mol/sec and put it into passing array dum_sfxrok(:,:,:) = loc_force_flux_weather_o_ocean(:,:,:)/conv_yr_s diff --git a/genie-rokgem/src/fortran/rokgem_data.f90 b/genie-rokgem/src/fortran/rokgem_data.f90 index 296c15b27..3772392b3 100644 --- a/genie-rokgem/src/fortran/rokgem_data.f90 +++ b/genie-rokgem/src/fortran/rokgem_data.f90 @@ -89,7 +89,8 @@ SUBROUTINE sub_load_goin_rokgem() print*,' routes to : ',max_drain_cells !--- WEATHERING PARAMETERS -------------------------------------------------------------------------------------------------- ! print*,'--- WEATHERING PARAMETERS ---' - print*,'short circuit atmosphere : ',opt_short_circuit_atm + print*,'short circuit atmosphere : ',opt_short_circuit_atm + print*,'scale (global) weathering with runoff : ',opt_weather_runoff print*,'weathering scheme ID string : ',trim(par_weathopt) print*,'global CaCO3 weathering scheme ID string : ',trim(opt_weather_CaCO3) print*,'global CaSiO3 weathering scheme ID string : ',trim(opt_weather_CaSiO3) diff --git a/genie-rokgem/src/fortran/rokgem_lib.f90 b/genie-rokgem/src/fortran/rokgem_lib.f90 index 1f9a97b4f..961a6fa50 100644 --- a/genie-rokgem/src/fortran/rokgem_lib.f90 +++ b/genie-rokgem/src/fortran/rokgem_lib.f90 @@ -57,8 +57,10 @@ MODULE rokgem_lib INTEGER::max_drain_cells ! maximum number of ocean cells a single land cell routes to NAMELIST /ini_rokgem_nml/max_drain_cells !--- WEATHERING PARAMETERS ----------------------------------------------------------------------------------------------------- ! - LOGICAL:: opt_short_circuit_atm !short circuit atmosphere by not taking CO2 directly, instead have less DIC put into ocean. + LOGICAL:: opt_short_circuit_atm ! short circuit atmosphere by not taking CO2 directly, instead have less DIC put into ocean. NAMELIST /ini_rokgem_nml/opt_short_circuit_atm + LOGICAL:: opt_weather_runoff ! scale (global) weathering with runoff + NAMELIST /ini_rokgem_nml/opt_weather_runoff CHARACTER(len=63)::par_weathopt ! weathering scheme ID string ('Global_avg','GKWM',or 'GEM_CO2') NAMELIST /ini_rokgem_nml/par_weathopt CHARACTER(len=63)::opt_weather_CaCO3 ! global CaCO3 weathering scheme ID diff --git a/genie-userconfigs/MS/wilsonetal.2018/readme.txt b/genie-userconfigs/MS/wilsonetal.2018/readme.txt new file mode 100644 index 000000000..642b2eade --- /dev/null +++ b/genie-userconfigs/MS/wilsonetal.2018/readme.txt @@ -0,0 +1,61 @@ +################################################################ +### readme.txt ################################################# +################################################################ + +For: ‘Linking marine plankton ecosystems and climate: A new modelling approach to the warm Early Eocene climate.’ +Jamie D. Wilson, Fanny M. Monteiro, Daniela N. Schmidt, Ben A. Ward and Andy Ridgwell + +################################################################ +08/11/2018 -- README.txt file creation (J.D.W) +################################################################ + +Provided is the code used to create the model experiments presented in the paper. +Also given are the configuration files necessary to run the model experiments. + +All experiments are run from:$HOME/cgenie.muffin/genie-main +(unless a different installation directory has been used) + +################## Model Experiments ########################### + +The commands to run the model configurations as listed in the Methods are: + +(1) Modern + +./runmuffin.sh muffin.CBE.worjh2.BASES MS/wilsonetal.2018 wilsonetal.worjh2.8P8Z.pal.1x 10000 + +(2) Late Paleocene Early Eocene Paleogeography + +./runmuffin.sh muffin.CBE.p0055c.BASES MS/wilsonetal.2018 wilsonetal.p0055c.8P8Z.pal.1x 10000 + +(3) Early Eocene CO2 and Climate + +./runmuffin.sh muffin.CBE.p0055c.BASES MS/wilsonetal.2018 wilsonetal.p0055c.8P8Z.pal.3x 10000 + +(4) Temperature Effects of Circulation and Ecology + +./runmuffin.sh muffin.CBE.p0055c.BASES MS/wilsonetal.2018 wilsonetal.p0055c.8P8Z.pal.1x.3xSST 5000 wilsonetal.p0055c.8P8Z.pal.1x + +./runmuffin.sh muffin.CBE.p0055c.BASES MS/wilsonetal.2018 wilsonetal.p0055c.8P8Z.pal.3x.1xSST 5000 wilsonetal.p0055c.8P8Z.pal.3x + +################## Supplementary Model Experiments ############## + +The commands to run the model configurations as listed in the Supplementary Material are: + +(S1) Modern with 3 x CO2 + +./runmuffin.sh muffin.CBE.worjh2.BASES MS/wilsonetal.2018 wilsonetal.worjh2.8P8Z.pal.3x 10000 + +(S2) Temperature Effects of Circulation and Ecology for Modern + +./runmuffin.sh muffin.CBE.worjh2.BASES MS/wilsonetal.2018 wilsonetal.p0055c.8P8Z.pal.1x.3xSST 5000 wilsonetal.p0055c.8P8Z.pal.1x + +./runmuffin.sh muffin.CBE.worjh2.BASES MS/wilsonetal.2018 wilsonetal.worjh2.8P8Z.pal.3x.1xSST 5000 wilsonetal.worjh2.8P8Z.pal.3x + +(S3) Modern with Iron Limitation + +./runmuffin.sh muffin.CBE.worjh2.BASESFeTDTL MS/wilsonetal.2018 wilsonetal.worjh2.8P8Z.pal.1x_Fe 10000 + +################################################################ +################################################################ +################################################################ + diff --git a/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.p0055c.8P8Z.pal.1x b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.p0055c.8P8Z.pal.1x new file mode 100644 index 000000000..51788c797 --- /dev/null +++ b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.p0055c.8P8Z.pal.1x @@ -0,0 +1,115 @@ +# +# +# --- CLIMATE -------------------------------------------------- +# +# set climate feedback +ea_36=y +# +# --- BIOLOGICAL NEW PRODUCTION ------------------------------------- +# +# biological scheme ID string +bg_par_bio_prodopt="NONE" +# +# --- ORGANIC MATTER EXPORT RATIOS ----------------------------- +# +#production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS --------------------------- +# +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +#exponent for modifier of CaCO3:POC export ratio +bg_par_bio_red_POC_CaCO3_pP=0.0 +# +# --- REMINERALIZATION ----------------------------------------- +# +#DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +#initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +#remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +#initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +#remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# + +# --- ECOGEM ---------------------------------------------------------- +# +gm_ctrl_debug_init =0 +eg_ctrl_debug_eco_init =.false. + +eg_par_ecogem_plankton_file ='8P8Z.eco' + +#################### Nitrogen ############################ +eg_useNO3 =.false. +eg_nquota =.false. +################## Phosphorus ############################ +eg_usePO4 =.true. +eg_pquota =.true. +######################## Iron ############################ +eg_useFe =.false. # JDW +eg_fquota =.false. # JDW +################# Chlorophyll ############################ +eg_chlquota =.true. +################# d13C ################################### +eg_useDIC_13C =.true. + +# Tuned +eg_qminP_a = 3.33e-3 +eg_qminFe_a = 1.00e-6 +eg_qmaxFe_a = 4.00e-6 + +eg_ns = 2 + +eg_respir_a = 0 + +eg_beta_mort_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_mort_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_mort_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) +eg_beta_graz_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_graz_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_graz_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) + +eg_par_bio_red_POC_CaCO3 = 0.0285 # underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +eg_par_bio_red_POC_CaCO3_pP = 0.7440 # exponent for modifier of CaCO3:POC export ratio + +eg_nsubtime =25 +eg_n_keco =1 + +# disable ECOGEM restarts (as not currently coded up / used) +eg_ctrl_ncrst =.false. +eg_ctrl_continuing =.false. +# + +# --- DATA SAVING ---------------------------------------------- +# +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +# +# --- FORCINGS ------------------------------------------------- +# +# specify forcings +bg_par_forcing_name="pyyyyz.RpCO2_Rp13CO2" +bg_par_atm_force_scale_val_3=834.0E-06 +bg_par_atm_force_scale_val_4=-4.9 +bg_par_windspeed_file="p0055c_windspeed.dat" +#re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.5196 +# +# --- MISC ----------------------------------------------------- +# +# Set initial ocean alkalinity +bg_ocn_init_12=1.975E-03 +# +bg_par_atm_force_scale_val_3=278.0E-06 # JDW +go_imld = 1 # JDW +# --- END ------------------------------------------------------ +# diff --git a/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.p0055c.8P8Z.pal.1x.3xSST b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.p0055c.8P8Z.pal.1x.3xSST new file mode 100644 index 000000000..b9d54cae5 --- /dev/null +++ b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.p0055c.8P8Z.pal.1x.3xSST @@ -0,0 +1,119 @@ +# +# +# --- CLIMATE -------------------------------------------------- +# +# set climate feedback +ea_36=y +# +# --- BIOLOGICAL NEW PRODUCTION ------------------------------------- +# +# biological scheme ID string +bg_par_bio_prodopt="NONE" +# +# --- ORGANIC MATTER EXPORT RATIOS ----------------------------- +# +#production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS --------------------------- +# +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +#exponent for modifier of CaCO3:POC export ratio +bg_par_bio_red_POC_CaCO3_pP=0.0 +# +# --- REMINERALIZATION ----------------------------------------- +# +#DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +#initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +#remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +#initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +#remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# + +# --- ECOGEM ---------------------------------------------------------- +# +gm_ctrl_debug_init =0 +eg_ctrl_debug_eco_init =.false. + +eg_par_ecogem_plankton_file ='8P8Z.eco' + +#################### Nitrogen ############################ +eg_useNO3 =.false. +eg_nquota =.false. +################## Phosphorus ############################ +eg_usePO4 =.true. +eg_pquota =.true. +######################## Iron ############################ +eg_useFe =.false. # JDW +eg_fquota =.false. # JDW +################# Chlorophyll ############################ +eg_chlquota =.true. +################# d13C ################################### +eg_useDIC_13C =.true. + +# Tuned +eg_qminP_a = 3.33e-3 +eg_qminFe_a = 1.00e-6 +eg_qmaxFe_a = 4.00e-6 + +eg_ns = 2 + +eg_respir_a = 0 + +eg_beta_mort_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_mort_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_mort_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) +eg_beta_graz_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_graz_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_graz_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) + +eg_par_bio_red_POC_CaCO3 = 0.0285 # underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +eg_par_bio_red_POC_CaCO3_pP = 0.7440 # exponent for modifier of CaCO3:POC export ratio + +eg_nsubtime =25 +eg_n_keco =1 + +# disable ECOGEM restarts (as not currently coded up / used) +eg_ctrl_ncrst =.false. +eg_ctrl_continuing =.false. +# + +# --- DATA SAVING ---------------------------------------------- +# +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +# +# --- FORCINGS ------------------------------------------------- +# +# specify forcings +bg_par_forcing_name="pyyyyz.RpCO2_Rp13CO2" +bg_par_atm_force_scale_val_3=834.0E-06 +bg_par_atm_force_scale_val_4=-4.9 +bg_par_windspeed_file="p0055c_windspeed.dat" +#re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.5196 +# +# --- MISC ----------------------------------------------------- +# +# Set initial ocean alkalinity +bg_ocn_init_12=1.975E-03 +# +bg_par_atm_force_scale_val_3=278.0E-06 # JDW +go_imld = 1 # JDW +eg_ctrl_force_T=.true. +# limit negative biomass +eg_ctrl_limit_neg_biomass=.true. +eg_par_ecogem_force_T_file="20170811.p0055c.8P8Z.pal.3x_SST.dat" +# --- END ------------------------------------------------------ +# diff --git a/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.p0055c.8P8Z.pal.3x b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.p0055c.8P8Z.pal.3x new file mode 100644 index 000000000..b008e312b --- /dev/null +++ b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.p0055c.8P8Z.pal.3x @@ -0,0 +1,114 @@ +# +# +# --- CLIMATE -------------------------------------------------- +# +# set climate feedback +ea_36=y +# +# --- BIOLOGICAL NEW PRODUCTION ------------------------------------- +# +# biological scheme ID string +bg_par_bio_prodopt="NONE" +# +# --- ORGANIC MATTER EXPORT RATIOS ----------------------------- +# +#production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS --------------------------- +# +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +#exponent for modifier of CaCO3:POC export ratio +bg_par_bio_red_POC_CaCO3_pP=0.0 +# +# --- REMINERALIZATION ----------------------------------------- +# +#DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +#initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +#remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +#initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +#remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# + +# --- ECOGEM ---------------------------------------------------------- +# +gm_ctrl_debug_init =0 +eg_ctrl_debug_eco_init =.false. + +eg_par_ecogem_plankton_file ='8P8Z.eco' + +#################### Nitrogen ############################ +eg_useNO3 =.false. +eg_nquota =.false. +################## Phosphorus ############################ +eg_usePO4 =.true. +eg_pquota =.true. +######################## Iron ############################ +eg_useFe =.false. # JDW +eg_fquota =.false. # JDW +################# Chlorophyll ############################ +eg_chlquota =.true. +################# d13C ################################### +eg_useDIC_13C =.true. + +# Tuned +eg_qminP_a = 3.33e-3 +eg_qminFe_a = 1.00e-6 +eg_qmaxFe_a = 4.00e-6 + +eg_ns = 2 + +eg_respir_a = 0 + +eg_beta_mort_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_mort_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_mort_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) +eg_beta_graz_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_graz_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_graz_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) + +eg_par_bio_red_POC_CaCO3 = 0.0285 # underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +eg_par_bio_red_POC_CaCO3_pP = 0.7440 # exponent for modifier of CaCO3:POC export ratio + +eg_nsubtime =25 +eg_n_keco =1 + +# disable ECOGEM restarts (as not currently coded up / used) +eg_ctrl_ncrst =.false. +eg_ctrl_continuing =.false. +# +# --- DATA SAVING ---------------------------------------------- +# +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +# +# --- FORCINGS ------------------------------------------------- +# +# specify forcings +bg_par_forcing_name="pyyyyz.RpCO2_Rp13CO2" +bg_par_atm_force_scale_val_3=834.0E-06 +bg_par_atm_force_scale_val_4=-4.9 +bg_par_windspeed_file="p0055c_windspeed.dat" +#re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.5196 +# +# --- MISC ----------------------------------------------------- +# +# Set initial ocean alkalinity +bg_ocn_init_12=1.975E-03 +# +bg_par_atm_force_scale_val_3=834.0E-06 +go_imld=1 # JDW +# --- END ------------------------------------------------------ +# diff --git a/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.p0055c.8P8Z.pal.3x.1xSST b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.p0055c.8P8Z.pal.3x.1xSST new file mode 100644 index 000000000..544cef82c --- /dev/null +++ b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.p0055c.8P8Z.pal.3x.1xSST @@ -0,0 +1,118 @@ +# +# +# --- CLIMATE -------------------------------------------------- +# +# set climate feedback +ea_36=y +# +# --- BIOLOGICAL NEW PRODUCTION ------------------------------------- +# +# biological scheme ID string +bg_par_bio_prodopt="NONE" +# +# --- ORGANIC MATTER EXPORT RATIOS ----------------------------- +# +#production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS --------------------------- +# +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +#exponent for modifier of CaCO3:POC export ratio +bg_par_bio_red_POC_CaCO3_pP=0.0 +# +# --- REMINERALIZATION ----------------------------------------- +# +#DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +#initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +#remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +#initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +#remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# + +# --- ECOGEM ---------------------------------------------------------- +# +# +gm_ctrl_debug_init =0 +eg_ctrl_debug_eco_init =.false. + +eg_par_ecogem_plankton_file ='8P8Z.eco' + +#################### Nitrogen ############################ +eg_useNO3 =.false. +eg_nquota =.false. +################## Phosphorus ############################ +eg_usePO4 =.true. +eg_pquota =.true. +######################## Iron ############################ +eg_useFe =.false. # JDW +eg_fquota =.false. # JDW +################# Chlorophyll ############################ +eg_chlquota =.true. +################# d13C ################################### +eg_useDIC_13C =.true. + +# Tuned +eg_qminP_a = 3.33e-3 +eg_qminFe_a = 1.00e-6 +eg_qmaxFe_a = 4.00e-6 + +eg_ns = 2 + +eg_respir_a = 0 + +eg_beta_mort_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_mort_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_mort_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) +eg_beta_graz_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_graz_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_graz_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) + +eg_par_bio_red_POC_CaCO3 = 0.0285 # underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +eg_par_bio_red_POC_CaCO3_pP = 0.7440 # exponent for modifier of CaCO3:POC export ratio + +eg_nsubtime =25 +eg_n_keco =1 + +# disable ECOGEM restarts (as not currently coded up / used) +eg_ctrl_ncrst =.false. +eg_ctrl_continuing =.false. +# + +# --- DATA SAVING ---------------------------------------------- +# +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +# +# --- FORCINGS ------------------------------------------------- +# +# specify forcings +bg_par_forcing_name="pyyyyz.RpCO2_Rp13CO2" +bg_par_atm_force_scale_val_3=834.0E-06 +bg_par_atm_force_scale_val_4=-4.9 +bg_par_windspeed_file="p0055c_windspeed.dat" +#re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.5196 +# +# --- MISC ----------------------------------------------------- +# +# Set initial ocean alkalinity +bg_ocn_init_12=1.975E-03 +# +bg_par_atm_force_scale_val_3=834.0E-06 +go_imld = 1 # JDW +eg_ctrl_force_T=.true. +eg_par_ecogem_force_T_file="20170811.p0055c.8P8Z.pal.1x_SST.dat" +# --- END ------------------------------------------------------ +# diff --git a/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.worjh2.8P8Z.pal.1x b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.worjh2.8P8Z.pal.1x new file mode 100644 index 000000000..267a807b5 --- /dev/null +++ b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.worjh2.8P8Z.pal.1x @@ -0,0 +1,111 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set no climate feedback +ea_36=n +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# biological scheme ID string +bg_par_bio_prodopt="NONE" +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.0485 +# exponent for modifier of CaCO3:POC export ratio +bg_par_bio_red_POC_CaCO3_pP=0.7440 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- IRON ---------------------------------------------------------- +# +# +# --- ECOGEM ---------------------------------------------------------- +# +gm_ctrl_debug_init =0 +eg_ctrl_debug_eco_init =.false. + +eg_par_ecogem_plankton_file ='8P8Z.eco' + +#################### Nitrogen ############################ +eg_useNO3 =.false. +eg_nquota =.false. +################## Phosphorus ############################ +eg_usePO4 =.true. +eg_pquota =.true. +######################## Iron ############################ +eg_useFe =.false. # JDW +eg_fquota =.false. # JDW +################# Chlorophyll ############################ +eg_chlquota =.true. +################# d13C ################################### +eg_useDIC_13C =.true. + +# Tuned +eg_qminP_a = 3.33e-3 +eg_qminFe_a = 1.00e-6 +eg_qmaxFe_a = 4.00e-6 + +eg_ns = 2 + +eg_respir_a = 0 + +eg_beta_mort_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_mort_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_mort_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) +eg_beta_graz_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_graz_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_graz_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) + +eg_par_bio_red_POC_CaCO3 = 0.0285 # underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +eg_par_bio_red_POC_CaCO3_pP = 0.7440 # exponent for modifier of CaCO3:POC export ratio + +eg_nsubtime =25 +eg_n_keco =1 + +# disable ECOGEM restarts (as not currently coded up / used) +eg_ctrl_ncrst =.false. +eg_ctrl_continuing =.false. +go_imld = 1 # JDW +# +# --- DATA SAVING ----------------------------------------------------- +# +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +# +# --- FORCINGS -------------------------------------------------------- +# +# specify forcings +bg_par_forcing_name="pyyyyz.RpCO2_Rp13CO2" +bg_par_atm_force_scale_val_3=278.0E-06 +bg_par_atm_force_scale_val_4=-6.5 +# +bg_par_atm_force_scale_val_3=278.0E-06 +bg_par_data_save_level=8 +ea_36=y +# --- END ------------------------------------------------------------- +# +# diff --git a/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.worjh2.8P8Z.pal.1x.3xSST b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.worjh2.8P8Z.pal.1x.3xSST new file mode 100644 index 000000000..552eeea21 --- /dev/null +++ b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.worjh2.8P8Z.pal.1x.3xSST @@ -0,0 +1,113 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set no climate feedback +ea_36=n +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# biological scheme ID string +bg_par_bio_prodopt="NONE" +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.0485 +# exponent for modifier of CaCO3:POC export ratio +bg_par_bio_red_POC_CaCO3_pP=0.7440 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- IRON ---------------------------------------------------------- +# +# +# --- ECOGEM ---------------------------------------------------------- +# +gm_ctrl_debug_init =0 +eg_ctrl_debug_eco_init =.false. + +eg_par_ecogem_plankton_file ='8P8Z.eco' + +#################### Nitrogen ############################ +eg_useNO3 =.false. +eg_nquota =.false. +################## Phosphorus ############################ +eg_usePO4 =.true. +eg_pquota =.true. +######################## Iron ############################ +eg_useFe =.false. # JDW +eg_fquota =.false. # JDW +################# Chlorophyll ############################ +eg_chlquota =.true. +################# d13C ################################### +eg_useDIC_13C =.true. + +# Tuned +eg_qminP_a = 3.33e-3 +eg_qminFe_a = 1.00e-6 +eg_qmaxFe_a = 4.00e-6 + +eg_ns = 2 + +eg_respir_a = 0 + +eg_beta_mort_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_mort_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_mort_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) +eg_beta_graz_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_graz_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_graz_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) + +eg_par_bio_red_POC_CaCO3 = 0.0285 # underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +eg_par_bio_red_POC_CaCO3_pP = 0.7440 # exponent for modifier of CaCO3:POC export ratio + +eg_nsubtime =25 +eg_n_keco =1 + +# disable ECOGEM restarts (as not currently coded up / used) +eg_ctrl_ncrst =.false. +eg_ctrl_continuing =.false. +go_imld = 1 # JDW +# +# --- DATA SAVING ----------------------------------------------------- +# +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +# +# --- FORCINGS -------------------------------------------------------- +# +# specify forcings +bg_par_forcing_name="pyyyyz.RpCO2_Rp13CO2" +bg_par_atm_force_scale_val_3=278.0E-06 +bg_par_atm_force_scale_val_4=-6.5 +# +bg_par_atm_force_scale_val_3=278.0E-06 +bg_par_data_save_level=8 +ea_36=y +eg_ctrl_force_T=.true. +eg_par_ecogem_force_T_file="20170811.worjh2.PO4.8P8Z.pal.3x_SST.dat" +# --- END ------------------------------------------------------------- +# +# diff --git a/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.worjh2.8P8Z.pal.1x_Fe b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.worjh2.8P8Z.pal.1x_Fe new file mode 100644 index 000000000..17836faf2 --- /dev/null +++ b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.worjh2.8P8Z.pal.1x_Fe @@ -0,0 +1,137 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set no climate feedback +ea_36=n +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# biological scheme ID string +bg_par_bio_prodopt="NONE" +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.0485 +# exponent for modifier of CaCO3:POC export ratio +bg_par_bio_red_POC_CaCO3_pP=0.7440 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- IRON ---------------------------------------------------------- +# +bg_opt_geochem_Fe='hybrid' +# aeolian Fe solubility +bg_par_det_Fe_sol=0.002014275 +#exponent for aeolian Fe solubility [use 1.0 for uniform solubility] +bg_par_det_Fe_sol_exp=0.500 +# modifier of the scavenging rate of dissolved Fe +bg_par_scav_Fe_sf_POC=0.344104 +bg_par_scav_Fe_sf_CaCO3=0.000 +bg_par_scav_Fe_sf_opal=0.000 +bg_par_scav_Fe_sf_det=0.000 +# no scavenged regeneration +bg_par_scav_fremin=0.0 +# return POFe +bg_ctrl_bio_NO_fsedFe=.false. +# Variable Fe:C +bg_ctrl_bio_red_fixedFetoC=.false. +# adjust pK'(FeL) +bg_par_K_FeL_pP=11.0 +#(max) C/Fe organic matter ratio +bg_par_bio_red_POFe_POC=250000.0 +#[FeT] dependent Fe:C ratio [Ridgwell, 2001] -- power +bg_par_bio_FetoC_pP=-0.4225 +#[FeT] dependent Fe:C ratio [Ridgwell, 2001] -- scaling +bg_par_bio_FetoC_K=103684.0 +#[FeT] dependent Fe:C ratio [Ridgwell, 2001] -- constant +bg_par_bio_FetoC_C=0.0 +# +# --- ECOGEM ---------------------------------------------------------- +# +gm_ctrl_debug_init =0 +eg_ctrl_debug_eco_init =.false. + +eg_par_ecogem_plankton_file ='8P8Z.eco' + +#################### Nitrogen ############################ +eg_useNO3 =.false. +eg_nquota =.false. +################## Phosphorus ############################ +eg_usePO4 =.true. +eg_pquota =.true. +######################## Iron ############################ +eg_useFe =.true. # JDW +eg_fquota =.true. # JDW +################# Chlorophyll ############################ +eg_chlquota =.true. +################# d13C ################################### +eg_useDIC_13C =.true. + +# Tuned +eg_qminP_a = 3.33e-3 +eg_qminFe_a = 1.00e-6 +eg_qmaxFe_a = 4.00e-6 + +eg_ns = 2 + +eg_respir_a = 0 + +eg_beta_mort_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_mort_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_mort_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) +eg_beta_graz_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_graz_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_graz_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) + +eg_par_bio_red_POC_CaCO3 = 0.0285 # underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +eg_par_bio_red_POC_CaCO3_pP = 0.7440 # exponent for modifier of CaCO3:POC export ratio + +eg_nsubtime =25 +eg_n_keco =1 + +# disable ECOGEM restarts (as not currently coded up / used) +eg_ctrl_ncrst =.false. +eg_ctrl_continuing =.false. +go_imld = 1 # JDW +# +# --- DATA SAVING ----------------------------------------------------- +# +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +# +# --- FORCINGS -------------------------------------------------------- +# +# specify forcings +bg_par_forcing_name="worjh2.RpCO2_Rp13CO2.FeMahowald2006" +bg_par_atm_force_scale_val_3=278.0E-06 +bg_par_atm_force_scale_val_4=-6.5 +# +bg_par_atm_force_scale_val_3=278.0E-06 +bg_par_data_save_level=8 +ea_36=y +# --- END ------------------------------------------------------------- +# +# diff --git a/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.worjh2.8P8Z.pal.3x b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.worjh2.8P8Z.pal.3x new file mode 100644 index 000000000..7b0696f1c --- /dev/null +++ b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.worjh2.8P8Z.pal.3x @@ -0,0 +1,111 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set no climate feedback +ea_36=n +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# biological scheme ID string +bg_par_bio_prodopt="NONE" +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.0485 +# exponent for modifier of CaCO3:POC export ratio +bg_par_bio_red_POC_CaCO3_pP=0.7440 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- IRON ---------------------------------------------------------- +# +# +# --- ECOGEM ---------------------------------------------------------- +# +gm_ctrl_debug_init =0 +eg_ctrl_debug_eco_init =.false. + +eg_par_ecogem_plankton_file ='8P8Z.eco' + +#################### Nitrogen ############################ +eg_useNO3 =.false. +eg_nquota =.false. +################## Phosphorus ############################ +eg_usePO4 =.true. +eg_pquota =.true. +######################## Iron ############################ +eg_useFe =.false. # JDW +eg_fquota =.false. # JDW +################# Chlorophyll ############################ +eg_chlquota =.true. +################# d13C ################################### +eg_useDIC_13C =.true. + +# Tuned +eg_qminP_a = 3.33e-3 +eg_qminFe_a = 1.00e-6 +eg_qmaxFe_a = 4.00e-6 + +eg_ns = 2 + +eg_respir_a = 0 + +eg_beta_mort_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_mort_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_mort_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) +eg_beta_graz_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_graz_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_graz_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) + +eg_par_bio_red_POC_CaCO3 = 0.0285 # underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +eg_par_bio_red_POC_CaCO3_pP = 0.7440 # exponent for modifier of CaCO3:POC export ratio + +eg_nsubtime =25 +eg_n_keco =1 + +# disable ECOGEM restarts (as not currently coded up / used) +eg_ctrl_ncrst =.false. +eg_ctrl_continuing =.false. +go_imld = 1 # JDW +# +# --- DATA SAVING ----------------------------------------------------- +# +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +# +# --- FORCINGS -------------------------------------------------------- +# +# specify forcings +bg_par_forcing_name="pyyyyz.RpCO2_Rp13CO2" +bg_par_atm_force_scale_val_3=278.0E-06 +bg_par_atm_force_scale_val_4=-6.5 +# +bg_par_atm_force_scale_val_3=834.0E-06 +bg_par_data_save_level=8 +ea_36=y +# --- END ------------------------------------------------------------- +# +# diff --git a/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.worjh2.8P8Z.pal.3x.1xSST b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.worjh2.8P8Z.pal.3x.1xSST new file mode 100644 index 000000000..674c1c781 --- /dev/null +++ b/genie-userconfigs/MS/wilsonetal.2018/wilsonetal.worjh2.8P8Z.pal.3x.1xSST @@ -0,0 +1,139 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set no climate feedback +ea_36=n +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# biological scheme ID string +bg_par_bio_prodopt="NONE" +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.0485 +# exponent for modifier of CaCO3:POC export ratio +bg_par_bio_red_POC_CaCO3_pP=0.7440 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- IRON ---------------------------------------------------------- +# +bg_opt_geochem_Fe='hybrid' +# aeolian Fe solubility +bg_par_det_Fe_sol=0.002014275 +#exponent for aeolian Fe solubility [use 1.0 for uniform solubility] +bg_par_det_Fe_sol_exp=0.500 +# modifier of the scavenging rate of dissolved Fe +bg_par_scav_Fe_sf_POC=0.344104 +bg_par_scav_Fe_sf_CaCO3=0.000 +bg_par_scav_Fe_sf_opal=0.000 +bg_par_scav_Fe_sf_det=0.000 +# no scavenged regeneration +bg_par_scav_fremin=0.0 +# return POFe +bg_ctrl_bio_NO_fsedFe=.false. +# Variable Fe:C +bg_ctrl_bio_red_fixedFetoC=.false. +# adjust pK'(FeL) +bg_par_K_FeL_pP=11.0 +#(max) C/Fe organic matter ratio +bg_par_bio_red_POFe_POC=250000.0 +#[FeT] dependent Fe:C ratio [Ridgwell, 2001] -- power +bg_par_bio_FetoC_pP=-0.4225 +#[FeT] dependent Fe:C ratio [Ridgwell, 2001] -- scaling +bg_par_bio_FetoC_K=103684.0 +#[FeT] dependent Fe:C ratio [Ridgwell, 2001] -- constant +bg_par_bio_FetoC_C=0.0 +# +# --- ECOGEM ---------------------------------------------------------- +# +gm_ctrl_debug_init =0 +eg_ctrl_debug_eco_init =.false. + +eg_par_ecogem_plankton_file ='8P8Z.eco' + +#################### Nitrogen ############################ +eg_useNO3 =.false. +eg_nquota =.false. +################## Phosphorus ############################ +eg_usePO4 =.true. +eg_pquota =.true. +######################## Iron ############################ +eg_useFe =.false. # JDW +eg_fquota =.false. # JDW +################# Chlorophyll ############################ +eg_chlquota =.true. +################# d13C ################################### +eg_useDIC_13C =.true. + +# Tuned +eg_qminP_a = 3.33e-3 +eg_qminFe_a = 1.00e-6 +eg_qmaxFe_a = 4.00e-6 + +eg_ns = 2 + +eg_respir_a = 0 + +eg_beta_mort_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_mort_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_mort_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) +eg_beta_graz_a = 0.8 # Maximum fraction to DOM as ESD --> zero +eg_beta_graz_b = 0.4 # Minimum fraction to DOM as ESD --> infinity +eg_beta_graz_c = 100.0 # Size at 50:50 partition (default = 100 µm^3) + +eg_par_bio_red_POC_CaCO3 = 0.0285 # underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +eg_par_bio_red_POC_CaCO3_pP = 0.7440 # exponent for modifier of CaCO3:POC export ratio + +eg_nsubtime =25 +eg_n_keco =1 + +# disable ECOGEM restarts (as not currently coded up / used) +eg_ctrl_ncrst =.false. +eg_ctrl_continuing =.false. +go_imld = 1 # JDW +# +# --- DATA SAVING ----------------------------------------------------- +# +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +# +# --- FORCINGS -------------------------------------------------------- +# +# specify forcings +bg_par_forcing_name="pyyyyz.RpCO2_Rp13CO2" +bg_par_atm_force_scale_val_3=278.0E-06 +bg_par_atm_force_scale_val_4=-6.5 +# +bg_par_atm_force_scale_val_3=834.0E-06 +bg_par_data_save_level=8 +ea_36=y +eg_ctrl_force_T=.true. +eg_par_ecogem_force_T_file="20170811.worjh2.PO4.8P8Z.pal.1x_SST.dat" +# --- END ------------------------------------------------------------- +# +#