updates to ed relevant namelist build options

cime/machines/env_mach_specific.lawrencium-lr2

@@ -23,9 +23,9 @@ module load cmake
 module load perl xml-libxml switch
 
 if ( $COMPILER == "intel" ) then
-    module load intel/2015.0.090
+    module load intel/2015.6.233
     module load openmpi
-    module load netcdf/4.3.2-intel-p
+    module load netcdf/4.4.0-intel-p
     module load mkl
 endif
 

cime/machines/env_mach_specific.lawrencium-lr3

@@ -0,0 +1,37 @@
+#! /bin/csh -f
+
+set CESM_REPO = `./xmlquery CCSM_REPOTAG -value`
+if($status == 0) then
+  set COMPILER            = `./xmlquery  COMPILER          -value`
+  set MPILIB              = `./xmlquery  MPILIB            -value`
+  set DEBUG               = `./xmlquery  DEBUG             -value`
+  set OS                  = `./xmlquery  OS                -value`
+  set PROFILE_PAPI_ENABLE = `./xmlquery  PROFILE_PAPI_ENABLE -value`
+  set PROJECT             = `./xmlquery  PROJECT           -value`
+else
+  echo $0 using settings from environment:
+endif
+echo "COMPILER=$COMPILER"
+echo "MPILIB=$MPILIB"
+echo "DEBUG=$DEBUG"
+echo "OS=$OS"
+echo "PROJECT=$PROJECT"
+
+source /etc/profile.d/modules.csh
+module purge
+module load cmake
+module load perl xml-libxml switch
+
+if ( $COMPILER == "intel" ) then
+    module load intel/2015.6.233
+    module load openmpi
+    module load netcdf/4.4.0-intel-p
+    module load mkl
+endif
+
+setenv NETCDF_PATH $NETCDF_DIR
+#-------------------------------------------------------------------------------
+
+limit coredumpsize 1000000
+limit stacksize unlimited
+

cime/scripts/Tools/config_compsets.xml

@@ -845,7 +845,7 @@ CAM[45]%L60  => CAM with 60 layers and full gravity wave spectrum:
 <CLM_BLDNML_OPTS compset="_CLM45%[^_]*BGCDV"     >-bgc bgc -dynamic_vegetation</CLM_BLDNML_OPTS>
 <CLM_BLDNML_OPTS compset="_CLM45%[^_]*CNDV-CROP" >-bgc cn -dynamic_vegetation -crop</CLM_BLDNML_OPTS>
 <CLM_BLDNML_OPTS compset="_CLM45%[^_]*BGCDV-CROP">-bgc bgc -dynamic_vegetation -crop</CLM_BLDNML_OPTS>
-<CLM_BLDNML_OPTS compset="_CLM45%[^_]*ED"        >-ed_mode</CLM_BLDNML_OPTS>
+<CLM_BLDNML_OPTS compset="_CLM45%[^_]*ED"        >-ed_mode -no-megan</CLM_BLDNML_OPTS>
 
 <CLM_BLDNML_OPTS compset="_CLM50%[^_]*SP"        >-bgc sp</CLM_BLDNML_OPTS>
 <CLM_BLDNML_OPTS compset="_CLM50%[^_]*BGC"       >-bgc bgc</CLM_BLDNML_OPTS>

components/clm/bld/CLMBuildNamelist.pm

@@ -159,7 +159,9 @@ OPTIONS
                               (can ONLY be turned on when BGC type is 'cn' or 'bgc')
                               This turns on the namelist variable: use_cndv
      -ed_mode                 Turn ED (Ecosystem Demography) : [on | off] (default is off)
-                              Sets the namelist variable use_ed and use_spit_fire.
+                              Sets the namelist variable use_ed, checks for use_spit_fire, 
+                              use_vertsoilc, use_century_decomp and use_lch4
+			      Also specifies defaults for these
      -glc_present             Set to true if the glc model is present (not sglc).
                               This is used for error-checking, to make sure other options are
                               set appropriately.
