You can find instructions on how to install the coupled LMDZOR model here. This page assumes you already have the model checked out, and simply want to activate the 11 hydrology, stomate, and/or the new physics. This all comes from tests done by Fuxing Wang without stomate. If you are not using stomate, you need to use an LAI map.

River routing is also deactivated here.

Most of the changes are done in the configuration files after the model is compiled, with the exception of the following.

Source code changes

New Physics

Note that you must use revision 1792 of LMD! The default revision checked out from svn is revision 1729. To update, go to modeles/LMDZ and type

svn update -r 1792

After that, you need to change the following lines:

LMDZ/libf/phylmd/fisrtilp.F90

Original

              zcond(i) = MAX(0.0,zqn(i)-zqs(i))/(1+zdqs(i))
              rhcl(i,k)=1.0
           ELSE
              zcond(i) = MAX(0.0,zqn(i)-zqs(i))*rneb(i,k)/(1+zdqs(i))

Modified:

              zcond(i) = MAX(0.0,zqn(i)-zqs(i))
              rhcl(i,k)=1.0
           ELSE
              zcond(i) = MAX(0.0,zqn(i)-zqs(i))*rneb(i,k)

11 layer hydrology (CWRR)

In modeles/LMDZ/libf/phylmd/surf_land_orchidee_mod.F90:

Original

    zlev(1:knon) = (100.*plev(1:knon))/((ps(1:knon)/RD*temp_air(1:knon))*RG)

Modified

    zlev(1:knon) = plev(1:knon)*RD*temp_air(1:knon)/((ps(1:knon)*100.0)*RG)

In modeles/LMDZ/libf/phylmd/coefcdrag.F90:

Original

       ztsolv(i) = ts(i)
       ztvd(i) = t(i) * (psol(i)/pref(i))**RKAPPA
       trm0(i) = 1. + RETV * max(qsurf(i),0.0)
       trm1(i) = 1. + RETV * max(q(i),0.0)
       ztsolv(i) = ztsolv(i) * trm0(i)
       ztvd(i) = ztvd(i) * trm1(i)
       zri1(i) = zdphi(i)*(ztvd(i)-ztsolv(i))/(zdu2(i)*ztvd(i))
 

Modified

       ztsolv(i) = ts(i)
!       ztvd(i) = t(i) * (psol(i)/pref(i))**RKAPPA
       trm0(i) = 1. + RETV * max(qsurf(i),0.0)
       trm1(i) = 1. + RETV * max(q(i),0.0)
       ztsolv(i) = ztsolv(i) * trm0(i)
!       ztvd(i) = ztvd(i) * trm1(i)
       ztvd(i) = (t(i)+zdphi(i)/RCPD/(1.+RVTMP2*q(i))) &
          *(1.+RETV*q(i))
       zri1(i) = zdphi(i)*(ztvd(i)-ztsolv(i))/(zdu2(i)*ztvd(i))

Once you have done this, you can create the makefile and compile as normal (found elsewhere on the Wiki).

Configuration Files

These changes are done in the directory where you will be submitting the job (config/LMDZOR_v5.2/JOBNAME, for example).

11 layer hydrology (CWRR)

Confirm the following lines are present. If a line exists but with something different after the = sign, then change it.

In PARAM/orchidee.def_CWRR:

   RIVER_ROUTING = n

In COMP/orchidee.card:

   DefSuffix=CWRR

New physics

Confirm the following lines are present.

In PARAM/orchidee.def_CWRR:

   evapot_corr=y

In COMP/lmdz.card:

   LMDZ_Physics=NPv3.1

In PARAM/physiq.def_L39_NPv3.1 (assuming you are using 39 levels, i.e. you compile with gmake LMD????-L39):

 iflag_thermals=17
 iflag_pbl=11 

STOMATE

In PARAM/orchidee.def_CWRR:

STOMATE_OK_STOMATE = y

In config.card:

SBG= (stomate, ORCHIDEE_1_9_5)

In COMP/stomate.card (useful mostly for debugging):

[OutputText]
List=   (out_orchidee, run.def)

No STOMATE

In PARAM/orchidee.def_CWRR:

LAI_MAP = y
STOMATE_OK_STOMATE = n

In COMP/orchidee.card:

ListNonDel= (${R_BC}/SRF/${config_UserChoices_TagName}/lai2D.nc, .)

In config.card, remove the following line:

SBG= (stomate, ORCHIDEE_1_9_5)