Skip to content

Commit

Permalink
Merge pull request #1501 from ESMCI/mvertens/somfix
Browse files Browse the repository at this point in the history 
Fix for SOM bug introduced in PR #1451 This PR fixes the SOM bug introduced in PR#1451

Test suite: scripts_regression_tests and

verified that the following two tests passed functionality and were bfb with d4dfbf1 (the CIME PR
before the aquaplanet SOM changes were introduced)
ERS_D_Ld5.f19_g16.ETEST.yellowstone_intel.cice-default
ERS_Ld5.f19_g16.ETEST.yellowstone_intel.cice-default
verified that the following two tests passed functionality (the QSC4 compset is a new CAM compset that has aquaplanet SOM functionality)
ERS_D_Ld5.f09_f09.QSC4.yellowstone_intel
ERS_Ld5.f09_f09.QSC4.yellowstone_intel
Test baseline: d4dbf1
Test namelist changes: None
Test status: bit for bit
Fixes #1497
User interface changes?: None
Code review:jedwards
  • Loading branch information
@jedwards4b
jedwards4b authored May 9, 2017
2 parents 3beac4a + f6a257a commit 10fbc43
Show file tree
Hide file tree
Showing 2 changed files with 61 additions and 26 deletions.
31 changes: 17 additions & 14 deletions src/components/data_comps/docn/cime_config/config_component.xml
View file
Original file line number Diff line number Diff line change
Expand Up @@ -41,7 +41,7 @@
driver. The atmosphere/ocean fluxes are computed in the
coupler. Therefore, the data ocean model does not compute fluxes like
the data ice model. DOCN has two distinct modes of operation. It can
arun as a pure data model, reading in ocean SSTs (normally
run as a pure data model, reading in ocean SSTs (normally
climatological) from input datasets, performing time/spatial
interpolations, and passing these to the coupler. Alternatively,
DOCN can compute updated SSTs by running as a slab ocean model where
Expand All @@ -50,18 +50,14 @@
from the driver.
--- A setting of prescribed assumes the only field in the input stream is SST.
It also assumes the SST is in Celsius and must be converted to Kelvin.
All other fields are set to zero except for ocean salinity, which
is set to a constant reference salinity value.
Normally the ice fraction data is found in the same data files that
provide SST data to the data ocean model. They are normally found in
the same file because the SST and ice fraction data are derived from
the same observational data sets and are consistent with each other.
to the data ocean model. They are normally found in the same file
because the SST and ice fraction data are derived from the same
observational data sets and are consistent with each other.
--- A setting of som (slab ocean model) mode is a prognostic mode. This mode
computes a prognostic sea surface temperature and a freeze/melt
potential (surface Q-flux) used by the sea ice model. This
All other fields are set to zero except for ocean salinity, which is set to a
constant reference salinity value. Normally the ice fraction data is found in
the same data files that provide SST data to the data ocean model. They are
normally found in the same file because the SST and ice fraction data are derived
from the same observational data sets and are consistent with each other.
--- Settings of som (slab ocean model) or som_aquap (aquaplanet slab ocean) are
prognostic modes which compute a prognostic sea surface temperature and a
freeze/melt potential (surface Q-flux) used by the sea ice model. This
calculation requires an external SOM forcing data file that includes
ocean mixed layer depths and bottom-of-the-slab Q-fluxes.
Scientifically appropriate bottom-of-the-slab Q-fluxes are normally
Expand All @@ -71,7 +67,14 @@
appropriate and is provided for testing and development purposes only.
Users must create scientifically appropriate data for their particular
application. A tool is available to derive valid SOM forcing.
Default is prescribed.</desc>
--- A setting of sst_aquapN (where “N” is an integer index value) is a
type of prescribed SST mode used specifically for an aquaplanet setup in
which global SSTs correspond to an analytic form set by the index value.
Currently, indices for 10 SST profiles are supported [e.g., index 3 corresponds
to the “QOBS” profile of Neale and Hoskins (2001, Atmos. Sci. Lett.)].
With source code modifications, it is possible for users to create their own
analytic SST distributions and match them to indices 11 or greater.
</desc>
</entry>

<entry id="DOCN_SOM_FILENAME">
Expand Down
56 changes: 44 additions & 12 deletions src/components/data_comps/docn/docn_comp_mod.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -79,6 +79,7 @@ module docn_comp_mod
integer(IN) :: kt,ks,ku,kv,kdhdx,kdhdy,kq,kswp ! field indices
integer(IN) :: kswnet,klwup,klwdn,ksen,klat,kmelth,ksnow,krofi
integer(IN) :: kh,kqbot
integer(IN) :: index_lat, index_lon

type(shr_strdata_type) :: SDOCN
type(mct_rearr) :: rearr
Expand Down Expand Up @@ -148,8 +149,6 @@ subroutine docn_comp_init( EClock, cdata, x2o, o2x, NLFilename )
integer(IN) :: nunit ! unit number
integer(IN) :: kmask ! field reference
integer(IN) :: kfrac ! field reference
integer(IN) :: klat ! field reference
integer(IN) :: klon ! field reference
logical :: ocn_present ! flag
logical :: ocn_prognostic ! flag
logical :: ocnrof_prognostic ! flag
Expand Down Expand Up @@ -360,7 +359,7 @@ subroutine docn_comp_init( EClock, cdata, x2o, o2x, NLFilename )
ocnrof_prognostic = .true.
endif

if (trim(ocn_mode) == 'SOM') then
if (trim(ocn_mode) == 'SOM' .or. trim(ocn_mode) == 'SOM_AQUAP') then
ocn_prognostic = .true.
endif

Expand Down Expand Up @@ -466,11 +465,11 @@ subroutine docn_comp_init( EClock, cdata, x2o, o2x, NLFilename )
kmask = mct_aVect_indexRA(ggrid%data,'mask')
imask(:) = nint(ggrid%data%rAttr(kmask,:))

klon = mct_aVect_indexRA(ggrid%data,'lon')
xc(:) = ggrid%data%rAttr(klon,:)
index_lon = mct_aVect_indexRA(ggrid%data,'lon')
xc(:) = ggrid%data%rAttr(index_lon,:)

klat = mct_aVect_indexRA(ggrid%data,'lat')
yc(:) = ggrid%data%rAttr(klat,:)
index_lat = mct_aVect_indexRA(ggrid%data,'lat')
yc(:) = ggrid%data%rAttr(index_lat,:)

call t_stopf('docn_initmctavs')

Expand Down Expand Up @@ -508,7 +507,7 @@ subroutine docn_comp_init( EClock, cdata, x2o, o2x, NLFilename )
endif
endif
call shr_mpi_bcast(exists,mpicom,'exists')
if (trim(ocn_mode) == 'SOM') then
if (trim(ocn_mode) == 'SOM' .or. trim(ocn_mode) == 'SOM_AQUAP') then
if (my_task == master_task) write(logunit,F00) ' reading ',trim(rest_file)
call shr_pcdf_readwrite('read',iosystem,SDOCN%io_type,trim(rest_file),mpicom,gsmap,rf1=somtp,rf1n='somtp')
endif
Expand Down Expand Up @@ -705,7 +704,7 @@ subroutine docn_comp_run( EClock, cdata, x2o, o2x)
o2x%rAttr(kswp ,n) = swp
enddo

case('SOM','SOM_AQUAP')
case('SOM')
lsize = mct_avect_lsize(o2x)
do n = 1,SDOCN%nstreams
call shr_dmodel_translateAV(SDOCN%avs(n),avstrm,avifld,avofld,rearr)
Expand Down Expand Up @@ -743,6 +742,41 @@ subroutine docn_comp_run( EClock, cdata, x2o, o2x)
enddo
endif ! firstcall

case('SOM_AQUAP')
lsize = mct_avect_lsize(o2x)
do n = 1,SDOCN%nstreams
call shr_dmodel_translateAV(SDOCN%avs(n),avstrm,avifld,avofld,rearr)
enddo
if (firstcall) then
do n = 1,lsize
if (.not. read_restart) then
somtp(n) = o2x%rAttr(kt,n) + TkFrz
endif
o2x%rAttr(kt,n) = somtp(n)
o2x%rAttr(kq,n) = 0.0_R8
enddo
else ! firstcall
tfreeze = shr_frz_freezetemp(o2x%rAttr(ks,:)) + TkFrz
do n = 1,lsize
!--- pull out h from av for resuse below ---
hn = avstrm%rAttr(kh,n)
!--- compute new temp ---
o2x%rAttr(kt,n) = somtp(n) + &
(x2o%rAttr(kswnet,n) + & ! shortwave
x2o%rAttr(klwup ,n) + & ! longwave
x2o%rAttr(klwdn ,n) + & ! longwave
x2o%rAttr(ksen ,n) + & ! sensible
x2o%rAttr(klat ,n) + & ! latent
x2o%rAttr(kmelth,n) - & ! ice melt
avstrm%rAttr(kqbot ,n) - & ! flux at bottom
(x2o%rAttr(ksnow,n)+x2o%rAttr(krofi,n))*latice) * & ! latent by prec and roff
dt/(cpsw*rhosw*hn)
!--- compute ice formed or melt potential ---
o2x%rAttr(kq,n) = (tfreeze(n) - o2x%rAttr(kt,n))*(cpsw*rhosw*hn)/dt ! ice formed q>0
somtp(n) = o2x%rAttr(kt,n) ! save temp
enddo
endif ! firstcall

end select

call t_stopf('docn_mode')
Expand All @@ -764,7 +798,7 @@ subroutine docn_comp_run( EClock, cdata, x2o, o2x)
close(nu)
call shr_file_freeUnit(nu)
endif
if (trim(ocn_mode) == 'SOM') then
if (trim(ocn_mode) == 'SOM' .or. trim(ocn_mode) == 'SOM_AQUAP') then
if (my_task == master_task) write(logunit,F04) ' writing ',trim(rest_file),currentYMD,currentTOD
call shr_pcdf_readwrite('write',iosystem,SDOCN%io_type,trim(rest_file),mpicom,gsmap,clobber=.true., &
rf1=somtp,rf1n='somtp')
Expand Down Expand Up @@ -879,8 +913,6 @@ subroutine prescribed_sst(xc, yc, lsize, sst_option, sst)
rlon(:) = xc(:) * pio180
rlat(:) = yc(:) * pio180

write(6,*)"DEBUG: sst_option is ",sst_option

! Control

if (sst_option < 1 .or. sst_option > 10) then
Expand Down

0 comments on commit 10fbc43

Please sign in to comment.