@@ -730,7 +732,6 @@ sub setup_cmdl_ed_mode {
     $var = "use_ed";
     $nl_flags->{$var} = ".false.";
     if ($nl_flags->{'ed_mode'} eq 1) {
-      message("Using ED (Ecosystem Demography).");
       $val = ".true.";
       $nl_flags->{$var} = $val;
     }
@@ -745,13 +746,42 @@ sub setup_cmdl_ed_mode {
         fatal_error("$var has a value ($val) that is NOT valid. Valid values are: @valid_values\n");
       }
 
-      $var = "use_ed_spit_fire";
-      $nl->set_variable_value($group, $var, $val);
-      if ( ! $definition->is_valid_value($var, $val) ) {
-        my @valid_values   = $definition->get_valid_values( $var );
-        fatal_error("$var has a value ($val) that is NOT valid. Valid values are: @valid_values\n");
+      # This section is a place-holder to test for modules that are not allowed with ED
+      # the defaults which are set in the logic section of the namelist builder will
+      # automatically set these correctly (well that is the assumption), but here we
+      # want to set a catch to fail and warn users if they explicitly set incompatible user namelist
+      # options
+
+#      my $var = "use_somevar";
+#      $val = $nl_flags->{$var};
+#      if ( defined($nl->get_value($var))  ) {
+#	  if ( $nl->get_value($var) == ".true." ) { 
+#	      fatal_error("$var was set to .true., which is incompatible when -ed_mode option is used.\n");
+#	  }
+#      }
+
+
+      # The following variables may be set by the user and are compatible with use_ed
+      # no need to set defaults, covered in a different routine
+      my @list  = (  "use_ed_spit_fire", "use_vertsoilc", "use_century_decomp", "use_lch4" );
+      foreach my $var ( @list ) {
+	  if ( defined($nl->get_value($var))  ) {
+	      $nl_flags->{$var} = $nl->get_value($var);
+	      $val = $nl_flags->{$var};
+	      my $group = $definition->get_group_name($var);
+	      $nl->set_variable_value($group, $var, $val);
+	      if (  ! $definition->is_valid_value( $var, $val ) ) {
+		  my @valid_values   = $definition->get_valid_values( $var );
+		  fatal_error("$var has a value ($val) that is NOT valid. Valid values are: @valid_values\n");
+	      }
+	  }
       }
+
+#      add_default($test_files, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'use_vertsoilc', 'use_ed'=>$nl_flags->{'use_ed'} );
+
+
     } else {
+	# we only dis-allow ed_spit_fire with non-ed runs
        $var = "use_ed_spit_fire";
        if ( defined($nl->get_value($var)) ) {
            fatal_error("$var is being set, but can ONLY be set when -ed_mode option is used.\n");
@@ -771,6 +801,7 @@ sub setup_cmdl_bgc {
 
   $val = $opts->{$var};
   $nl_flags->{'bgc_mode'} = $val;
+
   if ( $physv->as_long() == $physv->as_long("clm4_0") ) {
     if ( $nl_flags->{'bgc_mode'} ne "default" ) {
       fatal_error("-bgc option used with clm4_0 physics. -bgc can ONLY be used with clm4_5/clm5_0 physics");
@@ -802,29 +833,32 @@ sub setup_cmdl_bgc {
     if ( defined($nl->get_value("use_cn")) && ($nl_flags->{'use_cn'} ne $nl->get_value("use_cn")) ) {
       fatal_error("The namelist variable use_cn is inconsistent with the -bgc option");
     }
-    # If the variable has already been set use it, if not set to the value defined by the bgc_mode
-    my @list  = (  "use_lch4", "use_nitrif_denitrif", "use_vertsoilc", "use_century_decomp" );
-    my $ndiff = 0;
-    foreach my $var ( @list ) {
-       if ( ! defined($nl->get_value($var))  ) {
-          $nl_flags->{$var} = $setting;
-       } else {
-          if ( $nl->get_value($var) ne $setting ) {
-             $ndiff += 1;
-          }
-          $nl_flags->{$var} = $nl->get_value($var);
-       }
-       $val = $nl_flags->{$var};
-       my $group = $definition->get_group_name($var);
-       $nl->set_variable_value($group, $var, $val);
-       if (  ! $definition->is_valid_value( $var, $val ) ) {
-         my @valid_values   = $definition->get_valid_values( $var );
-         fatal_error("$var has a value ($val) that is NOT valid. Valid values are: @valid_values\n");
-       }
-    }
-    # If all the variables are different report it as an error
-    if ( $ndiff == ($#list + 1) ) {
-       fatal_error("You are contradicting the -bgc setting with the namelist variables: @list" );
+
+    if ( $opts->{"ed_mode"} != 1) {
+	# If the variable has already been set use it, if not set to the value defined by the bgc_mode
+	my @list  = (  "use_lch4", "use_nitrif_denitrif", "use_vertsoilc", "use_century_decomp" );
+	my $ndiff = 0;
+	foreach my $var ( @list ) {
+	    if ( ! defined($nl->get_value($var))  ) {
+		$nl_flags->{$var} = $setting;
+	    } else {
+		if ( $nl->get_value($var) ne $setting ) {
+		    $ndiff += 1;
+		}
+		$nl_flags->{$var} = $nl->get_value($var);
+	    }
+	    $val = $nl_flags->{$var};
+	    my $group = $definition->get_group_name($var);
+	    $nl->set_variable_value($group, $var, $val);
+	    if (  ! $definition->is_valid_value( $var, $val ) ) {
+		my @valid_values   = $definition->get_valid_values( $var );
+		fatal_error("$var has a value ($val) that is NOT valid. Valid values are: @valid_values\n");
+	    }
+	}
+	# If all the variables are different report it as an error
+	if ( $ndiff == ($#list + 1) ) {
+	    fatal_error("You are contradicting the -bgc setting with the namelist variables: @list" );
+	}
     }
 
     # Now set use_cn
@@ -1427,6 +1461,7 @@ sub process_namelist_inline_logic {
   setup_logic_supplemental_nitrogen($opts->{'test'}, $nl_flags, $definition, $defaults, $nl, $physv);
   setup_logic_snowpack($opts->{'test'}, $nl_flags, $definition, $defaults, $nl, $physv);
   setup_logic_atm_forcing($opts->{'test'}, $nl_flags, $definition, $defaults, $nl, $physv);
+  setup_logic_ed($opts->{'test'}, $nl_flags, $definition, $defaults, $nl, $physv);
 
   #########################################
   # namelist group: clm_humanindex_inparm #
@@ -2736,6 +2771,23 @@ sub setup_logic_atm_forcing {
 
 #-------------------------------------------------------------------------------
 
+sub setup_logic_ed {
+    #
+    # Set some default options related to Ecosystem Demography
+    #                                                                                                                                                            
+    my ($test_files, $nl_flags, $definition, $defaults, $nl, $physv) = @_;
+
+    if ($physv->as_long() >= $physv->as_long("clm4_5") && value_is_true( $nl_flags->{'use_ed'})  ) {
+	add_default($test_files, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'use_vertsoilc', 'use_ed'=>$nl_flags->{'use_ed'} );
+	add_default($test_files, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'use_century_decomp', 'use_ed'=>$nl_flags->{'use_ed'} );
+ 	add_default($test_files, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'use_ed_spit_fire', 'use_ed'=>$nl_flags->{'use_ed'} );
+	add_default($test_files, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'use_lch4', 'use_ed'=>$nl_flags->{'use_ed'} );
+	add_default($test_files, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'use_nitrif_denitrif', 'use_ed'=>$nl_flags->{'use_ed'} );
+    }
+}
+
+#-------------------------------------------------------------------------------                     
+
 sub write_output_files {
   my ($opts, $nl_flags, $defaults, $nl, $physv) = @_;
 
@@ -2840,6 +2892,7 @@ sub add_default {
     $var = $1;
   }
   # Query the definition to find which group the variable belongs to.  Exit if not found.
+
   my $group = $definition->get_group_name($var);
   unless ($group) {
     my $fname = $definition->get_file_name();
@@ -2860,6 +2913,7 @@ sub add_default {
     # in %settings to the get_value method to be used as attributes that are matched
     # when looking for default values.
     else {
+
       $val = $defaults->get_value($var, \%settings);
 
       # Truncate model_version appropriately

components/clm/bld/namelist_files/namelist_defaults_clm4_5.xml

@@ -1612,4 +1612,15 @@ lnd/clm2/surfdata_map/surfdata_ne30np4_simyr1850_c141219.nc</fsurdat>
 
 <use_aereoxid_prog use_cn=".true.">.true.</use_aereoxid_prog>
 
+<!-- =========================================  -->
+<!-- Defaults for use_ed = .true.               -->
+<!-- =========================================  -->
+
+<use_vertsoilc use_ed=".true.">.true.</use_vertsoilc>
+<use_century_decomp use_ed=".true.">.true.</use_century_decomp>
+<use_ed_spit_fire use_ed=".true.">.true.</use_ed_spit_fire>
+<use_lch4 use_ed=".true.">.false.</use_lch4>
+<use_nitrif_denitrif use_ed=".true.">.false.</use_nitrif_denitrif>
+
+
 </namelist_defaults>

components/clm/bld/namelist_files/namelist_definition_clm4_5.xml

@@ -427,13 +427,13 @@ Requires the CN model to work (either CN or CNDV).
 <entry id="use_vertsoilc" type="logical" category="bgc"
        group="clm_inparm" valid_values="" value=".false.">
 Turn on vertical soil carbon.
-Requires the CN model to work (either CN or CNDV).
+Requires the CN or ED model to work (either CN or CNDV).
 </entry>
 
 <entry id="use_century_decomp" type="logical" category="bgc"
        group="clm_inparm" valid_values="" value=".false.">
 Use parameters for decomposition from the CENTURY Carbon model
-Requires the CN model to work (either CN or CNDV).
+Requires the CN or ED model to work (either CN or CNDV).
 </entry>
 
 <entry id="use_extralakelayers" type="logical" category="physics"

components/clm/bld/test_build_namelist/t/input/namelist_defaults_clm4_5_test.xml

@@ -75,6 +75,18 @@ attributes from the config_cache.xml file (with keys converted to upper-case).
 <!-- Plant function types (relative to {csmdata}) -->
 <paramfile       >lnd/clm2/paramdata/clm_params.c130821.nc</paramfile>
 
+
+<!--- ========================================  -->
+<!--- Defaults for use_ed = .true.              -->
+<!--  ========================================  -->
+
+<use_vertsoilc use_ed=".true.">.true.</use_vertsoilc>
+<use_century_decomp use_ed=".true."  >.true.</use_century_decomp>
+<use_ed_spit_fire use_ed=".true.">.true.</use_ed_spit_fire>
+<use_lch4 use_ed=".true.">.false.</use_lch4>
+<use_nitrif_denitrif use_ed=".true.">.false.</use_nitrif_denitrif>
+
+
 <!--
 
      Original initial condition file at 1 degree resolution
@@ -1478,4 +1490,5 @@ lnd/clm2/surfdata_map/surfdata_ne30np4_simyr1850_c130927.nc</fsurdat>
 
 <use_aereoxid_prog use_cn=".true.">.true.</use_aereoxid_prog>
 
+
 </namelist_defaults>

components/clm/bld/test_build_namelist/t/template_test_XXX.pm

@@ -23,8 +23,14 @@ sub startup : Test(startup => 3) {
   $self->{definition} = Build::NamelistDefinition->new("t/input/namelist_definition_clm4_5_test.xml");
   isnt($self->{definition}, undef, (caller(0))[3] . " : namelist_definition object created.");
 
-  $self->{defaults} = Build::NamelistDefaults->new("t/input/namelist_defaults_clm4_5_test.xml");
+  $self->{defaults} = Build::NamelistDefaults->new("t/input/namelist_defaults_clm4_5_test.xml",$self->{cfg});
   isnt($self->{defaults}, undef,  (caller(0))[3] . " : namelist_defaults object created.");
+
+#  The next line may be usefull (set startup => 4) in arg call above:
+#  $self->{physv} = config_files::clm_phys_vers->new( $self->{cfg}->get('phys') );
+#  isnt($self->{physv}, undef,  (caller(0))[3] . " : phys_vers object created.");
+
+
 }
 
 sub shutdown : Test(shutdown) {