new nearest_lat mapalgo in cdeps - does not requre a mesh

Author: mvertens

bld/namelist_files/namelist_defaults_cam.xml

@@ -2012,6 +2012,9 @@
 
 <!-- co2 surface source (overwrite co2_diag sent to land and ocean) -->
 <co2_surface_source>.false.</co2_surface_source>
+<co2_surface_source_mapalgo>nearest_lat</co2_surface_source_mapalgo>
+<co2_surface_source_mesh_filename>none</co2_surface_source_mesh_filename>
+<co2_surface_source_lat_dimname>lat</co2_surface_source_lat_dimname>
 
 <!-- nitrogen deposition to surface components -->
 <stream_ndep_varlist>NDEP_NHx_month:NDEP_NOy_month</stream_ndep_varlist>

bld/namelist_files/namelist_definition.xml

@@ -10297,10 +10297,22 @@ Default: path
   If TRUE, set co2_diag sent to land and ocean to value read from co2_surface_source_data_filename
   Default: FALSE
 </entry>
+<entry id="co2_surface_source_mapalgo" type="char*256" category="CO2 surface source"
+       group="co2_surface_source_stream_nl" valid_values="nearest_lat,bilinear,consf" >
+  Mapping algorithm to map co2 data to model grid. Nearest lat assumes
+  that a mesh file is not needed and that the forcing data is lat/time or lat/lev/time.
+  In this case the mesh is set to none.
+  Default: nearest_lat
+</entry>
+<entry id="stream_co2_surface_source_lat_dimname" type="char*256" input_pathname="abs" category="CO2 surface source"
+       group="co2_surface_source_stream_nl" valid_values="" >
+  Latitude dimension name on stream_co2_surface_source_data_filename. Only used if co2_surface_source_mapalgo is 'nearest_lat'
+  Default: lat
+</entry>
 <entry id="stream_co2_surface_source_mesh_filename" type="char*256" input_pathname="abs" category="CO2 surface source"
        group="co2_surface_source_stream_nl" valid_values="" >
-  Grid mesh file corresponding to stream_co2_surface_source_data_filename.
-  Must be set by user.
+  Grid mesh file corresponding to stream_co2_surface_source_data_filename. 
+  Default: none
 </entry>
 <entry id="stream_co2_surface_source_data_filename" type="char*256" input_pathname="abs" category="CO2 surface source"
        group="co2_surface_source_stream_nl" valid_values="" >

src/cpl/nuopc/atm_stream_co2.F90

@@ -27,6 +27,8 @@ module atm_stream_co2
   type(shr_strdata_type) :: sdat_co2_surface_source     ! input data stream
 
   ! namelist variables
+  character(len=CL) :: stream_co2_surface_source_mapalgo
+  character(len=CL) :: stream_co2_surface_source_lat_dimname
   character(len=CL) :: stream_co2_surface_source_mesh_filename
   character(len=CL) :: stream_co2_surface_source_data_filename
   character(len=CL) :: stream_co2_surface_source_data_varname ! variable name for co2_surface_source on stream file(s)
@@ -72,6 +74,8 @@ subroutine stream_co2_surface_source_readnl(nlfile)
 
     ! Default values for namelist
     co2_surface_source = .false.
+    stream_co2_surface_source_mapalgo       = 'none'
+    stream_co2_surface_source_lat_dimname   = 'lat'
     stream_co2_surface_source_data_filename = ' '
     stream_co2_surface_source_mesh_filename = ' '
     stream_co2_surface_source_data_varname  = ' '
@@ -111,6 +115,12 @@ subroutine stream_co2_surface_source_readnl(nlfile)
     if (ierr /= 0) call endrun(trim(subname)//": FATAL: mpi_bcast: co2_surface_source")
 
     if (co2_surface_source) then
+       call mpi_bcast(stream_co2_surface_source_mapalgo, &
+            len(stream_co2_surface_source_mapalgo), mpi_character, masterprocid, mpicom, ierr)
+       if (ierr /= 0) call endrun(trim(subname)//": FATAL: mpi_bcast: stream_co2_surface_source_mapalgo")
+       call mpi_bcast(stream_co2_surface_source_lat_dimname, &
+            len(stream_co2_surface_source_lat_dimname), mpi_character, masterprocid, mpicom, ierr)
+       if (ierr /= 0) call endrun(trim(subname)//": FATAL: mpi_bcast: stream_co2_surface_source_lat_dimname")
        call mpi_bcast(stream_co2_surface_source_mesh_filename, &
             len(stream_co2_surface_source_mesh_filename), mpi_character, masterprocid, mpicom, ierr)
        if (ierr /= 0) call endrun(trim(subname)//": FATAL: mpi_bcast: stream_co2_surface_source_mesh_filename")
@@ -137,6 +147,10 @@ subroutine stream_co2_surface_source_readnl(nlfile)
        write(iulog,'(a)')  ' '
        if (co2_surface_source) then
           write(iulog,'(2a)')    subname,' co2 surface source override settings::'
+          write(iulog,'(3a)')    subname,'  stream_co2_surface_source_mapalgo       = ',&
+               trim(stream_co2_surface_source_mapalgo)
+          write(iulog,'(3a)')    subname,'  stream_co2_surface_source_lat_dimname   = ',&
+               trim(stream_co2_surface_source_lat_dimname)
           write(iulog,'(3a)')    subname,'  stream_co2_surface_source_data_filename = ',&
                trim(stream_co2_surface_source_data_filename)
           write(iulog,'(3a)')    subname,'  stream_co2_surface_source_mesh_filename = ',&
@@ -188,12 +202,13 @@ subroutine stream_co2_surface_source_init(rc)
             stream_fldlistFile  = (/stream_co2_surface_source_data_varname/),        &
             stream_fldListModel = (/stream_co2_surface_source_data_varname/),        &
             stream_lev_dimname  = 'null',                                            &
-            stream_mapalgo      = 'bilinear',                                        &
+            stream_mapalgo      = trim(stream_co2_surface_source_mapalgo),           &
             stream_offset       = 0,                                                 &
             stream_taxmode      = trim(stream_co2_surface_source_taxmode),           &
             stream_dtlimit      = 1.0e30_r8,                                         &
             stream_tintalgo     = 'linear',                                          &
             stream_name         = 'CO2_SURFACE_SOURCE data ',                        &
+            stream_lat_dimname  = trim(stream_co2_surface_source_lat_dimname),       & 
             rc                  = rc)
        if (ChkErr(rc,__LINE__,u_FILE_u)) return