Main menu

Navigation

Decrease precipitation by 10% in CLM5 + DATM

3 posts / 0 new
Last post
thiagods@...
Decrease precipitation by 10% in CLM5 + DATM

I am running the "I2000Clm50BgcCruGs" CLM5 case with Data Atmosphere (DATM). Specifically, I am attempting to run a control case with default settings and two experiment cases with 10% and 20% less precipitation.

Rather than changing the forcing data itself, I am trying to temporarily change the code excerpt responsible for reading the atmospheric data. My first approach was to change lines 765 to 769 of the file `cime/src/components/data_comps/datm/datm_comp_mod.F90` from

!-------------------------------------------------------------------------

! PRECIPITATION DATA

!-------------------------------------------------------------------------

avstrm%rAttr(sprec,n) = avstrm%rAttr(sprec,n)/86400.0_R8 ! convert mm/day to kg/m^2/s


to


!-------------------------------------------------------------------------

! PRECIPITATION DATA

!-------------------------------------------------------------------------

avstrm%rAttr(sprec,n) = avstrm%rAttr(sprec,n)/86400.0_R8 ! convert mm/day to kg/m^2/s

avstrm%rAttr(sprec,n) = avstrm%rAttr(sprec,n)*0.9_R8 ! decrease precip by 10%


However, that change caused no impact in the model output. Am I looking at the right file here? Should I look at the "cpl" or "lnd" components instead?

Any hints to help me find the right file/routine are highly appreciated. Thanks in advance.

oleson

I think you are changing the wrong lines.  You are changing the lines within the case:

case('CORE2_NYF','CORE2_IAF')

 

You need to change lines within this case (for CRUNCEP forcing):

case('CLMNCEP')

thiagods@...
Thank you very much Keith, that worked.
Log in or register to post comments

Who's new

  • 408251463@...
  • 1161561668@...
  • 645505358@...
  • lbqnj@...
  • q1034908455@...